Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 3,609,070 +G 100% 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.→G100.0% 15.6 235intergenic (+9/‑162)clpB/dinBchaperone ClpB/DNA polymerase IV
Reads supporting (aligned to +/- strand):  new base (74/153):  ref base (2/6):  total (76/159)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.12e-03

Only 100 of 269 total aligned reads displayed.

TTCCGGTGACCGCCGGTCCCGACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GT.GG.GG.GCTCTG.CCCCCGCTTCGCTCCCCCCGGGATATTT.GGGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATT  >  NC_014034/3608970‑3609163
                                                                                                      |                                                                                                 
ttCCGGTGCCCGCCGGTCCCGACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.Gt                                                                                                    <  2:5219535/101‑1 (MQ=255)
ttCCGGTGACCGCCGGTCCCGACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.Gt                                                                                                    <  1:2745483/101‑1 (MQ=255)
tcccGGTCCCCGCCGGTCCCGACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.Gt                                                                                                    <  2:5010455/99‑1 (MQ=255)
  cgggTGACCGCCGGTCCCGACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GT.gg                                                                                                 <  2:2044626/99‑1 (MQ=255)
   cGGTGCCCGCCGGTCCCGACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GT.gg.g                                                                                               <  2:4274166/101‑1 (MQ=255)
    ggTGACCGCCGGTCCCGACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GT.gg.gg                                                                                              >  1:2008629/1‑101 (MQ=255)
      ggcccgccgGTCCCGACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GT.gg.gg.gg                                                                                           <  2:194599/98‑2 (MQ=255)
       gACCGCCGGTCCCGACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GT.gg.gg.ggc                                                                                          <  2:2751368/101‑3 (MQ=255)
          cgccgGTCCCGACGGGCTGGGGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GT.gg.gg.ggcgct                                                                                       >  2:3783570/1‑96 (MQ=255)
           gccgGTCCCGACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GT.gg.gg.ggccctg                                                                                      >  1:1303961/1‑95 (MQ=255)
           gccgGTCCCGACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GG.gg.gg.gggggtg                                                                                      >  2:4566435/1‑95 (MQ=255)
             cgGTCCCGACGGGCTGGGGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTG.cc                                                                                    <  2:1130660/101‑1 (MQ=255)
              ggTCCCGACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTG.ccc                                                                                   <  1:5777369/101‑1 (MQ=255)
               gTCCCGACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTG.cccc                                                                                  >  2:3844151/1‑101 (MQ=255)
                tCCCGACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTG.ccccc                                                                                 >  2:2460255/1‑101 (MQ=255)
                tCCCGACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTG.ccccc                                                                                 <  1:3658808/101‑1 (MQ=255)
                 cccGACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTG.cccccc                                                                                <  2:876111/101‑2 (MQ=255)
                  ccGACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTG.ccccccg                                                                               <  2:3030994/101‑3 (MQ=255)
                      cGGGCTGGGGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTc                                                                            <  1:4980643/101‑1 (MQ=255)
                       gggCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCg                                                                           >  2:3884725/1‑101 (MQ=255)
                       gggCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCg                                                                           <  1:1043988/101‑1 (MQ=255)
                         gCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCt                                                                         <  1:2643518/101‑1 (MQ=255)
                           tggtggGCGGGGGAAAGGTGTGGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTcc                                                                       <  1:443250/101‑1 (MQ=255)
                            ggtggtCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTccc                                                                      >  1:2197239/1‑101 (MQ=255)
                            gatggtCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTccc                                                                      >  2:3151493/3‑101 (MQ=255)
                             gtggtCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTcccc                                                                     >  2:5286566/1‑101 (MQ=255)
                              tggtCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTccccc                                                                    <  1:3434800/101‑1 (MQ=255)
                              tggtCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTccccc                                                                    <  2:704047/101‑1 (MQ=255)
                               ggtCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTcccccc                                                                   <  2:2419420/101‑1 (MQ=255)
                               ggtCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTcccccc                                                                   >  2:1699856/1‑101 (MQ=255)
                                gtCGGCGGCAAGGTCTCGCCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACAGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCg                                                                  <  1:1596271/101‑1 (MQ=255)
                                 tCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCgg                                                                 <  1:1690506/101‑1 (MQ=255)
                                   ggcggcAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGa                                                               <  2:2247459/101‑1 (MQ=255)
                                   ggcggcAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGa                                                               <  1:4182457/101‑1 (MQ=255)
                                    gcggcAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGat                                                              <  1:2373302/101‑1 (MQ=255)
                                     cggcAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGata                                                             <  1:1520832/101‑1 (MQ=255)
                                      ggcAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGatat                                                            <  1:3768932/101‑1 (MQ=255)
                                       gcAAGGTTTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAtt                                                           <  1:589199/101‑1 (MQ=38)
                                       gcAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATAtt                                                           >  2:437380/1‑101 (MQ=39)
                                         aaGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTT.g                                                        <  1:3095817/101‑1 (MQ=39)
                                         aaGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GT.GGTGG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATACTT.g                                                        >  1:1847579/1‑101 (MQ=255)
                                           ggTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTT.ggg                                                      <  2:1655186/101‑1 (MQ=39)
                                              ctcGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTT.ggggca                                                   <  1:3619426/101‑4 (MQ=38)
                                              cccGGCCAGCCAGCGCGTCCAGCCCGAACCGCGCGTGCTGAACTGATACGGGAGGG.GG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTT.ggggca                                                   <  1:4872862/99‑4 (MQ=37)
                                               tggtccagccagCGCGTCCAGACCGAACCGCGGGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTT.ggggcaa                                                  <  1:3145629/97‑5 (MQ=37)
                                                 ggCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAgg                                                 >  2:4453060/1‑101 (MQ=38)
                                                   ccagccagCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTa                                               <  2:2660349/101‑1 (MQ=38)
                                                   ccagccagCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTa                                               <  2:2217769/101‑1 (MQ=38)
                                                    cagccagcGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTaa                                              <  2:413042/101‑1 (MQ=38)
                                                    cagccagcGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTaa                                              >  1:5045095/1‑101 (MQ=38)
                                                    cagccagcGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTaa                                              >  1:4066219/1‑101 (MQ=38)
                                                     agccagcGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTaaa                                             >  1:2485238/1‑101 (MQ=38)
                                                     agccagcGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTaaa                                             >  1:463264/1‑101 (MQ=38)
                                                      gccagcGCGTCCAGACCGACCCGCGCGTGCTGAACGGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTaaaa                                            <  2:2024464/101‑1 (MQ=37)
                                                       ccagcGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAg                                           <  2:4843204/101‑1 (MQ=38)
                                                        cagcGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGc                                          <  1:3013560/101‑1 (MQ=38)
                                                        cagcGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGCAAAAGc                                          >  2:5523675/1‑101 (MQ=37)
                                                         agcGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCa                                         <  1:1610542/101‑1 (MQ=38)
                                                          gcgcGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCAt                                        <  2:1403241/101‑1 (MQ=38)
                                                          gcgcGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCAt                                        <  1:5267639/101‑1 (MQ=38)
                                                          gcgcGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCAt                                        <  1:4725153/101‑1 (MQ=38)
                                                          gcgcGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCAt                                        <  1:1325106/101‑1 (MQ=38)
                                                           cgcgTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCAtt                                       <  2:5757381/101‑1 (MQ=38)
                                                            gcgTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTg                                      <  2:1311393/101‑1 (MQ=38)
                                                            gcgTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTg                                      <  2:3655886/101‑1 (MQ=38)
                                                             cgTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGc                                     <  1:4777008/101‑1 (MQ=38)
                                                              gTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGAATTGcc                                    >  2:1787806/1‑101 (MQ=37)
                                                                ccAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCt                                  <  2:4764811/101‑1 (MQ=38)
                                                                ccAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCt                                  >  1:5351596/1‑101 (MQ=37)
                                                                 cAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCtt                                 <  2:3526540/101‑1 (MQ=38)
                                                                   gACCGAGCCGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTa                               <  1:3833097/101‑1 (MQ=37)
                                                                     ccGAACCGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATc                             <  1:4566711/101‑1 (MQ=38)
                                                                      cGAACCGCGCTGGCTGAACTGATACGGGA.GT.GG.GGTGCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCt                            <  2:10525/101‑1 (MQ=37)
                                                                      cGAACCGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCt                            >  1:1122188/1‑101 (MQ=37)
                                                                       gAACCGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCtt                           <  1:5671068/101‑1 (MQ=38)
                                                                        aaCCGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCTTg                          <  2:5127171/101‑1 (MQ=38)
                                                                        aaCCGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCTTg                          >  1:5554068/1‑101 (MQ=37)
                                                                         aCCGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCTTGc                         <  2:4920956/101‑1 (MQ=38)
                                                                          ccGCGCGTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCTTGcc                        <  1:801199/101‑1 (MQ=38)
                                                                           cgcgcgTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCTTGCCt                       >  2:4069310/1‑101 (MQ=37)
                                                                            gcgcgTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTc                      <  1:2243104/101‑1 (MQ=38)
                                                                              gcgTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAg                    >  1:4275871/1‑101 (MQ=37)
                                                                              gcgTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAg                    <  2:5348003/101‑1 (MQ=38)
                                                                               cgTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGa                   >  1:2320418/1‑101 (MQ=37)
                                                                                gTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGat                  >  1:423668/1‑101 (MQ=37)
                                                                                gTGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGat                  >  2:5043566/1‑101 (MQ=37)
                                                                                 tGCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGata                 >  2:3175836/1‑101 (MQ=37)
                                                                                  gCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGatat                <  1:5256286/101‑1 (MQ=37)
                                                                                  gCTGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGatat                <  2:2979200/101‑1 (MQ=37)
                                                                                    tGAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATcc              >  2:5658744/1‑101 (MQ=37)
                                                                                     gAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATccc             >  1:1475993/1‑101 (MQ=37)
                                                                                     gAACTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATccc             >  2:1211525/1‑101 (MQ=37)
                                                                                      aaCTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCg            >  2:3259006/1‑101 (MQ=37)
                                                                                       aCTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCgg           <  1:1070938/101‑1 (MQ=37)
                                                                                        cTGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCggg          >  1:3976820/1‑101 (MQ=37)
                                                                                         tGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCgggg         <  1:737469/101‑1 (MQ=37)
                                                                                         tGATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCgggg         >  1:638774/1‑101 (MQ=37)
                                                                                          gATACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGAATTGCCCTTTATCTTGCCTCAGATATCCCggggg        >  1:1619922/1‑101 (MQ=25)
                                                                                            tACGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTg      <  1:2805230/101‑1 (MQ=37)
                                                                                              cGGGA.GTGGG.GG.GCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGaa    <  1:437424/101‑1 (MQ=37)
                                                                                                      |                                                                                                 
TTCCGGTGACCGCCGGTCCCGACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GT.GG.GG.GCTCTG.CCCCCGCTTCGCTCCCCCCGGGATATTT.GGGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATT  >  NC_014034/3608970‑3609163

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