Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,640,184 G→A 14.9% 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→A14.9% 69.4 ‑3.5 0.0193intergenic (‑160/+116)thrS/RCAP_rcc02462threonyl‑tRNA synthetase/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (15/14):  ref base (48/116):  total (63/130)
Fisher's exact test strand distribution p-value = 3.00e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 207 total aligned reads displayed.

TATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGC.TTTTGCGCAGGCAAAA.GGGGGGCGCTGCCCCCCTC.TTGG.CTTC.GGCCAATTCGCCCCCCGAGGTATTTGGAACAAGATGAAAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTTC  >  NC_014034/2640086‑2640285
                                                                                                       |                                                                                                     
tATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGC.TTTTGCGCAGGCAAAA.GGGGGGCGCTGCCCCCCTC.TTGGCCTTC.GGCCAATTCACCCCCCGAGGt                                                                                                      <  1:3111059/101‑1 (MQ=255)
tATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGC.TTTTGCGCAGGCAAAA.GGGGGGCGCTGCCCCCCTC.TTGGCCTTC.GGCCAATTCACCCCCCGAGGt                                                                                                      <  2:2212560/101‑1 (MQ=255)
 aTGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGC.TTTTGCGCAGGCAAAA.GGGGGGCGCTGCCCCCCTC.TTGGCCTTC.GGCCAATTCACCCCCCGAGGTa                                                                                                     <  2:1149663/101‑1 (MQ=255)
  tGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGC.TTTTGCGCAGGCAAAA.GGGGGGCGCTGCCCCCCTC.TTGGCCTTC.GGCCAATTCACCCCCCGAGGTAt                                                                                                    <  2:4576241/101‑1 (MQ=255)
    cGACGCGGCGGCAGGCGTGTCAACCGTGGCGC.TTTTGCGCAGGCAAAA.GGGGGGCGCTGCCCCCCTC.TTGGCCTTC.GGCCAATTCACCCCCCGAGGTAttt                                                                                                  <  1:2462549/101‑1 (MQ=255)
     gACGCGGCGGCAGGCGTGTCAACCGTGGCGC.TTTTGCGCAGGCAAAA.GGGGGGCGCTGCCCCCCTC.TTGGCCTTC.GGCCAATTCACCCCCCGAGGTATTTg                                                                                                 <  2:596595/101‑1 (MQ=255)
              gcAGGCGTGTCAACCGTGGCGC.TTTTGCGCAGGCAAAA.GGGGGGCGCTGCCCCCCTC.TTGGCCTTC.GGCCAATTCACCCCCCGAGGTATTTGGAGCAAGAt                                                                                        <  1:4822526/101‑1 (MQ=37)
               cAGGCGTGTCAACCGTGGCGC.TTTTGCGCAGGCAAAA.GGGGGGCGCTGCCCCCCTC.TTGGCCTTC.GGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATg                                                                                       <  1:3322436/101‑1 (MQ=37)
                 ggCGTGTCAACCGTGGCGC.TTTGGCGCAGGCAAAA.GGGGGGCGCTGCCCCCCTC.TTGGCCTTC.GGCCAATTCCCCCCCCGAGGTATTTGGAGCAAGATGaa                                                                                     <  2:3870815/101‑1 (MQ=25)
                 gccGTGCCAACCGT.GCGCTTTTTGCG.AGGCAAAAGGGGGGGCCCTCCCCCCCCC.TGGC.CTTCGGGCAAATTCCCCCCCCGAGGTATTTGGAGCAAGATGaa                                                                                     <  2:1571107/99‑1 (MQ=255)
                  gCGTGTCAACCGTGGCGC.TTTTGCGCAGGCAAAA.GGGGGGCGCTGCCCCCCTC.TTGGCCTTC.GGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATGaaa                                                                                    <  2:1311014/101‑1 (MQ=25)
                   cGTGTCACCCGTGGCGC.TTTTGCGCAGGCAAAA.GGGGGGCGCTGCCCCCCTC.TTGGCCTTC.GGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATGaaag                                                                                   <  2:1998669/101‑2 (MQ=25)
                       tccaCCGTGGCGC.TTTTGCGCAGGCAAAA.GGGGGGCGCTGCCCCCCTC.TTGGCCTTC.GGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCa                                                                              <  2:4055416/98‑1 (MQ=21)
                       tCAACCGTGGCGC.TTTTGCGCAGGCAAAA.GGGGGGCGCTGCGCCCCTC.TTGGCCTTC.GGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCa                                                                              >  2:1060519/1‑101 (MQ=21)
                       tCAACCGTGGCGC.TTTTGCGCAGGCAAAA.GGGGGGCGCTGCCCCCCTC.TTGGCCTTC.GGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCa                                                                              <  2:3592964/101‑1 (MQ=25)
                        cAACCGTGGCGC.TTTTGCGCAGGCAAAA.GGGGGGCGCTGCCCCCCTC.TTGGCCTTC.GGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCAt                                                                             <  2:4332313/101‑1 (MQ=25)
                         aaCCGTGGCGC.TTTTGCGCAGGCAAAA.AGGGGGGGCTGCCCCCCCCTTTGC.CCTT.CGCCCATTCACCCCCCCAGGGATTTTGGGCAAAATGaaagacccag                                                                             >  1:5053521/1‑94 (MQ=255)
                          aCCGTGGCGC.TTTTGCGCAGGCAAAA.GGGGGGCGCTGCCCCCCTC.TTGGCCTTC.GGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATgg                                                                           <  1:3436251/101‑1 (MQ=25)
                            cGTGGCGC.TTTTGCGCAGGCAAAA.GGGGGGCGCTGCCCCCCTC.TTGGCCTTC.GGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGt                                                                         <  1:2715080/101‑1 (MQ=25)
                              tGGCGC.TTTTGCGCAGGCAAAA.GGGGGGCGCTGCCCCCCTC.TTGGCCTTC.GGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTc                                                                       <  1:2159917/101‑1 (MQ=25)
                              tGGCGC.TTTTGCGCAGGCAAAA.GGGGGGCGCTGCCCCCCTC.TTGGCCTTC.GGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTc                                                                       <  1:1426972/101‑1 (MQ=25)
                              tGGCGC.TTTTGCGCAGGCAAAA.GGGGGGCGCTGCCCCCCTC.TTGGCCTTC.GGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTc                                                                       <  1:1368349/101‑1 (MQ=25)
                              gggCGC.TTTTGCGGAGGCAAAA.GGGGGGCGCTGCCCCCCTC.TTGGCCTTC.GGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTc                                                                       <  1:1756064/100‑1 (MQ=21)
                               ggCGC.TTTTGCGCCGGCAAAA.GGGGGGCGCTGCCCCCCTC.TTGGCCTTC.GGCCAATTCACCCCCCAAGGAATTTGGAGCAAGATG.AAAGACCATGGGTTCa                                                                      <  2:879748/101‑1 (MQ=21)
                               ggCGC.TTTTGCGCAGGCAAAA.GGGGGGCGCTGCCCCCCTC.TTGGCCTTC.GGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCa                                                                      <  2:2736335/101‑1 (MQ=25)
                                gcgc.TTTTGCGCAGGCAAAA.GGGGGGCGCTGCCCCCCTC.TTGGCCTTC.GGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCAt                                                                     <  2:1596138/101‑1 (MQ=25)
                                   c.TTTTGCGCAGGCAAAA.GGGGGGCGCTGCCCCCCTC.TTGGCCTTC.GGCCAATTCCCCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCtt                                                                  <  1:3124062/101‑1 (MQ=25)
                                     ttttGCGCAGGCAAAA.GGGGGGCGCTGCCCCCCTC.TTGGCCTTC.GGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTg                                                                 <  2:4296686/101‑1 (MQ=25)
                                       ttGCGCAGGCAAAA.GGGGGGCGCTGCCCCCCTC.TTGGCCTTC.GGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGcc                                                               <  1:4656310/101‑1 (MQ=25)
                                        tGCGCAGGCAAAA.GGGGGGCGCTGCCCCCCTC.TTGGCCTTC.GGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCt                                                              <  1:2189074/101‑1 (MQ=25)
                                        tGCGCAGGCAAAA.GGGGGGCGCTGCCCCCCTC.TTGGCCTTC.GGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCt                                                              <  1:4611862/101‑1 (MQ=25)
                                         gcgcAGGCAAAA.GGGGGGCGCTGCCCCCCTC.TTGGCCTTC.GGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTc                                                             <  2:4826725/101‑1 (MQ=25)
                                         gcgcAGGCAAAA.GGGGGGCGCTGCCCCCCTC.TTGGCCTTC.GGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTc                                                             <  2:2550142/101‑1 (MQ=25)
                                          cgcAGGCAAAA.GGGGGGCGCTGCCCCCCTC.TTGGCCTTC.GGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCa                                                            <  1:1089099/101‑1 (MQ=25)
                                           gcAGGCAAAA.GGGGGGCGCTGCCCCCCTC.TTGGCCTTC.GGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCaa                                                           <  1:4774475/101‑1 (MQ=25)
                                            gccccgaaa.AGGGGGCGCTGCCCCCCTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAat                                                         <  1:3887285/95‑1 (MQ=14)
                                            gccccgaaa.AGGGGGCGCTGCCCCCCTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAat                                                         <  1:3195352/95‑1 (MQ=14)
                                             aGGCAAAA.GGGGGGCGCTGCCACCCTC.TTGGCCTTC.GGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAat                                                         <  1:1852141/101‑1 (MQ=25)
                                              ggCAAAA.GGGGGGCGCTGCCCCCCTC.TTGGCCTTC.GGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAata                                                        <  2:2932073/101‑1 (MQ=21)
                                              cccgaaa.AGGGGGCGCTGCCCCCCTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAatat                                                       <  1:2965209/97‑1 (MQ=17)
                                                cgaaa.AGGGGGCGCTGCCCCCCTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATcc                                                     <  2:439416/99‑1 (MQ=17)
                                                 gaaa.AGGGGGCGCTGCCCCCCTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA.TGGGGGTT.ATCTTGCCCCAAATATccc                                                    >  1:3780954/2‑101 (MQ=17)
                                                  aaa.AGGGGGCGCTGCCCCCCTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCg                                                   >  1:2578103/1‑101 (MQ=17)
                                                  aaa.AGGGGGCGCTGCCCCCCTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCg                                                   >  1:621196/1‑101 (MQ=17)
                                                  aaa.AGGGGGCGCTGCCCCCCTC.TTGGCCTTC.CGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCg                                                   >  1:4487786/1‑101 (MQ=17)
                                                    aaagggggCGCTGCCCCCCTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCgg                                                  <  1:4141951/98‑1 (MQ=18)
                                                    aaagggggCGCTGCCCCCCTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCgg                                                  <  1:531251/98‑1 (MQ=17)
                                                      ggggggCGCTGCCCCCCTC.TTGGCCTTC.GGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCggg                                                 <  1:289084/101‑1 (MQ=21)
                                                      agggggCGCTGCCCCCCTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCgggg                                                <  2:3252288/100‑1 (MQ=18)
                                                      agggggCGCTGCCCCCCTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCgggg                                                <  2:3534039/100‑1 (MQ=18)
                                                       gggggCGCTGCCCCCCTC.TTGGCCTTC.GGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCgggg                                                <  2:3521188/101‑1 (MQ=21)
                                                        ggggCGCTGCCCCCCTC.TTGGCCTTC.GGCCAATTCCCCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCgggggt                                              <  1:106036/101‑2 (MQ=18)
                                                        ggggCGCTGCCCCCCTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCgggggt                                              <  1:4633495/101‑2 (MQ=18)
                                                        ggggCGCTGCCCCCCTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCgggggg                                              >  2:3196178/1‑101 (MQ=18)
                                                          ggCGCTGCCCCCCTC.TTGGCCTTC.GGCCAATTCCCCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGa                                             <  1:3928763/101‑1 (MQ=25)
                                                            cgcTGCCCCCCTC.TTGGCCTTC.GGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGt                                           >  2:3075701/1‑101 (MQ=25)
                                                             gcTGCCCCCCTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCgggggtgaggc                                         >  2:1413667/1‑95 (MQ=21)
                                                             gcTGCCCCCCTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCgggggtgaggc                                         >  2:3701700/1‑95 (MQ=21)
                                                              cTGCCCCCCTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCgggggtgaggcg                                        >  2:1513103/1‑94 (MQ=25)
                                                                gCCCCCCTC.TTGGCCTTC.GGCCACTTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTg                                       <  1:4928576/101‑1 (MQ=25)
                                                                gCCCCCCTC.TTGGCCTTC.GGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTg                                       <  1:1181267/101‑1 (MQ=37)
                                                                gCCCCCCTC.TTGGCCTTC.GGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTg                                       <  2:4038199/101‑1 (MQ=37)
                                                                  cccccTC.TTGGCCTTC.GGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCgggggggaggcgaagc                                    >  1:2604683/1‑91 (MQ=25)
                                                                  cccccTC.TTGGCCTTC.GGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGaa                                     <  2:5053191/101‑1 (MQ=37)
                                                                  cccccTC.TTGGCCTTC.GGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGaa                                     >  1:1640045/1‑101 (MQ=37)
                                                                    cccTC.TTGGCCTTC.GGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAg                                   <  1:2337068/101‑1 (MQ=37)
                                                                     ccTC.TTGGCCTTC.GGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAgg                                  >  2:4735993/1‑101 (MQ=37)
                                                                     ccTC.TTGGCCTTC.GGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAgg                                  <  2:4643736/101‑1 (MQ=37)
                                                                     ccTC.TTGGCCTTC.GGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAgg                                  >  2:3890968/1‑101 (MQ=37)
                                                                     ccTC.TTGGCCTTC.GGCCAATTCACCCCCCGAGGTATTAGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGCGGGAGTCCTGAAAgg                                  <  1:4260532/101‑1 (MQ=37)
                                                                     ccTC.TTGGCCTGC.GGACAATTCACCCCCCGAGGTATCTGGAGCAAGATG.AAAGACCATTGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAgg                                  <  1:3021335/101‑1 (MQ=21)
                                                                      ctc.tTGGCCTTC.GGCCATTTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGa                                 >  2:1866864/1‑101 (MQ=25)
                                                                      ctc.tTGGCCTTC.GGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGa                                 >  1:2784059/1‑101 (MQ=37)
                                                                        c.tTGGCCTTC.GGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACg                               >  2:711154/1‑101 (MQ=37)
                                                                        c.tTGGCCTTC.GGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACg                               <  2:1247148/101‑1 (MQ=37)
                                                                           tggccTTC.GGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACggg                             >  2:3724532/5‑101 (MQ=37)
                                                                           tggccTTC.GGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACggg                             >  1:1493249/5‑101 (MQ=37)
                                                                           tggccTTC.GGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACggg                             >  1:3565189/5‑101 (MQ=37)
                                                                            ggccTTC.GGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACgggg                            >  1:2542152/4‑101 (MQ=37)
                                                                              ccTTC.GGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACgggggg                          >  2:2102751/2‑101 (MQ=37)
                                                                               cTTC.GGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGc                         <  2:1544919/101‑1 (MQ=255)
                                                                               cTTC.GGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGc                         >  2:2881398/1‑101 (MQ=37)
                                                                                ttC.GGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCa                        >  1:4739035/1‑101 (MQ=37)
                                                                                  c.GGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCaga                      >  2:2433240/1‑101 (MQ=255)
                                                                                    ggCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCagag                     >  1:4385269/1‑101 (MQ=255)
                                                                                    cgCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCagag                     <  2:4922051/100‑1 (MQ=37)
                                                                                      ccAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGcc                   >  1:4082888/1‑101 (MQ=37)
                                                                                      ccAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGcc                   <  1:2759494/101‑1 (MQ=255)
                                                                                        aaTTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGcccc                 >  1:3685505/1‑101 (MQ=37)
                                                                                         aTTCCCCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGccccc                <  2:4227003/101‑1 (MQ=255)
                                                                                          ttCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGcccccc               >  1:4000454/1‑101 (MQ=25)
                                                                                            caccccccGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTc             >  2:2186221/3‑101 (MQ=255)
                                                                                            caccccccGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTc             >  1:2237304/3‑101 (MQ=255)
                                                                                            caccccccGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTc             >  2:1731516/3‑101 (MQ=255)
                                                                                              ccccccGGCGTATTTGGAGCAAAGTG.AAAGACCCTGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTccc           <  2:3626203/101‑1 (MQ=21)
                                                                                              ccccccGATGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTccc           >  1:2272521/1‑101 (MQ=25)
                                                                                                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)
                                                                                                       |                                                                                                     
TATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGC.TTTTGCGCAGGCAAAA.GGGGGGCGCTGCCCCCCTC.TTGG.CTTC.GGCCAATTCGCCCCCCGAGGTATTTGGAACAAGATGAAAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTTC  >  NC_014034/2640086‑2640285

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