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.5 ‑1.2 ‑0.0202intergenic (+19/‑152)clpB/dinBchaperone ClpB/DNA polymerase IV
Reads supporting (aligned to +/- strand):  new base (76/122):  ref base (3/1):  total (79/123)
Fisher's exact test strand distribution p-value = 3.02e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.78e-01

Only 100 of 220 total aligned reads displayed.

GCCGGTCCCGACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGT.GGG.GGCTCT.G.CC.CCC..GCTTCGCTCCCCCCGGGATATTT.G.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAA  >  NC_014034/3608981‑3609176
                                                                                                       |                                                                                                     
gccgGTCCCGACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGG.GGCTCT.g                                                                                                        >  1:2659091/1‑101 (MQ=255)
      cccGACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGG.GGCTCT.G.cc.cccc                                                                                                <  1:4410060/101‑2 (MQ=255)
       ccGACGGGCTGGTGGTCGGCGGCAAGGTCTCGCCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGG.GGCTCT.G.cc.ccccg                                                                                               <  1:1781415/101‑3 (MQ=255)
        cGACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGG.GGCTCT.G.cc.ccccgc                                                                                              <  2:837553/101‑4 (MQ=255)
         gACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGG.GGCTCT.G.cc.ccccgct                                                                                             <  1:4583727/101‑5 (MQ=255)
          aCGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGG.GGCTCT.G.cc.ccccgctt                                                                                            >  2:216057/1‑96 (MQ=255)
           cGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTc                                                                                          <  1:1007139/101‑1 (MQ=255)
              gCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCt                                                                                       <  2:1732452/101‑1 (MQ=255)
               cTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTc                                                                                      <  1:4476786/101‑1 (MQ=255)
               cTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCNTCGCTc                                                                                      >  1:1832935/1‑101 (MQ=255)
               cTGGTGGTCGGCGGCAAGGTCTCGCCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTc                                                                                      <  2:4316253/101‑1 (MQ=255)
                tggtggTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTcc                                                                                     <  1:799465/101‑1 (MQ=255)
                  gtggtCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTcccc                                                                                   >  2:3940711/1‑101 (MQ=255)
                    ggtCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTcccccc                                                                                 <  2:2220079/101‑1 (MQ=255)
                    ggtCGGCGCCAAGGCCTCGGCCAGCCAGCGCGTTCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTcccccc                                                                                 <  1:4025101/101‑1 (MQ=255)
                     gtCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCg                                                                                >  1:68618/1‑101 (MQ=255)
                     gtCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCg                                                                                >  2:1647254/1‑101 (MQ=255)
                      ggggcggcAAGGTCTGGGGCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATGCGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCgg                                                                               <  2:372267/99‑1 (MQ=255)
                       cggcggCAAGGTCTCGGCCAGCCAGCGCGTCCAGCCCGAACCGCGCGTGCTGAACTGATACGGGAGT.GGGNNNNNNT.G.CCNCNC..CCNTCGCTCCCCCCggg                                                                              >  1:3614007/1‑101 (MQ=255)
                       cggcggCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCggg                                                                              <  1:3463148/101‑1 (MQ=255)
                        ggcggcAGGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGa                                                                             <  2:2377244/101‑1 (MQ=255)
                        ggcggcAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGa                                                                             <  1:1903825/101‑1 (MQ=255)
                          cggcAGGTTCCCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGata                                                                           <  2:1996225/101‑1 (MQ=255)
                            gcAAGGTCTCGCCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATAtt                                                                         <  1:4220030/101‑1 (MQ=38)
                              aaGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTT.g                                                                      >  2:58720/1‑101 (MQ=39)
                                ggTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTT.g.gg                                                                   >  2:1240711/1‑101 (MQ=39)
                                  tctcGTCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTT.g.gggc                                                                 <  1:2872786/101‑3 (MQ=38)
                                  tctcGGCCAGCCAGCGCGTCCATACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTT.g.gggc                                                                 <  2:538337/101‑3 (MQ=38)
                                  tctcGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTT.g.gggc                                                                 >  2:2523327/1‑99 (MQ=38)
                                   ctcGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTT.g.gggca                                                                >  2:297648/1‑98 (MQ=38)
                                   ctcGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTT.g.gggca                                                                >  2:383148/1‑98 (MQ=38)
                                   ctcGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCCGGATATTT.g.gggca                                                                >  1:1750346/1‑98 (MQ=255)
                                    tcGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTT.g.gggcaa                                                               >  2:2162669/1‑97 (MQ=38)
                                      ggCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAgg                                                              <  2:4477113/101‑1 (MQ=38)
                                       gccagccaGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGt                                                             >  1:293478/1‑101 (MQ=38)
                                        ccagccagCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTa                                                            <  1:818751/101‑1 (MQ=38)
                                        ccagccagCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTa                                                            >  2:3098128/1‑101 (MQ=38)
                                          agccagcGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTaaa                                                          <  2:1897923/101‑1 (MQ=38)
                                           gccagcGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTaaaa                                                         <  2:3350307/101‑1 (MQ=38)
                                            ccagcGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAg                                                        <  2:2291976/101‑1 (MQ=38)
                                            ccagcGCGTCCAGACCGAACCGCGCGCGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGACATTTGG.GGCAAGGTAAAAg                                                        >  1:2476690/1‑101 (MQ=37)
                                             cagcGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGc                                                       <  2:1135134/101‑1 (MQ=38)
                                             cagcGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGc                                                       >  2:970461/1‑101 (MQ=37)
                                             cagcGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGc                                                       >  1:1178691/1‑101 (MQ=37)
                                               gcgcGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCAt                                                     <  2:4157545/101‑1 (MQ=38)
                                                ggcgTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCAtt                                                    <  1:3427233/100‑1 (MQ=37)
                                                cgcgTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCAtt                                                    <  1:922991/101‑1 (MQ=38)
                                                  cgTCCGGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGc                                                  <  1:2180538/101‑1 (MQ=37)
                                                   gTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGcc                                                 <  2:725227/101‑1 (MQ=38)
                                                   gTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGcc                                                 >  2:3406401/1‑101 (MQ=37)
                                                     ccAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCt                                               <  2:3358788/101‑1 (MQ=38)
                                                     ccAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCt                                               <  1:2322047/101‑1 (MQ=38)
                                                      cAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCtt                                              <  2:2871115/101‑1 (MQ=38)
                                                      cAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCtt                                              <  1:4541748/101‑1 (MQ=38)
                                                       aGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCttt                                             >  1:3782884/1‑101 (MQ=37)
                                                         aCCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTAt                                           <  1:1460277/101‑1 (MQ=38)
                                                          ccGAACCGCGCGTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATc                                          <  1:3872431/101‑1 (MQ=38)
                                                          ccGAACCGCGCGTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATc                                          <  2:1194396/101‑1 (MQ=38)
                                                           cGAACCGCGCGTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCt                                         >  1:489218/1‑101 (MQ=37)
                                                           cGAACCGCGCGTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCt                                         <  1:1986370/101‑1 (MQ=38)
                                                            gAACCGCGCGTTCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCtt                                        <  1:3659215/101‑1 (MQ=37)
                                                            gAACCGCGCGTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCtt                                        >  1:1238085/1‑101 (MQ=37)
                                                            gAACCGCGCGTGCTGAACTGATACAGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCtt                                        <  2:2638405/101‑1 (MQ=37)
                                                             aaCCGCGCGTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTg                                       <  1:3084701/101‑1 (MQ=38)
                                                             aaCCGCGCGTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTg                                       >  1:2829437/1‑101 (MQ=37)
                                                                cgcgcgTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCt                                    <  2:1380301/101‑1 (MQ=38)
                                                                cgcgcgTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCt                                    >  1:4143935/1‑101 (MQ=37)
                                                                  cgcgTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCa                                  <  2:684306/101‑1 (MQ=38)
                                                                   gcgTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCat                                 <  1:2984505/101‑2 (MQ=38)
                                                                    cgTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGa                                >  2:1422770/1‑101 (MQ=37)
                                                                    cgTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGa                                >  2:2403080/1‑101 (MQ=37)
                                                                     gTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGat                               <  2:4240665/101‑1 (MQ=38)
                                                                     gTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGat                               <  2:268744/101‑1 (MQ=38)
                                                                     gTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGat                               <  2:4661249/101‑1 (MQ=38)
                                                                     gTGCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGat                               <  2:3626616/101‑1 (MQ=38)
                                                                       gCTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGatat                             >  2:1572396/1‑101 (MQ=37)
                                                                        cTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATc                            <  2:3553045/101‑1 (MQ=38)
                                                                        cTGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATc                            <  2:3281427/101‑1 (MQ=38)
                                                                         tGAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATcc                           >  2:4366117/1‑101 (MQ=37)
                                                                          gAACTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATccc                          <  2:2818562/101‑1 (MQ=37)
                                                                             cTGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCggg                       >  1:4561799/1‑101 (MQ=37)
                                                                              tGATACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCgggg                      <  1:2630122/101‑1 (MQ=37)
                                                                                aTACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGt                    <  1:58726/101‑1 (MQ=37)
                                                                                 tACGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTg                   >  1:3266244/1‑101 (MQ=37)
                                                                                  aCGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAACAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGGGa                  >  2:765278/1‑101 (MQ=25)
                                                                                  aCGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGa                  >  2:2593619/1‑101 (MQ=25)
                                                                                   cGGGAGTGGGG.GGCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAGGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGaa                 >  2:3153688/1‑101 (MQ=25)
                                                                                     ggggggggg.ggCTCTCC.CC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAAtt               <  1:3236462/95‑1 (MQ=25)
                                                                                       agtgggg.ggCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAAtttt             <  2:2473939/97‑1 (MQ=37)
                                                                                        gtgggg.ggCTCT.GCCC.CCC..GCTTCGCTCCCCCCGTGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTg            <  2:495864/98‑1 (MQ=37)
                                                                                        gtgggg.ggCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTg            >  2:4264269/4‑101 (MQ=37)
                                                                                          gggg.ggCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGgc          <  1:2052997/100‑1 (MQ=38)
                                                                                          gggg.ggCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGgc          <  1:1315293/100‑1 (MQ=37)
                                                                                            gg.ggCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGgcgc        >  1:2133330/1‑101 (MQ=37)
                                                                                            gg.ggCTCT.G.CC.CCCCGGCTTCGCTCCCCCCGGGATATTT.GCGGCAAGATGAAGGGGCATGATTTCATCTTGCCCGAAATATCCCAGGGGTGAAttggccaa         >  2:3794856/1‑95 (MQ=2)
                                                                                            gg.ggCTCG.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGgcgc        <  1:917776/101‑1 (MQ=37)
                                                                                             g.ggCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCa       <  1:3406655/101‑1 (MQ=38)
                                                                                             g.ggCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCa       <  2:2149278/101‑1 (MQ=38)
                                                                                               ggCTCT.GCCC.CCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAg      >  1:1941671/1‑101 (MQ=37)
                                                                                                gCTCT.G.CC.CCCCGGCTTCGCTCCCCCCGGGATATTT.GCGGCAAGATGAAGGGGCATGATTTCATCTTGCCCGAAATATCCCAGGGGTGAAttggccaaagg      >  2:2840926/1‑92 (MQ=2)
                                                                                                       |                                                                                                     
GCCGGTCCCGACGGGCTGGTGGTCGGCGGCAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGT.GGG.GGCTCT.G.CC.CCC..GCTTCGCTCCCCCCGGGATATTT.G.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAA  >  NC_014034/3608981‑3609176

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