Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 3,609,070 +G 97.2% intergenic (+9/‑162) clpB → / → dinB chaperone ClpB/DNA polymerase IV

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140343,609,0701.→G97.2% 12.2 0.0 ‑1.3250intergenic (+9/‑162)clpB/dinBchaperone ClpB/DNA polymerase IV
Reads supporting (aligned to +/- strand):  new base (95/148):  ref base (3/4):  total (98/152)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.59e-02

Only 100 of 282 total aligned reads displayed.

TCCGGTGACCGCCGGTCCCGACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GT.GGGG.GCT.CTG.CCCCCGCTTCGCTCCCCCCGGGATATTT.G.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATT  >  NC_014034/3608971‑3609163
                                                                                                     |                                                                                                  
tCCGGTGACCGCCGGTCCCGACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GT.g                                                                                                   >  2:2309183/1‑101 (MQ=255)
  gggTGCCCGCCGGTCCCGACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GT.ggg                                                                                                 <  2:942014/100‑1 (MQ=255)
    gTGACCGCCGGTCCCGACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GT.gggg.g                                                                                              <  2:2620285/101‑1 (MQ=255)
      ggccgccgGTCCCGACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GT.gggg.ggc                                                                                            <  1:1128614/99‑3 (MQ=255)
      gACCGCCGGTCCCGACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GT.gggg.ggc                                                                                            >  1:1230843/1‑99 (MQ=255)
       aCCGCCNGTCCCGACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GT.gggg.ggct                                                                                           <  1:2497173/101‑4 (MQ=255)
        ccgccgGTCCCGACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GT.gggg.ggctc                                                                                          <  1:3312360/101‑5 (MQ=255)
         cgccgGTCCCGACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GT.gggg.ggctct                                                                                         <  1:3565430/101‑6 (MQ=255)
           ccgGGCCCGACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTG.c                                                                                      <  2:1342178/101‑1 (MQ=255)
            cgGTCCCGACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTG.cc                                                                                     <  2:1870779/101‑1 (MQ=255)
             ggTCCCGACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTG.ccc                                                                                    >  1:3943097/1‑101 (MQ=255)
              gTCCCGACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTG.cccc                                                                                   <  1:762650/101‑1 (MQ=255)
               tCCCGACGGGCGGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTG.ccccc                                                                                  <  1:466296/101‑1 (MQ=255)
                cccGACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGCGCTGAACTGATACGGGA.GTCGGGG.GCT.CGG.cccgcc                                                                                 >  2:1923855/1‑98 (MQ=255)
                   gACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTG.ccccccgct                                                                              <  2:3673119/101‑5 (MQ=255)
                    aCGGGCTGGGGGTCGGCGGCAAGGTCTCGGTCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTG.ccccccgctt                                                                             <  1:4586181/101‑6 (MQ=255)
                     cGGGTTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTc                                                                             <  1:794429/101‑1 (MQ=255)
                     cGGGCTGGTGGTCGGGGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTc                                                                             <  2:1143231/101‑1 (MQ=255)
                     cGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTc                                                                             <  1:4635373/101‑1 (MQ=255)
                     cGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GT.GGGGCGCTCCTC.CCCCCGCTTc                                                                             >  2:2731815/1‑101 (MQ=255)
                        gcgggtgggCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCt                                                                          <  2:1772468/98‑1 (MQ=255)
                           ggtggtCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTccc                                                                       <  2:1422952/101‑1 (MQ=255)
                           ggtggtCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTccc                                                                       >  2:2849662/1‑101 (MQ=255)
                            gtggtCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTcccc                                                                      >  2:2856883/1‑101 (MQ=255)
                             tggtCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTccccc                                                                     >  1:1446638/1‑101 (MQ=255)
                             tggtCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTccccc                                                                     <  1:1493637/101‑1 (MQ=255)
                              ggtCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTcccccc                                                                    >  1:1324496/1‑101 (MQ=255)
                              ggtCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTcccccc                                                                    <  1:1609397/101‑1 (MQ=255)
                               gccggcggCAAGGTCCCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCg                                                                   <  1:1815670/99‑1 (MQ=255)
                               gccggcggCAAGGCCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCg                                                                   <  1:3865828/99‑1 (MQ=255)
                                gcggcggCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCgg                                                                  <  2:333526/100‑1 (MQ=255)
                                 gggcggcAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCggg                                                                 <  1:3526870/100‑1 (MQ=255)
                                 cggcggCAAGGTCTGGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCggg                                                                 <  1:1191575/101‑1 (MQ=255)
                                  ggcggcAAGTTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCTCGGGa                                                                <  2:3740672/101‑1 (MQ=255)
                                   gcggcAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGat                                                               <  2:70478/101‑1 (MQ=255)
                                    gggcAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGata                                                              <  2:3912487/100‑1 (MQ=255)
                                     ggcAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGatat                                                             <  2:622107/101‑1 (MQ=255)
                                      gcAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGGGCTGAACTGATACGGGA.GTGGGGG.GCC.CTGCCCCCCGCTTCGCTCCCCCCGGGacttt                                                            >  1:3842967/1‑97 (MQ=255)
                                      gcAAGGTCTCGCCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGATAtt                                                            <  1:3371083/101‑1 (MQ=38)
                                       cAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGGG.GGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGATAttt                                                           <  2:1482605/101‑1 (MQ=38)
                                       cAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGATAttt                                                           <  1:2911586/101‑1 (MQ=39)
                                       cAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGATAttt                                                           <  1:3638091/101‑1 (MQ=39)
                                        gaGGTCTCGGCCAGCCAGCTCGTCTAGACCGAACCGCGCGTGCTGAACCGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCGCCCCCCGGGATAttgg                                                          <  2:1297134/100‑3 (MQ=255)
                                          ggTCGCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGATATTT.g.gg                                                      <  1:999799/101‑1 (MQ=38)
                                            tctcGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGATATTT.g.gggc                                                    <  1:546868/101‑3 (MQ=38)
                                            tctcGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGATATTT.g.gggc                                                    <  2:2154309/101‑3 (MQ=38)
                                            tctcGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGCGATATTT.g.gggc                                                    >  2:1225075/1‑99 (MQ=255)
                                             ctcGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGATATTT.g.gggca                                                   >  1:643913/1‑98 (MQ=38)
                                                 gccagccaGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGt                                                <  1:715052/101‑1 (MQ=38)
                                                 gccagccaGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGt                                                <  1:1399967/101‑1 (MQ=38)
                                                    agccagcGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTaaa                                             >  1:2031873/1‑101 (MQ=38)
                                                     gccagcGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTaaaa                                            >  2:4321741/1‑101 (MQ=38)
                                                      ccagcGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAg                                           <  1:422230/101‑1 (MQ=38)
                                                       cagcGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCACCCGGGATACTT.GAGGCAAGGCAAAAGc                                          >  2:3350642/1‑101 (MQ=255)
                                                        agcGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCa                                         >  1:1690956/1‑101 (MQ=37)
                                                         gcgcGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCAt                                        <  1:4882604/101‑1 (MQ=38)
                                                          ggcgGCCAGACCGAACCGCGCGTGCTGACCTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCAtt                                       <  2:439586/100‑1 (MQ=25)
                                                          cgcgTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCAtt                                       <  1:946540/101‑1 (MQ=38)
                                                          cgcgTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCAtt                                       <  2:4668367/101‑1 (MQ=38)
                                                           gcgTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTg                                      <  2:4908692/101‑1 (MQ=37)
                                                             gTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGcc                                    >  2:4791393/1‑101 (MQ=37)
                                                             gTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGcc                                    >  1:1094674/1‑101 (MQ=37)
                                                             gTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGcc                                    <  1:2519924/101‑1 (MQ=38)
                                                                 aGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCttt                                >  2:1732153/1‑101 (MQ=37)
                                                                 aGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCttt                                <  1:4283575/101‑1 (MQ=38)
                                                                  gACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTa                               >  1:595332/1‑101 (MQ=37)
                                                                  gACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTa                               <  2:1685249/101‑1 (MQ=38)
                                                                   aCCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTAt                              <  2:3648804/101‑1 (MQ=38)
                                                                      gAACCGCGCGTGCTGAACTGATGCGGGAGGG.GGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCtt                           <  2:4241743/101‑1 (MQ=37)
                                                                      gAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCtt                           <  1:3210001/101‑1 (MQ=38)
                                                                       aaCCGCGCGTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTg                          <  2:2068686/101‑1 (MQ=38)
                                                                       aaCCGCGCGTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCCTTATCTTg                          >  1:4200711/1‑101 (MQ=37)
                                                                          cgcgcgTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCt                       <  2:87254/101‑1 (MQ=38)
                                                                            cgcgTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCa                     <  1:4536235/101‑1 (MQ=38)
                                                                            cgcgTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCa                     <  2:1163217/101‑1 (MQ=38)
                                                                            cgcgTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCa                     <  1:411245/101‑1 (MQ=38)
                                                                               gTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGat                  <  2:4747563/101‑1 (MQ=38)
                                                                               gTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGat                  >  2:2936478/1‑101 (MQ=37)
                                                                               gTGCTGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGat                  <  1:3222121/101‑1 (MQ=38)
                                                                                 gCTGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGatat                <  1:1565943/101‑1 (MQ=37)
                                                                                 gCTGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGatat                <  2:3868214/101‑1 (MQ=37)
                                                                                 gCTGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGatat                >  2:3176437/1‑101 (MQ=37)
                                                                                  cTGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATc               <  2:1980283/101‑1 (MQ=37)
                                                                                   tGAACTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATcc              >  1:1625083/1‑101 (MQ=37)
                                                                                    taaCTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATccc             <  2:4204090/100‑1 (MQ=37)
                                                                                    gAACTGATACGGAA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATccc             >  1:2740295/1‑101 (MQ=25)
                                                                                       cTGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCggg          >  1:882296/1‑101 (MQ=37)
                                                                                        tGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCggtg         >  2:3836113/1‑99 (MQ=25)
                                                                                        tGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCgggg         >  2:4108474/1‑101 (MQ=37)
                                                                                        tGATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCgggg         >  1:4047478/1‑101 (MQ=37)
                                                                                         gATACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCggggg        >  1:3129004/1‑101 (MQ=37)
                                                                                          aTACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGGATTGCCCTTTATCTTGCCTAAGATATCCCgggggg       >  2:772797/1‑100 (MQ=25)
                                                                                          aTACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGt       >  2:158506/1‑101 (MQ=37)
                                                                                           tACGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTg      <  1:2643576/101‑1 (MQ=38)
                                                                                            aCGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGa     >  2:4523729/1‑101 (MQ=25)
                                                                                            aCGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGa     >  1:1299142/1‑101 (MQ=25)
                                                                                             cGGGA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGaa    <  1:423510/101‑1 (MQ=37)
                                                                                              gggA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAAt   <  1:4480307/101‑1 (MQ=37)
                                                                                               ggA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAAtt  >  1:2325443/1‑101 (MQ=25)
                                                                                               ggA.GTGGGGG.GCT.CTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAAtt  <  2:2233188/101‑1 (MQ=37)
                                                                                                     |                                                                                                  
TCCGGTGACCGCCGGTCCCGACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GT.GGGG.GCT.CTG.CCCCCGCTTCGCTCCCCCCGGGATATTT.G.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATT  >  NC_014034/3608971‑3609163

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 26 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG