Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,732,696 +C 100% intergenic (‑120/‑283) RCAP_rcc01601 ← / → tgt FeS assembly SUF system protein/queuine tRNA‑ribosyltransferase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,732,6961.→C100.0% 75.6 234intergenic (‑120/‑283)RCAP_rcc01601/tgtFeS assembly SUF system protein/queuine tRNA‑ribosyltransferase
Reads supporting (aligned to +/- strand):  new base (39/168):  ref base (21/6):  total (60/174)
Fisher's exact test strand distribution p-value = 1.67e-09
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.19e-12

Only 100 of 289 total aligned reads displayed.

CTAGAGTTTATATAGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.A.CC..GG.G.C..GGG..CGG..G...GCCCCCGCACCCCCGGCGTATTTGC.GCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCAGCCGCATGCCGTCC  >  NC_014034/1732596‑1732784
                                                                                                       |                                                                                                      
 tAGAGTTTATATAGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.A.c                                                                                                       <  2:4500392/101‑1 (MQ=255)
 tAGAGTTTATATAGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.A.c                                                                                                       <  2:4549614/101‑1 (MQ=255)
   gagTTTATATAGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.A.ccc                                                                                                     <  1:2668154/101‑2 (MQ=255)
    agTTTATATAGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.A.cccg                                                                                                    <  2:3129685/101‑3 (MQ=255)
    agTTTATATAGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.A.cccg                                                                                                    <  1:7973257/101‑3 (MQ=255)
     gTTTATATAGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.A.Cc..cg.g                                                                                                >  1:1643816/1‑101 (MQ=255)
       ttATATAGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.A.Cc..cg.ggc                                                                                              >  2:5100299/1‑99 (MQ=255)
       ttATATAGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCCCGGCGGGGCCCCGGC.G.A.Cc..cg.ggg                                                                                              >  1:223/1‑99 (MQ=255)
         atataGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.A.Cc..cg.g.g..ggg                                                                                         >  1:2331098/1‑101 (MQ=255)
          tataGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGTCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.A.Cc..cg.g.g..gggg                                                                                        >  1:7329098/1‑100 (MQ=255)
           ataGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.A.Cc..cg.gggcggg                                                                                          >  2:7885396/1‑95 (MQ=255)
            taGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.ACCc..gg.g.c..ggggg                                                                                       <  2:5726826/101‑3 (MQ=255)
            taGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.ACCc..gg.g.c..ggggg                                                                                       >  1:4791151/1‑99 (MQ=255)
             aGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGCTTCGCCGCGGCGGGGCCGCGGC.G.A.Cc..cg.g.c..gcggggc                                                                                     >  1:9001388/1‑95 (MQ=255)
              gggTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.ACCc..gg.g.c..gggggct                                                                                     <  2:5247596/101‑5 (MQ=255)
                gTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.ACCc..gg.g.c..gggggctcc                                                                                   <  1:1511753/101‑7 (MQ=255)
                gTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.ACCc..gg.g.c..gggggctcc                                                                                   <  1:859805/101‑7 (MQ=255)
                gTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.ACCc..gg.g.c..gggggctcc                                                                                   <  2:4146603/101‑7 (MQ=255)
                gTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.ACCC..GG.G.C..GGG...GG..G...Gcc                                                                           >  1:1158653/1‑101 (MQ=255)
                 tATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGCGG.C.CC..GG.G.C..GGG..GGG..C...Gcg                                                                           >  2:6414545/1‑100 (MQ=255)
                 tATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.ACCC..GG.G.C..GGG...GG..G...Gccg                                                                          >  1:3569021/1‑100 (MQ=255)
                 caTCGGGGGGAGAGGTAAAGGCCGTCGGTGCGCGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.GCC.Cc..gg.g.c..gggggctccg                                                                                  <  2:2888033/100‑8 (MQ=255)
                  aTCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.A.CC..CG.G.G..GGG..Ggg..g...ggcgc                                                                         >  1:2528202/1‑97 (MQ=255)
                  aTCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCAGCGTGCCAAGCAGGGTTTCGACGCGGCGGGGACGCAGG.G.A.CC..CG.G.T..CGG..GGG..G...Gccaa                                                                         >  2:4884481/1‑99 (MQ=255)
                    cGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.ACCC..GG.G.G..GGG..GGG..G...Gcgccc                                                                        >  1:6770678/1‑101 (MQ=255)
                    cGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.ACCC..GG.G.C..AGG...GG..G...Gccggcc                                                                       >  2:832119/1‑97 (MQ=255)
                    cGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.A.CC..CG.G.G..GGG..GGG..G...Gccgccc                                                                       >  2:2328319/1‑100 (MQ=255)
                     ggggggAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.ACCC..GG.GCC..GGG..GGG..G...Ggcccc                                                                        >  1:5881036/1‑101 (MQ=255)
                     ggggggAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.A.Cc..cg.g.g..ggg..gggcgggggccc                                                                           >  1:5714952/1‑92 (MQ=255)
                      gggggAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGG.G.ACCC..GG.G.C..GGG..GGG..G...Cccccccc                                                                      >  2:7477722/1‑99 (MQ=255)
                      gggggAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.ACCC..GG.G.C..GGG..GGG..G...Gccccccc                                                                      >  1:4534886/1‑99 (MQ=255)
                      gggggAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.ACCC..GG.G.C..GGG..GGG..G...Cccccccc                                                                      >  1:9144233/1‑99 (MQ=255)
                      gggggAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.A.CCCGGG.G.C..GGG..Ggg..g...gggcccc                                                                       >  1:3523465/1‑95 (MQ=255)
                       ggggAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.ACCC..GG.G.C..GGG..GGG..T...Gcccccccc                                                                     >  1:2965753/1‑98 (MQ=255)
                        gggAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.ACCC..GG.G.C..GGG...GG..G...Gccggcccccc                                                                   >  2:3736332/1‑93 (MQ=255)
                          gagaGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.ACCC..GG.G.C..GGG..GGG..G...GCGCCCCCCccc                                                                  >  2:67034/1‑101 (MQ=255)
                             agGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.ACCC..GG.G.CGGGGG..CTC..C...GCCCCCGCAcccc                                                                 <  1:1938919/101‑1 (MQ=25)
                             agGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.ACCC..GG.G.C..GGG..GGG..G...CCGCCCCCCcccccc                                                               >  2:5470271/1‑100 (MQ=255)
                              ggTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.ACCC..GG.G.CGGGGG..CTC..C...GCCCCCGCAccccc                                                                <  1:6026382/101‑1 (MQ=25)
                               gTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.ACCC..GG.G.CGGGGG..CTC..C...GCCCCCGCACCCCCg                                                               <  1:832132/101‑1 (MQ=25)
                               gTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.ACCC..GG.G.C..GGG..GGG..GCC.GCCCCCGCCCCACCg                                                               >  2:6015676/1‑101 (MQ=255)
                                tAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.ACCC..GG.G.CGGGGG..CTC..C...GCCCCCGCACCCCCgg                                                              <  2:4827655/101‑1 (MQ=25)
                                 aaaGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.ACCC..GG.G.CGGGGG..CTC..C...GCCCCCGCACCCCCGGc                                                             <  1:2103293/101‑1 (MQ=25)
                                 aaaGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.ACCC..GG.G.C..GGG..GGGTCC...GGCCCCGCACCCCCGGc                                                             >  2:4200446/1‑101 (MQ=255)
                                     gCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.ACCC..GG.G.CGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTAt                                                         <  1:7183287/101‑1 (MQ=25)
                                     gCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.ACCC..GG.G.CGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTAt                                                         <  1:5192050/101‑1 (MQ=25)
                                     gCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.ACCC..GG.G.CGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTAt                                                         <  1:5592399/101‑1 (MQ=25)
                                      ccGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.ACCC..GG.G.CGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTAtt                                                        <  2:8831134/101‑1 (MQ=25)
                                       cggcGGTGCGGGGCAGAAAGCCGCGTGCCAAGAAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.ACCC..GG.G.CGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTAttt                                                       <  2:1032939/98‑1 (MQ=18)
                                       cGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.ACCC..GG.G.CGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTAttt                                                       <  1:5546090/101‑1 (MQ=25)
                                       cGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.ACCC..GG.G.CGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTAttt                                                       <  2:8936504/101‑1 (MQ=25)
                                       cGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.ACCC..GG.G.CGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTAttt                                                       <  2:2372562/101‑1 (MQ=25)
                                          cGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.ACCC..GG.G.CGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTgc.g                                                   <  2:8541410/101‑1 (MQ=21)
                                           ggTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.ACCC..GG.G.CGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTgc.gc                                                  <  1:3214200/101‑1 (MQ=25)
                                            gTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.ACCC..GG.G.CGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGC.Gcc                                                 <  1:1786473/101‑1 (MQ=21)
                                             tGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.ACCC..GG.G.CGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGC.GCCa                                                <  2:3395766/101‑1 (MQ=21)
                                             ggCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.ACCC..GG.G.CGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGC.GCCa                                                <  2:2996180/100‑1 (MQ=21)
                                              gCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.ACCC..GG.G.CGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGC.GCCaa                                               <  2:3297784/101‑1 (MQ=21)
                                               cGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.ACCC..GG.G.CGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGC.GCCaag                                              <  2:9051217/101‑1 (MQ=21)
                                               cGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.ACCC..GG.G.CGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGC.GCCaag                                              <  1:3996934/101‑1 (MQ=21)
                                                ggggCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.ACCC..GG.G.C..GGG..GGGCCGC..GCCCCCCCCCCCCCCGCG.ACTTGC.GCCaaca                                             >  2:2505423/1‑99 (MQ=255)
                                                    cGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.ACCC..GG.G.CGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGC.GCCAAgaagaa                                         <  2:3568913/101‑1 (MQ=18)
                                                    cGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.ACCC..GG.G.CGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGC.GCCAAgaagaa                                         <  2:2612523/101‑1 (MQ=18)
                                                     ggAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.ACCC..GG.G.CGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGC.GCCAAGAAGaaa                                        <  1:9007783/101‑1 (MQ=18)
                                                      gAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.ACCC..GG.G.CGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGC.GCCAAGAAGAAAc                                       <  2:7035830/101‑1 (MQ=18)
                                                      gAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.ACCC..GG.G.CGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGC.GCCAAGAAGAAAc                                       <  1:2471666/101‑1 (MQ=18)
                                                       gaaGCCGCGTGCCAAGCAGGGGTTCGCCGCGGCGGGGCCGCGGC.GCC.CC..GG.G.CGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGC.GCCAAGAAGAAACa                                      <  2:4216946/100‑1 (MQ=18)
                                                       aaaGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.ACCC..GG.G.CGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGC.GCCAAGAAGAAACa                                      <  1:940709/101‑1 (MQ=18)
                                                       aaaGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.ACCC..GG.G.CGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGC.GCCAAGAAGAAACa                                      <  1:313337/101‑1 (MQ=18)
                                                        aaGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.ACCC..GG.G.CGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGC.GCCAAGAAGAAACaa                                     <  2:875586/101‑1 (MQ=18)
                                                        aaGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.ACCC..GG.G.CGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGC.GCCAAGAAGAAACaa                                     <  1:6354158/101‑1 (MQ=18)
                                                        aaGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.ACCC..GG.G.CGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGC.GCCAAGAAGAAACaa                                     <  1:6193790/101‑1 (MQ=18)
                                                        aaGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.ACCC..GG.G.CGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGC.GCCAAGAAGAAACaa                                     <  2:8220967/101‑1 (MQ=18)
                                                          gCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCTCGGC.G.ACCC..GG.G.CGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGC.GCCAAGAAGAAACAAGc                                   <  1:8161434/101‑1 (MQ=18)
                                                          gCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.ACCC..GG.G.CGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGC.GCCAAGAAGAAACAAGc                                   <  1:7696852/101‑1 (MQ=18)
                                                            cgcgTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.ACCC..GG.G.CGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGC.GCCAAGAAGAAACAAGCgg                                 <  1:2893085/101‑1 (MQ=18)
                                                                tGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.ACCC..GG.G.CGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGC.GCCAAGAAGAAACAAGCGGGAcc                             <  1:5140001/101‑1 (MQ=18)
                                                                  ccAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.ACCC..GG.G.CGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGC.GCCAAGAAGAAACAAGCGGGACCAg                           <  2:771328/101‑1 (MQ=18)
                                                                   cAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.ACCC..GG.G.CGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGC.GCCAAGAAGAAACAAGCGGGACCAgg                          <  2:3766147/101‑1 (MQ=18)
                                                                   cAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.ACCC..GG.G.CGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGC.GCCAAGAAGAAACAAGCGGGACCAgg                          <  2:1716176/101‑1 (MQ=18)
                                                                    aaGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.A.CC..CG.G.G..GGG..GGG..GTCCCCCCCCCCACCCCCCGCGGATTTGC.GCCCAGAAGAAAAAAGCGGGACCCGGc                         >  1:3895827/1‑101 (MQ=255)
                                                                     aGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.ACCC..GG.G.C..GGG..GGG..GCC.CCCCCCCCCCCCCCCGCGCATTTGCGGCCAAGAAG.AACAAACGGGAcccagcc                        >  2:7293894/1‑96 (MQ=255)
                                                                        aGGGTTTCGCCGCGGCGGGGCCGCGGC.G.ACCC..GG.G.CGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGC.GCCAAGAAGAAACAAGCGGGACCAGGccccc                     <  2:926006/101‑1 (MQ=18)
                                                                         gggTTTCGCCGCGGCGGGGCCGCGGC.G.ACCC..GG.G.CGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGC.GCCAAGAAGAAACAAGCGGGACCAGGcccccc                    <  1:1221964/101‑1 (MQ=18)
                                                                         gggTTTCGCCGCGGCGGGGCCGCGGC.G.ACCC..GG.G.CGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGC.GCCAAGAAGAAACAAGCGGGACCAGGcccccc                    <  2:1343097/101‑1 (MQ=18)
                                                                         gggTTTCGCCGCGGCGGGGCCGCGGC.G.A.CCCGGGTG.G..GGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGC.GCCAAGAAGAAACAAGCGGGACCAGGcccccc                    <  2:1236038/101‑1 (MQ=18)
                                                                           gTTTCGCCGCGGCGGGGCCGCGGC.G.ACCC..GG.G.CGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGC.GCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTg                  <  2:6888674/101‑1 (MQ=18)
                                                                             ttCGCCGCGGCGGGGCCGCGGC.G.ACCC..GG.G.CGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGC.GCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCa                <  1:3675185/101‑1 (MQ=18)
                                                                             ttCGCCGCGGCGGGGCCGCGGC.G.ACCC..GG.G.CGGGGG..CTC..C...CCCCCCGCACCCCCGGCGTATTTGC.GCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCa                <  2:9000966/101‑1 (MQ=25)
                                                                             ttCGCCGCGGCGGGGCCGCGGC.G.ACCC..GG.G.C..GGG..GGG..GCC.GCCCCCCCCCCCCCCGCGGATTTTC.GGCAACCAGAAAAAAACGGGGCCCGGgccccccccc                >  2:1814001/1‑97 (MQ=255)
                                                                                gccgcGGCGGGGCCGCGGC.G.ACCC..GG.G.CGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGC.GCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCAGcc             <  2:8624204/101‑1 (MQ=18)
                                                                                gccgcGGCGGGGCCGCGGC.G.ACCC..GG.G.CGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGC.GCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCAGcc             <  1:4210938/101‑1 (MQ=18)
                                                                                  cgcgGGGGGGCCGCGGC.G.ACCC..GG.G.C..GGGGCCTC..C...GCCCCCGCACTCCCGGCGTATTTGC.GCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCAGCCGc           <  1:1395764/101‑1 (MQ=25)
                                                                                    cggcggGGCCGCGGC.G.ACCC..GG.G.CGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGC.GCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCAGCCGCAt         <  2:2769005/101‑1 (MQ=18)
                                                                                     ggcggGGCCGCGGC.G.ACCC..GG.G.CGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGC.GCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCAGCCGCATg        <  2:8387048/101‑1 (MQ=18)
                                                                                     ggcggGGCCGCGGC.G.ACCC..GG.G.CGGGGG..CTC..C...GCCCCCGCACCCCCCGCGTATTTGC.GCCAAGAAGAAAAAAGCGGGACCAGGCCCCCCCGCAGCCGCCTg        >  1:7318501/1‑101 (MQ=17)
                                                                                       cggGGCCGCGGC.G.ACCC..GG.G.CGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGC.GCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCAGCCGCATGcc      <  1:8814534/101‑1 (MQ=18)
                                                                                          ggCCGCGGC.G.ACCC..GG.G.CGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGC.GCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCAGCCGCATGCCGTc   <  2:2308389/101‑1 (MQ=18)
                                                                                          ggCCGCGGC.G.ACCC..GG.G.CGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGC.GCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCAGCCGCATGCCGTc   <  2:5356495/101‑1 (MQ=18)
                                                                                           gCCGCGGC.G.ACCC..GG.G.CGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGC.GCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCAGCCGCATGCCGTcc  <  1:7573093/101‑1 (MQ=18)
                                                                                                       |                                                                                                      
CTAGAGTTTATATAGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.A.CC..GG.G.C..GGG..CGG..G...GCCCCCGCACCCCCGGCGTATTTGC.GCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCAGCCGCATGCCGTCC  >  NC_014034/1732596‑1732784

Base quality scores:  ATCG  < 3 ≤  ATCG  < 16 ≤  ATCG  < 24 ≤  ATCG  < 33 ≤  ATCG  < 39 ≤  ATCG