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

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140343,609,1290C→A 1.2% 2.7 ‑1.0 ‑0.2420intergenic (+68/‑103)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 = 5.83e-01

Only 100 of 441 total aligned reads displayed.

GTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGT.GGGGGCTCTG.CC.CCCGCTTCGCTCCCCCCGGGATATTT.G.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGA.GGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGTTAGAATTGATC  >  NC_014034/3609032‑3609228
                                                                                                      |                                                                                                    
 tCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCC.CCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGccc                                                                                                     <  2:2189342/101‑1 (MQ=38)
 tCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCC.CCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGccc                                                                                                     >  2:5970030/1‑101 (MQ=37)
 tCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCC.CCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGccc                                                                                                     >  2:5130165/1‑101 (MQ=37)
 tCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCC.CCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTCccc                                                                                                     <  1:5597096/101‑1 (MQ=37)
  ccAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCC.CCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCt                                                                                                    >  2:4141161/1‑101 (MQ=37)
   cAGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCC.CCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCtt                                                                                                   <  2:1976320/101‑1 (MQ=38)
    aGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCC.CCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCttt                                                                                                  <  2:573471/101‑1 (MQ=38)
    aGACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCC.CCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCttt                                                                                                  <  2:212721/101‑1 (MQ=38)
     gACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCC.CCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTa                                                                                                 >  2:6580168/1‑101 (MQ=37)
     gACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCC.CCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTa                                                                                                 <  2:538812/101‑1 (MQ=38)
     gACCGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCC.CCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTa                                                                                                 <  1:3731959/101‑1 (MQ=38)
        cGAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCC.CCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCt                                                                                              <  1:6720885/101‑1 (MQ=38)
         gAACCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCC.CCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCtt                                                                                             >  1:3092540/1‑101 (MQ=37)
          aaCCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCC.CCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTg                                                                                            <  1:1236688/101‑1 (MQ=38)
          aaCCGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCC.CCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTg                                                                                            >  2:126943/1‑101 (MQ=37)
            ccGCGCGTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCC.CCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGcc                                                                                          <  2:3603023/101‑1 (MQ=38)
                gcgTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCC.CCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAg                                                                                      <  1:5960285/101‑1 (MQ=38)
                  gTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCC.CCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGat                                                                                    <  2:6312089/101‑1 (MQ=38)
                  gTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCC.CCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGat                                                                                    <  2:3907768/101‑1 (MQ=38)
                   tGCTGAACTGATACGGGAGTGGGGGGCTCTGCCC.CCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGata                                                                                   <  2:5335168/101‑1 (MQ=38)
                   tGCTGAACTGATACGGGAGTGGGGGGCTCTGCCC.CCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGata                                                                                   <  1:2928869/101‑1 (MQ=38)
                    gCTGAACTGATACGGGAGTGGGGGGCTCTGCCC.CCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGatat                                                                                  <  2:5489854/101‑1 (MQ=38)
                      tGAACTGATACGGGAGTGGGGGGCTCTGCCC.CCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGAATTGCCCTTTTTCTGTCCTCCGATATcc                                                                                >  2:6948141/1‑101 (MQ=25)
                       gAACTGATACGGGAGTGGGGGGCTCTGCCC.CCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATccc                                                                               >  2:6163919/1‑101 (MQ=37)
                        aaCTGATACGGGAGTGGGGGGCTCTGCCC.CCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCg                                                                              <  2:63870/101‑1 (MQ=37)
                        aaCTGATACGGGAGTGGGGGGCTCTGCCC.CCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCg                                                                              >  1:1883687/1‑101 (MQ=37)
                        aaCTGATACGGGAGTGGGGGGCTCTGCCC.CCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCg                                                                              >  2:3417805/1‑101 (MQ=37)
                         aCTGATACGGGAGTGGGGGGCTCTGCCC.CCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCgg                                                                             <  2:4276784/101‑1 (MQ=37)
                           tGATACGGGAGTGGGGGGCTCTGCCC.CCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCgggg                                                                           >  1:3190135/1‑101 (MQ=37)
                            gATACGGGAGTGGGGGGCTCTGCCC.CCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCggggg                                                                          >  2:1692590/1‑101 (MQ=37)
                               aCGGGAGTGGGGGGCTCTGCCC.CCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGa                                                                       <  1:4409267/101‑1 (MQ=37)
                               aCGGGAGTGGGGGGCTCTGCCC.CCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGa                                                                       <  1:5536689/101‑1 (MQ=37)
                                cGGGAGTGGGGGGCTCTGCCC.CCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGaa                                                                      <  1:1533653/101‑1 (MQ=37)
                                cGGGAGTGGGGGGCTCTGCCC.CCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGaa                                                                      <  2:4762946/101‑1 (MQ=37)
                                 gggAGTGGGGGGCTCTGCCC.CCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAAt                                                                     <  2:2372058/101‑1 (MQ=37)
                                 gggAGTGGGGGGCTCTGCCC.CCCGCTTCGCTCCACCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAAt                                                                     >  1:5021957/1‑101 (MQ=25)
                                 gggAGTGGGGGGCTCTG.CCGCCGGCTTCGCGCCCCCCGGGGGATTTGG.GGCAAGGTAAAAGCATTGCGCTTTAGCTTGCCTCAGAGATCCCgggggcgtgt                                                                     >  2:2254919/1‑96 (MQ=255)
                                  ggAGTGGGGGGCTCTGCCC.CCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAAtt                                                                    >  2:2825769/1‑101 (MQ=25)
                                   gagtggggggCTCTGCCC.CCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAAttt                                                                   <  1:6164238/96‑1 (MQ=37)
                                    agtggggggCTCTGCCC.CCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAAtttt                                                                  <  1:3431142/97‑1 (MQ=37)
                                     gtggggggCTCTGCCC.CCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTg                                                                 >  1:6441932/4‑101 (MQ=37)
                                      tggggggCTCTGCCC.CCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTgg                                                                >  1:947206/3‑101 (MQ=37)
                                        gggggCTCTGCCC.CCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCNGGGGGTGAATTTTGgcg                                                              >  2:2726294/1‑101 (MQ=37)
                                         gggtctgtgccc.cccGCTTCGCTCCCCCCGTGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGTATTTTGgcgc                                                             <  2:681704/91‑1 (MQ=37)
                                             ctctGCCC.CCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGcc                                                         <  2:2056284/101‑1 (MQ=37)
                                              tctgccc.cccGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCa                                                        >  1:1794536/6‑101 (MQ=255)
                                              tctgccc.cccGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCa                                                        >  1:2577277/6‑101 (MQ=255)
                                               ctgccc.cccGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCaa                                                       >  1:2001603/5‑101 (MQ=255)
                                                tgccc.cccGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCaaa                                                      <  2:3142860/98‑1 (MQ=255)
                                                 gccc.cccGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCaaaa                                                     <  2:3343235/99‑1 (MQ=255)
                                                 gccc.cccGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCaaaa                                                     <  2:4526553/99‑1 (MQ=255)
                                                  ccc.cccGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGAATTGCCCATTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAag                                                    <  1:2800169/100‑1 (MQ=37)
                                                    c.cccGCTTCGCTCCCCCCGGGATATTT.GAGGCGAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAaga.g                                                 >  1:6272088/1‑101 (MQ=255)
                                                        cGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGA.gggg                                              >  2:3989611/1‑101 (MQ=255)
                                                        cGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGA.gggg                                              >  1:1449694/1‑101 (MQ=255)
                                                          cTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGA.gggggc                                            >  1:3467647/1‑99 (MQ=255)
                                                            tCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGA.gggggcag                                          >  2:1072397/1‑97 (MQ=255)
                                                             cGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGA.gggggcagc                                         <  1:6948514/101‑6 (MQ=255)
                                                             cGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGA.gggggcagc                                         <  2:2577181/101‑6 (MQ=255)
                                                              gCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAgcg                                         >  2:2986459/1‑101 (MQ=25)
                                                              gCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAgcg                                         <  1:2330827/101‑1 (MQ=25)
                                                               cTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAgcgc                                        >  2:2803805/1‑101 (MQ=21)
                                                                tCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGcc                                       >  1:5421856/1‑101 (MQ=21)
                                                                 ccccccGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGccc                                      >  2:6720381/1‑101 (MQ=21)
                                                                 ccccccGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGccc                                      >  1:3132701/1‑101 (MQ=21)
                                                                  cccccGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGcccc                                     >  1:4097975/1‑101 (MQ=21)
                                                                   ccccGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGccccc                                    >  2:4921046/1‑101 (MQ=21)
                                                                      cGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTcc                                 >  1:3560701/1‑101 (MQ=21)
                                                                        ggATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTcccc                               >  1:753979/1‑101 (MQ=21)
                                                                         gATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTccccc                              <  1:4141356/101‑1 (MQ=18)
                                                                         gATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTccccc                              <  1:568616/101‑1 (MQ=18)
                                                                          atatTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCa                             >  1:1054010/1‑101 (MQ=21)
                                                                           tatTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCAc                            >  2:4662338/1‑101 (MQ=21)
                                                                            atttgg.ggCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACg                           >  1:1127044/6‑101 (MQ=255)
                                                                            atttgg.ggCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACg                           <  1:5941104/96‑1 (MQ=255)
                                                                            atttgg.ggCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACg                           >  1:2187288/6‑101 (MQ=255)
                                                                             tttgg.ggCAAGGTAACAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACgg                          <  1:3278816/97‑1 (MQ=18)
                                                                             tttgg.ggCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACgg                          >  1:6514744/5‑101 (MQ=255)
                                                                              ttgg.ggCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCTCCCACGGc                         >  2:3751944/4‑101 (MQ=18)
                                                                               tgg.ggCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGcc                        >  2:2667270/3‑101 (MQ=255)
                                                                               tgg.ggCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGcc                        >  1:2708988/3‑101 (MQ=255)
                                                                                 g.ggCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCgg                      >  1:5270108/1‑101 (MQ=255)
                                                                                     cAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTc                   >  2:880871/1‑101 (MQ=255)
                                                                                     cAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTc                   <  1:1677278/101‑1 (MQ=255)
                                                                                     cAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTc                   >  1:4202569/1‑101 (MQ=255)
                                                                                     cAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTc                   <  2:1782931/101‑1 (MQ=255)
                                                                                       aGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCt                 >  1:2468526/1‑101 (MQ=255)
                                                                                        ggTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTccgc                >  2:1540088/1‑99 (MQ=255)
                                                                                        ggTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTg                >  1:4646291/1‑101 (MQ=255)
                                                                                          tAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCg              >  1:4871270/1‑101 (MQ=255)
                                                                                           aaaaGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGt             >  1:3877919/1‑101 (MQ=255)
                                                                                            aaaGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGtt            >  2:1365348/1‑101 (MQ=255)
                                                                                             aaGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGTTa           >  1:7117077/1‑101 (MQ=255)
                                                                                               gCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGTTAGa         <  2:5641084/101‑1 (MQ=255)
                                                                                                cATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTCCGTTAGat        >  2:5835734/1‑100 (MQ=255)
                                                                                                 aTTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGTTAGAAt       >  1:5786794/1‑101 (MQ=255)
                                                                                                 aTTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGTTAGAAt       <  2:6496053/101‑1 (MQ=255)
                                                                                                   tGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGTCCGGGTCCTGCGGTAGAATTg     >  1:1271265/1‑101 (MQ=18)
                                                                                                   tGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGTTAGAATTg     <  1:1754669/101‑1 (MQ=255)
                                                                                                   tGCCCTTTATCTTGCCTCAGATAGCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGTTAGAATTg     <  2:7033611/101‑1 (MQ=255)
                                                                                                      |                                                                                                    
GTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAGT.GGGGGCTCTG.CC.CCCGCTTCGCTCCCCCCGGGATATTT.G.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGA.GGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGTTAGAATTGATC  >  NC_014034/3609032‑3609228

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 30 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG