Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 3,609,080 +C 98.0% 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.→C98.0% 4.2 ‑0.4 ‑0.2196intergenic (+19/‑152)clpB/dinBchaperone ClpB/DNA polymerase IV
Reads supporting (aligned to +/- strand):  new base (75/117):  ref base (2/2):  total (77/119)
Fisher's exact test strand distribution p-value = 6.46e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.13e-01

Only 100 of 214 total aligned reads displayed.

GCCGGTCCCGACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GT.GGGGGCTCT.G.CCCCC..GCTTCGCTCCCCCCGGGATATTT.G.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCC  >  NC_014034/3608981‑3609174
                                                                                                       |                                                                                                  
gccgGTCCCGGCGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.g                                                                                                     <  2:1652622/101‑1 (MQ=255)
gccgGTCCCGACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.g                                                                                                     >  1:2631448/1‑101 (MQ=255)
  cgGTCCCGACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.G.cc                                                                                                  <  2:2145730/101‑1 (MQ=255)
   ggTCCCGACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.G.ccc                                                                                                 >  1:1747369/1‑101 (MQ=255)
    gTCCCGACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.G.cccc                                                                                                <  2:4425876/101‑1 (MQ=255)
     tCCCGACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.G.ccccc                                                                                               <  2:3423900/101‑1 (MQ=255)
      cccGACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.G.cccccc                                                                                              <  2:771256/101‑2 (MQ=255)
       ccGACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.G.ccccccg                                                                                             <  1:1075519/101‑3 (MQ=255)
       ccGACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.G.ccccccg                                                                                             <  1:1227075/101‑3 (MQ=255)
         gACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.G.ccccccgct                                                                                           <  1:3528142/101‑5 (MQ=255)
         gACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTCGGGGGCTCG.G.tccccagct                                                                                           >  2:3762445/1‑97 (MQ=255)
           cGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTc                                                                                        <  1:3635189/101‑1 (MQ=255)
             ggCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGc                                                                                      >  1:2018946/1‑101 (MQ=255)
              gCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCt                                                                                     >  1:282969/1‑101 (MQ=255)
              cgtccgatctggcggcAAGGTCACGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCt                                                                                     <  1:378783/91‑1 (MQ=255)
               cTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTc                                                                                    >  2:3915070/1‑101 (MQ=255)
               cTGGTGGTCGGCGGCAAGGTCTCGGCCACCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCCCTc                                                                                    >  2:3143635/1‑101 (MQ=255)
                  gtggtCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTcccc                                                                                 >  2:1026606/1‑101 (MQ=255)
                    ggtGGGCGGCAAGGTCTCGGCCAGCCAGCGCGCCCAGACCGACCCGCGCGTGCTGAACTGATAGGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTcccccc                                                                               <  1:3718731/101‑1 (MQ=255)
                    ggtCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTcccccc                                                                               <  1:1936682/101‑1 (MQ=255)
                     gtCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCg                                                                              <  1:1991801/101‑1 (MQ=255)
                     gtCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCg                                                                              >  1:2190949/1‑101 (MQ=255)
                       tggcggcAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCggg                                                                            <  2:3103768/100‑1 (MQ=255)
                       cggcggCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCggg                                                                            <  2:425649/101‑1 (MQ=255)
                       cggcggCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCggg                                                                            <  1:1727954/101‑1 (MQ=255)
                        ggcggcAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGa                                                                           <  1:3495843/101‑1 (MQ=255)
                         gcggcAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGat                                                                          <  1:1025951/101‑1 (MQ=255)
                           ggcAAGGTGTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGatat                                                                        <  1:3584164/101‑1 (MQ=255)
                             cAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATAttt                                                                      >  1:2287168/1‑101 (MQ=39)
                               aGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTT.g.g                                                                  <  2:1694808/101‑1 (MQ=39)
                                ggTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTT.g.gg                                                                 <  2:2634840/101‑1 (MQ=39)
                                ggTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTT.g.gg                                                                 >  1:3881015/1‑101 (MQ=39)
                                 gTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTT.g.ggg                                                                <  2:2126892/101‑2 (MQ=39)
                                   ctcGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTT.g.gggca                                                              <  2:3847652/101‑4 (MQ=38)
                                     cgcccagccagCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTT.g.gggcaag                                                            <  1:1654315/98‑6 (MQ=37)
                                     cGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTT.g.gggcaag                                                            <  1:3423519/101‑6 (MQ=38)
                                     cGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTT.g.gggcaag                                                            >  2:3399609/1‑96 (MQ=38)
                                      ggCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAgg                                                            >  2:4180702/1‑101 (MQ=38)
                                      ggCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAgg                                                            <  1:2980109/101‑1 (MQ=38)
                                       gccagccaGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGt                                                           >  2:2228520/1‑101 (MQ=38)
                                       gccagccaGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGt                                                           <  2:3619225/101‑1 (MQ=38)
                                        ccagccagCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTa                                                          <  1:3231126/101‑1 (MQ=38)
                                        ccagccagCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTa                                                          <  1:1395456/101‑1 (MQ=38)
                                         cagccagcGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTaa                                                         >  1:4203077/1‑101 (MQ=38)
                                           gccagcGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGNTCCCCCCGGGATATTTGG.GGCAAGGCaaaa                                                       >  1:86435/1‑101 (MQ=37)
                                           gccagcGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTaaaa                                                       <  2:563123/101‑1 (MQ=38)
                                           gccagcGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTaaaa                                                       <  1:1626014/101‑1 (MQ=38)
                                            ccagcGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAg                                                      <  2:1846820/101‑1 (MQ=38)
                                            ccagcGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAg                                                      >  1:4418620/1‑101 (MQ=38)
                                             cagcGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGc                                                     <  2:3524822/101‑1 (MQ=38)
                                             cagcGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGc                                                     >  2:2967461/1‑101 (MQ=37)
                                              agcGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCa                                                    <  1:813248/101‑1 (MQ=38)
                                                cgcgTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCAtt                                                  <  1:3447018/101‑1 (MQ=38)
                                                cgcgTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCAtt                                                  <  2:4171542/101‑1 (MQ=38)
                                                cgcgTCCAGACCGAACCGCGCGTGCTGAACTGAGACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCAtt                                                  <  1:3861478/101‑1 (MQ=37)
                                                 gcgTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTg                                                 <  2:1711778/101‑1 (MQ=38)
                                                  cgTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGc                                                >  1:3544223/1‑101 (MQ=37)
                                                   gTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGcc                                               >  2:2762017/1‑101 (MQ=37)
                                                   gTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGcc                                               >  2:1215952/1‑101 (MQ=37)
                                                    tCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGccc                                              >  1:2876897/1‑101 (MQ=37)
                                                       aGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCttt                                           <  2:3236929/101‑1 (MQ=37)
                                                       aGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCttt                                           >  1:2198289/1‑101 (MQ=37)
                                                       aGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCttt                                           >  1:3891116/1‑101 (MQ=37)
                                                       aGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCttt                                           >  2:2287520/1‑101 (MQ=37)
                                                        gACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTa                                          <  2:3435713/101‑1 (MQ=38)
                                                        gACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTa                                          >  1:548422/1‑101 (MQ=37)
                                                         aCCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTAt                                         <  1:1330792/101‑1 (MQ=38)
                                                         aCCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCCGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTAt                                         >  1:250743/1‑101 (MQ=37)
                                                          ccGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATc                                        <  2:1330961/101‑1 (MQ=38)
                                                            gAACCGCGCGTGCTGAACTGATACGGGAGGG.GGGGGCTCTCC.CCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCtt                                      <  1:956319/101‑1 (MQ=37)
                                                              aCCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGc                                    >  1:2406278/1‑101 (MQ=37)
                                                              aCCGCGCGTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTCCGCTCCCCCCGGGGTATTTGG.GGCAAGGTAAAAGCAGTGCCCTTTATCTTGc                                    >  1:1021765/1‑101 (MQ=37)
                                                                cgcgcgTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCt                                  >  2:35445/1‑101 (MQ=37)
                                                                  cgcgTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCa                                >  2:66982/1‑101 (MQ=37)
                                                                  cgcgTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCa                                <  2:1058500/101‑1 (MQ=38)
                                                                     gTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGat                             <  1:2218326/101‑1 (MQ=38)
                                                                     gTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGat                             >  2:1400797/1‑101 (MQ=37)
                                                                     gTGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGat                             <  1:1525196/101‑1 (MQ=38)
                                                                      tGCTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGata                            >  2:3980723/1‑101 (MQ=37)
                                                                        cTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATc                          <  2:463126/101‑1 (MQ=38)
                                                                        cTGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATc                          >  1:701666/1‑101 (MQ=37)
                                                                         tGAACTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATcc                         >  2:4371917/1‑101 (MQ=37)
                                                                             cTGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCACTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCggg                     >  2:2540584/1‑101 (MQ=25)
                                                                              tGATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCgggg                    >  2:1120571/1‑101 (MQ=37)
                                                                               gATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCggggg                   >  1:977205/1‑101 (MQ=37)
                                                                               gATACGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCggggg                   >  1:2154761/1‑101 (MQ=37)
                                                                                  aCGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCACGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGa                >  1:3000806/1‑101 (MQ=25)
                                                                                   cGGGA.GTGGGGGGCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGaa               >  1:71607/1‑101 (MQ=25)
                                                                                       gagtggggggCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAAttt            <  1:2116529/96‑1 (MQ=25)
                                                                                        agtggggggCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAAtttt           <  2:855971/97‑1 (MQ=37)
                                                                                        agtggggggCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAAtttt           <  1:4166529/97‑1 (MQ=37)
                                                                                        agtggggggCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAAtttt           <  2:2198150/97‑1 (MQ=37)
                                                                                         gtggggggCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTg          >  2:4212324/4‑101 (MQ=37)
                                                                                            gggggCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGgcg       <  1:4180963/101‑1 (MQ=37)
                                                                                            gggggCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGgcg       >  2:1285665/1‑101 (MQ=37)
                                                                                             ggggCGCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGgcgc      <  2:2018817/101‑1 (MQ=37)
                                                                                              gggCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCa     <  1:4024239/101‑1 (MQ=38)
                                                                                               ggCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAg    <  1:2967682/101‑1 (MQ=37)
                                                                                               ggCTCT.GCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAg    <  1:1830033/101‑1 (MQ=38)
                                                                                               ggCTCT.G.CCCCCCGGCTTCGCTCCCCCCGGGATATTT.GCGGCAAGATGAAGGGGCATGATTTCATCTTGCCCGAAATATCCCAGGGGTGAAttggccaaag     >  2:280510/1‑93 (MQ=2)
                                                                                                       |                                                                                                  
GCCGGTCCCGACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GT.GGGGGCTCT.G.CCCCC..GCTTCGCTCCCCCCGGGATATTT.G.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCC  >  NC_014034/3608981‑3609174

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