Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 3,609,080 +C 97.5% intergenic (+19/‑152) clpB → / → dinB chaperone ClpB/DNA polymerase IV

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140343,609,0801.→C97.5% 6.8 ‑1.1 ‑0.5202intergenic (+19/‑152)clpB/dinBchaperone ClpB/DNA polymerase IV
Reads supporting (aligned to +/- strand):  new base (66/130):  ref base (3/2):  total (70/132)
Fisher's exact test strand distribution p-value = 3.41e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 3.27e-01

Only 100 of 226 total aligned reads displayed.

CCGGTCCCGACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGT.GGGGG.CTCTG.C.CC.CC..GCTTCGCTCCCCCCGGGATATTT.G.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCC  >  NC_014034/3608982‑3609174
                                                                                                     |                                                                                                    
ccgGTCCCGACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGG.CTCTG.c                                                                                                     <  2:3041122/101‑1 (MQ=255)
 cgGTCCCGACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGG.CTCTG.c.c                                                                                                   >  2:485857/1‑101 (MQ=255)
 cgGTCCCGACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGG.CTCTG.c.c                                                                                                   <  2:3149393/101‑1 (MQ=255)
   gTCCCGACGGGCTGGTGGTGGGGGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGG.CTCTG.c.cc.c                                                                                                <  2:564517/101‑1 (MQ=255)
    tCCCGACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGG.CTCTG.c.cc.cc                                                                                               <  1:2987774/101‑1 (MQ=255)
      ccGACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGG.CTCTG.c.cc.cccg                                                                                             <  2:934516/101‑3 (MQ=255)
      ccGACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGG.CTCTG.c.cc.cccg                                                                                             <  1:1912528/101‑3 (MQ=255)
      ccGACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGCCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGG.CTCTG.t.cc.cccg                                                                                             >  2:1607959/1‑99 (MQ=255)
       cGACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGG.CTCTG.c.cc.cccgc                                                                                            >  2:3327619/1‑98 (MQ=255)
         aCGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGG.CTCTG.c.cc.cccgctt                                                                                          <  1:217467/101‑6 (MQ=255)
          cGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTc                                                                                        <  1:1221929/101‑1 (MQ=255)
           gggCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCg                                                                                       >  1:2360310/1‑101 (MQ=255)
             gCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCt                                                                                     <  1:1086479/101‑1 (MQ=255)
             gCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGG.CTCTG.C.CCTCC..GCTTCGcc                                                                                     >  1:1529764/1‑100 (MQ=255)
                    gtCGGCGGGATGTTCTTGACCATCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATAGGGGAGT.GGGGGACTCTGCC.CC.CC..GCTGCGTTCCCCCCg                                                                              <  2:2060110/101‑1 (MQ=255)
                     tCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCgg                                                                             <  1:620047/101‑1 (MQ=255)
                     tCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCgg                                                                             <  2:371333/101‑1 (MQ=255)
                      cggcggCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCggg                                                                            >  2:3475485/1‑101 (MQ=255)
                      cggcggCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCggg                                                                            >  2:3780658/1‑101 (MQ=255)
                      cggcggCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGG.CTCTG.C.CC.CCA.GCTTCGCTCCCCCCggg                                                                            >  2:3507351/1‑101 (MQ=255)
                         gggcAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGata                                                                         <  2:24822/100‑1 (MQ=255)
                         cggcAAGTTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGata                                                                         <  2:939520/101‑1 (MQ=255)
                         cggcAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGata                                                                         <  1:4137611/101‑1 (MQ=255)
                          ggcAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGatat                                                                        <  2:157847/101‑1 (MQ=255)
                           gcAAGGTCTCGGACAGCCAGCGAGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATAtt                                                                       <  2:686361/101‑1 (MQ=38)
                            cAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATAttt                                                                      >  1:1954027/1‑101 (MQ=39)
                            cAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATAttt                                                                      <  2:2307651/101‑1 (MQ=39)
                               ggTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTT.g.gg                                                                 >  2:13873/1‑101 (MQ=39)
                                gTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGGTATTT.g.ggg                                                                >  2:2715043/1‑100 (MQ=255)
                                gTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTT.g.ggg                                                                <  2:1798172/101‑2 (MQ=39)
                                   tcGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTT.g.gggcaa                                                             <  2:390278/101‑5 (MQ=38)
                                    cGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTT.g.gggcaag                                                            >  1:3365667/1‑96 (MQ=38)
                                     ggCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAgg                                                            >  2:2261447/1‑101 (MQ=38)
                                     ggCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAgg                                                            <  1:3464317/101‑1 (MQ=38)
                                      gccagccaGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGt                                                           <  1:45147/101‑1 (MQ=38)
                                       ccagccagCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTa                                                          <  2:595088/101‑1 (MQ=38)
                                       ccagccagCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTa                                                          <  1:3820453/101‑1 (MQ=38)
                                         agccagcGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTaaa                                                        <  1:1894487/101‑1 (MQ=38)
                                         agccagcGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTaaa                                                        <  2:2019359/101‑1 (MQ=38)
                                         agccagcGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTaaa                                                        <  2:833525/101‑1 (MQ=38)
                                          gccagcGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTaaaa                                                       <  2:4181559/101‑1 (MQ=38)
                                          gccagcGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTaaaa                                                       <  2:3922983/101‑1 (MQ=38)
                                            cagcGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGc                                                     >  2:1219191/1‑101 (MQ=37)
                                            cagcGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGc                                                     <  2:2964739/101‑1 (MQ=38)
                                              gcgcGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCAt                                                   <  2:3633770/101‑1 (MQ=37)
                                               cgcgTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCAtt                                                  <  1:636511/101‑1 (MQ=38)
                                               cgcgTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGG.CTCTG.CTCC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCAtt                                                  <  2:3658930/101‑1 (MQ=38)
                                                gcgTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTg                                                 <  1:2679810/101‑1 (MQ=38)
                                                 cggccAGACGGAATCGCGGGTGCTGATCTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTgg                                                <  2:1250934/98‑2 (MQ=21)
                                                 cgTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGc                                                <  2:3423334/101‑1 (MQ=38)
                                                 cgTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGc                                                <  2:4020654/101‑1 (MQ=38)
                                                  gTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGcc                                               >  2:947762/1‑101 (MQ=37)
                                                   tCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGccc                                              >  2:268771/1‑101 (MQ=37)
                                                     cAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCtt                                            <  1:1301628/101‑1 (MQ=38)
                                                     cAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCtt                                            <  2:2578148/101‑1 (MQ=38)
                                                      aGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCttt                                           >  1:1738050/1‑101 (MQ=37)
                                                         ccGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATc                                        <  1:1546955/101‑1 (MQ=38)
                                                          cGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCt                                       <  2:3743359/101‑1 (MQ=38)
                                                           gAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGTTATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCtt                                      >  1:4009239/1‑101 (MQ=37)
                                                           gAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCtt                                      <  2:3018973/101‑1 (MQ=38)
                                                            aaCCGCGCGTGCTGAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTg                                     <  1:2102559/101‑1 (MQ=38)
                                                            aaCCGCGCGTGCTGAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTg                                     <  1:3296614/101‑1 (MQ=38)
                                                              ccGCGCGTGCTGAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGcc                                   <  1:2245367/101‑1 (MQ=38)
                                                               cgcgcgTGCTGAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCt                                  <  2:3375864/101‑1 (MQ=38)
                                                               cgcgcgTGCTGAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCACTGCCCTTTATCTTGCCt                                  >  1:1132770/1‑101 (MQ=37)
                                                                gcgcgTGCTGAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTc                                 <  1:3100585/101‑1 (MQ=38)
                                                                 cgcgTGCTGAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCa                                <  1:2606972/101‑1 (MQ=38)
                                                                  gcgTGCTGAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAg                               <  2:1048147/101‑1 (MQ=38)
                                                                   cgTGCTGAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGa                              >  1:371677/1‑101 (MQ=37)
                                                                    gTGCTGAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGat                             <  2:1877566/101‑1 (MQ=38)
                                                                    gTGCTGAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGat                             <  1:3536293/101‑1 (MQ=38)
                                                                    gTGCTGAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGat                             <  2:1529660/101‑1 (MQ=38)
                                                                    gTGCTGAACTGAGACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGat                             >  1:2296960/1‑101 (MQ=37)
                                                                     tGCTGAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGata                            <  2:1131397/101‑1 (MQ=38)
                                                                      gCTGAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGatat                           >  1:3003522/1‑101 (MQ=37)
                                                                        tGAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATcc                         >  1:3216813/1‑101 (MQ=37)
                                                                         gAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATccc                        <  1:1389715/101‑1 (MQ=37)
                                                                         gAACTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATccc                        <  1:2962721/101‑1 (MQ=37)
                                                                           aCTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCgg                      >  1:4124772/1‑101 (MQ=37)
                                                                           aCTGATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCgg                      <  1:902674/101‑1 (MQ=37)
                                                                              gATACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCGCggggg                   >  1:699697/1‑101 (MQ=37)
                                                                               aTACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGt                  <  1:3324378/101‑1 (MQ=37)
                                                                               aTACGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGt                  >  1:1155267/1‑101 (MQ=37)
                                                                                 gcGGGAGT.GGGGGTCTCTGCC.CC.CC..GCTTCGCTCCCCCCGAGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGa                <  2:699620/100‑1 (MQ=37)
                                                                                 aCGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGa                <  1:2852003/101‑1 (MQ=37)
                                                                                 aCGGGAGTGGGGGG.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGa                >  2:2323555/1‑101 (MQ=25)
                                                                                     gagtgggggg.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATCGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAAttt            >  2:4226862/6‑101 (MQ=25)
                                                                                        tgggggg.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTgg         <  1:2409894/99‑1 (MQ=38)
                                                                                        tgggggg.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTgg         >  2:1445531/3‑101 (MQ=37)
                                                                                         gggggg.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGgc        >  1:2692011/2‑101 (MQ=37)
                                                                                          ggggg.CTCTG.C.CC.CCCGGCTTCGCTCCCCCCGGGATATTT.GCGGCAAGATGAAGGGGCATGATTTCATCTTGCCCGAAATATCCCAGGGGTGAAttggcca        >  2:3218750/1‑96 (MQ=2)
                                                                                           gggg.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGgcgc      <  1:2153247/101‑1 (MQ=38)
                                                                                           gggg.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGgcgc      <  1:36532/101‑1 (MQ=37)
                                                                                           gggg.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGgcgc      <  1:3760605/101‑1 (MQ=38)
                                                                                            ggg.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCa     >  2:3278513/1‑101 (MQ=37)
                                                                                             gg.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAg    >  1:969536/1‑101 (MQ=37)
                                                                                             gg.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAg    <  2:1028312/101‑1 (MQ=38)
                                                                                             gg.CTCTGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAg    <  2:215133/101‑1 (MQ=38)
                                                                                                ctctGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGcc  <  1:2731527/101‑1 (MQ=37)
                                                                                                ctctGCC.CC.CC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGcc  <  1:343987/101‑1 (MQ=37)
                                                                                                     |                                                                                                    
CCGGTCCCGACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGT.GGGGG.CTCTG.C.CC.CC..GCTTCGCTCCCCCCGGGATATTT.G.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCC  >  NC_014034/3608982‑3609174

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 29 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG