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

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,640,1980A→G15.2% 45.1 ‑2.4 ‑0.0131intergenic (‑174/+102)thrS/RCAP_rcc02462threonyl‑tRNA synthetase/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (14/6):  ref base (53/58):  total (67/64)
Fisher's exact test strand distribution p-value = 8.92e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.89e-01

Only 100 of 139 total aligned reads displayed.

GGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGG..CTTCGGCCAATT.CGCCCCCCGAGGTATTTGGAACAAGATGAAAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTTCTCAAGCGTTCC  >  NC_014034/2640099‑2640296
                                                                                                       |                                                                                                  
   aGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGC.CTTCGGCCAATT.CCCCCCCCGAGGTATTTGGAGCAATATGa                                                                                                 <  2:4096019/101‑1 (MQ=25)
    ggCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGC.CTTCGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATGaa                                                                                                <  1:1203772/101‑1 (MQ=25)
     gCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGC.CTTCGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATGaaa                                                                                               <  2:3662876/101‑1 (MQ=25)
         gtCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGC.CTTCGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATGaaagacc                                                                                           <  2:4056931/101‑5 (MQ=21)
            aaCCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGC.CTTCGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATg                                                                                       <  1:3137833/101‑1 (MQ=25)
            aaCCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGC.CTTCGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATg                                                                                       <  2:2487968/101‑1 (MQ=21)
              ccGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGC.CTTCGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATggg                                                                                     <  1:1451997/101‑1 (MQ=25)
                   gcgcTTCTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGC.CTTCGGCCAATT.CCCCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCAt                                                                                <  1:685582/101‑1 (MQ=21)
                     gcTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGC.CTTCGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCt                                                                              <  1:2672922/101‑1 (MQ=25)
                      cTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGC.CTTCGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCtt                                                                             <  1:3089224/101‑1 (MQ=21)
                       ttttGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGC.CTTCGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTg                                                                            <  2:3042280/101‑1 (MQ=25)
                           gcgcAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGC.CTTCGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTc                                                                        <  1:251102/101‑1 (MQ=25)
                            cgcAGGCAAAAGGGGGGGGCGGCCCCCCCCTTTGGC.CTTCGGCCAATTACCCCCCCCGGGGT.TTTGG.GGAAGAAGAAAAAACAAGGGGTT.ATCTTTCCTCa                                                                       >  2:1067304/1‑101 (MQ=255)
                            cgcAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGC.CTTCGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCa                                                                       <  1:925952/101‑1 (MQ=25)
                              gccccgaaaAGGGGGCGCTGCCCCCCTC.TTGGC.CTTCGGCCAATT.CACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAat                                                                    >  1:1398393/7‑101 (MQ=14)
                              cAGGCAAAAGGGGGGGGCCGCCCCCCCC.TTTGCCCTTGGGCCAATT.ACCCCCCCG.GGGATTTTGGGGAAAATGAAAA.ACCATGGGTTTATCTTTCCTCaaa                                                                     >  2:2451290/1‑101 (MQ=255)
                               ccccgaaaAGGGGGCGCTGCCCCCCTC.TTGGC.CTTCGGCCAATT.CACCCCCCGGGATATTTGCGGCAAGGCGAAAGGA..TGGGGTTTATCTTGCCCCAAata                                                                   >  2:723395/6‑101 (MQ=9)
                                ggCAAAAGGGGGGCGCTGCCCCCCTC.TTGGC.CTTCGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAata                                                                   <  1:210318/101‑1 (MQ=21)
                                cccgaaaAGGGGGCGCTGCCCCCCTC.TTGGC.CTTCGGCCAATT.CACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAatat                                                                  <  1:536724/97‑1 (MQ=17)
                                cccgaaaAGGGGGCGCTGCCCCCCTC.TTGGC.CTTCGGCCAATT.CACCCCCCGGGACATTTGCGGCAAGGTGAAAGGA..TGGGGCTTATCTTGCCCCAAatat                                                                  >  2:1809181/5‑101 (MQ=9)
                                 gCAAAAGGGGGGCGCTGCCCCCCTC.TTGGC.CTTCGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAatat                                                                  <  2:1115952/101‑1 (MQ=25)
                                  cAAAAGGGGGGCGCTGCCCCCCTC.TTGGC.CTTCGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATc                                                                 >  2:3948178/1‑101 (MQ=21)
                                   gaaaAGGGGGCGCTGCCCCCCTC.TTGGC.CTTCGGCCAATT.CACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATccc                                                               <  2:2171381/100‑1 (MQ=17)
                                    aaaGGGGGGAGCTGCCCCCCTC.TTGGC.CTTTGGCCAATT.CACCCCCCGCGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATacc                                                               >  2:551273/1‑101 (MQ=18)
                                    aaaAGGGGGCGCTGCCCCCCTC.TTGGC.CTTCGGCCAATT.CACCCCCCGGGATATTTGCGGCAAGGGGAAAGGA.TGGGGGTT.ATCTTGCCCCAAATATCCCg                                                              >  1:131524/1‑101 (MQ=17)
                                      aagggggCGCTGCCCCCCTC.TTGGC.CTTCGGCCAATT.CACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCggg                                                            >  2:1863998/3‑101 (MQ=18)
                                       ggggggCGCTGCCCCCCTC.TTGGC.CTTCGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCggg                                                            <  2:211406/101‑1 (MQ=25)
                                       ggggggCCCTGCCCCCCTC.TTGGC.CTTCGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCggg                                                            <  2:2338473/101‑1 (MQ=25)
                                       agggggCGCTGCCCCCCTC.TTGGC.CTTCGGCCAATT.CACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCgggg                                                           >  2:4098878/2‑101 (MQ=18)
                                       agggggCGCTGCCCCCCTC.TTGGC.CTTCGGCCAATT.CACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCgggg                                                           >  1:529732/2‑101 (MQ=18)
                                        gggggCGCTGCCCCCCTC.TTGGC.CTTCGGCCAATT.CACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCggggg                                                          >  2:1643891/1‑101 (MQ=18)
                                        gggggCGCTGCCCCCCTC.TTGGC.CTTCGGCCAATT.CACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCggggg                                                          <  1:1490286/101‑1 (MQ=18)
                                        gggggCGCTGCCCCCCTC.TTGGC.CTTCGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCgggg                                                           <  2:1443836/101‑1 (MQ=25)
                                         ggggCGCTGCCCCCCTC.TTGGC.CTTCGGCCAATT.CACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCgggggt                                                         >  2:3483856/1‑100 (MQ=18)
                                         ggggCGCTGCCCCCCTC.TTGGC.CTTCGGCCAATT.CACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCgggggg                                                         >  2:4027452/1‑101 (MQ=18)
                                         ggggCGCTGCCCCCCTC.TTGGC.CTTCGGCCAATT.CACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCgggggg                                                         >  2:2878552/1‑101 (MQ=18)
                                          gggCGCTGCCCCCCTC.TTGGC.CTTCGGCCAATT.CACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA.TGGGGGTT.ATCTTGCCCCAAATATCCCgggggtg                                                        >  1:4120356/1‑99 (MQ=18)
                                          gggCGCTGCCCCCCTC.TTGGC.CTTCGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCgggggg                                                         <  1:3920565/101‑1 (MQ=25)
                                           ggCGCTGCCCCCCTC.TTGGC.CTTCGGCCAATT.CACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCgggggtga                                                       <  1:439035/101‑4 (MQ=21)
                                           ggCGCTGCCCCCCTC.TTGGC.CTTCGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGa                                                        <  2:2075963/101‑1 (MQ=25)
                                            gcgcTGCCCCCCTC.TTGGC.CTTCGGCCAATT.CACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCgggggtgag                                                      <  2:3270389/101‑5 (MQ=21)
                                            gcgcTGCCCCCCTC.TTGGC.CTTCGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAg                                                       <  2:3909766/101‑1 (MQ=25)
                                              gcTGCCCCCCTC.TTGGC.CTTCGGCCAATT.CACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCgggggtgaggc                                                    >  2:3670734/1‑95 (MQ=21)
                                              gcTGCCCCCCTC.TTGGC.CTTCGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGNGNGtn                                                     >  1:1363739/1‑100 (MQ=25)
                                               cTGCCCCCCTC.TTGGC.CTTCGGCCAATT.CACCCCCCGGGAAATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCgggggtgaggcg                                                   >  1:960724/1‑94 (MQ=25)
                                                 gCCCCCCTC..TGG..CTTC..CCCATT.CCCCCCCCGGGGTATTTCGGGCAAGATGAAAGGCATCCGGTTTCATCTTGCTGCAAATATCCCGGGGGGAGTccgcttg                                               <  1:3647441/101‑6 (MQ=17)
                                                    ccccTC.TTGGC.CTTCGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGaaa                                               <  1:3978114/101‑1 (MQ=37)
                                                     tcttc.tTGGC.CTTCGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAg                                              <  1:337064/98‑1 (MQ=25)
                                                     cccTC.TTGGC.CTTCGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAg                                              >  1:1208728/1‑101 (MQ=37)
                                                     cccTC.TTGGC.CTTCGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAg                                              >  1:1352283/1‑101 (MQ=37)
                                                      ccTC.TTGGC.CTTCGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAgg                                             >  2:2793688/1‑101 (MQ=37)
                                                       ctc.tTGGC.CTTCGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGa                                            >  1:3750789/1‑101 (MQ=37)
                                                       ctc.tTGGC.CTTCGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGa                                            >  1:336164/1‑101 (MQ=37)
                                                        tc.tTGGC.CTTCGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGAc                                           >  1:2479405/1‑101 (MQ=37)
                                                        tc.tTGGC.CTTCGGCCAATT.CAACCCCCGAGGAATTTGGAGCAAGATG.AAAGACCATGGGTTAATCTTGAATCAAATATCCCGGGGGGAGTCCTGAAAgaga                                           >  2:337013/1‑98 (MQ=25)
                                                            ttggccTTCGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACgg                                         <  2:1208634/96‑1 (MQ=37)
                                                             tggccTTCGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACggg                                        >  1:3496281/5‑101 (MQ=37)
                                                             tggccTTCGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACagg                                        >  2:3866125/5‑101 (MQ=25)
                                                              ggccTTCGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACgggg                                       <  1:1280080/98‑1 (MQ=37)
                                                              ggccTTCGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACgggg                                       >  1:2898838/4‑101 (MQ=37)
                                                                 cTTCGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGc                                    >  1:824355/1‑101 (MQ=37)
                                                                 cTTCGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGc                                    >  2:1764107/1‑101 (MQ=37)
                                                                  ttCGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCa                                   >  2:4017796/1‑101 (MQ=37)
                                                                    cGGCCAATT.CCCCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCaga                                 <  1:1922506/101‑1 (MQ=255)
                                                                    cGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCaga                                 >  1:907509/1‑101 (MQ=37)
                                                                    cGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCaga                                 >  2:3534400/1‑101 (MQ=255)
                                                                    cGGCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCaga                                 >  1:3003586/1‑101 (MQ=255)
                                                                     ggCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCagag                                <  1:3158304/101‑1 (MQ=255)
                                                                      gCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGc                               >  1:1952819/1‑101 (MQ=37)
                                                                      gCCAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGc                               >  2:490034/1‑101 (MQ=37)
                                                                       ccAATT.CGCCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGcc                              <  2:3002284/101‑1 (MQ=255)
                                                                        cAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGccc                             >  2:1709144/1‑101 (MQ=255)
                                                                        cAATT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGccc                             >  1:3660143/1‑101 (MQ=37)
                                                                         aaTT.CACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGcccc                            <  2:67344/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.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGAGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTc                        >  1:1174196/3‑101 (MQ=37)
                                                                               accccccGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTcc                       >  1:4120075/2‑101 (MQ=255)
                                                                                ccccccGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTccc                      >  1:2928867/1‑101 (MQ=255)
                                                                                  ccccGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCgg                    <  1:1709259/101‑1 (MQ=255)
                                                                                  ccccGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCgg                    <  1:3093020/101‑1 (MQ=255)
                                                                                   cccGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCggg                   <  2:225005/101‑1 (MQ=255)
                                                                                     cGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTc                 <  1:3948438/101‑1 (MQ=255)
                                                                                      gAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCt                >  1:2761305/1‑101 (MQ=255)
                                                                                      gAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCt                >  2:2863901/1‑101 (MQ=255)
                                                                                      gAGGTAGTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCt                <  1:2254934/101‑1 (MQ=255)
                                                                                       aGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCtt               <  1:2793907/101‑1 (MQ=255)
                                                                                        ggTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCttt              >  2:1466929/1‑101 (MQ=255)
                                                                                         gTATTTGGAGCAAGATG.AAAGACCATGCGTTCATCTGGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTtc             <  1:1785083/101‑1 (MQ=25)
                                                                                           aTTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTtctc           >  2:4159715/1‑101 (MQ=255)
                                                                                            tttGGAGCACGATG.AAAGACCATGGGATCATCTTGCCTCAAATAGCCCGGGGGAAGTCCTGAAAGGACGCGGTGCTGACCCCCCCACCCCACTCTTTTTCa          >  2:752225/1‑101 (MQ=16)
                                                                                             ttGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTTCTCaa         >  1:3191625/1‑101 (MQ=255)
                                                                                               ggAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTTCTCAAGc       >  1:1490583/1‑101 (MQ=255)
                                                                                                gagcAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTTCTCAAGCg      >  2:1188532/4‑101 (MQ=255)
                                                                                                gagcAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTTCTCAAGCg      <  2:1952693/98‑1 (MQ=255)
                                                                                                   cAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGTAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTTCTCAAGCGTTc   >  1:3591916/1‑101 (MQ=255)
                                                                                                   cAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTTCTCAAGCGTTc   >  2:1060009/1‑101 (MQ=255)
                                                                                                   cAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTTCTCAAGCGTTc   <  1:161938/101‑1 (MQ=255)
                                                                                                       |                                                                                                  
GGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGG..CTTCGGCCAATT.CGCCCCCCGAGGTATTTGGAACAAGATGAAAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTTCTCAAGCGTTCC  >  NC_014034/2640099‑2640296

Base quality scores:  ATCG  < 3 ≤  ATCG  < 16 ≤  ATCG  < 27 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG