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

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,640,1840G→A19.5% 74.7 0.0 0.0153intergenic (‑160/+116)thrS/RCAP_rcc02462threonyl‑tRNA synthetase/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (15/15):  ref base (59/64):  total (74/79)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 168 total aligned reads displayed.

GTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TT.GG.CTTCGGCCAATTCGCCCCCCGAGGTATTTGGAACAAGATGAAAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTTC  >  NC_014034/2640085‑2640285
                                                                                                      |                                                                                                     
   tGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TT.GGCCTTCGGCCAATTCACCCCCCGAGGTAt                                                                                                    <  2:3945141/101‑1 (MQ=255)
   tGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TT.GGCCTTCGGCCAATTCACCCCCCGAGGTAt                                                                                                    <  2:3665347/101‑1 (MQ=255)
     cGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.CT.GGCCCTCTGCCAATTCCCCCCCCGCGGGAttt                                                                                                  >  2:2264664/1‑101 (MQ=25)
      gACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TT.GGCCTTCGGCCAATTCCCCCCCCGAGGTCTTTg                                                                                                 >  2:1754583/1‑101 (MQ=37)
            gcggcAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TT.GGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCaa                                                                                           <  2:5268482/101‑1 (MQ=37)
            gcggcAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TT.GGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCaa                                                                                           <  2:2526534/101‑1 (MQ=37)
              ggcAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TT.GGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGa                                                                                         <  1:1830338/101‑1 (MQ=37)
               gcAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TT.GGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGAt                                                                                        <  1:285942/101‑1 (MQ=37)
               gcAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TT.GGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGAt                                                                                        <  2:1535203/101‑1 (MQ=37)
                cAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TT.GGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATg                                                                                       <  1:2467963/101‑1 (MQ=25)
                 aGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCACCTC.TT.GGCCTTCGGCCAATTCACCCCCCCCGGTATTTGGAGCAAGATGa                                                                                      >  2:1223828/1‑101 (MQ=25)
                  ggCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TT.GGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATGaa                                                                                     <  2:2036113/101‑1 (MQ=25)
                     gtgtCAGCCGGGGCCCTTTTGCGCAGGGAAAAGGGGGGCGCTGCCCCCCTC.CT.GGCCTTCGGGCAATTCCCCCCCCGAGGTATTTGGAGCAAGATGaaaga                                                                                  <  1:4690161/101‑3 (MQ=18)
                     gtgtCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TT.GGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATGaaaga                                                                                  <  1:154107/101‑3 (MQ=25)
                      tgtCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TT.GGCCTTCGGCCAATTCACCCCCCGAGGGATTTGGAGCAAGATGaaagac                                                                                 >  2:3232790/1‑98 (MQ=25)
                         cAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTACCCCCCTC.TT.GGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCAt                                                                             <  2:1437485/101‑1 (MQ=25)
                          aaCCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TT.GGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATg                                                                            <  1:4438297/101‑1 (MQ=21)
                          aaCCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TT.GGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATg                                                                            <  1:4009758/101‑1 (MQ=25)
                            ccGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TT.GGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATggg                                                                          <  2:3086863/101‑1 (MQ=25)
                                   gcTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TT.GGCCTTCGGCCAATTCCCCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCt                                                                   <  1:1050008/101‑1 (MQ=25)
                                   gcTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TT.GGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCt                                                                   <  1:1376150/101‑1 (MQ=21)
                                        tGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TT.GGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCt                                                              <  2:2800121/101‑1 (MQ=25)
                                         gcgcAGGCAAAAGGGGGGCGCTGCCCCCCTC.TT.GGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTc                                                             <  2:4815659/101‑1 (MQ=25)
                                             aGGCAAAAGGGGGGCGCTGCCCCCCTC.TT.GGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAat                                                         <  2:4445298/101‑1 (MQ=21)
                                              ggCAAAAGGGGGGCGCTGCCCCCCTCTTT.GC.CTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.GAAGACCATGGGTTCATTTTGCCTCAAata                                                        >  2:4396765/1‑101 (MQ=21)
                                              cccgaaaAGGGGGCGCTGCCCCCCTC.TT.GGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAatat                                                       >  1:3613185/5‑101 (MQ=17)
                                               gCAAAAGGGGCGCGCCTCCCCCCTC.TT.GGCCTTCGGCCAATTCGCCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAatat                                                       <  2:2183357/101‑1 (MQ=25)
                                                 gaaaAGGGGGCGCTGCCCCCCTC.TT.GGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATccc                                                    <  1:3834882/100‑1 (MQ=17)
                                                 gaaaAGGGGGCGCTGCCCCCCTC.TT.GGCCTTCGGCCAATTCACCCCCCCGGATATTTGCGGCAAGGGGAAAGGA.TGGGGGTT.ATCTTGCCCCCAATATccc                                                    >  1:2075324/2‑101 (MQ=9)
                                                  aaaAGGGGGCGCTGCCCCCCTC.TT.GGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCg                                                   <  1:4119471/101‑1 (MQ=17)
                                                  aaaAGGGGGCGCTGCCCCCCTC.TT.GGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCg                                                   >  2:4791669/1‑101 (MQ=17)
                                                   aaagggggCGCTGCCCCCCTC.TT.GGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCgg                                                  >  1:1855996/4‑101 (MQ=17)
                                                   aaGGGGGGCGCTGCCCCCCTC.TT.GGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCg                                                   <  2:4554574/101‑1 (MQ=25)
                                                    aagggggCGCTGCCCCCCTC.TT.GGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCggg                                                 >  1:1167496/3‑101 (MQ=18)
                                                    aagggggCGCTGCCCCCCTC.TT.GGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCggg                                                 >  1:1867500/3‑101 (MQ=18)
                                                    aagggggCGCTGCCCCCCTC.TT.GGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCggg                                                 <  1:2114397/99‑1 (MQ=18)
                                                    aagggggCGCTGCCCCCCTC.TT.GGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGACCCAAATATCCCggg                                                 >  1:3353817/3‑101 (MQ=17)
                                                     agggggCGCTGCCCCCCTC.TT.GGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCgggg                                                <  2:5208637/100‑1 (MQ=18)
                                                      gggggCGCTGCCCCCCTC.TT.GGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGTAAGGA..TGGGGTTTATCTTGCCCCAAATATTCCggggg                                               >  2:5542638/1‑101 (MQ=17)
                                                      gggggCGCTGCCCCCCTC.TT.GGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCggggg                                               >  2:3177246/1‑101 (MQ=18)
                                                      gggggCGCTGCCCCCCTC.TT.GGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCggggg                                               <  1:2358007/101‑1 (MQ=18)
                                                        gggCGCTGCCCCCCTC.TT.GGCCTTCGGCCAATTCTCCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCcggggtg                                             <  2:164463/101‑3 (MQ=18)
                                                        gggCGCTGCCCCCCTC.TT.GGCCTTCGGCCAATTCCCCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCgggggg                                              <  1:1754682/101‑1 (MQ=25)
                                                        gggCGCTGCCCCCCTC.TT.GGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCgggggtg                                             <  1:3393296/101‑3 (MQ=18)
                                                        gggCCCCGCCCCCCTC.TT.GGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCgggggg                                              <  1:3942700/101‑1 (MQ=21)
                                                         ggCGCTGCCCCCCTC.TTGGG.CTTCGGCCAAGTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGa                                             <  1:4327505/101‑1 (MQ=25)
                                                         ggCGCTGCCCCCCTC.TT.GGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCgggggtga                                            >  2:1178424/1‑98 (MQ=18)
                                                          gcgcTGCCCCCCTC.TT.GGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCgggggtgag                                           <  2:3114982/101‑5 (MQ=21)
                                                           cgcTGCCCCCCTC.TT.GGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCgggggtgagg                                          <  2:4543097/101‑6 (MQ=21)
                                                            gcTGCCCCCCTC.TT.GGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCgggggtgaggc                                         >  2:341214/1‑95 (MQ=21)
                                                            gcTGCCCCCCTC.TT.GGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCgcgggtgaggc                                         >  2:4322320/1‑95 (MQ=21)
                                                            gcTGCCCCCCTC.TT.GGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTc                                          >  1:3601608/1‑101 (MQ=25)
                                                            gcTGCCCCCCTC.TT.GGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTc                                          >  2:2350610/1‑101 (MQ=25)
                                                              tGCCCCCCTC.TT.GGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCt                                        >  2:684030/1‑101 (MQ=25)
                                                               gCCCCCCTC.TT.GGCCTTCGGCCAATTCACCCCCCCCGGTATTTTGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGCGGGGAGTCCtc                                       >  2:3333491/1‑100 (MQ=25)
                                                               gCCCCCCTC.TT.GGCCTTCGGCCAATTCACCCCCCCAGGGATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCCGGGGGAGGCCTg                                       >  2:609878/1‑101 (MQ=25)
                                                                ccccccTC.TT.GGCCTTCGGCCAATTCACCCCCCGAGGGATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGa                                      >  2:4108986/1‑101 (MQ=25)
                                                                   cccTC.TT.GGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAg                                   >  2:4566733/1‑101 (MQ=37)
                                                                    ccTC.TT.GGCCTTCGGAAAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATTCCGGGGGGAGTCCTGAAAgg                                  <  1:134978/101‑1 (MQ=25)
                                                                     ctc.tT.GGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTTAAAGGa                                 >  1:952612/1‑101 (MQ=25)
                                                                     ctc.tT.GGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGa                                 <  1:5235355/101‑1 (MQ=37)
                                                                      tc.tT.GGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGAc                                <  2:5167840/101‑1 (MQ=37)
                                                                       c.tT.GGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACg                               <  2:2122091/101‑1 (MQ=37)
                                                                          ttggccTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACgg                              <  1:1039159/96‑1 (MQ=37)
                                                                          ttggccTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACgg                              >  1:2715115/6‑101 (MQ=37)
                                                                           tggccTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACggg                             >  2:4207044/5‑101 (MQ=37)
                                                                           tggccTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACggg                             >  1:1281741/5‑101 (MQ=37)
                                                                           tggccTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACggg                             >  2:85464/5‑101 (MQ=37)
                                                                           tggccTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACggg                             >  1:1488014/5‑101 (MQ=37)
                                                                             gccTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACggggg                           >  2:1477467/3‑101 (MQ=37)
                                                                              ccTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACgggggg                          <  2:5782305/100‑1 (MQ=37)
                                                                                  cGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGTGAGTCCTGAAAGGACGGGGGGCaga                      >  1:3523901/1‑101 (MQ=37)
                                                                                  cGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCaga                      >  1:83945/1‑101 (MQ=255)
                                                                                  cGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCaga                      <  1:2264790/101‑1 (MQ=255)
                                                                                   ggCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCagag                     >  1:5661350/1‑101 (MQ=255)
                                                                                   ggCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCagag                     >  2:1950199/1‑101 (MQ=37)
                                                                                      cAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGccc                  >  1:617288/1‑101 (MQ=255)
                                                                                       aaTTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGcccc                 >  1:3717657/1‑101 (MQ=37)
                                                                                       aaTTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGcccc                 >  1:4386708/1‑101 (MQ=37)
                                                                                        aTTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAATGACGGGGGGCAGAGccccc                >  1:3248765/1‑101 (MQ=37)
                                                                                        aTTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGccccc                >  2:4712850/1‑101 (MQ=37)
                                                                                         ttCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGcccccc               >  2:1132129/1‑101 (MQ=25)
                                                                                         ttCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGcccccc               >  1:4570900/1‑101 (MQ=25)
                                                                                          tcaccccccGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCt              <  1:4815591/98‑1 (MQ=255)
                                                                                            accccccGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTcc            >  2:5580576/2‑101 (MQ=255)
                                                                                            accccccGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTcc            >  2:2696365/2‑101 (MQ=255)
                                                                                              cccccGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCg          <  2:952536/101‑1 (MQ=255)
                                                                                              cccccGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCg          <  1:2926825/101‑1 (MQ=255)
                                                                                              cccccGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCg          <  1:4054548/101‑1 (MQ=255)
                                                                                               ccccGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCgg         >  2:3234337/1‑101 (MQ=255)
                                                                                                 ccGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCGCGAAAGGACGGGGGGCAGAGGCCCCCTCCCGGGt       >  1:4970158/1‑101 (MQ=25)
                                                                                                 ccGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGt       <  1:587100/101‑1 (MQ=255)
                                                                                                 ccGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCCCCCgggc       >  1:2370244/1‑100 (MQ=25)
                                                                                                  cGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTc      >  2:1638340/1‑101 (MQ=255)
                                                                                                   gAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCt     >  1:3826016/1‑101 (MQ=255)
                                                                                                   gAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCt     <  1:1950310/101‑1 (MQ=255)
                                                                                                    aGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCtt    >  1:160650/1‑101 (MQ=255)
                                                                                                     ggTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCttt   <  2:3857431/101‑1 (MQ=255)
                                                                                                      gTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTtc  >  1:1453986/1‑101 (MQ=255)
                                                                                                      gTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTtc  >  2:4066177/1‑101 (MQ=255)
                                                                                                      |                                                                                                     
GTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TT.GG.CTTCGGCCAATTCGCCCCCCGAGGTATTTGGAACAAGATGAAAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTTC  >  NC_014034/2640085‑2640285

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 27 ≤  ATCG  < 33 ≤  ATCG  < 41 ≤  ATCG