Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 3,609,108 +G 98.3% intergenic (+47/‑124) clpB → / → dinB chaperone ClpB/DNA polymerase IV

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140343,609,1081.→G98.3% 8.0 ‑2.4 ‑0.7358intergenic (+47/‑124)clpB/dinBchaperone ClpB/DNA polymerase IV
Reads supporting (aligned to +/- strand):  new base (150/200):  ref base (5/1):  total (155/203)
Fisher's exact test strand distribution p-value = 8.96e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.24e-01

Only 100 of 393 total aligned reads displayed.

CAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGT.GGGGGCTCTG.CCCCC..GCTTCGCTCCCCCCGGGATATTT.G.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGA.GGGGCAGCGCCCCCTCCCCCAC  >  NC_014034/3609010‑3609202
                                                                                                       |                                                                                                
 aaGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTT.g                                                                                                 >  2:1787768/1‑101 (MQ=39)
  aGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTT.g.g                                                                                               <  1:3431926/101‑1 (MQ=39)
     tctcGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTT.g.gggc                                                                                            <  2:4896584/101‑3 (MQ=38)
      ctcGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTT.g.gggca                                                                                           >  1:6860614/1‑98 (MQ=38)
       tcGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTT.g.gggcaa                                                                                          <  1:641370/101‑5 (MQ=38)
       tcGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTT.g.gggcaa                                                                                          <  2:5536852/101‑5 (MQ=38)
           ccagccagCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTa                                                                                       <  2:5741225/101‑1 (MQ=38)
           ccagccagCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTa                                                                                       <  2:894734/101‑1 (MQ=38)
             agccagcGCGTCCAGACCGAACCGCGCGTGCTGAACTGCTACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTAG.GGCAAGGTaaa                                                                                     <  1:5093489/101‑1 (MQ=37)
             agccagcGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTaaa                                                                                     <  2:2883036/101‑1 (MQ=38)
             agccagcGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTaaa                                                                                     <  2:3093726/101‑1 (MQ=38)
              gccagcGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTaaaa                                                                                    <  1:807040/101‑1 (MQ=38)
                cagcGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGc                                                                                  <  1:6058371/101‑1 (MQ=38)
                 agcGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCa                                                                                 <  2:5015439/101‑1 (MQ=38)
                   cgcgTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCAtt                                                                               <  2:6229631/101‑1 (MQ=38)
                       tCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGccc                                                                           <  2:2189342/101‑1 (MQ=38)
                       tCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGccc                                                                           >  2:296226/1‑101 (MQ=37)
                       tCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGccc                                                                           <  2:1557993/101‑1 (MQ=38)
                       tCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGccc                                                                           >  1:5650200/1‑101 (MQ=37)
                        ccAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCt                                                                          <  2:6776465/101‑1 (MQ=38)
                         cAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCtt                                                                         <  2:1976320/101‑1 (MQ=38)
                          aGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCttt                                                                        <  2:459906/101‑1 (MQ=38)
                          aGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCttt                                                                        <  2:5364560/101‑1 (MQ=38)
                          aGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCttt                                                                        <  2:212721/101‑1 (MQ=38)
                           gACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTa                                                                       >  2:6580168/1‑101 (MQ=37)
                             ccGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATc                                                                     <  1:1885919/101‑1 (MQ=38)
                                aaCCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTg                                                                  >  2:126943/1‑101 (MQ=37)
                                aaCCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGAATTGCCCTTTATCTTg                                                                  >  2:2800067/1‑101 (MQ=37)
                                 aCCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGc                                                                 <  1:1382435/101‑1 (MQ=38)
                                 aCCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGc                                                                 <  1:110556/101‑1 (MQ=38)
                                       cgTGCTGAACTGATACGGGAGTGGGGGGCTTTGCCCCCC..GCTTCGCTCCCCCCCGGATATTTGG.GGCAAGGGAAAAGCGTTGCCCTTTATCTTGCCTCTGa                                                           >  1:3376358/1‑101 (MQ=255)
                                       cgTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGa                                                           >  1:1496114/1‑101 (MQ=37)
                                        gTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGat                                                          <  2:3907768/101‑1 (MQ=38)
                                        gTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGat                                                          <  2:6312089/101‑1 (MQ=38)
                                         tGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGata                                                         <  1:4960106/101‑1 (MQ=38)
                                           cTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATc                                                       <  1:5944074/101‑1 (MQ=38)
                                            tGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGAATTGCCCTTTTTCTGTCCTCCGATATcc                                                      >  2:6948141/1‑101 (MQ=25)
                                             gAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATccc                                                     >  2:1417602/1‑101 (MQ=37)
                                             gAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATccc                                                     >  2:5651371/1‑101 (MQ=37)
                                             gAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATccc                                                     >  1:5338381/1‑101 (MQ=37)
                                              aaCTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCg                                                    <  2:63870/101‑1 (MQ=37)
                                              aaCTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCg                                                    <  2:6490113/101‑1 (MQ=37)
                                                 tGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCgggg                                                 >  1:3190135/1‑101 (MQ=37)
                                                 tGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCgggg                                                 >  1:21757/1‑101 (MQ=37)
                                                 tGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCgggg                                                 >  1:71250/1‑101 (MQ=37)
                                                  gATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCggggg                                                >  1:2342032/1‑101 (MQ=37)
                                                  gATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCggggg                                                >  2:1692590/1‑101 (MQ=37)
                                                  gATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCGTGCCTCAGATATCCCggggg                                                >  1:4369180/1‑101 (MQ=37)
                                                   aTACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGt                                               <  2:2714473/101‑1 (MQ=37)
                                                    tACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTg                                              >  2:3349670/1‑101 (MQ=37)
                                                    tACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTg                                              <  1:3756773/101‑1 (MQ=37)
                                                     aCGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCCGGATATTTGG.GGCAAGGTAAACGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGa                                             >  1:299732/1‑101 (MQ=37)
                                                      cGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGaa                                            <  2:4762946/101‑1 (MQ=37)
                                                      cGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGaa                                            >  2:1058802/1‑101 (MQ=25)
                                                       gggAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAAt                                           >  2:6674003/1‑101 (MQ=25)
                                                         gagtggggggCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAAttt                                         <  1:6164238/96‑1 (MQ=37)
                                                          agtggggggCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAAtttt                                        <  2:7025453/97‑1 (MQ=37)
                                                           gtggggggCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTg                                       >  1:4764543/4‑101 (MQ=37)
                                                           gtggggggCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTg                                       <  1:3561066/98‑1 (MQ=37)
                                                              gggggCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGgcg                                    <  1:3138438/101‑1 (MQ=37)
                                                               gggtctgtgcccccc..GCTTCGCTCCCCCCGTGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGTATTTTGgcgc                                   <  2:681704/91‑1 (MQ=37)
                                                                gggCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCa                                  >  2:4105621/1‑101 (MQ=37)
                                                                 ggCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAg                                 <  1:6715915/101‑1 (MQ=37)
                                                                 ggCTCTG.CCCCCCGGCTTCGCTCCCCCCGGGATATTT.GCGGCAAGATGAAGGGGCATGATTTCATCTTGCCCGAAATATCCCAGGGGTGAAttggccaaag                                  >  2:6064346/1‑93 (MQ=2)
                                                                  gCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGc                                <  1:2885509/101‑1 (MQ=37)
                                                                  gCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGc                                <  1:6537520/101‑1 (MQ=37)
                                                                   ctctGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGcc                               <  2:2056284/101‑1 (MQ=37)
                                                                    tctgcccccc..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCa                              >  1:1794536/6‑101 (MQ=255)
                                                                    tctgcccccc..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCa                              <  2:5169106/96‑1 (MQ=255)
                                                                     ctgcccccc..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCaa                             <  2:5933695/97‑1 (MQ=255)
                                                                     ctgcccccc..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCaa                             >  1:326681/5‑101 (MQ=255)
                                                                     ctgcccccc..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCaa                             >  1:5575702/5‑101 (MQ=255)
                                                                      tgcccccc..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCaaa                            <  2:3142860/98‑1 (MQ=255)
                                                                      tgcccccc..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCaaa                            >  2:5812703/4‑101 (MQ=255)
                                                                       gcccccc..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCaaaa                           >  1:5702382/3‑101 (MQ=255)
                                                                       gcccccc..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCaaaa                           <  1:6274761/99‑1 (MQ=255)
                                                                       gcccccc..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCaaaa                           >  1:5443994/3‑101 (MQ=255)
                                                                       gcccccc..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCaaaa                           <  1:4273671/99‑1 (MQ=255)
                                                                        cccccc..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAag                          <  2:3998698/100‑1 (MQ=255)
                                                                         ccccc..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAaga                         >  2:2237457/1‑101 (MQ=255)
                                                                         ccccc..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAaga                         <  1:2480788/101‑1 (MQ=255)
                                                                          cccc..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAaga.g                       <  1:271870/101‑1 (MQ=255)
                                                                             c..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGCAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGAGAATTTTGCCGCAGCCAAAAGC.gggg                    >  2:2653384/1‑101 (MQ=37)
                                                                                 cTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGA.gggggc                  <  2:3437945/101‑3 (MQ=255)
                                                                                 cTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGA.gggggc                  >  1:3467647/1‑99 (MQ=255)
                                                                                  ttCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGA.gggggca                 <  2:3700948/101‑4 (MQ=255)
                                                                                    cGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGA.gggggcagc               <  2:2399405/101‑6 (MQ=255)
                                                                                    cGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGA.gggggcagc               <  1:6948514/101‑6 (MQ=255)
                                                                                        ccccccGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGTGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGccc            >  1:4600557/1‑101 (MQ=21)
                                                                                         cccccGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGcccc           >  1:4097975/1‑101 (MQ=21)
                                                                                          ccccGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGccccc          <  2:5414077/101‑1 (MQ=21)
                                                                                          ccccGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGccccc          <  1:126944/101‑1 (MQ=21)
                                                                                            ccGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTc        <  1:1675362/101‑1 (MQ=21)
                                                                                              gggATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTccc      >  1:3932691/1‑101 (MQ=21)
                                                                                                gATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTccccc    <  1:3123230/101‑1 (MQ=18)
                                                                                                gATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTccccc    >  2:4208961/1‑101 (MQ=21)
                                                                                                gATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTccccc    <  1:4141356/101‑1 (MQ=18)
                                                                                                 atatTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCa   <  2:5495741/101‑1 (MQ=18)
                                                                                                 atatTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCa   >  1:1054010/1‑101 (MQ=21)
                                                                                                  tatTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCAc  <  2:437404/101‑1 (MQ=18)
                                                                                                       |                                                                                                
CAAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGT.GGGGGCTCTG.CCCCC..GCTTCGCTCCCCCCGGGATATTT.G.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGA.GGGGCAGCGCCCCCTCCCCCAC  >  NC_014034/3609010‑3609202

Base quality scores:  ATCG  < 3 ≤  ATCG  < 24 ≤  ATCG  < 31 ≤  ATCG  < 35 ≤  ATCG  < 41 ≤  ATCG