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

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,640,1920A→G15.6% 60.6 ‑0.7 ‑0.0160intergenic (‑168/+108)thrS/RCAP_rcc02462threonyl‑tRNA synthetase/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (13/12):  ref base (59/76):  total (72/88)
Fisher's exact test strand distribution p-value = 5.14e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.93e-01

Only 100 of 168 total aligned reads displayed.

CGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTG.CCCCCCTC.TTGG.CTTCGGCCAATTCGCCCCCCGAGGTATTTGGAACAAGATGAAAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTTCTCAAGC  >  NC_014034/2640095‑2640291
                                                                                                    |                                                                                                   
 ggcgggAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCa                                                                                                  <  2:3254120/101‑1 (MQ=25)
 ggcggcAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCa                                                                                                  <  1:3301104/101‑1 (MQ=25)
  gcggcAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCaa                                                                                                 <  1:3237643/101‑1 (MQ=37)
   gggcAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAg                                                                                                <  2:3368113/100‑1 (MQ=37)
   cggcAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGACGGATTTGGAGCAAg                                                                                                >  2:2352871/1‑101 (MQ=25)
    ggcAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGa                                                                                               <  2:3924693/101‑1 (MQ=37)
    ggcAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGa                                                                                               <  2:3240426/101‑1 (MQ=37)
     gcAGGCGTGTCAACCGTGGCGCTTCTGCGCAGGCAAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGAt                                                                                              <  1:2344927/101‑1 (MQ=25)
     gcAGGCGGGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTAGAGCAAGAt                                                                                              <  1:3463949/101‑1 (MQ=25)
      cAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATg                                                                                             <  1:3759899/101‑1 (MQ=37)
        ggCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATGaa                                                                                           <  1:4281870/101‑1 (MQ=25)
         gCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATGaaa                                                                                          <  2:2911592/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:3277733/101‑5 (MQ=25)
                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)
                    gTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGtt                                                                              >  1:2352669/1‑101 (MQ=21)
                    ctGGTGCTTTTGCGCAGGCAAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGtt                                                                              <  1:1475094/100‑1 (MQ=21)
                     tGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTc                                                                             <  1:4109557/101‑1 (MQ=25)
                      ggCGCTTTTGCGCAGGCAAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCa                                                                            <  1:3999301/101‑1 (MQ=25)
                       gcgcTTTTGCGCAGGCAAAAGGGGTGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCAt                                                                           <  1:526971/101‑1 (MQ=21)
                           ttttGCGCAGGCAAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTg                                                                       <  2:4416073/101‑1 (MQ=25)
                                 cgacccgaaaAGGGGGCGCTG.CCCCCCGC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCaaa                                                                <  1:77379/94‑1 (MQ=17)
                                  gccccgaaaAGGGGGCGCTG.CCCCCCTC.TTGGTCTTCGGGCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAat                                                               <  1:3862236/95‑1 (MQ=9)
                                  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)
                                    ggCAAAAGGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAata                                                              <  2:2394440/101‑1 (MQ=25)
                                      cgaaaAGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGGGAAAGGA.TGGGGGTT.ATCTTGCCCCAAATATcc                                                           >  1:3545823/3‑101 (MQ=9)
                                        aaaGGGGGGCGCTC.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGTGGTATTTGGAGCAAGATG.AAAGACTATGGGTTCATCTTGCCTCAAATATccc                                                          <  1:3142659/101‑1 (MQ=18)
                                        aaaAGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCg                                                         >  2:1113414/1‑101 (MQ=17)
                                        aaaAGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCg                                                         >  1:3071039/1‑101 (MQ=17)
                                        aaaAGGGGGCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCg                                                         <  1:4557588/101‑1 (MQ=17)
                                         aaagggggCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCgg                                                        >  2:490676/4‑101 (MQ=17)
                                          aagggggCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCggg                                                       >  2:187717/3‑101 (MQ=18)
                                           agggggCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCgggg                                                      <  1:2454564/100‑1 (MQ=18)
                                           agggggCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCgggg                                                      >  1:4582159/2‑101 (MQ=18)
                                            gggggCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCggggg                                                     >  2:2963111/1‑101 (MQ=18)
                                             ggggCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCggggg                                                     <  1:4196744/101‑1 (MQ=25)
                                              gggCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCgggggtg                                                   <  1:2542500/101‑3 (MQ=18)
                                              gggCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGCGGGATTTGGAGAAAGATGAAAA.ACCATGGGTTCATCTTTCCTCAAATATCCCgggggg                                                    >  2:4276616/1‑101 (MQ=21)
                                              gggCGCTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCCCGGCATTTTGAGCAAGATG.AAAGACCATGGGGTCATCTTTCCTCAAATATCCCcggggg                                                    >  2:1127926/1‑101 (MQ=16)
                                                gcgcTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCgggggtgag                                                 <  2:2658265/101‑5 (MQ=21)
                                                 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)
                                                   cTG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTcc                                               <  1:1923838/101‑1 (MQ=25)
                                                    tG.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCt                                              >  1:3944849/1‑101 (MQ=25)
                                                     g.CCCCCCTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGGATTTTGAGCAAGATG.AAAGACCATGGGGTCATCTTGCCTCAAATATCCCGGGGGGGGTCCTg                                             >  1:2907301/1‑101 (MQ=25)
                                                         ccccTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGaaa                                          <  2:430373/101‑1 (MQ=37)
                                                          cccTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAg                                         <  1:3399039/101‑1 (MQ=37)
                                                          cccTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAg                                         >  2:3294340/1‑101 (MQ=37)
                                                          cccTC.TTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAg                                         <  2:2289244/101‑1 (MQ=37)
                                                            ctc.tTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGa                                       >  1:1723413/1‑101 (MQ=37)
                                                              c.tTGGCCTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACg                                     <  1:1567318/101‑1 (MQ=37)
                                                                 tggccTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACggg                                   >  1:4490279/5‑101 (MQ=37)
                                                                 tggccTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACggg                                   >  1:726826/5‑101 (MQ=37)
                                                                 tggccTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGACCTGAAAGGACggg                                   <  2:2008335/97‑1 (MQ=25)
                                                                  ggccTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACgggg                                  >  1:2561448/4‑101 (MQ=37)
                                                                     cTTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGc                               >  2:2452919/1‑101 (MQ=37)
                                                                        cGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCaga                            <  2:3464292/101‑1 (MQ=255)
                                                                        cGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCaga                            >  1:2789889/1‑101 (MQ=37)
                                                                        cGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGAGGGCaga                            >  2:4177219/1‑101 (MQ=37)
                                                                           ccAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGcc                         <  2:2352536/101‑1 (MQ=255)
                                                                            cAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGccc                        <  2:1048494/101‑1 (MQ=255)
                                                                            cAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGAGGGGGCCCTGAAAGGACGGGGGGCAGAGccc                        >  2:629234/1‑101 (MQ=25)
                                                                             aaTTCACCCCCCGAGGTATTTGGGGCAAGATG.AAAGACCGTGGGATCATCTTGCCTCAAATATCGCGCCGGGAGTCCTGAAAGGAGGGGGGGCACAGcgcc                       >  2:3818712/1‑101 (MQ=19)
                                                                             aaTTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGcccc                       >  2:2727196/1‑101 (MQ=37)
                                                                             aaTTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGcccc                       >  2:1247573/1‑101 (MQ=255)
                                                                              aTTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGccccc                      >  1:1362666/1‑101 (MQ=255)
                                                                              aTTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGccccc                      >  2:16762/1‑101 (MQ=37)
                                                                              aTTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGccccc                      >  2:1213101/1‑101 (MQ=37)
                                                                              aTTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCACGGGGGCAGTCCTGATAGGACGGGGGACAGACccccc                      >  2:1900101/1‑101 (MQ=25)
                                                                               ttCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGTGGCAGAGcccccc                     >  2:2976755/1‑101 (MQ=25)
                                                                                tcaccccccGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCACAGCCCCCCt                    >  2:459148/4‑101 (MQ=25)
                                                                                 caccccccGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTc                   >  2:203479/3‑101 (MQ=255)
                                                                                 caccccccGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTc                   <  2:3146606/99‑1 (MQ=255)
                                                                                    cccccGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCg                >  2:2975898/1‑101 (MQ=255)
                                                                                     ccccGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCgg               <  2:1454458/101‑1 (MQ=255)
                                                                                      cccGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCggg              >  1:3533448/1‑101 (MQ=255)
                                                                                      cccGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCgag              >  1:665174/1‑99 (MQ=255)
                                                                                       ccGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCGGAGCCCCCCTCCCGGGt             >  1:20988/1‑101 (MQ=37)
                                                                                        cGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTc            >  1:4041019/1‑101 (MQ=255)
                                                                                         gAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCt           >  1:137667/1‑101 (MQ=255)
                                                                                           ggTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCttt         <  2:1766876/101‑1 (MQ=255)
                                                                                            gTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTtc        >  2:3407244/1‑101 (MQ=255)
                                                                                             tATTTGGAGCAAGATG.AAAGTCCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTtct       <  1:380306/101‑1 (MQ=255)
                                                                                             tATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTtct       >  1:3176031/1‑101 (MQ=255)
                                                                                             tATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTtct       >  1:88993/1‑101 (MQ=255)
                                                                                             tATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTtct       >  2:2754765/1‑101 (MQ=255)
                                                                                              aTTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTtctc      >  2:705936/1‑101 (MQ=255)
                                                                                               tttGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTTCTCa     <  2:3686629/101‑1 (MQ=255)
                                                                                               tttGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTTCTCa     <  1:1042070/101‑1 (MQ=255)
                                                                                                ttGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTTCTCaa    >  2:809106/1‑101 (MQ=255)
                                                                                                  ggAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTTCTCAAGc  >  2:4621315/1‑101 (MQ=255)
                                                                                                  ggAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTTCTCAAGc  <  1:1041480/101‑1 (MQ=255)
                                                                                                    |                                                                                                   
CGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTG.CCCCCCTC.TTGG.CTTCGGCCAATTCGCCCCCCGAGGTATTTGGAACAAGATGAAAAGACCATGGGTTCATCTTGCCTCAAATATCCCGGGGGGAGTCCTGAAAGGACGGGGGGCAGAGCCCCCCTCCCGGGTCTTTCTCAAGC  >  NC_014034/2640095‑2640291

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