Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,494,015 C→G 3.0% intergenic (+63/‑179) RCAP_rcc02319 → / → RCAP_rcc02320 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,494,0150C→G 3.0% 5.3 ‑1.8 ‑0.3166intergenic (+63/‑179)RCAP_rcc02319/RCAP_rcc02320hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (4/1):  ref base (69/91):  total (74/92)
Fisher's exact test strand distribution p-value = 1.71e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 4.49e-01

Only 100 of 184 total aligned reads displayed.

GCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAA.GGGGGGG.CG..CTGCCCCCGTTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.AGCAGAAGGGGGGCAAGCG.CC.CCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGA  >  NC_014034/2493916‑2494114
                                                                                                       |                                                                                                         
gcgcCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGGGG.CG..CTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAAc                                                                                                           <  1:1815448/101‑1 (MQ=21)
 cgccgcGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGGGG.CG..CTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAAcc                                                                                                          >  2:2541273/1‑101 (MQ=21)
  gccgcGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGGGG.CG..CTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCg                                                                                                         <  1:269230/101‑1 (MQ=21)
  gccgcGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGGGG.CG..CTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCg                                                                                                         <  1:4489113/101‑1 (MQ=21)
   ccgcGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGGGG.CG..CTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCgg                                                                                                        <  2:4804522/101‑1 (MQ=21)
   ccgcGGCGGGCGCGCGGCACGACAGGCCGAGTAAGGGGGGGG.CG..CTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCgg                                                                                                        <  1:3951055/101‑1 (MQ=21)
     gcggcgGGCGCGCGGCACGACAGGCCGAGGAA.GGGGGGG.GGGCCGGCCCCCCGCCCCCAACGGGACCCCCCCCCGGCGATTTCCACCAAAATAAAAACGGa                                                                                                       >  1:902978/1‑101 (MQ=255)
      cggcggGCGCGCGGCACGACAGGCCGAGTAAGGGGGGGG.CG..CTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATg                                                                                                     <  1:610261/101‑1 (MQ=21)
        gcgggcgCGCGGCACGACAGGCCGAGTAAGGGGGGGG.CG..CTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCCGCAAGATGAAACCGGATGcc                                                                                                   >  1:2245161/1‑101 (MQ=18)
        gcgggcgCGCGGCACGACAGGCCGAGTAAGGGGGGGG.CG..CTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGcc                                                                                                   <  1:1338150/101‑1 (MQ=21)
         cgggcgCGCGGCACGACAGGCCGAGTAAGGGGGGGG.CG..CTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCt                                                                                                  <  2:3568878/101‑1 (MQ=21)
          gggcgCGCGGCACGACAGGCCGAGTAAGGGGGGGG.CG..CTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCtt                                                                                                 <  1:4471302/101‑1 (MQ=21)
          gggcgCGCGGCACGACAGGCCGAGGGA.GGGGGGGACG..GTGCCCCCCGTCCGCAAGCGGCCTCTCACCCGGTTTTTGGCAGCGAGAGGGAACCTGAttgcct                                                                                                 >  2:162605/1‑96 (MQ=9)
           ggcgCGCGGCACGACAGGCCGAGTAAGGGGGGGG.CG..CTGCCCCCCGTCCGCAAGCGGGCTCGCCCCCGGCTAATTGCAGCAAGATTAAAGCGGGTGGCttt                                                                                                >  2:578350/1‑101 (MQ=12)
           ggcgCGCGGCACGACAGGCCGAGTAAGGGGGGGG.CG..CTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCttt                                                                                                <  2:1945107/101‑1 (MQ=21)
           ggcgCGCGGCACGACAGGCCGAGTAAGGGGGGGG.CG..CTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCttt                                                                                                <  2:3009561/101‑1 (MQ=21)
            gcgcgcGGCACGACAGGCCGAGTAAGGGGGGGG.CG..CTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTc                                                                                               <  1:4292234/101‑1 (MQ=21)
                 cgGCACGACAGGCCGAGTAAGGGGGGGG.CG..CTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGCTATTTGCAGCAAGATGAAACCGGATGCCTTTCATCtt                                                                                          >  2:3262291/1‑101 (MQ=21)
                 cagCACGACAGGCCGAGTAAGGGGGGGG.CG..CTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCtt                                                                                          <  2:5014979/99‑1 (MQ=18)
                   gCACGACAGGCCGAGTAAGGGGGGGG.CG..CTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGc                                                                                        >  2:4969158/1‑101 (MQ=21)
                     ccGACAGGCCGAGTAAGGGGGGGG.CG..CTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGccc                                                                                      <  1:276658/100‑1 (MQ=18)
                     aCGACAGGCCGAGTAAGGGGGGGG.CG..CTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGccc                                                                                      <  1:3164033/101‑1 (MQ=21)
                     aCGACAGGCCGAGTAAGGGGGGGG.CG..CTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGccc                                                                                      <  1:1085005/101‑1 (MQ=21)
                      cGACAGGCCGAGTAAGGGGGGGG.CG..CTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCg                                                                                     <  1:389726/101‑1 (MQ=21)
                         cAGGCCGAGTAAGGGGGGGG.CG..CTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGaaa                                                                                  >  2:2443809/1‑101 (MQ=21)
                          aGGCCGAGTAAGGGGGGGG.CG..CTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAat                                                                                 <  2:3888637/101‑1 (MQ=21)
                          aGGCCGAGTAAGGGGGGGG.CC..CTTCCCCCCCTCCGCAAGGGGGCCCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAat                                                                                 <  2:3696304/101‑1 (MQ=21)
                           ggCCGAGTAAGGGGGGGG.CG..CTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAata                                                                                <  1:739234/101‑1 (MQ=21)
                            gCCGAGTAAGGGGGGGG.CG..CTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAatat                                                                               <  2:1716759/101‑1 (MQ=21)
                             ccGAGTAAGGGGGGGG.CG..CTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATc                                                                              <  1:168796/101‑1 (MQ=21)
                              cGAGTAAGGGGGGGG.CG..CTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATcc                                                                             <  2:2103101/101‑1 (MQ=18)
                              cGAGTAAGGGGGGGG.CG..CTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATcc                                                                             <  2:1492206/101‑1 (MQ=18)
                              cGAGTAAGGGGGGGG.CG..CTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATcc                                                                             >  1:1926281/1‑101 (MQ=18)
                              cGAGTAAGGGGGGGG.CG..CTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATcc                                                                             <  2:4482938/101‑1 (MQ=21)
                                aGTAAGGGGGGGG.CG..CTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTTCCCGAAATATccca                                                                           >  1:2444096/1‑100 (MQ=18)
                                 gtaagggggggg.CG..CTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCgg                                                                          >  2:50858/6‑101 (MQ=255)
                                 gtaagggggggg.CG..CTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCgg                                                                          >  1:551448/6‑101 (MQ=255)
                                  taagggggggg.CG..CTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCCAGATGAAACCGGGTGCCTTTCATCTTGCCCGAAATATTCCggg                                                                         >  2:1893391/5‑101 (MQ=18)
                                  taaggaggggg.CG..CTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCCAGATGAAACCGGGTGCCTTTCATCTTGCCCGAAATATCCCggg                                                                         >  1:333415/7‑101 (MQ=17)
                                   aagggggggg.CG..CTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCgggg                                                                        >  1:3722445/4‑101 (MQ=255)
                                    agggggggg.CG..CTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCCGCAAGATGAAACCCGATGCCTTTCATCTTGCCCGAACTATCCCggggg                                                                       >  2:1216163/3‑101 (MQ=18)
                                    agggggggg.CG..CTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCggggg                                                                       >  1:4942366/3‑101 (MQ=255)
                                     gggggggg.CG..CTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCgggggg                                                                      >  2:2513676/2‑101 (MQ=255)
                                        ggggg.CG..CTGCCCCCCGTGCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATCAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGa                                                                   <  1:5026068/101‑1 (MQ=18)
                                          ggg.CG..CTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.t                                                                <  1:231201/101‑1 (MQ=255)
                                          ggg.CG..CTGCCCCCCGTCCGCAAGCGGACTCCCCCCCGGCTATTTGCAGCAAGATGAAACCGGAAGCATTTCCTCTTTCCCGAAAAATCCCggggggggcat                                                                 >  1:1162576/1‑96 (MQ=17)
                                           gg.CG..CTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.Tg                                                               <  1:2366847/101‑1 (MQ=255)
                                           gg.CG..CTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.Tg                                                               <  2:4312072/101‑1 (MQ=255)
                                            g.cg..cTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGCAATATCCCGGGGGGTGAA.Tgg                                                              <  1:2876840/101‑1 (MQ=255)
                                            g.cg..cTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.Tgg                                                              >  2:4718543/1‑101 (MQ=255)
                                               g..cTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..g                                                          >  1:944133/1‑101 (MQ=255)
                                               g..cTGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..g                                                          >  2:2348482/1‑101 (MQ=255)
                                                   tGCCCCCCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aa                                                       >  1:2942809/1‑100 (MQ=255)
                                                    ccccccCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aag                                                      <  2:295632/100‑3 (MQ=255)
                                                     cccccG.TCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagcc                                                    <  1:3490442/101‑5 (MQ=255)
                                                     cccccG.TCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagcc                                                    <  1:1753572/101‑5 (MQ=255)
                                                     cccccCGTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagc                                                     <  1:779015/101‑4 (MQ=255)
                                                      ccccG.TCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagcca                                                   >  2:1580692/1‑96 (MQ=255)
                                                        cccgtCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAAAATCCCGGGGGGTGAA.TGGC..G.aagccag                                                  >  1:2081755/5‑95 (MQ=21)
                                                           gtCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGCAATATCCCGGGGGGCGAA.TGGC..G.aagccagagg                                               >  2:4054418/2‑92 (MQ=18)
                                                           gtCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.aagccagagg                                               <  2:4234184/100‑10 (MQ=255)
                                                                 aaGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGCa                                          >  2:2286492/1‑101 (MQ=18)
                                                                  aGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGCag                                         >  2:2861193/1‑100 (MQ=18)
                                                                  aGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGCag                                         <  2:1010559/101‑2 (MQ=18)
                                                                   gCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGCagc                                        <  2:2245006/101‑3 (MQ=18)
                                                                    cGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGCaggg                                       >  2:2108998/1‑98 (MQ=18)
                                                                     ggACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGCagcgc                                      >  1:4922635/1‑97 (MQ=17)
                                                                      gACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCG.cc                                   >  1:1658484/1‑101 (MQ=17)
                                                                      gACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCG.cc                                   >  1:1954739/1‑101 (MQ=17)
                                                                       aCTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCG.cc.c                                 >  2:380330/1‑101 (MQ=18)
                                                                         tCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCG.cc.ccc                               >  2:3924517/1‑101 (MQ=17)
                                                                           cccccGGAATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCGGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCG.cc.cccct                             <  2:4240711/101‑3 (MQ=17)
                                                                               cGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCCTCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCC.CCCTGccc                          >  1:1996975/1‑101 (MQ=17)
                                                                               cGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCC.CCCTGccc                          <  2:1778835/101‑1 (MQ=17)
                                                                                 ggATATTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AGGCCAAGAGGGGGCAGCG.CC.CCCTGCCCcg                        >  1:414303/1‑101 (MQ=2)
                                                                                 ggATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCC.CCCTGCCCcg                        <  1:1580824/101‑1 (MQ=17)
                                                                                 ggATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCC.CCCTGCCCcg                        >  1:4557965/1‑101 (MQ=17)
                                                                                  gATATTTGCGGCAAGATGAAGGGGCATG.ATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGGCCAA.AGGCCAAGAGGGGGCAGCG.CC.CCCTGCCCcgc                       >  1:1106325/1‑100 (MQ=2)
                                                                                  gATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGCGAA.TGGC..GAAGCCAGAGGGGGGCAGCGG.CC.CCCTGCCCcgg                       >  2:3780334/1‑101 (MQ=17)
                                                                                   atatTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGGGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCC.CCCTGCCCcgcg                      >  1:2509077/1‑99 (MQ=17)
                                                                                    tatTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCC.CCCTGCCCcggcg                     >  2:4221845/1‑101 (MQ=18)
                                                                                     atTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCC.CCCTGCCCcggcgg                    >  2:715411/1‑101 (MQ=18)
                                                                                      tttGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGt                   <  2:1530092/101‑1 (MQ=17)
                                                                                       ttGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTc                  <  1:1445728/101‑1 (MQ=17)
                                                                                        tGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTcc                 >  1:3887436/1‑101 (MQ=17)
                                                                                        tGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTcc                 >  2:1798383/1‑101 (MQ=17)
                                                                                        tGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGGcc                 >  1:3779656/1‑101 (MQ=17)
                                                                                        tGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGCcc                 >  1:3021628/1‑101 (MQ=17)
                                                                                        tGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGG.GAGCG.CCGCCCTGCCCCGGCGGTcc                 >  2:1902392/1‑101 (MQ=18)
                                                                                        tGCACNNAGATGCAACCGGATTCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTcc                 <  1:1976333/101‑1 (MQ=14)
                                                                                            gcaAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGGGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTCCGTgg             >  1:4272173/1‑101 (MQ=17)
                                                                                             caAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCCGAGGGGGGC.AGCGCCC.CCCTCCCCCAGCCGTCCCTggc            >  2:887342/1‑101 (MQ=9)
                                                                                             caAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTCCGTggc            <  2:256841/101‑1 (MQ=17)
                                                                                               aGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTCCGTggcgg          >  2:3403828/1‑101 (MQ=18)
                                                                                               aGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGGGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCC.CCCTGCCCCCGCGGTCCGTggcgg          >  2:3144664/1‑101 (MQ=17)
                                                                                                    aaaCCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTCCGTGGCGGCGCCt     <  2:5000912/101‑1 (MQ=17)
                                                                                                    aaaCCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTCCGTGGCGGCGCCt     >  2:29184/1‑101 (MQ=17)
                                                                                                    aaaCCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGGGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTCCGTGGCGGCGCCt     >  1:3099186/1‑101 (MQ=18)
                                                                                                      aCCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTCCGTGGCGGCGCCTTg   <  1:4969488/101‑1 (MQ=17)
                                                                                                      aCCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..GAAGCCAGAGGGGGGC.AGCGCCC.CCCTGCCCCGGCGGTCCGTGGCGGCGCCTTg   >  2:3661741/1‑101 (MQ=18)
                                                                                                       |                                                                                                         
GCGCCGCGGCGGGCGCGCGGCACGACAGGCCGAGTAA.GGGGGGG.CG..CTGCCCCCGTTCCGCAAGCGGACTCCCCCCGGGATATTTGCAGCAAGATGAAACCGGATGCCTTTCATCTTGCCCGAAATATCCCGGGGGGTGAA.TGGC..G.AGCAGAAGGGGGGCAAGCG.CC.CCCTGCCCCGGCGGTCCGTGGCGGCGCCTTGA  >  NC_014034/2493916‑2494114

Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 24 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG