Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,732,696 +C 91.9% 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.→C91.9% 24.6 ‑3.9 ‑1.9123intergenic (‑120/‑283)RCAP_rcc01601/tgtFeS assembly SUF system protein/queuine tRNA‑ribosyltransferase
Reads supporting (aligned to +/- strand):  new base (18/95):  ref base (5/5):  total (23/100)
Fisher's exact test strand distribution p-value = 1.98e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.21e-02

Only 100 of 158 total aligned reads displayed.

TAGAGTTTATATAGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.A..CC..GG.GC..GGG..CGG..G...GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAAC.AAGCGGGACCAGGCCCCCCTGCAGCCGCATGCCGTCC  >  NC_014034/1732597‑1732784
                                                                                                      |                                                                                                      
tAGAGTTTATATAGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.A..c                                                                                                      >  2:863953/1‑101 (MQ=255)
 agagTTTATATAGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGGCGCGGC.G.A..cc                                                                                                     >  2:1427672/1‑101 (MQ=255)
 agagTTTATATAGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.A..cc                                                                                                     >  2:786624/1‑101 (MQ=255)
 agagTTTATATAGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.A..cc                                                                                                     <  2:2394280/101‑1 (MQ=255)
 agagTTTATATAGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.A..cc                                                                                                     >  1:4977044/1‑101 (MQ=255)
  gagTTTATATAGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.C..ccc                                                                                                    >  1:6638051/1‑100 (MQ=255)
  gagTTTATATAGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.A..ccc                                                                                                    <  2:3028388/101‑2 (MQ=255)
  gagTTTATATAGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.A..ccc                                                                                                    >  1:4037539/1‑100 (MQ=255)
   agTTTATATAGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.A..cccg                                                                                                   >  2:2188361/1‑99 (MQ=255)
   agTTTATATAGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.A..cccg                                                                                                   <  2:1716404/101‑3 (MQ=255)
     tttATATAGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.A..Cc..cg.gg                                                                                              >  2:3915445/1‑100 (MQ=255)
      ttATATAGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCGGCGGC.G.A..Cc..cg.ggg                                                                                             >  1:4645419/1‑99 (MQ=255)
      ttATATAGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.A..Cc..cg.ggc                                                                                             <  1:3342181/101‑3 (MQ=255)
      ttATATAGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.A..Cc..cg.ggc                                                                                             >  1:3981787/1‑99 (MQ=255)
       tatgtaGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGGTTCGCCGCGGCCGGGCCGCGGC.G.A..ccccggcg                                                                                               >  1:6840622/1‑95 (MQ=255)
       tatataGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.AC.Cc..gg.gc..g                                                                                           >  1:1577619/1‑101 (MQ=255)
       tatataGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.AC.Cc..gg.gc..g                                                                                           >  2:6237839/1‑101 (MQ=255)
        atataGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCTGCGGGGCCGCGGC.G.AC.Cc..gg.gc..gg                                                                                          >  1:6000359/1‑101 (MQ=255)
        atataGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.AC.Cc..gg.gc..gg                                                                                          >  2:1532159/1‑101 (MQ=255)
        atataGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.AC.Cc..gg.gc..gg                                                                                          <  1:1618774/101‑1 (MQ=255)
        atataGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.A..Cc..cg.gg..ggg                                                                                         >  1:2462177/1‑101 (MQ=255)
        atataGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGAAGGGTTTCGCCGCGGCGGGGCGGCGGG.G.A..Cc..cg.gg..ggg                                                                                         >  1:2277572/1‑101 (MQ=255)
          ataGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.AC.Cc..gg.gc..gggg                                                                                        >  2:3933041/1‑100 (MQ=255)
          ataGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.AC.Cc..gg.gc..gggg                                                                                        >  2:5511221/1‑100 (MQ=255)
           taGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.AC.Cc..gg.gc..ggggg                                                                                       <  2:253174/101‑3 (MQ=255)
            aGGGTATGGGGGGGAGAGGTAAAGGCGGTCGGTGCGGGGGGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.AC.Cc..gg.gc..gggggc                                                                                      <  2:1052322/101‑4 (MQ=255)
            aGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.AC.Cc..gg.gc..gggggc                                                                                      <  2:6084043/101‑4 (MQ=255)
            aGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.AC.Cc..gg.gc..ggg..ggg                                                                                    >  1:191976/1‑101 (MQ=255)
             gggTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.AC.Cc..gg.gc..gggggct                                                                                     <  2:5532921/101‑5 (MQ=255)
             gggTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.AC.Cc..gg.gc..gggggct                                                                                     <  2:5952168/101‑5 (MQ=255)
              ggTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.AC.Cc..gg.gc..gggggctc                                                                                    <  1:5591348/101‑6 (MQ=255)
              ggTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.AC.Cc..gg.gc..gggggctc                                                                                    <  2:3759192/101‑6 (MQ=255)
               gTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCTGCGGGGCCGCGGC.G.AC.Cc..gg.gc..gggggctcc                                                                                   <  2:6000054/101‑7 (MQ=255)
               gTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.AC.Cc..gg.gc..gggggctcc                                                                                   <  2:1407920/101‑7 (MQ=255)
               gTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.A..CC..CG.GG..GGG..Ggg..g...ggc                                                                           >  2:2824573/1‑99 (MQ=255)
               gTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGAGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.AC.Cc..gg.gc..gggggctcc                                                                                   <  2:4063004/101‑7 (MQ=255)
                tATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.AC.Cc..gg.gc..gggggctccg                                                                                  <  2:4721922/101‑8 (MQ=39)
                tATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGACGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.AC.CC..GG.GC..GGG..GGG..G...Gcg                                                                           >  2:1013727/1‑100 (MQ=255)
                tATCGGGGGGAGAGGTAAAGGCCGGCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCACGGTTTCGCCGCGGCGGGGCCGCAGC.G.C..CC..CG.CG..CGG..CGG..G...Gccg                                                                          >  2:4010455/1‑100 (MQ=255)
                 aTCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.AC.Cc..gg.gc..gggggctccgc                                                                                 <  2:2695812/101‑9 (MQ=38)
                 aTCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.A..CC..CG.GG..GGG..GGG..G...Gccgc                                                                         >  1:6424062/1‑99 (MQ=255)
                  tCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.AC.Cc..gg.gc..gggggctccgcc                                                                                <  1:1145715/101‑10 (MQ=38)
                   cGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.AC.Cc..gg.gc..gggggctccgccc                                                                               <  2:6487519/101‑11 (MQ=38)
                   cGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.A..CC..CG.GG..GGG..GGG..G...Gccgccc                                                                       >  2:6700189/1‑100 (MQ=255)
                     gggggAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.AC.CC..GG.GC..GGG..GGG..G...Gcgccccc                                                                      >  2:4465567/1‑99 (MQ=255)
                      ggggAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.ACCCc..gg.gc..ggg..ggg..gcggccccc                                                                         >  2:1460091/1‑93 (MQ=255)
                         gagaGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.AC.CC..GG.GC..GGG...GG..G...GCCGCCCCCcccc                                                                 >  1:1261984/1‑101 (MQ=255)
                          agagGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.AC.CC..GG.GC..GGG...GG..G...GCCCGCCCCccccc                                                                >  2:637309/1‑101 (MQ=255)
                             ggTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.AC.CC..GG.GCGGGGG..CTC..C...GCCCCCGCAccccc                                                                <  2:5363609/101‑1 (MQ=25)
                             ggTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.AC.CC..GG.GCGGGGG..CTC..C...GCCCCCGCAccccc                                                                <  1:1877549/101‑1 (MQ=25)
                             ggTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.AC.CC..GG.GCGGGGG..CTC..C...GCCCCCGCAccccc                                                                <  1:2961222/101‑1 (MQ=25)
                              gTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.AC.CC..GG.GCGGGGG..CTC..C...GCCCCCGCACCCCCg                                                               <  1:5485686/101‑1 (MQ=25)
                                aaaGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.AC.CC..GG.GC..GGG..GGC..G...GCCGCCCCCccccccccgc                                                           >  2:1716318/1‑96 (MQ=255)
                                aaaGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.A..CCCGGG.GC..GGG..GGG..GC..GCCCCCCCCcccccccc                                                             >  2:6132646/1‑98 (MQ=255)
                                  aGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.AC.CC..GG.GCGGGGG..CTC..C...GCCCCCGCACCCCCGGCGt                                                           <  1:3227348/101‑1 (MQ=25)
                                    gCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.AC.CC..GG.GCGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTAt                                                         <  2:3526686/101‑1 (MQ=25)
                                     ccGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.AC.CC..GG.GCGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTAtt                                                        <  1:393330/101‑1 (MQ=25)
                                          ggTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.AC.CC..GG.GCGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTgcgc                                                   <  1:818411/101‑1 (MQ=21)
                                           gTGCGGGGCGGAAAGCCGCGTGCCAAGGAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.AC.CC..GG.GCGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGCGcc                                                  <  2:4233433/101‑1 (MQ=21)
                                           gTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.AC.CC..GG.GCGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGCGcc                                                  <  2:1978579/101‑1 (MQ=21)
                                              gggggCGGAAAGCCGAGTGCCCAGCAGGGTTTCTCCGCGGCGTGGCCGCGGC.GCC..CC..GG.GC..GGGGCCTC..C...GCTCCGGCACCCCCGGCGTATTTGCGCCaag                                               <  2:3029426/100‑1 (MQ=255)
                                                    ggAAAGCCGCGTGCCAAGCAGGGGTTCGCCGCGGCGGGGCCGCGGC.GCC..CC..GG.GCGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGaaa                                         <  1:5308600/101‑1 (MQ=18)
                                                     gAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.AC.CC..GG.GCGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGCGCCAAGACGAAAc                                        <  1:5635746/101‑1 (MQ=18)
                                                       aaGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.AC.CC..GG.GC..GGG..GGG..GCCCGCCCCCCCCCCCCCCGCGCTTTTGGCCCAAGAAGaaaaa                                       >  2:6763712/1‑99 (MQ=255)
                                                        aGCCGCGTGCGAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.AC.CC..GG.GCGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAAC.AAg                                    <  2:780283/101‑1 (MQ=18)
                                                        aGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.AC.CC..GG.GC..GGG..GGG..GCC.GCCCCCCCCCCCCCCGCGCATTTTGGGCCAGAAGGAAC.aaa                                    >  2:4999290/1‑100 (MQ=255)
                                                           cgcgTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.AC.CC..GG.GCGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAAC.AAGCgg                                 <  2:6597525/101‑1 (MQ=18)
                                                           cgcgTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.AC.CC..GG.GCGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAAC.AAGCgg                                 <  1:3567709/101‑1 (MQ=18)
                                                            gcgTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.AC.CC..GG.GCGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAAC.AAGCggg                                <  2:244602/101‑1 (MQ=18)
                                                             cgTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.AC.CC..GG.GCGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAAC.AAGCGGGa                               <  2:2826463/101‑1 (MQ=18)
                                                              gTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.AC.CC..GG.GCGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAAC.AAGCGGGAc                              <  2:1826402/101‑1 (MQ=18)
                                                               tGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.AC.CC..GG.GCGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAAC.AAGCGGGAcc                             <  2:6977936/101‑1 (MQ=18)
                                                               tGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.AC.CC..GG.GCGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAAC.AAGCGGGAcc                             <  1:1716382/101‑1 (MQ=18)
                                                                 ccAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.AC.CC..GG.GCGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAAC.AAGCGGGACCAg                           <  2:3836160/101‑1 (MQ=18)
                                                                    aGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.AC.CC..GG.GCGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAAC.AAGCGGGACCAGGcc                        <  1:6493003/101‑1 (MQ=18)
                                                                       aGGGTTTCGCCGCGGCGGGGCCGCGGC.G.AC.CC..GG.GCGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAAC.AAGCGGGACCAGGccccc                     <  1:4887241/101‑1 (MQ=18)
                                                                       aGGGTTTCGCCGCGGCGGGGCCGCGGC.G.AC.CC..GG.GCGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAAC.AAGCGGGACCAGGccccc                     <  1:6132958/101‑1 (MQ=18)
                                                                        gggTTTCGCCGCGGGGGGGCCGCGGCGG.C..CC..GG.GCGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAAC.AAGCGGGACCAGGcccccc                    <  1:3196958/101‑1 (MQ=18)
                                                                        gggTTTCGCCGCGGCGGGGCCGCGGC.G.AC.CC..GG.GCGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAAC.AAGCGGGACCAGGcccccc                    <  2:1742346/101‑1 (MQ=18)
                                                                        gggTTTCGCCGCGGCGGGGCCGCGGC.G.AC.CC..GG.GCGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAAC.AAGCGGGACCAGGcccccc                    <  2:5195435/101‑1 (MQ=18)
                                                                        gggTTTCGCCGCGGCGGGGCCGCGGC.G.AC.CC..GG.GCGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAAC.AAGCGGGACCAGGcccccc                    <  1:1230755/101‑1 (MQ=18)
                                                                        gggTTTCGCCGCGGCGGGGCCGCGGC.G.AC.CC..GG.GCGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAAC.AAGCGGGACCAGGcccccc                    <  1:4542202/101‑1 (MQ=18)
                                                                        gggTTTCGCCGCGGCGGGGCCGCGGC.G.A..CCCGGGCGC..GGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAAC.AAGCGGGACCAGGcccccc                    <  2:1577558/101‑1 (MQ=18)
                                                                         ggTTTCGCCGCGGCGGGGCCGCGGC.G.AC.CC..GG.GCGGGGG..CTC..C...GCCCCAGCACCCCCGGCGTATTTGCGCCAAGAAGAAAC.AAGCGGGACCAGGCCCCCCt                   <  1:374638/101‑1 (MQ=21)
                                                                          gTTTCGCCGCGGCGGGGCCGCGGC.G.AC.CC..GG.GCGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAAC.AAGCGGGACCAGGCCCCCCTg                  <  1:534902/101‑1 (MQ=18)
                                                                          gTTTCGCCGCGGCGGGGCCGCGGC.G.AC.CC..GG.GCGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAAC.AAGCGGGACCAGGCCCCCCTg                  <  2:2469620/101‑1 (MQ=18)
                                                                           tttCGCCGCGGCGGGGCCGCGGC.G.AC.CC..GG.GCGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGCGCCAAGCAGAAAC.AAGCGGGACTAGGCCCCCCTGc                 <  1:6279734/101‑1 (MQ=18)
                                                                            ttCGCCGCGGCGGGGCCGCGGC.G.AC.CC..GG.GCGGGGG..CCC..C...GCCCCCCCACCCCCCGCGGATTTTTGCCCACAAGAAACAAAAGGGGGCCAGGccccccccc                 >  2:5771149/1‑98 (MQ=9)
                                                                             tCGCCGCGGCGGGGCCGCGGC.G.AC.CC..GG.GCGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAAC.AAGCGGGACCAGGCCCCCCTGCAg               <  1:4631786/101‑1 (MQ=18)
                                                                              cgccgcGGCGGGGCCGCGGC.G.AC.CC..GG.GCGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAAC.AAGCGGGACCAGGCCCCCCTGCAGc              <  2:5125247/101‑1 (MQ=18)
                                                                                ccgcGGCGGGGCCGCGGC.G.AC.CC..GG.GC..GGG..GGGTCC...CCCCCCCCCCCCCCCGCGGATTTTCCCCAAAAAGAAAC.AAAAGGGGCCCGGccccccccccgccc            >  1:2132679/1‑93 (MQ=255)
                                                                                 ctcggcgtGGCCGCGGC.G.AC.CC..GG.GCAGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAAC.AAGCGGGACCAGGCCCCCCTGCAGCCGc           <  2:4626178/99‑1 (MQ=18)
                                                                                   cggcggGGCCGCGGC.G.AC.CC..GG.GCGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAAC.AAGCGGGACCAGGCCCCCCTGCAGCCGCAt         <  1:4081539/101‑1 (MQ=18)
                                                                                   cggcggGGCCGCGGC.G.AC.CC..GG.GCGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAAC.AAGCGGGACCAGGCCCCCCTGCAGCCGCAt         <  1:2580262/101‑1 (MQ=18)
                                                                                    ggcggGGCCGCGGC.G.AC.CC..GG.GCGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAAC.AAGCGGGACCAGGCCCCCCTGCAGCCGCATg        <  1:1013762/101‑1 (MQ=18)
                                                                                    ggcggGGCCGCGGC.G.AC.CC..GG.GCGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAAC.AAGCGGGACCAGGCCCCCCTGCAGCCGCATg        <  1:6764066/101‑1 (MQ=18)
                                                                                      cggGGCCGCGGC.G.AC.CC..GG.GCGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAAC.AAGCGGGACCAGGCCCCCCTGCAGCCGCATGcc      <  1:3200648/101‑1 (MQ=18)
                                                                                       ggggCCGCGGC.G.AC.CC..GG.GCGGGGG..CTC..C...GCCCCCCCACCCCCCGCGTATTTGCGCCCAGAAGAAAC.AAGCGGGGCCCAGCCCCCCCGCACCCGcccgccc     >  2:2165089/1‑95 (MQ=12)
                                                                                        gggCCGCGGC.G.AC.CC..GG.GCGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAAC.AAGCGGGACCAGGCCCCCCTGCAGCCGCATGCCGt    <  1:2875940/101‑1 (MQ=18)
                                                                                          gCCGCGGC.G.AC.CC..GG.GCGGGGG..CTC..C...GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAAC.AAGCGGGACCAGGCCCCCCTGCAGCCGCATGCCGTcc  <  1:6593760/101‑1 (MQ=18)
                                                                                                      |                                                                                                      
TAGAGTTTATATAGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGCGGC.G.A..CC..GG.GC..GGG..CGG..G...GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAAC.AAGCGGGACCAGGCCCCCCTGCAGCCGCATGCCGTCC  >  NC_014034/1732597‑1732784

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 27 ≤  ATCG  < 33 ≤  ATCG  < 39 ≤  ATCG