Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,640,184 G→A 15.1% 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→A15.1% 58.8 ‑1.8 ‑0.0159intergenic (‑160/+116)thrS/RCAP_rcc02462threonyl‑tRNA synthetase/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (12/12):  ref base (46/89):  total (58/101)
Fisher's exact test strand distribution p-value = 1.68e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.98e-01

Only 100 of 169 total aligned reads displayed.

TATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTG.CCCCCCTC.TTGG.CTTCGGCCAATTCGCCCCCCGAGGTATTTGGAACAAGATGAAAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTTC  >  NC_014034/2640086‑2640285
                                                                                                     |                                                                                                     
tATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGt                                                                                                      <  2:1228405/101‑1 (MQ=255)
 aTGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTa                                                                                                     <  2:708431/101‑1 (MQ=255)
  tGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTAt                                                                                                    <  2:4474337/101‑1 (MQ=255)
   gCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTAtt                                                                                                   <  1:1056033/101‑1 (MQ=255)
   gCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTAtt                                                                                                   <  2:200493/101‑1 (MQ=255)
   gCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTAtt                                                                                                   <  1:809415/101‑1 (MQ=255)
    cGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTAttt                                                                                                  <  2:761649/101‑1 (MQ=255)
    cGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTAttt                                                                                                  <  2:1958958/101‑1 (MQ=255)
    cGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTAttt                                                                                                  <  2:1265487/101‑1 (MQ=255)
          gggggcCGGCGTGTCAACCGTTGCGCTTTTGCGCAGGCAAAGGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCa                                                                                            <  1:1490855/98‑1 (MQ=21)
          ggcgggAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCa                                                                                            <  2:3254120/101‑1 (MQ=25)
           gcggcAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCaa                                                                                           <  1:2778093/101‑1 (MQ=25)
            gggcAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAg                                                                                          <  2:3368113/100‑1 (MQ=37)
            cggcAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGACGGATTTGGAGCAAg                                                                                          >  2:2352871/1‑101 (MQ=25)
             ggcAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGa                                                                                         <  2:3240426/101‑1 (MQ=37)
             ggcAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGa                                                                                         <  2:3924693/101‑1 (MQ=37)
             ggcAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGa                                                                                         <  1:1966229/101‑1 (MQ=37)
              gcAGGCGTGTCAACCGTGGCGCTTCTGCGCAGGCAAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGAt                                                                                        <  1:2344927/101‑1 (MQ=25)
               cAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATg                                                                                       <  1:3759899/101‑1 (MQ=37)
               cAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATg                                                                                       <  1:3683988/101‑1 (MQ=25)
                 ggCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATGaa                                                                                     <  1:3058003/101‑1 (MQ=25)
                  gCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATGaaa                                                                                    <  2:2911592/101‑1 (MQ=25)
                  gCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATGaaa                                                                                    <  1:2814169/101‑1 (MQ=25)
                   ggtgtCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTG.CCCCCCTCTTTGC.CTTCGGCCAATTCCCCCCCCGAGGTATTTGGAGCAAGATGaaag                                                                                   <  1:361249/100‑2 (MQ=25)
                    gtgtCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATGaaaga                                                                                  >  2:607453/1‑99 (MQ=25)
                      gtCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATGaaagacc                                                                                <  2:1861871/101‑5 (MQ=21)
                         aaCCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATg                                                                            <  2:1804605/101‑1 (MQ=21)
                          aCCGGGGCGCTTTGGCGCAGGCAAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCCCCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATgg                                                                           <  2:188200/101‑1 (MQ=21)
                            cGTGGCGCTTTTGCGCAGGCAAAAGGGGGGGGCTGCCCCCCCCC.TTGG.CTTCTGCCAATTCCCACCCCCCGGGATTTTGGGCAAGATGAAAAAACAAAggg                                                                          >  1:567806/1‑101 (MQ=255)
                              tGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTc                                                                       <  1:1837944/101‑1 (MQ=25)
                              tGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTc                                                                       <  1:4109557/101‑1 (MQ=25)
                               ggCGCTTTTGCGCAGGCAAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCa                                                                      <  1:3999301/101‑1 (MQ=25)
                                gcgcTTTTGCGCAGGCAAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCAt                                                                     <  2:2220102/101‑1 (MQ=25)
                                    ttttGCGCAGGCAAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTg                                                                 <  2:4416073/101‑1 (MQ=25)
                                        gcgcAGGCAAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTc                                                             <  2:1669204/101‑1 (MQ=25)
                                         cgcAGGCAAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCa                                                            <  2:4058540/101‑1 (MQ=25)
                                          cgacccgaaaAGGGGGCGCTG.CCCCCCGC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCaaa                                                          <  1:77379/94‑1 (MQ=17)
                                           cAGGCAAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCaaa                                                          <  2:1112464/101‑1 (MQ=25)
                                            ccccgaaaAGGGGGCTCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAata                                                        <  1:1989155/96‑1 (MQ=17)
                                            aGGCAAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAat                                                         <  1:4694603/101‑1 (MQ=21)
                                            aGGCAAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAat                                                         <  2:2473578/101‑1 (MQ=25)
                                              ccgaaaAGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATc                                                      >  2:4203224/4‑101 (MQ=17)
                                              ccgaaaAGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATc                                                      <  2:2346504/98‑1 (MQ=17)
                                               cAAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATc                                                      <  2:3398727/101‑1 (MQ=25)
                                                caaaGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATcc                                                     <  2:2358587/100‑1 (MQ=21)
                                                 aaaGGGGGGCGCTC.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGTGGTATTTGGAGCAAGATG.AAAGACTATGGGTTCATCTTGCCTCAAATATccc                                                    <  1:3142659/101‑1 (MQ=18)
                                                 aaaAGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCg                                                   <  1:3545633/101‑1 (MQ=17)
                                                 aaaAGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCg                                                   >  2:1113414/1‑101 (MQ=17)
                                                 aaaAGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCg                                                   <  1:4557588/101‑1 (MQ=17)
                                                  aaagggggCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCgg                                                  >  2:490676/4‑101 (MQ=17)
                                                    ggggggCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCggg                                                 >  1:3464504/1‑101 (MQ=21)
                                                    agggggCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCgggg                                                >  2:2484948/2‑101 (MQ=18)
                                                    agggggCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCgggg                                                <  1:2454564/100‑1 (MQ=18)
                                                     gggggCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCggggg                                               >  2:2963111/1‑101 (MQ=18)
                                                      ggggCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCgggggg                                              >  1:920075/1‑101 (MQ=18)
                                                      ggggCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCggggg                                               <  1:4196744/101‑1 (MQ=25)
                                                       gggCGCTG.CCCCCCTC.TTGGCCTTCGGGCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCgggggg                                              <  1:3107933/101‑1 (MQ=25)
                                                       gggCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCgggggtg                                             <  1:2542500/101‑3 (MQ=18)
                                                       gggCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCCCGGCATTTTGAGCAAGATG.AAAGACCATGGGGTCATCTTTCCTCAAATATCCCcggggg                                              >  2:1127926/1‑101 (MQ=16)
                                                        ggCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGa                                             <  2:3708193/101‑1 (MQ=25)
                                                         gcgcTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCgggggtgag                                           <  2:2658265/101‑5 (MQ=21)
                                                         gcgcTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAg                                            <  1:3925227/101‑1 (MQ=25)
                                                          cgcTG.CCCCCCTC.TTGGCCTTCGGCCAATTCTCCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCgggggtgagg                                          <  2:2803561/101‑6 (MQ=21)
                                                          cgcTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCgggggtgagg                                          <  1:1167898/101‑6 (MQ=21)
                                                          cgcTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCgggggtgagg                                          <  1:2825796/101‑6 (MQ=21)
                                                          cgcTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGt                                           <  1:2353004/101‑1 (MQ=25)
                                                             tG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCt                                        <  1:3350898/101‑1 (MQ=25)
                                                             tG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCt                                        >  1:3944849/1‑101 (MQ=25)
                                                                  ccccTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGaaa                                    <  2:430373/101‑1 (MQ=37)
                                                                  ccccTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGaaa                                    >  1:3286893/1‑101 (MQ=37)
                                                                   cccTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAg                                   <  1:3399039/101‑1 (MQ=37)
                                                                   cccTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAg                                   >  2:3294340/1‑101 (MQ=37)
                                                                      tc.tTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGAc                                <  1:1459164/101‑1 (MQ=37)
                                                                       c.tTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACg                               >  1:4085158/1‑101 (MQ=25)
                                                                       c.tTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACg                               <  1:1567318/101‑1 (MQ=37)
                                                                          tggccTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACggg                             >  2:4381375/5‑101 (MQ=37)
                                                                          tggccTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACggg                             >  1:4490279/5‑101 (MQ=37)
                                                                           ggccTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACgggg                            >  1:2561448/4‑101 (MQ=37)
                                                                                 cGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCaga                      <  2:3464292/101‑1 (MQ=255)
                                                                                   gCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGc                    <  1:4382493/101‑1 (MQ=255)
                                                                                    ccAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGcc                   <  2:3257122/101‑1 (MQ=255)
                                                                                     cAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGAGGGGGCCCTGAAAGGACGGGGGGCAGAGccc                  >  2:629234/1‑101 (MQ=25)
                                                                                      aaTTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGcccc                 >  2:2727196/1‑101 (MQ=37)
                                                                                       aTTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGccccc                >  1:4674788/1‑101 (MQ=255)
                                                                                       aTTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGccccc                >  1:1362666/1‑101 (MQ=255)
                                                                                       aTTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGccccc                >  2:1213101/1‑101 (MQ=37)
                                                                                       aTTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCACGGGGGCAGTCCTGATAGGACGGGGGACAGACccccc                >  2:1900101/1‑101 (MQ=25)
                                                                                        ttCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGTGGCAGAGcccccc               >  2:2976755/1‑101 (MQ=25)
                                                                                        ttCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGcccccc               >  2:2976754/1‑101 (MQ=25)
                                                                                         tcaccccccGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCACAGCCCCCCt              >  2:459148/4‑101 (MQ=25)
                                                                                          caccccccGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTc             >  2:203479/3‑101 (MQ=255)
                                                                                           accccccGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTcc            >  1:1145054/2‑101 (MQ=255)
                                                                                           accccccGAGGGATTTGGAGCAAGATG.AAAGACCATGGGTGCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTcc            <  1:2169463/100‑1 (MQ=25)
                                                                                             cccccGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCg          <  2:2421040/101‑1 (MQ=255)
                                                                                              ccccGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCgg         <  2:1454458/101‑1 (MQ=255)
                                                                                               cccGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCgag        >  1:665174/1‑99 (MQ=255)
                                                                                                  gAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCt     >  1:137667/1‑101 (MQ=255)
                                                                                                   aGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCtt    >  2:1393269/1‑101 (MQ=255)
                                                                                                   aGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCtt    >  1:2350304/1‑101 (MQ=255)
                                                                                                    ggTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCttt   <  2:1766876/101‑1 (MQ=255)
                                                                                                     |                                                                                                     
TATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTG.CCCCCCTC.TTGG.CTTCGGCCAATTCGCCCCCCGAGGTATTTGGAACAAGATGAAAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTTC  >  NC_014034/2640086‑2640285

Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 26 ≤  ATCG  < 33 ≤  ATCG  < 41 ≤  ATCG