Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 3,609,085 +G 1.5% intergenic (+24/‑147) clpB → / → dinB chaperone ClpB/DNA polymerase IV

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140343,609,0852.→G 1.5% 2.1 ‑2.9 ‑1.1201intergenic (+24/‑147)clpB/dinBchaperone ClpB/DNA polymerase IV
Reads supporting (aligned to +/- strand):  new base (3/0):  ref base (75/123):  total (78/123)
Fisher's exact test strand distribution p-value = 5.71e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.82e-02

Only 100 of 231 total aligned reads displayed.

CCCGACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GT.GGG.GGCTCT.G.CC.CCCGCTTCGCTCCCCCCGGGATATTT.GGGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGG  >  NC_014034/3608987‑3609184
                                                                                                                                                                                                             
 ccGACGGGCTGGTGGTCGGCGGCAAGGTCTCGCCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.G.cc.ccccg                                                                                                    <  1:1781415/101‑3 (MQ=255)
  cGACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.G.cc.ccccgc                                                                                                   <  2:837553/101‑4 (MQ=255)
   gACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.G.cc.ccccgct                                                                                                  <  1:4583727/101‑5 (MQ=255)
    aCGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.G.cc.ccccgctt                                                                                                 >  2:216057/1‑96 (MQ=255)
    aCGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.G.cc.ccccgctt                                                                                                 <  2:4126360/101‑6 (MQ=255)
    aCGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.G.cc.ccccgctt                                                                                                 <  2:4380207/101‑6 (MQ=255)
     ggggCTGGTGGGCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTc                                                                                                 <  2:1024156/100‑1 (MQ=255)
     cGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTTCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTc                                                                                                 <  2:511493/101‑1 (MQ=255)
     cGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTc                                                                                                 <  2:1374909/101‑1 (MQ=255)
     cGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTc                                                                                                 <  1:1007139/101‑1 (MQ=255)
      gggCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCg                                                                                                >  2:4302492/1‑101 (MQ=255)
      gggCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCg                                                                                                <  2:3828259/101‑1 (MQ=255)
      gggCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCg                                                                                                >  2:2885367/1‑101 (MQ=255)
          tggtggTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTcc                                                                                            >  1:372131/1‑101 (MQ=255)
          tggtggTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTcc                                                                                            >  2:719232/1‑101 (MQ=255)
            gtggtCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTcccc                                                                                          >  2:3940711/1‑101 (MQ=255)
             tggtCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTccccc                                                                                         >  2:175067/1‑101 (MQ=255)
              ggtCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTcccccc                                                                                        <  2:2820922/101‑1 (MQ=255)
               gtCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCg                                                                                       >  1:68618/1‑101 (MQ=255)
               gtCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCg                                                                                       >  2:1647254/1‑101 (MQ=255)
                gcggcggCAGGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCgg                                                                                      <  1:482854/100‑1 (MQ=255)
                gcggcggCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCgg                                                                                      <  2:3018284/100‑1 (MQ=255)
                 cggcggCAAGGTCTCGGCCAGCCAGCGCGTCCAGCCCGAACCGCGCGTGCTGAACTGATACGGGA.GT.GGGNNNNNNT.G.CCNCNCCCNTCGCTCCCCCCggg                                                                                     >  1:3614007/1‑101 (MQ=255)
                 cggcggCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCggg                                                                                     <  1:3463148/101‑1 (MQ=255)
                  ggcggcAGGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGa                                                                                    <  2:2377244/101‑1 (MQ=255)
                  ggcggcAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGa                                                                                    <  1:1903825/101‑1 (MQ=255)
                   gcggcAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGat                                                                                   <  1:2741911/101‑1 (MQ=255)
                    cggcAGGTTCCCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGata                                                                                  <  2:1996225/101‑1 (MQ=255)
                     ggcAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGatat                                                                                 <  2:2147733/101‑1 (MQ=255)
                      gcAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATAtt                                                                                <  1:4540599/101‑1 (MQ=39)
                      gcAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GT.GGGTGGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATAtt                                                                                >  1:150947/1‑101 (MQ=255)
                      gcAAGGTCTCGCCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATAtt                                                                                <  1:4220030/101‑1 (MQ=38)
                       cAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATAttt                                                                               <  1:4457739/101‑1 (MQ=39)
                        aaGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATATTT.g                                                                             <  2:2169869/101‑1 (MQ=39)
                        aaGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATATTT.g                                                                             <  2:4614979/101‑1 (MQ=39)
                        aaGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATATTT.g                                                                             >  2:58720/1‑101 (MQ=39)
                          ggTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATATTT.ggg                                                                           >  2:1240711/1‑101 (MQ=39)
                            tctcGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATATTT.ggggc                                                                         <  2:3559380/101‑3 (MQ=38)
                            tctcGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATATTT.ggggc                                                                         >  2:2523327/1‑99 (MQ=38)
                             ctcGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCCGGATATTT.ggggca                                                                        >  1:1750346/1‑98 (MQ=255)
                              tcGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATATTT.ggggcaa                                                                       >  2:2162669/1‑97 (MQ=38)
                                  ccagccagCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTa                                                                    >  2:2742224/1‑101 (MQ=38)
                                   cagccagcGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTaa                                                                   <  2:46780/101‑1 (MQ=38)
                                   cagccagcGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTaa                                                                   <  2:2644329/101‑1 (MQ=38)
                                    agccagcGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTaaa                                                                  <  2:1897923/101‑1 (MQ=38)
                                      ccagcGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAg                                                                <  2:2291976/101‑1 (MQ=38)
                                      ccagcGCGTCCAGACCGAACCGCGCGCGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGACATTTGGGGCAAGGTAAAAg                                                                >  1:2476690/1‑101 (MQ=37)
                                       cagcGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGc                                                               >  2:970461/1‑101 (MQ=37)
                                       cagcGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGc                                                               >  1:1178691/1‑101 (MQ=37)
                                         gcgcGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCAt                                                             <  2:667746/101‑1 (MQ=38)
                                          ggcgTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCAtt                                                            <  1:3427233/100‑1 (MQ=37)
                                          cgcgTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCAtt                                                            <  1:4038556/101‑1 (MQ=38)
                                            cgTCCGGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGc                                                          <  1:2180538/101‑1 (MQ=37)
                                            cgTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGc                                                          <  1:2977908/101‑1 (MQ=38)
                                             gTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGcc                                                         >  2:3406401/1‑101 (MQ=37)
                                             gTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGcc                                                         >  2:3767750/1‑101 (MQ=37)
                                               ccAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCt                                                       <  1:2322047/101‑1 (MQ=38)
                                                cAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCtt                                                      <  1:2926221/101‑1 (MQ=38)
                                                 aGACCGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCttt                                                     >  1:3782884/1‑101 (MQ=37)
                                                   actgAACCGCGATTGATGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTAt                                                   <  1:746292/98‑1 (MQ=25)
                                                    ccGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATc                                                  <  2:2841184/101‑1 (MQ=38)
                                                    ccGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATc                                                  <  1:3872431/101‑1 (MQ=38)
                                                    ccGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATc                                                  <  2:1194396/101‑1 (MQ=38)
                                                     cGAACCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCt                                                 <  1:1986370/101‑1 (MQ=38)
                                                      gAACCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCtt                                                >  1:1238085/1‑101 (MQ=37)
                                                       aaCCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCTTg                                               >  1:2829437/1‑101 (MQ=37)
                                                       aaCCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCTTg                                               >  1:918349/1‑101 (MQ=37)
                                                       aaCCGCGCGTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCTTg                                               <  1:3084701/101‑1 (MQ=38)
                                                          cgcgcgTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCTTGCCt                                            >  1:4143935/1‑101 (MQ=37)
                                                           gcgcgTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTc                                           >  2:3286620/1‑101 (MQ=37)
                                                            cgcgTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCa                                          <  1:811993/101‑1 (MQ=38)
                                                             gcgTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAg                                         <  1:4403668/101‑1 (MQ=38)
                                                              cgTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGa                                        >  2:1422770/1‑101 (MQ=37)
                                                              cgTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGa                                        >  2:2403080/1‑101 (MQ=37)
                                                               gTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGat                                       <  2:4240665/101‑1 (MQ=38)
                                                               gTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGat                                       <  2:3275346/101‑1 (MQ=38)
                                                               gTGCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGat                                       >  2:3353813/1‑101 (MQ=37)
                                                                 gCTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGatat                                     >  2:1572396/1‑101 (MQ=37)
                                                                  cTGAACTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATc                                    <  2:3553045/101‑1 (MQ=38)
                                                                     aaCTGATCCGGGAGGG.GGG.GGCTCTCC.CC.CCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCg                                 <  1:856495/101‑1 (MQ=25)
                                                                       cTGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCggg                               >  1:4561799/1‑101 (MQ=37)
                                                                        tGATACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCgggg                              <  1:2630122/101‑1 (MQ=37)
                                                                          aTACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGt                            <  1:58726/101‑1 (MQ=37)
                                                                           tACGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTg                           >  1:3266244/1‑101 (MQ=37)
                                                                            aCGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAACAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGGGa                          >  2:765278/1‑101 (MQ=25)
                                                                            aCGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGa                          >  2:2593619/1‑101 (MQ=25)
                                                                            aCGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGa                          >  2:4335012/1‑101 (MQ=25)
                                                                             cGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGaa                         <  2:2476524/101‑1 (MQ=25)
                                                                             cGGGA.GTGGGG.GGCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGaa                         >  2:314293/1‑101 (MQ=25)
                                                                                ggggggggg.ggCTCTCC.CC.CCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAAtt                       <  1:3236462/95‑1 (MQ=25)
                                                                                     gggg.ggCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGgc                  <  1:1315293/100‑1 (MQ=37)
                                                                                     gggg.ggCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGgc                  <  1:2052997/100‑1 (MQ=38)
                                                                                       gg.ggCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGgcgc                >  1:2133330/1‑101 (MQ=37)
                                                                                        g.ggCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCa               <  1:3406655/101‑1 (MQ=38)
                                                                                          ggCTCT.GCCC.CCCGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAg              >  1:1941671/1‑101 (MQ=37)
                                                                                               ctgccc.cccGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCaa          <  2:3636428/97‑1 (MQ=255)
                                                                                                tgccc.cccGATTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCTTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCaaa         <  1:2903508/98‑1 (MQ=37)
                                                                                                 gccc.cccGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCaaaa        <  2:4158121/99‑1 (MQ=255)
                                                                                                       ccGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAggg   <  1:2378180/101‑1 (MQ=255)
                                                                                                        cGCTTCGCTCCCCCCGGGATATTTGGGGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAgggg  >  1:1844953/1‑101 (MQ=255)
                                                                                                                                                                                                             
CCCGACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGA.GT.GGG.GGCTCT.G.CC.CCCGCTTCGCTCCCCCCGGGATATTT.GGGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGG  >  NC_014034/3608987‑3609184

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