Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 3,609,128 C→G 1.2% intergenic (+67/‑104) clpB → / → dinB chaperone ClpB/DNA polymerase IV

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140343,609,1280C→G 1.2% 2.0 ‑1.0 ‑0.2420intergenic (+67/‑104)clpB/dinBchaperone ClpB/DNA polymerase IV
Reads supporting (aligned to +/- strand):  new base (4/1):  ref base (213/202):  total (217/203)
Fisher's exact test strand distribution p-value = 3.73e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.61e-01

Only 100 of 441 total aligned reads displayed.

CGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGT.GGGGGCTCTG.CCCCC..GCTTCGCTCCCCCCGGGATATTT.G.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGA.GGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGTTAGAATTGAT  >  NC_014034/3609031‑3609227
                                                                                                       |                                                                                                    
  tCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGccc                                                                                                    >  2:5970030/1‑101 (MQ=37)
  tCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGccc                                                                                                    <  2:1557993/101‑1 (MQ=38)
  tCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGccc                                                                                                    >  2:5130165/1‑101 (MQ=37)
    cAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCtt                                                                                                  <  1:3692659/101‑1 (MQ=38)
    cAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCtt                                                                                                  >  1:3453518/1‑101 (MQ=37)
     tgACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCttt                                                                                                 <  1:3876800/100‑1 (MQ=38)
     aGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCttt                                                                                                 <  2:122037/101‑1 (MQ=38)
     aGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCttt                                                                                                 <  2:573471/101‑1 (MQ=38)
     aGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCttt                                                                                                 <  2:459906/101‑1 (MQ=38)
     aGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCttt                                                                                                 <  1:4338362/101‑1 (MQ=38)
     aGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCttt                                                                                                 <  1:4134385/101‑1 (MQ=38)
      gACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTa                                                                                                <  2:538812/101‑1 (MQ=38)
      gACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTa                                                                                                >  2:6580168/1‑101 (MQ=37)
          gAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATGGCCCTTTATCtt                                                                                            <  1:5124517/101‑1 (MQ=37)
           aaCCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTg                                                                                           <  1:1236688/101‑1 (MQ=38)
            aCCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGc                                                                                          <  2:1733129/101‑1 (MQ=38)
             ccGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGcc                                                                                         <  2:2954925/101‑1 (MQ=38)
              cgcgcgTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCt                                                                                        <  2:4076622/101‑1 (MQ=38)
                 gcgTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAg                                                                                     <  1:4034858/101‑1 (MQ=38)
                   gTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGat                                                                                   <  2:6312089/101‑1 (MQ=38)
                    tGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGata                                                                                  <  1:4462566/101‑1 (MQ=38)
                    tGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGata                                                                                  <  1:4960106/101‑1 (MQ=38)
                    tGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGata                                                                                  <  2:5335168/101‑1 (MQ=38)
                    tGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGata                                                                                  <  2:2164628/101‑1 (MQ=38)
                     gCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGatat                                                                                 <  2:5489854/101‑1 (MQ=38)
                     gCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGatat                                                                                 <  1:5677309/101‑1 (MQ=37)
                       tGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGAATTGCCCTTTTTCTGTCCTCCGATATcc                                                                               >  2:6948141/1‑101 (MQ=25)
                        gAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATccc                                                                              <  1:5944782/101‑1 (MQ=37)
                        gAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATccc                                                                              <  2:4203396/101‑1 (MQ=37)
                        gAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATccc                                                                              >  2:6163919/1‑101 (MQ=37)
                         aaCTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCg                                                                             <  2:63870/101‑1 (MQ=37)
                            tGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCgggg                                                                          >  1:21757/1‑101 (MQ=37)
                              aTACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGt                                                                        <  2:2821528/101‑1 (MQ=38)
                                aCGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGa                                                                      >  1:6909269/1‑101 (MQ=25)
                                 cGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGaa                                                                     >  2:1058802/1‑101 (MQ=25)
                                 cGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGaa                                                                     >  1:5584373/1‑101 (MQ=25)
                                       tggggggCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTgg                                                               >  1:3164575/3‑101 (MQ=37)
                                          gggtctgtgcccccc..GCTTCGCTCCCCCCGTGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGTATTTTGgcgc                                                            <  2:681704/91‑1 (MQ=37)
                                          ggggCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGgcgc                                                            >  1:5884698/1‑101 (MQ=37)
                                             gCTCTG.CCCCCCGGCTTCGCTCCCCCCGGGATATTT.GCGGCAAGATGAAGGGGCATGATTTCATCTTGCCCGAAATATCCCAGGGGTGAAttggccaaagg                                                          >  2:3115225/1‑92 (MQ=2)
                                              ctctGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGcc                                                        >  2:2652525/1‑101 (MQ=37)
                                               tctgcccccc..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCa                                                       >  1:609605/6‑101 (MQ=255)
                                                ctgcccccc..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCaa                                                      <  1:1417652/97‑1 (MQ=255)
                                                ctgcccccc..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCaa                                                      >  1:326681/5‑101 (MQ=255)
                                                   cccccc..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAag                                                   >  1:4197969/2‑101 (MQ=255)
                                                    ccccc..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAaga                                                  >  2:2237457/1‑101 (MQ=255)
                                                      ccc..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGA.gg                                               <  1:1880583/101‑1 (MQ=255)
                                                           gCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGA.ggggg                                            <  1:7080250/101‑2 (MQ=255)
                                                           gCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGA.ggggg                                            >  1:3419234/1‑100 (MQ=255)
                                                            cTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGA.gggggc                                           <  2:3437945/101‑3 (MQ=255)
                                                             ttCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGA.gggggca                                          <  2:3700948/101‑4 (MQ=255)
                                                               cGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGA.gggggcagc                                        >  1:1969760/1‑96 (MQ=255)
                                                                gCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAgcg                                        <  1:482084/101‑1 (MQ=25)
                                                                  tCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGcc                                      <  1:6231583/101‑1 (MQ=21)
                                                                   ccccccGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGccc                                     >  2:6720381/1‑101 (MQ=21)
                                                                     ccccGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGccccc                                   >  1:2439062/1‑101 (MQ=21)
                                                                     ccccGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGccccc                                   <  2:208535/101‑1 (MQ=21)
                                                                       ccGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTc                                 <  1:1675362/101‑1 (MQ=21)
                                                                       ccGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTc                                 >  2:1424650/1‑101 (MQ=21)
                                                                         gggATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTccc                               >  1:1520467/1‑101 (MQ=21)
                                                                           gATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTccccc                             >  1:2539407/1‑101 (MQ=21)
                                                                            atatTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCa                            >  1:3086442/1‑101 (MQ=21)
                                                                            atatTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCa                            >  1:1054010/1‑101 (MQ=21)
                                                                             tatTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCAc                           <  2:437404/101‑1 (MQ=18)
                                                                              atttgg.ggCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTTAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACg                          <  1:725672/96‑1 (MQ=21)
                                                                              atttgg.ggCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACg                          >  2:2459182/6‑101 (MQ=255)
                                                                              atttgg.ggCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACg                          >  1:3125631/6‑101 (MQ=255)
                                                                              atttgg.ggCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACg                          >  1:1127044/6‑101 (MQ=255)
                                                                               tttgg.ggCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACgg                         <  1:544611/97‑1 (MQ=255)
                                                                                 tgg.ggCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGcc                       >  1:6505339/3‑101 (MQ=255)
                                                                                      gCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCgggc                   >  1:2780733/1‑100 (MQ=255)
                                                                                       cAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTc                  >  2:880871/1‑101 (MQ=255)
                                                                                         aGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCt                >  1:1693383/1‑101 (MQ=255)
                                                                                         aGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCt                <  1:6638669/101‑1 (MQ=255)
                                                                                          ggTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTg               <  1:3549636/101‑1 (MQ=255)
                                                                                          ggTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTg               >  2:2710303/1‑101 (MQ=255)
                                                                                           gTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGc              <  1:1787832/101‑1 (MQ=255)
                                                                                             aaaaGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGt            >  1:2924968/1‑101 (MQ=255)
                                                                                              aaaGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGtt           >  2:3970555/1‑101 (MQ=255)
                                                                                              aaaGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGtt           <  2:2341937/101‑1 (MQ=255)
                                                                                              aaaGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGtt           >  2:3863565/1‑101 (MQ=255)
                                                                                              aaaGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGtt           <  2:325871/101‑1 (MQ=255)
                                                                                                aGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGTTAg         >  2:4909851/1‑101 (MQ=255)
                                                                                                aGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGTTAg         >  2:1890358/1‑101 (MQ=255)
                                                                                                 gCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGTTAGa        >  1:2672149/1‑101 (MQ=255)
                                                                                                 gCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGTTAGa        <  2:1335863/101‑1 (MQ=255)
                                                                                                 gCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGTTAGa        <  2:5641084/101‑1 (MQ=255)
                                                                                                  cATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTCCGTTAGat       >  2:5835734/1‑100 (MQ=255)
                                                                                                  cATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGGCCTGCGTTAGaa       <  1:6352093/101‑1 (MQ=255)
                                                                                                   aTTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGTTAGAAt      >  1:2328227/1‑101 (MQ=255)
                                                                                                   aTTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGTTAGAAt      <  2:6496053/101‑1 (MQ=255)
                                                                                                   aTTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGTTAGAAt      >  1:853153/1‑101 (MQ=255)
                                                                                                   aTTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGTTAGAAt      >  2:2028104/1‑101 (MQ=255)
                                                                                                   aTTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGTTAGAAt      <  2:3352453/101‑1 (MQ=255)
                                                                                                     tGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGTCCGGGTCCTGCGGTAGAATTg    >  1:1271265/1‑101 (MQ=18)
                                                                                                     tGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGTTAGAATTg    <  1:5377978/101‑1 (MQ=255)
                                                                                                     tGCCCTTTATCTTGCCTCAGATAGCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGTTAGAATTg    <  2:7033611/101‑1 (MQ=255)
                                                                                                      gCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGTTAGAATTGa   >  2:4761344/1‑101 (MQ=255)
                                                                                                      gCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGTTAGAATTGa   >  2:2758044/1‑101 (MQ=255)
                                                                                                      gCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGTTAGAATGGa   >  1:5716939/1‑101 (MQ=255)
                                                                                                       |                                                                                                    
CGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGT.GGGGGCTCTG.CCCCC..GCTTCGCTCCCCCCGGGATATTT.G.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGA.GGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGTTAGAATTGAT  >  NC_014034/3609031‑3609227

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