Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 3,609,109 +C 1.0% intergenic (+48/‑123) clpB → / → dinB chaperone ClpB/DNA polymerase IV

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140343,609,1091.→C 1.0% 4.6 ‑0.9 0.0518intergenic (+48/‑123)clpB/dinBchaperone ClpB/DNA polymerase IV
Reads supporting (aligned to +/- strand):  new base (3/2):  ref base (210/302):  total (213/305)
Fisher's exact test strand distribution p-value = 4.07e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 590 total aligned reads displayed.

AAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGT.GGGG.GCTCTG.CCCCCGCTTCGCTCCCCCCGGGATAT.TT.GG.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGA.GGGGCAGCGCCCCCTCCCCCACGGCCGG  >  NC_014034/3609011‑3609208
                                                                                                                                                                                                             
aaGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TT.g                                                                                                       >  2:4330598/1‑101 (MQ=39)
aaGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TT.g                                                                                                       >  2:6275881/1‑101 (MQ=39)
  ggTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TT.gg.g                                                                                                    <  2:4226091/101‑1 (MQ=39)
     ctcGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TT.gg.ggca                                                                                                 >  2:1142282/1‑98 (MQ=38)
      tcGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TT.gg.ggcaa                                                                                                <  2:7279705/101‑5 (MQ=38)
       cGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TT.gg.ggcaag                                                                                               <  1:749020/101‑6 (MQ=38)
       cGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGG.GCTCTGCCCCCCGCTTCGCTCCCCCCCGGATAT.TT.gg.ggcaag                                                                                               >  2:7369830/1‑96 (MQ=38)
          ccagccagCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTa                                                                                             <  1:5530876/101‑1 (MQ=38)
          ccagccagCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTa                                                                                             <  2:1687687/101‑1 (MQ=38)
           cagccagcGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTaa                                                                                            >  1:4023388/1‑101 (MQ=38)
           cagccagcGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTaa                                                                                            >  1:7093968/1‑101 (MQ=38)
           cagccagcGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTaa                                                                                            >  1:6544269/1‑101 (MQ=38)
             gccagcGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTaaaa                                                                                          <  2:5999936/101‑1 (MQ=38)
             gccagcGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTaaaa                                                                                          <  1:534378/101‑1 (MQ=38)
               cagcGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGc                                                                                        <  2:2072505/101‑1 (MQ=38)
               cagcGCGTCCAGACCGAACCACGCGTGCTGAACTGATACGGGAGTGGGGG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAACAGc                                                                                        <  1:7199332/101‑1 (MQ=37)
                     gTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGcc                                                                                  >  1:4199652/1‑101 (MQ=37)
                     gTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGcc                                                                                  >  1:1017086/1‑101 (MQ=37)
                      tCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGccc                                                                                 <  2:4723255/101‑1 (MQ=38)
                      tCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGccc                                                                                 <  2:8017179/101‑1 (MQ=38)
                      tCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGccc                                                                                 >  2:8809246/1‑101 (MQ=37)
                       ccAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCt                                                                                <  2:1811067/101‑1 (MQ=38)
                       ccAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCt                                                                                >  1:4586666/1‑101 (MQ=37)
                       ccAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCt                                                                                <  1:7966727/101‑1 (MQ=38)
                        cAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCtt                                                                               <  2:6809624/101‑1 (MQ=38)
                        cAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCtt                                                                               <  2:7200442/101‑1 (MQ=38)
                          gACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTa                                                                             <  2:2564765/101‑1 (MQ=38)
                          gACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTa                                                                             <  2:3686726/101‑1 (MQ=38)
                          gACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTa                                                                             <  1:2265702/101‑1 (MQ=38)
                           gccGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTAt                                                                            <  2:9054747/100‑1 (MQ=38)
                           aCCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTAt                                                                            <  1:6862664/101‑1 (MQ=38)
                            ccGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATc                                                                           <  1:839191/101‑1 (MQ=38)
                            ccGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATc                                                                           <  1:2771495/101‑1 (MQ=38)
                            ccGAACCGCGCGTGCTGAACTGATACGGGAGT.GGGGAGCTCTGCCCCCCGCTTCGCTTCCCCCTGGTTAT.TT.GGAGCATGGTAAAAGCATTGCCCTTTAtt                                                                           >  2:5603808/1‑100 (MQ=255)
                              gAACCGCGCGTGCTGAACTGATACGGGAGTGGGGG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATCtt                                                                         <  2:5703479/101‑1 (MQ=38)
                              gAACCGCGCGTGCTGAACTGATACGGGAGTGGGGG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATCtt                                                                         <  2:8562474/101‑1 (MQ=38)
                              gAACCGCGCGTGCTGAACTGATACGGGAGTGGGGG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATCtt                                                                         <  2:648927/101‑1 (MQ=38)
                               aaCCGCGCGTGCTGAACTGATACGGGAGTGGGGG.GCTCTGCCCCCCGCTTCCCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATCTTg                                                                        <  2:4969734/101‑1 (MQ=37)
                                aCCGCGCGTGCTGAACTGATACGGGAGTGGGGG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATCTTGc                                                                       <  1:517464/101‑1 (MQ=38)
                                aCCGCGCGTGCTGAACTGATACGGGAGTGGGGG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATCTTGc                                                                       >  2:7652442/1‑101 (MQ=37)
                                  cgcgcgTGCTGAACTGATACGGGAGTGGGGG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCt                                                                     <  1:3891746/101‑1 (MQ=38)
                                    cgcgTGCTGAACTGATACGGGAGTGGGGG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCa                                                                   <  1:1445865/101‑1 (MQ=38)
                                    cgcgTGCTGAACTGATACGGGAGTGGGGG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCa                                                                   <  2:2746890/101‑1 (MQ=38)
                                      cgTGCTGAACTGATACGGGAGTGGGGG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGa                                                                 <  2:2419371/101‑1 (MQ=38)
                                      cgTGCTGAACTGATACGGGAGTGGGGG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTGTATCTTGCCTCAGa                                                                 >  2:1355299/1‑101 (MQ=37)
                                       gTGCTGAACTGATACGGGAGTGGGGG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGat                                                                >  1:1525962/1‑101 (MQ=37)
                                       gTGCTGAACTGATACGGGAGTGGGGG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGat                                                                <  2:3298862/101‑1 (MQ=38)
                                       gTGCTGAACTGATACGGGAGTGGGGG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGat                                                                >  1:6314981/1‑101 (MQ=37)
                                       gTGCTGAACTGATACGGGAGTGGGGG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGat                                                                <  1:8888813/101‑1 (MQ=38)
                                         gCTGAACTGATACGGGAGTGGGGG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGatat                                                              <  1:1270883/101‑1 (MQ=37)
                                         gCTGAACTGATACGGGAGTGGGGG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGatat                                                              <  1:1173748/101‑1 (MQ=38)
                                              aCTGATACGGGAGTGGGGG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCgg                                                         >  1:20091/1‑101 (MQ=37)
                                                tGATACGGGAGTGGGGG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGAATTGCCCTTTATCTTGCCTCAGATATCCCgggg                                                       >  2:5483361/1‑101 (MQ=37)
                                                 gATACGGGAGTGGGGG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCggggg                                                      >  1:1774846/1‑101 (MQ=37)
                                                 gATACGGGAGTGGGGG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCggggg                                                      >  2:3895986/1‑101 (MQ=37)
                                                  aTACGGGAGTGGGGG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGt                                                     >  1:931059/1‑101 (MQ=37)
                                                  aTACGGGAGTCGGGG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGt                                                     <  1:8106379/101‑1 (MQ=37)
                                                   tACGGGAGTGGGGG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTg                                                    <  1:2392046/101‑1 (MQ=38)
                                                   tACGGGAGTGGGGG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTg                                                    >  2:5241882/1‑101 (MQ=37)
                                                       ggAGTGGGGG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAAtt                                                >  2:815096/1‑101 (MQ=25)
                                                        gagtggggg.gCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAAttt                                               <  2:5900328/96‑1 (MQ=37)
                                                         agtggggg.gCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGGGAAtttt                                              >  1:3290611/5‑101 (MQ=37)
                                                            ggggg.gCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGgc                                           <  1:1388375/100‑1 (MQ=38)
                                                            ggggg.gCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGgc                                           >  2:3007537/2‑101 (MQ=37)
                                                             gggg.gCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGgcg                                          >  2:5387020/1‑101 (MQ=37)
                                                               gg.gCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCa                                        <  2:2821093/101‑1 (MQ=37)
                                                                g.gCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAg                                       >  1:8746007/1‑101 (MQ=37)
                                                                   ctctGCCCCCCGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGcc                                     >  1:3692664/1‑101 (MQ=37)
                                                                   ctctGCCCCCCGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGcc                                     >  2:5932330/1‑101 (MQ=37)
                                                                    tctgccccccGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCa                                    <  1:1654837/96‑1 (MQ=255)
                                                                      tgccccccGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCaaa                                  <  1:4982537/98‑1 (MQ=255)
                                                                      tgccccccGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCaaa                                  <  2:408175/98‑1 (MQ=255)
                                                                      tgccccccGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCaaa                                  <  2:3417987/98‑1 (MQ=255)
                                                                       gccccccGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCaaaa                                 >  2:8399020/3‑101 (MQ=255)
                                                                       gccccccGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCaaaa                                 <  2:1499009/99‑1 (MQ=255)
                                                                        ccccccGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAag                                <  1:6117388/100‑1 (MQ=255)
                                                                        ccccccGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAag                                <  1:5756646/100‑1 (MQ=255)
                                                                         cccccGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAaga                               >  1:4457229/1‑101 (MQ=255)
                                                                           cccGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGA.gg                            <  1:3007635/101‑1 (MQ=255)
                                                                           cccGCTTCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGA.gg                            <  2:12568/101‑1 (MQ=255)
                                                                               ctgcGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGA.gggggc                        <  1:5921302/98‑3 (MQ=255)
                                                                               cTTCGCTCCCCCCGGGATATCTG.GG.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGA.ggcggc                        <  1:8817845/101‑5 (MQ=25)
                                                                                ttCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGA.gggggca                       >  1:2598956/1‑98 (MQ=255)
                                                                                 tCGCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGA.gggggcag                      >  1:7868723/1‑97 (MQ=255)
                                                                                   gCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAgcg                     >  1:4351874/1‑101 (MQ=25)
                                                                                   gCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAgcg                     <  1:7679811/101‑1 (MQ=25)
                                                                                   gCTCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAgcg                     <  2:7093683/101‑1 (MQ=25)
                                                                                     tCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGcc                   >  2:3552303/1‑101 (MQ=21)
                                                                                     tCCCCCCGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGcc                   >  1:4661208/1‑101 (MQ=21)
                                                                                      ccccccGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGccc                  >  1:8572907/1‑101 (MQ=21)
                                                                                      ccccccGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGccc                  <  1:9006366/101‑1 (MQ=21)
                                                                                       cccccGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGcccc                 <  2:1621666/101‑1 (MQ=21)
                                                                                        ccccGGGATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGCGGGGGCAGCGccccc                <  1:7261889/101‑1 (MQ=21)
                                                                                              gATAT.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTccccc          <  1:3547050/101‑1 (MQ=18)
                                                                                               atat.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTcccccc         >  1:8279426/1‑100 (MQ=18)
                                                                                                tat.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCAc        >  2:487263/1‑101 (MQ=21)
                                                                                                tat.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCAc        >  1:7393213/1‑101 (MQ=21)
                                                                                                tat.TTGGG.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCAc        <  2:2222034/101‑1 (MQ=18)
                                                                                                    ttggg.gCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGc     <  1:672816/98‑1 (MQ=255)
                                                                                                      ggg.gCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCg   >  2:7512615/2‑101 (MQ=255)
                                                                                                                                                                                                             
AAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGT.GGGG.GCTCTG.CCCCCGCTTCGCTCCCCCCGGGATAT.TT.GG.GCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGA.GGGGCAGCGCCCCCTCCCCCACGGCCGG  >  NC_014034/3609011‑3609208

Base quality scores:  ATCG  < 3 ≤  ATCG  < 24 ≤  ATCG  < 31 ≤  ATCG  < 35 ≤  ATCG  < 41 ≤  ATCG