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

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,640,1380G→A 8.7% 19.6 ‑3.6 ‑0.2112intergenic (‑114/+162)thrS/RCAP_rcc02462threonyl‑tRNA synthetase/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (4/6):  ref base (11/91):  total (15/97)
Fisher's exact test strand distribution p-value = 2.76e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.81e-01

Only 100 of 129 total aligned reads displayed.

CGGTTTGGCGGCTACAAAGCGCCCGGTTGCGGGTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCT.C.TT.G.G.C.TTCGGCCAATTCGCCCCCCGAGGTATTTGGAACAAGATGAAAAGACCATGGGTTCATCTTGCCTCAAATATCCCG  >  NC_014034/2640038‑2640236
                                                                                                    |                                                                                                        
cgGTTTGGCGGCTACAAAGCGCCCGGTTGCGGGTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAg                                                                                                          <  1:1818346/101‑1 (MQ=255)
cgGTTTGGCGGCTACAAAGCGACCGGTTGCGGGTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAg                                                                                                          >  1:3438702/1‑101 (MQ=255)
  gTTTGGCGGCTACAAAGCGCCCGGTTGCGGGTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAggg                                                                                                        <  1:945852/101‑1 (MQ=255)
    ttGGCGGCTACAAAGCGCCCGGTTGCGGGTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAggggg                                                                                                      <  2:1044531/101‑1 (MQ=255)
    ttGGCGGCTACAAAGCGCCCGGTTGCGGGTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAggggg                                                                                                      <  1:2752149/101‑1 (MQ=255)
     tGGCGGCTACAAAGCGCCCGGTTGCGGGTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAgaagga                                                                                                     >  1:2584528/1‑96 (MQ=255)
       gcggcTACAAAGCGCCCGGTTGCGGGTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGgcg                                                                                                   <  2:2910887/101‑1 (MQ=255)
         ggcTACAAAGCGCCCGGTTGCGGGTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCt                                                                                                 <  2:3204844/101‑1 (MQ=255)
         ggcTACAAAGCGCCCGGTTGCGGGTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCt                                                                                                 <  1:507283/101‑1 (MQ=255)
         ggcTACAAAGCGCCCGGTTGCGGGTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCt                                                                                                 <  1:4064872/101‑1 (MQ=255)
          gcTACAAAGCGCCCGGTTGCGGGTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTg                                                                                                <  2:1315371/101‑1 (MQ=255)
           cTACAAAGCGCCCGTTTGCGGGTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGc                                                                                               <  1:2855087/101‑1 (MQ=255)
           cTACAAAGCGCCCGGTTGCGGGTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGc                                                                                               <  1:625799/101‑1 (MQ=255)
           cTACAAAGCGCCCGGTTGCGGGTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGc                                                                                               <  1:338697/101‑1 (MQ=255)
            tACAAAGCGCCCGGTTGCGGGTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGcc                                                                                              <  1:1069630/101‑1 (MQ=255)
             aCAAAGCGCCCGGTTGCGGGTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGccc                                                                                             <  1:1096507/101‑1 (MQ=255)
              cAAAGCGCCCGGTTGCGGGTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGcccc                                                                                            <  2:2527727/101‑1 (MQ=255)
                aaGCGCCCGGTTGCGGGTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGcccccc                                                                                          <  2:3902901/101‑1 (MQ=255)
                 aGCGCCCGGTTGCGGGTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCct                                                                                         <  2:2103801/101‑1 (MQ=255)
                   cgcCCGGTTGCGGGTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCct.c.t                                                                                     <  2:1462303/101‑1 (MQ=255)
                    gcCCGGTTGCGGGTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAGCCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCT.C.tt                                                                                    <  1:160226/101‑1 (MQ=255)
                    gcCCGGTTGCGGGTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCT.C.tt                                                                                    <  1:2851745/101‑1 (MQ=255)
                    gcCCGGTTGCGGGTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCT.C.tt                                                                                    <  1:3962300/101‑1 (MQ=255)
                    gcCCGGTTGCGGGTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCT.C.tt                                                                                    <  1:1994126/101‑1 (MQ=255)
                    gcCCGGTTGCGGGTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCT.C.tt                                                                                    <  2:2564803/101‑1 (MQ=255)
                    gcCCGGTTGCGGGTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCT.C.tt                                                                                    <  2:4081743/101‑1 (MQ=255)
                     cccGGTTGCGGGTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCT.C.TT.g                                                                                  <  2:4655727/101‑1 (MQ=255)
                       cGGTTGCGGGTTATGATGCGCCGGGTATGCGACGCGGCGGTAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCT.C.TT.G.G.c                                                                              <  2:4214977/101‑1 (MQ=255)
                        ggTTGCGGGTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCT.C.TT.G.G.cc                                                                             <  1:2340572/101‑2 (MQ=255)
                          ttGCGGGTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCT.C.TT.G.G.cctt                                                                           <  2:13897/101‑4 (MQ=255)
                          ttGCGGGTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCT.C.TT.G.G.cctt                                                                           <  1:1021698/101‑4 (MQ=255)
                           tGCGGGTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCT.C.TT.G.G.ccttc                                                                          >  2:1242490/1‑97 (MQ=255)
                           tGCGGGTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCT.C.TT.G.G.ccttc                                                                          <  1:1678021/101‑5 (MQ=255)
                           tGCGGGTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCT.C.TT.G.G.ccttc                                                                          <  2:1880270/101‑5 (MQ=255)
                             cGGGTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGGGCTGCCCCCCC.CTTT.G.C.C.TTCgg                                                                        >  1:2241583/1‑101 (MQ=255)
                              gggTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCT.C.TT.G.GCC.TTCGGc                                                                       <  2:3452555/101‑1 (MQ=255)
                                gTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCT.C.TT.G.GCC.TTCGGCCa                                                                     <  2:1845453/101‑1 (MQ=255)
                                gTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCT.C.TT.G.GCC.TTCGGCCa                                                                     <  1:4242047/101‑1 (MQ=255)
                                  tATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCT.C.TT.G.GCC.TTCGGCCAAt                                                                   <  2:4224664/101‑1 (MQ=255)
                                     gatgCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGTCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCT.C.TT.G.GCC.TTCGGCCAATTca                                                                <  1:3874659/101‑2 (MQ=255)
                                        gcgcCGGGTATGCGACGCGGCTGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCT.C.TT.G.GCC.TTCGGCCAATTCAccc                                                             <  2:2196345/101‑1 (MQ=37)
                                        gcgcCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCT.C.TT.G.GCC.TTCGGCCAATTCAccc                                                             <  2:1043641/101‑1 (MQ=255)
                                         cgcCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCT.C.TT.G.GCC.TTCGGCCCATTCAcccc                                                            >  2:4459636/1‑101 (MQ=38)
                                         cgcCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCT.C.TT.G.GCC.TTCGGCCAATTCAcccc                                                            <  1:742829/101‑1 (MQ=255)
                                          gcCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGGTGCCCCCCT.CTTT.G.C.C.TTCGGCCAATTCCccccc                                                           >  1:2171036/1‑101 (MQ=255)
                                          gcCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCT.C.TT.G.GCC.TTCGGCCAATTCAccccc                                                           <  1:4625344/101‑1 (MQ=255)
                                           ccGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCT.C.TT.G.GCC.TTCGGCCAATTCAcccccc                                                          <  2:2761747/101‑1 (MQ=255)
                                           ccGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCT.C.TT.G.GCC.TTCGGCCAATTCAcccccc                                                          <  2:1145390/101‑1 (MQ=255)
                                            cGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCT.C.TT.G.GCC.TTCGGCCAATTCACCCCCCg                                                         >  1:1371832/1‑101 (MQ=255)
                                             gggTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCT.C.TT.G.GCC.TTCGGCCAATTCACCCCCCGa                                                        <  2:1414109/101‑1 (MQ=255)
                                              ggTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCT.C.TT.G.GCC.TTCGGCCAATTCACCCCCCGAg                                                       <  2:1554785/101‑1 (MQ=255)
                                               gTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCT.C.TT.G.GCC.TTCGGCCAATTCACCCCCCGAgg                                                      <  2:2171866/101‑1 (MQ=255)
                                               gTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCT.C.TT.G.GCC.TTCGGCCAATTCACCCCCCGAgg                                                      <  2:3356069/101‑1 (MQ=255)
                                                tATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCT.C.TT.G.GCC.TTCGGCCAATTCACCCCCCGAGGt                                                     <  2:820444/101‑1 (MQ=255)
                                                 aTGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCT.C.TT.G.GCC.TTCGGCCAATTCACCCCCCGAGGTa                                                    <  2:3438448/101‑1 (MQ=255)
                                                   gCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCT.C.TT.G.GCC.TTCGGCCAATTCACCCCCCGAGGTAtt                                                  <  2:756326/101‑1 (MQ=255)
                                                            cggcAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCT.C.TTGG.G.C.TTCGGCCCATTCACCCCCCCAGGGATTTTGaggaag                                         >  2:3682733/1‑96 (MQ=25)
                                                             ggcAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCT.C.TT.G.GCC.TTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGa                                        <  1:345314/101‑1 (MQ=37)
                                                              gccagCGTGTGAACCGTGGCGCTTTTGCTCAGGCGAAAGGGGGGCTGTGCCCCCCT.C.TT.G.GCC.TTCGGGCAATTCACCCCCCGAGGTATTTGGAGCAAGAt                                       <  2:513021/97‑1 (MQ=21)
                                                              gcAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCT.C.TT.G.GCC.TTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGAt                                       <  2:3801667/101‑1 (MQ=37)
                                                               cAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCT.C.TT.G.GCC.TTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATg                                      <  2:1805999/101‑1 (MQ=37)
                                                                aGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCT.C.TT.G.GCC.TTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATGa                                     <  2:2367902/101‑1 (MQ=37)
                                                                aGGCGTATCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCGCCCCT.C.TT.G.GCC.TTCGGCCATTTCACCCCCCGAGGTACTTGGAGCAAGACGa                                     <  1:2659316/101‑1 (MQ=25)
                                                                 ggCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCT.C.TT.G.GCC.TTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATGaa                                    <  1:2184471/101‑1 (MQ=25)
                                                                 ggCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCC.C.TT.G.G.CGCTTGGGCCATTTACCCCCCCAGGTATTTGGAGGAAGATGaa                                    >  1:2089534/1‑101 (MQ=255)
                                                                  gCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCT.C.TT.G.GCC.TTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATGaaa                                   <  1:1897055/101‑1 (MQ=25)
                                                                   cGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGGGCTGCCCCCCT.C.TTCT.C.C.TTCGGCCAATTCACCCCCCCCGGTATTTTGAGCAAGATGaaaa                                  >  2:1682965/1‑101 (MQ=25)
                                                                   cGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCT.C.TT.G.GCC.TTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATGaaag                                  <  1:58641/101‑2 (MQ=25)
                                                                      gtCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCT.C.TT.G.GCC.TTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATGaaagacc                               <  2:4104802/101‑5 (MQ=25)
                                                                      gtCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCT.C.TT.G.GCC.TTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATGaaagacc                               <  1:1109101/101‑5 (MQ=25)
                                                                        catccGTGGCGCTTTTGCGCACGCCATATCGGGGCGCAGCGCCCCT.C.TT.GCA.C.TTAGGCCAATTCGCCCCCCGAGGTATTTGGAGCAAGATGaaatgccat                             <  2:2492458/98‑7 (MQ=255)
                                                                          aCCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCT.C.TT.G.GCC.TTCGGCCAATTCCCCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATgg                          <  1:1652240/101‑1 (MQ=25)
                                                                           ccGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCT.C.TT.G.GCC.TTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATggg                         <  1:401648/101‑1 (MQ=25)
                                                                              tGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCT.C.TTCG.T.C.TTTGGCCAATTCAGCCTCCGAGGTATTTTTGCCAAGATG.CAAGAGCATGGGTTc                      >  2:3141611/1‑101 (MQ=16)
                                                                              tGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCT.C.TT.G.GCC.TTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTc                      <  2:553602/101‑1 (MQ=21)
                                                                              tGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCT.C.TT.G.GCC.TTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTc                      <  1:2114179/101‑1 (MQ=21)
                                                                              gggCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCT.C.TT.G.GCC.TTCGGCCAATTCCCCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTc                      <  2:2139839/100‑1 (MQ=25)
                                                                               ggCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCT.C.TT.G.GCC.TTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCa                     <  2:2789973/101‑1 (MQ=25)
                                                                                gcgcTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCT.C.TT.G.GCC.TTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCAt                    <  2:898430/101‑1 (MQ=25)
                                                                                 cgcTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCT.C.TT.G.GCC.TTCGGCCAATTCATCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATc                   >  1:1989832/1‑101 (MQ=21)
                                                                                   cTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCT.C.TT.G.GCC.TTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCtt                 <  1:3304586/101‑1 (MQ=21)
                                                                                       tGCGCAGGCAAAAGGGGGGCGCTGCCCCCCT.C.TT.G.GCC.TTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGcgt             <  2:65143/101‑3 (MQ=21)
                                                                                        gcgcAGGCAAAAGGGGGGCGCTGCCCCCCT.C.TTAG.C.C.TTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTc            <  1:4563587/101‑1 (MQ=25)
                                                                                        gcgcAGGCAAAAGGGGGGCGCTGCCCCCCT.C.TT.G.GCC.TTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTc            <  1:2311517/101‑1 (MQ=21)
                                                                                        gcgcAGGCAAAAGGGGGGCGCTGCCCCCCT.C.TT.G.GCC.TTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTc            <  2:2584354/101‑1 (MQ=25)
                                                                                         cgcAGGCAAAAGGGGGGCGCTGCCCCCCT.C.TT.G.GCC.TTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCa           <  2:3073057/101‑1 (MQ=25)
                                                                                          gcAGGCAAAAGGGGGGCGCTGCCCCCCT.C.TT.G.GCC.TTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCaa          <  2:843951/101‑1 (MQ=25)
                                                                                           gccccgaaaAGGGGGCGCTGCCCCCCT.C.TT.G.GCC.TTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAat        <  2:2641040/95‑1 (MQ=14)
                                                                                            agccAAAAGGGGGCCGCTGCCCCCCT.C.TT.G.GTC.TTCGGCCAATTCCCCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAat        <  2:746158/98‑1 (MQ=25)
                                                                                             ggCAAAAGGGGGGCGCTGCCCCCCT.C.TT.G.GCC.TTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAata       <  2:3575984/101‑1 (MQ=21)
                                                                                             cccgaaaAGGGGGCGCTGCCCCCCT.C.TT.G.GCC.TTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAatat      <  2:2473130/97‑1 (MQ=17)
                                                                                               cgaaaAGGGGGCGCTGCCCCCCT.C.TT.G.GCC.TTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATcc    >  1:1204884/3‑101 (MQ=17)
                                                                                               cgaaaAGGGGGCGCTGCCCCCCT.C.TT.G.GCC.TTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATcc    >  2:4119652/3‑101 (MQ=17)
                                                                                               cgaaaAGGGGGCGCTGCCCCCCT.C.TT.G.GCC.TTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATcc    <  2:575675/99‑1 (MQ=17)
                                                                                               cgaaaAGGGGGCGCTGCCCCCCT.C.TT.G.GCC.TTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATcc    <  1:4234149/99‑1 (MQ=17)
                                                                                                gaaaAGGGGGCGCTGCCCCCCT.C.TT.G.GCC.TTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATccc   <  1:1215644/100‑1 (MQ=17)
                                                                                                 aaaGGGGGGCGCTGCCCCCCT.C.TT.G.GCC.TTCGGCCAATTCACCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATccc   <  1:518442/101‑1 (MQ=25)
                                                                                                 aaaGGGGGGC.CCGCCCCCCTCC.TT.G.G.C.CTCCGGCAAATCCCCCCCCGAGGTATTTGGAGCAAGATG.AAAGACCATGGGTTCATCTTGCCTCAAATATCCCg  <  1:3006725/101‑1 (MQ=21)
                                                                                                 aaaAGGGGGCGCTGCCCCCCT.C.TT.G.GCC.TTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCg  <  2:2748614/101‑1 (MQ=17)
                                                                                                 aaaAGGGGGCCCTGCCCCCCT.C.TT.G.GCC.TTCGGCCAATTCACCCCCCGGGATATTTGCGGCAAGGTGAAAGGA..TGGGGTTTATCTTGCCCCAAATATCCCg  <  1:258985/101‑1 (MQ=18)
                                                                                                    |                                                                                                        
CGGTTTGGCGGCTACAAAGCGCCCGGTTGCGGGTTATGATGCGCCGGGTATGCGACGCGGCGGCAGGCGTGTCAACCGTGGCGCTTTTGCGCAGGCAAAAGGGGGGCGCTGCCCCCCT.C.TT.G.G.C.TTCGGCCAATTCGCCCCCCGAGGTATTTGGAACAAGATGAAAAGACCATGGGTTCATCTTGCCTCAAATATCCCG  >  NC_014034/2640038‑2640236

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