Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,640,182 A→G 15.9% intergenic (‑158/+118) thrS ← / ← RCAP_rcc02462 threonyl‑tRNA synthetase/hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,640,1820A→G15.9% 44.8 ‑3.6 ‑0.0132intergenic (‑158/+118)thrS/RCAP_rcc02462threonyl‑tRNA synthetase/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (14/7):  ref base (42/68):  total (57/75)
Fisher's exact test strand distribution p-value = 2.84e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.00e-01

Only 100 of 143 total aligned reads displayed.

GGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGG.CTT.CGGCCAATT.CGCCCCCCGAGGTATTTGGAACAAGATGAAAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTT  >  NC_014034/2640083‑2640283
                                                                                                       |                                                                                                     
gggTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGCCTT.CGGCCAATT.CACCCCCCGa                                                                                                       <  2:2484439/101‑1 (MQ=255)
 ggTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGCCTT.CGGCCAATT.CACCCCCCGAg                                                                                                      <  1:667616/101‑1 (MQ=255)
     tGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGCCTT.CGGCCCATT.CACCCCCCGAGGGAt                                                                                                  >  2:1279984/1‑101 (MQ=37)
     tGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGCCTT.CGGCCAATT.CACCCCCCGAGGTAt                                                                                                  <  2:3161143/101‑1 (MQ=255)
     tGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGCCTT.CGGCCAATT.CACCCCCCGAGGTAt                                                                                                  <  2:1988366/101‑1 (MQ=255)
      tcGACGAGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGCCTT.CGGCCAATT.CCCCCCCCGAGGTAtt                                                                                                 <  2:585320/100‑1 (MQ=25)
      gCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGCCTT.CGGCCAATT.CACCCCCCGAGGTAtt                                                                                                 <  2:3669146/101‑1 (MQ=255)
       cGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGCCTT.CGGCCAATT.CACCCCCCGAGGTAttt                                                                                                <  2:2418875/101‑1 (MQ=255)
       cGACGCGGCGGCAGGCGTGTCAACCGGGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGCCTT.CGGCCAATT.CACCCCCCGAGGTAttt                                                                                                <  1:820020/101‑1 (MQ=25)
        gACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGCCTT.CGGCCAATT.CACCCCCCGAGGTATTTg                                                                                               <  1:663839/101‑1 (MQ=255)
        gACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGCCTT.CGGCCAATT.CACCCCCCGAGGTATTTg                                                                                               <  2:2452333/101‑1 (MQ=255)
         aCGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGCCTT.CGGCCAATT.CACCCCCCGAGGTATTTgg                                                                                              <  2:2249054/101‑1 (MQ=255)
         aCGCGGCGGCAGGCGTGTCAACCGTGCCGCGTTGGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGCCTT.CGGCCAATT.CCCCCCCCGAGGTATTTgg                                                                                              <  1:3925569/101‑1 (MQ=25)
           gcggcgGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGCCTT.CGGCCAATT.CGCCCCCCGAGGTATTTGGag                                                                                            <  1:1966406/101‑2 (MQ=255)
            cggcggCAGGCGTGTCAACCGGGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGCCTT.CGGCCAATT.CACCCCCCGAGGTATTTGGagc                                                                                           <  1:758480/101‑3 (MQ=25)
             ggcggcAGGCGTGTCAACCGTGGCGCTTTTGCGAAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGCCTT.CGGCCAATT.CACCCCCCGAGGTATTTGGAGCa                                                                                          <  1:3662671/101‑1 (MQ=25)
                ggcAGGCGTGTCAACCGTGGCGCTTTGGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGCCTT.CGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGa                                                                                       <  2:694869/101‑1 (MQ=25)
                 gcAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCCGCCCCCCCC.CTTC.CTTCCTGCCAATC.CTCCCCCCCCGCGATTTGGAGCaaaaa                                                                                      >  1:1838174/1‑98 (MQ=255)
                  cAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGCCTT.CGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATg                                                                                     <  2:2116866/101‑1 (MQ=37)
                   aGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGCCTT.CGGCCAATT.CCCCCCCCGAGGTATTTGGAGCAATATGa                                                                                    <  2:4096019/101‑1 (MQ=25)
                    ggggttccacccGTGGCACTTTTGCCAAGGC.AAAGGGGGGCGCTGCCCCCCTC.TTGGCCTT.CGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATGaaa                                                                                  <  2:3804816/91‑1 (MQ=17)
                    ggCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGCCTT.CGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATGaa                                                                                   <  1:1203772/101‑1 (MQ=25)
                         gtCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGCCTT.CGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATGaaagacc                                                                              <  2:4056931/101‑5 (MQ=21)
                            aaCCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGCCTT.CGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATg                                                                          <  2:2487968/101‑1 (MQ=21)
                             aCCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGCCTT.CGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATgg                                                                         <  1:3180511/101‑1 (MQ=25)
                              ccGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGCCTT.CTGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATggg                                                                        <  2:3662767/101‑1 (MQ=25)
                              ccGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGCCTT.CGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATggg                                                                        <  1:1451997/101‑1 (MQ=25)
                                     gcTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGCCTT.CGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCt                                                                 <  1:2672922/101‑1 (MQ=25)
                                       ttttGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGCCTT.CGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTg                                                               <  2:3042280/101‑1 (MQ=25)
                                          tGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGCCTT.CGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCt                                                            <  1:3837973/101‑1 (MQ=21)
                                            cgcAGGCAAAAGGGGGGGGCGGCCCCCCCCTTTGGCCTT.CGGCCAATTACCCCCCCCGGGGT.TTTGG.GGAAGAAGAAAAAACAAGGGGTT.ATCTTTCCTCa                                                          >  2:1067304/1‑101 (MQ=255)
                                            cgcAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGCCTT.CGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCa                                                          <  1:925952/101‑1 (MQ=25)
                                              gccccgaaaAGGGGGCGCTGCCCCCCTC.TTGGCCTT.CGGCCAATT.CACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAat                                                       >  1:1398393/7‑101 (MQ=14)
                                               ccccgaaaAGGGGGCGCTGCCCCCCTC.TTGGCCTT.CGGCCAATT.CACCCCCCGGGATATTTGCGGCAAGGCGAAAGGA..TGGGGTTTATCTTGCCCCAAata                                                      >  2:723395/6‑101 (MQ=9)
                                                ggCAAAAGGGGGGCGCTGCCCCCCTC.TTGGCCTT.CGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAata                                                      <  1:210318/101‑1 (MQ=21)
                                                 gCAAAAGGGGGGCGCTGCCCCCCTC.TTGGCCTT.CGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAatat                                                     <  2:3949620/101‑1 (MQ=25)
                                                  cAAAAGGGGGGCGCTGCCCCCCTC.TTGGCCTT.CGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATc                                                    <  1:4044003/101‑1 (MQ=25)
                                                   gaaaAGGGGGCGCTGCCCCCCTC.TTGGCCTT.CGGCCAATT.CACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATccc                                                  >  1:2592976/2‑101 (MQ=17)
                                                   gaaaAGGGGGCGCTGCCCCCCTC.TTGGCCTT.CGGCCAATT.CACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATccc                                                  <  2:1086767/100‑1 (MQ=17)
                                                    aaaGGGGGGAGCTGCCCCCCTC.TTGGCCTT.TGGCCAATT.CACCCCCCGCGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATacc                                                  >  2:551273/1‑101 (MQ=18)
                                                    aaaAGGGGGCGCTGCCCCCCTC.TTGGCCTT.CGGCCAATT.CACCCCCCGGGATATTTGCGGCAAGGGGAAAGGA.TGGGGGTT.ATCTTGCCCCAAATATCCCg                                                 >  1:131524/1‑101 (MQ=17)
                                                      aagggggCGCTGCCCCCCTC.TTGGCCTT.CGGCCAATT.CACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCggg                                               >  2:1863998/3‑101 (MQ=18)
                                                       ggggggCGCTGCCCCCCTC.TTGGCCTT.CGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCggg                                               <  2:211406/101‑1 (MQ=25)
                                                       ggggggCGCTGCCCCCCTC.TTGGCCTT.CGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCggg                                               >  1:1356940/1‑101 (MQ=21)
                                                       agggggCGCTGCCCCCCTC.TTGGCCTT.CGGCCAATT.CACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCgggg                                              >  1:529732/2‑101 (MQ=18)
                                                        gggggCGCTGCCCCCCTC.TTGGCCTT.CGGCCAATT.CACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCggggg                                             >  2:1643891/1‑101 (MQ=18)
                                                        gggggCGCTGCCCCCCTC.TTGGCCTT.CGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCgggg                                              <  2:1443836/101‑1 (MQ=25)
                                                         ggggCGCTGCCCCCCTC.TTGGCCTT.CGGCCAATT.CACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCgggggt                                            >  2:3483856/1‑100 (MQ=18)
                                                           ggCGCTGCCCCCCTC.TTGGCCTT.CGGCCAATT.CACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCgggggtga                                          <  1:439035/101‑4 (MQ=21)
                                                           ggCGCTGCCCCCCTC.TTGGCCTT.CGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGa                                           <  2:2075963/101‑1 (MQ=25)
                                                            gcgcTGCCCCCCTC.TTGGCCTT.CGGCCAATT.CACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCgggggtgag                                         <  2:3270389/101‑5 (MQ=21)
                                                            gcgcTGCCCCCCTC.TTGGCCTT.CGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAg                                          <  2:3909766/101‑1 (MQ=25)
                                                             cgcTGCCCCCCTC.TTGGCCTT.CGGCCCATT.CACCCCCCCAGGTATTTTGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCCCAAATATCCCcgggggggg                                         >  1:3457657/1‑98 (MQ=21)
                                                             cgcTGCCCCCCTC.TTGGCCTT.CGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGt                                         <  2:3468420/101‑1 (MQ=25)
                                                              gcTGCCCCCCTC.TTGGCCTT.CGGCCAATT.CACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCgggggtgaggc                                       >  2:3670734/1‑95 (MQ=21)
                                                               cTGCCCCCCTC.TTGGCCTT.CGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTcc                                       >  2:236615/1‑101 (MQ=25)
                                                                 gCCCCCCTC..TGG.CTT.C..CCCATT.CCCCCCCCGGGGTATTTCGGGCAAGATGAAAGGCATCCGGTTTCATCTTGCTGCAAATATCCCGGGGGGAGTccgcttg                                  <  1:3647441/101‑6 (MQ=17)
                                                                   cccccTC.TTGGCCTT.CGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGaa                                   >  1:313449/1‑101 (MQ=37)
                                                                   cccccTC.TTGGCCTT.CGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGaa                                   <  1:14955/101‑1 (MQ=25)
                                                                     cccTC.TTGGCCTT.CGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAg                                 >  1:1208728/1‑101 (MQ=37)
                                                                     cccTC.TTGGCCTT.CGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAg                                 >  1:1352283/1‑101 (MQ=37)
                                                                      ccTC.TTGGCCTT.CGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAgg                                >  2:2793688/1‑101 (MQ=37)
                                                                       ctc.tTGGCCTT.CGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGa                               >  1:3750789/1‑101 (MQ=37)
                                                                       ctc.tTGGCCTT.CGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGa                               >  1:336164/1‑101 (MQ=37)
                                                                        tc.tTGGCCTT.CGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGAc                              >  1:2479405/1‑101 (MQ=37)
                                                                           ttggccTT.CGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACgg                            >  1:144686/6‑101 (MQ=37)
                                                                           ttggccTT.CGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACgg                            <  2:1208634/96‑1 (MQ=37)
                                                                            tggccTT.CGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACggg                           >  1:3496281/5‑101 (MQ=37)
                                                                            tggccTT.CGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACagg                           >  2:3866125/5‑101 (MQ=25)
                                                                            tggccTT.CGGCCAATT.CACCCCCCGAGGTATTTGGAGAAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACggg                           >  1:4197583/5‑101 (MQ=37)
                                                                             ggccTT.CGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACgggg                          <  1:1280080/98‑1 (MQ=37)
                                                                             ggccTT.CGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACgggg                          >  1:2898838/4‑101 (MQ=37)
                                                                              gtcTT.CGTCCAATT.CCCCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACggggg                         <  1:2980666/99‑1 (MQ=37)
                                                                               ccTT.CGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACgggggg                        >  2:468205/2‑101 (MQ=37)
                                                                                cTT.CGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGc                       >  2:1764107/1‑101 (MQ=37)
                                                                                 tt.CGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCa                      >  2:4017796/1‑101 (MQ=37)
                                                                                 ctgcGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCGTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGc                       <  1:886318/98‑1 (MQ=37)
                                                                                    cGGCCAATT.CCCCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCaga                    <  1:1922506/101‑1 (MQ=255)
                                                                                    cGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCaga                    >  1:3003586/1‑101 (MQ=255)
                                                                                    cGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCaga                    >  1:907509/1‑101 (MQ=37)
                                                                                     ggCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCagag                   >  1:2892209/1‑101 (MQ=37)
                                                                                     ggCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCagag                   <  1:3158304/101‑1 (MQ=255)
                                                                                       ccAATT.CGCCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGcc                 <  2:3002284/101‑1 (MQ=255)
                                                                                         aaTT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGcccc               >  2:3597345/1‑101 (MQ=37)
                                                                                          aTT.CGCCCCCCGAGGTATTTGGAGCAAGATG.AAAGCCCATGGGGTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGccccc              <  2:463192/101‑1 (MQ=25)
                                                                                          aTT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGccccc              >  1:2473321/1‑101 (MQ=37)
                                                                                             tcaccccccGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCt            <  2:3330312/98‑1 (MQ=255)
                                                                                              caccccccGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTc           >  2:4042349/3‑101 (MQ=255)
                                                                                              caccccccGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTc           >  2:1520481/3‑101 (MQ=255)
                                                                                              caccccccGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGAGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTc           >  1:1174196/3‑101 (MQ=37)
                                                                                               accccccGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTcc          >  1:4120075/2‑101 (MQ=255)
                                                                                                 cccccGAGGTATTTGGAGCAAGATG.AAAGACCATGGGCTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCg        <  1:3620756/101‑1 (MQ=255)
                                                                                                  ccccGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCgg       <  1:3093020/101‑1 (MQ=255)
                                                                                                  ccccGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCgg       <  1:1709259/101‑1 (MQ=255)
                                                                                                   cccGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCggg      <  2:225005/101‑1 (MQ=255)
                                                                                                     cGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTc    <  1:3948438/101‑1 (MQ=255)
                                                                                                      gAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCt   >  2:2863901/1‑101 (MQ=255)
                                                                                                      gAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCt   >  2:1804722/1‑101 (MQ=255)
                                                                                                       aGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCtt  <  1:2793907/101‑1 (MQ=255)
                                                                                                       aGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCtt  >  1:2040531/1‑101 (MQ=255)
                                                                                                       |                                                                                                     
GGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGG.CTT.CGGCCAATT.CGCCCCCCGAGGTATTTGGAACAAGATGAAAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTT  >  NC_014034/2640083‑2640283

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