Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,146,389 (T)6→7 coding (132/879 nt) pflC → putative pyruvate formate‑lyase 2 activating enzyme PflC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,146,3891.T93.8% 42.2 / ‑3.0 16coding (132/879 nt)pflCputative pyruvate formate‑lyase 2 activating enzyme PflC
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base T (9/6);  total (10/6)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.65e-01

GCGCAACAGCGCCCACAAGATTATGCCGGGCTGGTGGTGCGCGTTGCCGGATACAGCGCCTTCTTTGTCGAACTGTCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CTTTTTT‑AAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCG  >  NC_000913/4146167‑4146544
                                                                                                                                                                                                                                |                                                                                                                                                           
gcgcAACAGCGCCCACAAGATTATGCCGGGCTGGTGGTGCGCGTTGCCGGATACAGCGCCTTCTTTGTCGAACTGTCGAAGGAGATCCAGGATGACATCCTCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTCTTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CTTTTTTTAAAGGCTGTCCGCATCTTTGCCCgtgg                                                                                                                                  >  1:197514/1‑251 (MQ=255)
gcgcAACAGCGCCCACAAGATTATGCCGGGCTGGTGGTGCGCGTTGCCGGATACAGCGCCTTCTTTGTCGAACTGTCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CTTTTTTTAAAGGCTGTCCGCATCTTTGCCCgtgg                                                                                                                                  >  1:119923/1‑251 (MQ=255)
                                ggtggtGCGCGTTGCCGGATACAGCGCCTTCTTTGTCGCACTGTCGAAGGAGATCCAGGATGACATCATCCGCCGGCCAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCGTTTTCAACATTCAGCGTTATTCACTGACTGACGGCGAGGGCATTCGTACGGTGGT‑CTTTTTTTACAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAAt                                                                                                  >  1:121738/1‑251 (MQ=255)
                                   ggtgCGCGTTGCCGGATACAGCGCCTTCTTTGTCGAACTGTCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CTTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGc                                                                                                                             >  1:82641/1‑221 (MQ=255)
                                                ggATACAGCGCCTTCTTTGTCGAACTGTCGAAGGAGATCCAGGATGACATCATCCGCCGGAAAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGCCGGCGAGGGCATTCGTCCGGTGGT‑CTTTTTTTAAAGGc                                                                                                                                                       >  1:168759/1‑182 (MQ=255)
                                                  aTACAGCGCCTTCTTTGTCGAACTGTCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTAGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CTTTTTTTAAAGGCTGTCCGCATCTTTg                                                                                                                                         >  1:107882/1‑194 (MQ=255)
                                                    aCAGCGCCTTCTTTGTCGAACTGTCGAAGGAGATCCAGGATGACATCATCGGCCGGAAAGCGAATCAGCTGTAACGTTGCGGAAACGCGCCGCGATGATCTGGAGCGCATTTTCAACATTCAGCGATATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CTTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAgg                                                                              >  1:225510/1‑251 (MQ=255)
                                                     cAGCGCCTTCTTTGTCGAACTGTCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CTTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCg                                                                                                         <  1:87771/223‑1 (MQ=255)
                                                                    cGAACTGTCGAAGGAGATCCAGGATGACATCATCCGACGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGTCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTCCGGTGGT‑CTTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGCCGGTACGCAGAGAGGCGAAATGTCTGCACtg                                                              <  1:14224/251‑1 (MQ=255)
                                                                       aCTGTCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CTTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGc                                                           <  1:108403/251‑1 (MQ=255)
                                                                          ntCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CTTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGa                                                         >  1:153574/2‑250 (MQ=255)
                                                                                                       gccgGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CTTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAAt                                                                                                  >  1:111140/1‑180 (MQ=255)
                                                                                                             cAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CTTTTTTTAAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGATTCGATCTCCGGCAAAATCCAGACGGTACGCagagagn                                                                              <  1:49462/194‑2 (MQ=255)
                                                                                                                                 tgtgGAAACGCGCCGCGATGATGTGGCGCGCATTGTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGGATTCGTACGGTGGT‑CTTTTTTTAAAGGCTGCCCGCATCTTTGCCCGTGGTGTGTTAAGCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGGGTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTcg  <  1:43741/250‑1 (MQ=255)
                                                                                                                                                           cgcgcATTTTCAACATTCAGCGTTATTCACTGAATGAAGGCGAGGGCATTCGTACGGTGGTATTTTTTT‑AAAGGCTGTCCGCATCTTTGccc                                                                                                                                      >  1:167508/1‑92 (MQ=255)
                                                                                                                                                                                   nttNACTGANTGACNGCGAGGGCATTCGTANGGTGGT‑CTTTTTTTAAAGGNNGNCCGCANCTTTg                                                                                                                                         <  2:107882/64‑1 (MQ=255)
                                                                                                                                                                                                                        cttttttt‑AAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGGGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACgg          <  1:41540/153‑1 (MQ=255)
                                                                                                                                                                                                                                |                                                                                                                                                           
GCGCAACAGCGCCCACAAGATTATGCCGGGCTGGTGGTGCGCGTTGCCGGATACAGCGCCTTCTTTGTCGAACTGTCGAAGGAGATCCAGGATGACATCATCCGCCGGACAGCGCATCAGCTGTAACGTTGTGGAAACGCGCCGCGATGATGTGGCGCGCATTTTCAACATTCAGCGTTATTCACTGAATGACGGCGAGGGCATTCGTACGGTGGT‑CTTTTTT‑AAAGGCTGTCCGCATCTTTGCCCGTGGTGTGCTAATCCGGAGTCGATCTCCGGCAAAATCCAGACGGTACGCAGAGAGGCGAAATGTCTGCACTGTGCGAAATGTTTGCGTGATGCGGATGAATGCCCCTCCGGGGCGTTTGAACGGATTGGTCG  >  NC_000913/4146167‑4146544

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑