Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,397,495 C→G 20.5% T459T (ACC→TCG RCAP_rcc01307 → hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,397,4950C→G20.5% 109.4 ‑1.8 ‑0.9210T459T (ACC→ACGRCAP_rcc01307hypothetical protein
Reads supporting (aligned to +/- strand):  new base (20/23):  ref base (98/68):  total (119/91)
Fisher's exact test strand distribution p-value = 1.68e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.20e-01

Only 100 of 244 total aligned reads displayed.

GAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGA  >  NC_014034/1397395‑1397595
                                                                                                       |                                                                                                    
gAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGcacc                                                                                                      >  2:3605571/1‑101 (MQ=255)
 aaaCTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGcacca                                                                                                     <  2:3868046/101‑1 (MQ=255)
  aaCTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGcaccac                                                                                                    <  1:468067/101‑1 (MQ=255)
  aaCTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGcaccac                                                                                                    >  1:2127539/1‑101 (MQ=255)
   aCTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCaccacc                                                                                                   >  1:3465451/1‑101 (MQ=255)
    ctctTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCg                                                                                                  >  2:3104602/1‑101 (MQ=255)
           aTGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTc                                                                                           >  1:2532386/1‑101 (MQ=255)
             ggATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTccc                                                                                         <  2:866112/101‑1 (MQ=255)
              gATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTcccc                                                                                        <  2:4143128/101‑1 (MQ=255)
               agatcgCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTccccc                                                                                       >  2:2919982/6‑101 (MQ=255)
                tatGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTcccccc                                                                                      >  1:1908813/1‑101 (MQ=255)
                 atGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCcg                                                                                     >  1:3268135/1‑101 (MQ=255)
                 atGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCcg                                                                                     <  1:822844/101‑1 (MQ=255)
                  tGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACAACCGGTCGCTCCCCCcgc                                                                                    >  2:1450740/1‑101 (MQ=255)
                   ggCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCcgcg                                                                                   >  1:3661065/1‑101 (MQ=255)
                   ggCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TC.TC.AACCCGCAA.CCGCCCCGCACCCCCGGGCGCCTCCCCCCCgg                                                                                  >  1:2762949/1‑101 (MQ=21)
                    gCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCgg                                                                                  <  1:2338710/101‑1 (MQ=255)
                    gCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCgg                                                                                  >  2:774519/1‑101 (MQ=255)
                      aaaaCCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCgggg                                                                                >  2:1106069/1‑101 (MQ=255)
                      aaaaCCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCgggg                                                                                >  1:2849111/1‑101 (MQ=255)
                      aaaaCCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCgggg                                                                                >  2:4509547/1‑101 (MQ=255)
                        aaCCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGc                                                                              <  2:2449850/101‑1 (MQ=255)
                          ccGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGt                                                                            >  1:1713864/1‑101 (MQ=255)
                           cGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTg                                                                           >  2:2003457/1‑101 (MQ=255)
                              cgataccgcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTgg                                                                          <  2:520386/92‑1 (MQ=0)
                              ccGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGAt                                                                        <  1:1020711/101‑1 (MQ=255)
                              ccGCGACATGGCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGAt                                                                        <  2:1589370/101‑1 (MQ=255)
                               cgcgACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                       >  2:960085/1‑101 (MQ=255)
                                ataccgcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCTCCCGCGGGGGCGTGGAt                                                                        >  2:3783012/8‑101 (MQ=0)
                                 taccgcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGATc                                                                       >  1:567312/7‑101 (MQ=2)
                                 cgACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                                     >  1:1000454/1‑101 (MQ=255)
                                  accgcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                      <  2:2069985/96‑1 (MQ=2)
                                   ccgcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                                     >  1:2046598/5‑101 (MQ=2)
                                    cATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAaac                                                                  >  2:4058662/1‑101 (MQ=255)
                                    cATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAaac                                                                  <  1:243507/101‑1 (MQ=255)
                                    cATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAaac                                                                  >  2:2535641/1‑101 (MQ=255)
                                     gcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                                   >  1:3124721/3‑101 (MQ=2)
                                      tGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAaacaa                                                                >  2:4246922/1‑101 (MQ=255)
                                      cgtccgtCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAaac                                                                  <  2:79404/95‑1 (MQ=2)
                                      cgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAaac                                                                  >  1:4422136/2‑101 (MQ=2)
                                      cgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAaac                                                                  <  2:4048450/100‑1 (MQ=2)
                                       gTCGCTCCCCCCGCGGGGCCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAcacaac                                                               >  2:2290384/1‑101 (MQ=255)
                                         cGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGGGGATCGAaacggt                                                               >  2:4436945/1‑98 (MQ=2)
                                         cGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCAGCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAccc                                                             <  1:960162/101‑1 (MQ=255)
                                         cGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAccc                                                             >  1:1470355/1‑101 (MQ=255)
                                          gCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCg                                                            <  1:3104803/101‑1 (MQ=255)
                                           cTCCCCCCGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAaacggtgt                                                             >  2:1868245/1‑96 (MQ=2)
                                           cTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGa                                                           <  2:599689/101‑1 (MQ=255)
                                            tCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGAt                                                          >  1:2248958/1‑101 (MQ=255)
                                             ccccccGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAaacggtgtcc                                                           >  1:625758/1‑94 (MQ=2)
                                              cccccGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAaacggtgtcca                                                          <  2:2626139/101‑9 (MQ=2)
                                              cccccGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAaacggtgtcca                                                          <  2:2327641/101‑9 (MQ=2)
                                              cccccGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAaacggtgtcca                                                          <  1:1271523/101‑9 (MQ=2)
                                              cccccGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAaacggtgtcca                                                          <  1:3502971/101‑9 (MQ=2)
                                              cccccGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCa                                                        >  1:4256316/1‑101 (MQ=255)
                                                cccGCGGGGGCGTGGATCGAAACATCGATGATCAAATCCGCAATCGGCCGGGCTCGATC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAaacggtgtccagc                                                        <  2:4673500/101‑11 (MQ=2)
                                                cccGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAaacaaagatcggaag                                                      >  2:176320/1‑91 (MQ=255)
                                                cccGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATg                                                      <  1:2091823/101‑1 (MQ=255)
                                                   gcgGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGcc                                                   <  2:4405685/101‑1 (MQ=255)
                                                   gcgGGGGCGTGGATCGAAACACCG.CCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGcc                                                   <  2:2794916/101‑1 (MQ=255)
                                                    cgGGGGCGTGGATCGAAACGCCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGccc                                                  <  2:3634172/101‑1 (MQ=255)
                                                    cgGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGccc                                                  >  2:3392606/1‑101 (MQ=255)
                                                     gggggCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCg                                                 >  2:197033/1‑101 (MQ=255)
                                                        ggCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCt                                              <  2:2649471/101‑1 (MQ=255)
                                                        ggCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCt                                              >  1:1320637/1‑101 (MQ=255)
                                                        ggCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCt                                              >  2:1066458/1‑101 (MQ=255)
                                                         gCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCCCCGGTCCCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTc                                             <  1:3901285/101‑1 (MQ=255)
                                                           gTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAGACAACCCGATCATGGCCCGCCTCAg                                           >  2:3040432/1‑101 (MQ=255)
                                                           gTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAg                                           <  1:2726043/101‑1 (MQ=255)
                                                           gTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAg                                           >  1:3693404/1‑101 (MQ=255)
                                                            tGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAgc                                          >  2:240149/1‑101 (MQ=255)
                                                               aTCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGcgcc                                       <  1:1106155/101‑1 (MQ=255)
                                                                  gAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCcgcg                                    >  1:1068291/1‑101 (MQ=255)
                                                                  gAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCcgcg                                    >  1:2627197/1‑101 (MQ=255)
                                                                   aaaCACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCgcgc                                   >  2:3320444/1‑101 (MQ=255)
                                                                      cacCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGa                                >  1:2991073/1‑101 (MQ=255)
                                                                       acCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCACCGCCGCGCTGAc                               >  1:3786/1‑101 (MQ=255)
                                                                         cG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGt                             >  2:1574541/1‑101 (MQ=255)
                                                                          g.TCCTC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTc                            <  1:2453368/101‑1 (MQ=255)
                                                                            tcatc.aCCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCg                           >  1:1623952/1‑101 (MQ=255)
                                                                             catc.aCCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGc                          >  1:4741389/1‑101 (MQ=255)
                                                                             catc.aCCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGc                          <  1:115198/101‑1 (MQ=255)
                                                                             catc.aCCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGc                          >  2:4582071/1‑101 (MQ=255)
                                                                              atc.aCCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCt                         >  2:252274/1‑101 (MQ=255)
                                                                                c.accacGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTcc                       >  1:4667190/1‑101 (MQ=255)
                                                                                c.accacGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTcc                       <  1:4629726/101‑1 (MQ=255)
                                                                                  accacGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTccc                      <  1:329581/101‑1 (MQ=255)
                                                                                     ccGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTcccccc                   <  1:4735339/100‑1 (MQ=255)
                                                                                     acGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAAAACCCGATCATCGACCGCCTCAGCGGCGCGCTGATGTCGCTtccccc                   >  2:2453389/1‑101 (MQ=37)
                                                                                      cGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTCCCCCcg                  <  2:1696223/101‑1 (MQ=255)
                                                                                        cAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTCCCCCcgcg                <  2:3751349/101‑1 (MQ=255)
                                                                                         aa.CCGCCCTGCACCACCGGTCGCGCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCgg               <  2:1143827/101‑1 (MQ=255)
                                                                                              gCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGc           <  1:2570571/101‑1 (MQ=255)
                                                                                                ccTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGt         <  1:2920158/101‑1 (MQ=255)
                                                                                                     accaccGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATc    <  1:774593/101‑1 (MQ=255)
                                                                                                      ccaccGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCg   <  1:928142/101‑1 (MQ=255)
                                                                                                      ccaccGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCg   >  2:2239793/1‑101 (MQ=255)
                                                                                                      ccaccGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCg   >  2:1551603/1‑101 (MQ=255)
                                                                                                       caccGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGa  >  1:4591815/1‑101 (MQ=255)
                                                                                                       caccGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGa  >  2:3233188/1‑101 (MQ=255)
                                                                                                       |                                                                                                    
GAAACTCTTCAATGGATATGGCAAAACCGTCCGCGACATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCG.TCATC.ACCACGCAA.CCGCCCTGCACCACCGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAACCCGATCATGGCCCGCCTCAGCGCCGCGCTGACGTCGCTCCCCCCGCGGGGGCGTGGATCGA  >  NC_014034/1397395‑1397595

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 27 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG