Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,640,191 G→C 14.3% intergenic (‑167/+109) thrS ← / ← RCAP_rcc02462 threonyl‑tRNA synthetase/hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,640,1910G→C14.3% 66.8 ‑1.9 0.0202intergenic (‑167/+109)thrS/RCAP_rcc02462threonyl‑tRNA synthetase/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (15/14):  ref base (62/111):  total (77/125)
Fisher's exact test strand distribution p-value = 1.47e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 216 total aligned reads displayed.

ACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGG.CTTCGGCCAATTCGCCCCCCGAGGTATTTGGAACAAGATGAAAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTTCTCAAGC  >  NC_014034/2640092‑2640291
                                                                                                     |                                                                                                    
aCGCGGCGGGAGGCGTGTCAACCGTGGCGGTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTgg                                                                                                      <  2:1238442/101‑1 (MQ=25)
aCGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCTGAGGTATTTgg                                                                                                      <  1:4106284/101‑1 (MQ=37)
  gcggcgGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGag                                                                                                    <  2:124730/101‑2 (MQ=255)
   cggcggCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGagc                                                                                                   <  2:1682282/101‑3 (MQ=255)
    ggcggcAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGCGGCGCTGCCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCa                                                                                                  <  1:4246106/101‑1 (MQ=37)
      cggcAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAg                                                                                                <  1:2776809/101‑1 (MQ=25)
        gcAGGCGTGTCACCCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGAt                                                                                              <  2:419417/101‑1 (MQ=25)
        gcAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGAt                                                                                              <  1:4822526/101‑1 (MQ=37)
         cAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATg                                                                                             <  1:3322436/101‑1 (MQ=37)
           ggCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATGaa                                                                                           <  2:2929298/101‑1 (MQ=25)
           ggCGTGTCAACCGTGGCGCTTTGGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGCCTTCGGCCAATTCCCCCCCCGAGGTATTTGGAGCAAGATGaa                                                                                           <  2:3870815/101‑1 (MQ=25)
             cGTGTCACCCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATGaaag                                                                                         <  2:1998669/101‑2 (MQ=25)
             cGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATGaaag                                                                                         <  2:4678613/101‑2 (MQ=25)
                 tccaCCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCa                                                                                    <  2:4055416/98‑1 (MQ=21)
                 tCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCGCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCa                                                                                    >  2:1060519/1‑101 (MQ=21)
                 tCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCa                                                                                    <  2:3592964/101‑1 (MQ=25)
                  cAACCGTGGCGCTTTTGCGGAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCAt                                                                                   <  1:2839818/101‑1 (MQ=21)
                  cAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCAt                                                                                   <  2:2604540/101‑1 (MQ=25)
                   aaCCGTGGCGCTTTTGCGCAGGCAAAAAGGGGGGGCTGCCCCCCCCTTTGC.CCTTCGCCCATTCACCCCCCCAGGGATTTTGGGCAAAATGaaagacccag                                                                                   >  1:5053521/1‑94 (MQ=255)
                    aCCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATgg                                                                                 <  1:3436251/101‑1 (MQ=25)
                      cGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGt                                                                               <  1:2715080/101‑1 (MQ=25)
                       gTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGtt                                                                              <  2:1577555/101‑1 (MQ=25)
                        tGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTc                                                                             <  1:1426972/101‑1 (MQ=25)
                        gggCGCTTTTGCGGAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTc                                                                             <  1:1756064/100‑1 (MQ=21)
                         ggCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCa                                                                            <  1:2465438/101‑1 (MQ=25)
                         ggCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCa                                                                            <  2:2736335/101‑1 (MQ=25)
                             cTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGCCTTCGGCCAATTCCCCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCtt                                                                        <  1:3124062/101‑1 (MQ=25)
                                ttGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGcc                                                                     <  1:4656310/101‑1 (MQ=25)
                                 tGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCt                                                                    <  1:4611862/101‑1 (MQ=25)
                                  gcgcAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTc                                                                   <  2:4418752/101‑1 (MQ=25)
                                    gcAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCaa                                                                 <  1:4774475/101‑1 (MQ=25)
                                     gccccgaaaAGGGGGCGCTGCCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAat                                                               <  1:3887285/95‑1 (MQ=14)
                                     gccccgaaaAGGGGGCGCTGCCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAat                                                               <  1:3195352/95‑1 (MQ=14)
                                      aGGCAAAAGGGGGGCGCTGCCACCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAat                                                               <  1:1852141/101‑1 (MQ=25)
                                       cccgaaaAGGGGGCGCTGCCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAatat                                                             <  1:2965209/97‑1 (MQ=17)
                                          gaaaAGGGGGCGCTGCCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATccc                                                          >  2:1559102/2‑101 (MQ=17)
                                           aaaAGGGGGCGCTGCCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCg                                                         >  1:621196/1‑101 (MQ=17)
                                           aaaAGGGGGCGCTGCCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCg                                                         >  1:2578103/1‑101 (MQ=17)
                                           aaaAGGGGGCGCTGCCCCCCTC.TTGGCCTTCCGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCg                                                         >  1:4487786/1‑101 (MQ=17)
                                            aaagggggCGCTGCCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCgg                                                        <  1:531251/98‑1 (MQ=17)
                                              agggggCGCTGCCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCgggg                                                      <  2:3534039/100‑1 (MQ=18)
                                              agggggCGCTGCCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCgggg                                                      <  2:3252288/100‑1 (MQ=18)
                                               gggggCGCTGCCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCgggg                                                      <  2:3521188/101‑1 (MQ=21)
                                               gggggCGCTGCCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCCGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCCAAAATCCCcgggg                                                     >  1:4414993/1‑101 (MQ=12)
                                                ggggCGCTGCCCCCCTC.TTGGCCTTCGGCCAATTCCCCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCgggggt                                                    <  1:106036/101‑2 (MQ=18)
                                                ggggCGCTGCCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCgggggt                                                    <  1:4633495/101‑2 (MQ=18)
                                                  ggCGCTGCCCCCCTC.TTGGCCTTCGGCCAATTCCCCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGa                                                   <  1:3928763/101‑1 (MQ=25)
                                                    cgcTGCCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATTTTGGCCCAAATAATCCgggggtgagg                                                >  2:4307982/1‑96 (MQ=18)
                                                      cTGCCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCgggggtgaggcg                                              >  2:1513103/1‑94 (MQ=25)
                                                        gCCCCCCTC.TTGGCCTTCGGCCACTTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTg                                             <  1:4928576/101‑1 (MQ=25)
                                                        gCCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTg                                             <  2:4038199/101‑1 (MQ=37)
                                                          cccccTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCgggggggaggcgaagc                                          >  1:2604683/1‑91 (MQ=25)
                                                          cccccTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGaa                                           >  1:1640045/1‑101 (MQ=37)
                                                           ccccTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGaaa                                          <  2:4965158/101‑1 (MQ=37)
                                                            cccTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAg                                         <  1:2337068/101‑1 (MQ=37)
                                                             ccTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAgg                                        <  2:4643736/101‑1 (MQ=37)
                                                             ccTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAgg                                        >  2:3890968/1‑101 (MQ=37)
                                                             ccTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTAGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGCGGGAGTCCTGAAAgg                                        <  1:4260532/101‑1 (MQ=37)
                                                             ccTC.TTGGCCTGCGGACAATTCACCCCCCGAGGTATCTGGAGCAAGATG.AAAGACCATTGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAgg                                        <  1:3021335/101‑1 (MQ=21)
                                                              ctc.tTGGCCTTCGGCCATTTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGa                                       >  2:1866864/1‑101 (MQ=25)
                                                                c.tTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACg                                     >  2:711154/1‑101 (MQ=37)
                                                                  ttggccTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGCTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACgg                                    >  1:372153/6‑101 (MQ=37)
                                                                   tggccTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACggg                                   >  1:1493249/5‑101 (MQ=37)
                                                                    ggccTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACgggg                                  >  1:2542152/4‑101 (MQ=37)
                                                                      ccTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACgggggg                                >  2:2102751/2‑101 (MQ=37)
                                                                      ccTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACgggggg                                >  2:4868682/2‑101 (MQ=37)
                                                                       cTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGc                               <  1:2186958/101‑1 (MQ=255)
                                                                       cTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGc                               >  2:2881398/1‑101 (MQ=37)
                                                                       cTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGc                               >  2:3119214/1‑101 (MQ=37)
                                                                          cGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCaga                            >  2:2433240/1‑101 (MQ=255)
                                                                          cGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCaga                            <  2:4384895/101‑1 (MQ=255)
                                                                            gccacgtcaccccccGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGc                          <  1:1348731/92‑1 (MQ=37)
                                                                            gCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGc                          >  2:1312392/1‑101 (MQ=255)
                                                                              cAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGccc                        >  1:4019638/1‑101 (MQ=37)
                                                                               aaTTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGcccc                       >  2:551824/1‑101 (MQ=37)
                                                                               aaTTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGcccc                       >  2:3002397/1‑101 (MQ=255)
                                                                               aaTTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGcccc                       >  1:3586543/1‑101 (MQ=37)
                                                                                aTTCCCCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGccccc                      <  2:4227003/101‑1 (MQ=255)
                                                                                 ttCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGcccccc                     >  2:4803369/1‑101 (MQ=25)
                                                                                  tcaccccccGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCt                    >  1:4674168/4‑101 (MQ=25)
                                                                                   caccccccGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTc                   >  1:2237304/3‑101 (MQ=255)
                                                                                   caccccccGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTc                   >  2:2186221/3‑101 (MQ=255)
                                                                                     ccccccGGCGTATTTGGAGCAAAGTG.AAAGACCCTGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTccc                 <  2:3626203/101‑1 (MQ=21)
                                                                                     ccccccGAGGTATTTGGAGCAAGATG.AAAGACCAGGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTccc                 <  1:3561985/101‑1 (MQ=255)
                                                                                      cccccGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCg                >  2:3708087/1‑101 (MQ=255)
                                                                                       ccccGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCgg               <  2:1451870/101‑1 (MQ=255)
                                                                                           gAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCt           <  2:4602737/101‑1 (MQ=255)
                                                                                             ggTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCttt         >  2:1474203/1‑101 (MQ=255)
                                                                                              gTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTtc        >  1:969299/1‑101 (MQ=255)
                                                                                                aTTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTtctc      >  1:2916835/1‑101 (MQ=255)
                                                                                                 tttGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTTCTCa     <  2:3685227/101‑1 (MQ=255)
                                                                                                 tttGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTTCTCa     >  1:2893631/1‑101 (MQ=255)
                                                                                                  ttGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTTCTCaa    <  2:2298380/101‑1 (MQ=255)
                                                                                                  ttGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTTCTCaa    <  2:1639913/101‑1 (MQ=255)
                                                                                                  ttGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTTCTCaa    >  1:1101441/1‑101 (MQ=255)
                                                                                                   tGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTTCTCAAg   >  2:1540955/1‑101 (MQ=255)
                                                                                                   tGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTTCTCAAg   <  1:1084044/101‑1 (MQ=255)
                                                                                                    ggAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTTCTCAAGc  <  1:2020934/101‑1 (MQ=255)
                                                                                                    ggAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTTCTCAAGc  >  1:140193/1‑101 (MQ=255)
                                                                                                    ggAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTTCTCAAGc  >  2:834375/1‑101 (MQ=255)
                                                                                                     |                                                                                                    
ACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCTC.TTGG.CTTCGGCCAATTCGCCCCCCGAGGTATTTGGAACAAGATGAAAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTTCTCAAGC  >  NC_014034/2640092‑2640291

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