Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 570,448 A→C 39.8% D599A (GAC→GCC)  fdoG_2 → Formate dehydrogenase‑O major subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB122570,4480AC39.8% 16.2 / 4.9 33D599A (GAC→GCC) fdoG_2Formate dehydrogenase‑O major subunit
Reads supporting (aligned to +/- strand):  ref base A (9/10);  new base C (11/3);  total (20/13)
Fisher's exact test for biased strand distribution p-value = 8.73e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGGGCCTGGGCCTGGCCGGCCAACCGCCGCATCCTTTACAACCGTGCCTCCAGCGACCCGCAAGGCAAGCCGTGGGACCCGGAGAAAAAGCGCCTGGTATGGTGGAACGGCAAGGCCTGGACCGGCACCGACGTGCCGGACTTCAAGGTCGACTCGCCGCCCGAGGCGGGGATGAACCCTTTCATCATGAACCCCGAAGGCGTGGCGCGTTTCTTCGCCATCGACAAGATGGCCGAAGGCCCGTTCCCCGAGCACTACGA  >  ACB122/570317‑570577
                                                                                                                                   |                                                                                                                                 
ctgggcctgggcCTGGCCGGCCAACCGCCGCATCCTTTACAACCGTGCCTCCAGCGACCCGCAAGGCAAGCCGTGGGACCCGGAGAAAAAGCGCCTGGTATGGTGGAACGGCAAGGCCTGGACCGGCCCCGACGt                                                                                                                                >  6:108473/1‑135 (MQ=255)
ctgggcctgggcCTGGCCGGCCAACCGCCGCATCCTTTACAACCGTGCCTCCAGCGACCCGCAAGGCAAGCCGTGGGACCCGGAGAAAAAGCGCCTGGTATGGTGGAACGGCAAGGCCTGGACCCGCAccgccgg                                                                                                                                >  4:33587/1‑134 (MQ=255)
 ggggcgtgggcgtGGCGGGCCAACCGCCGAATCCTTTACCCCCGTGCCTCCCGCGACCCGCACGCAACGCGGGGGGACCCGGATAAAACGCGCCTGGTATGGTGGAACGGCAAGGCCTGGACCGGCACCGACGTg                                                                                                                               <  4:297378/134‑1 (MQ=255)
      ctgggcctggCCGGCCAACCGCCGCATCCTTTACAACCGTGCCTCCAGCGACCCGCAAGGCAAGCCGTGGGACCCGGAGAAAAAGCGCCTGGTATGGTGGAACGGCAAGGCCTGGACCGGCCCCGCCGTGGCggcc                                                                                                                         >  7:196003/1‑134 (MQ=255)
      ctgggcctggCCGGCCAACCGCCGCATCCTTTACAACCGTGCCTCCAGCGACCCGCAAGGCAAGCCGTGGGACCCGGAGAAAAAGCGCCTGGTATGGTGGAACGGCAAGGCCTGGACCCGCACCGCCGTGCCgggc                                                                                                                         >  7:42519/1‑134 (MQ=255)
           cctggCCGGCCAACCGCCGCATCCTTTACTACCGTGCCTCCAGCGACCCGCAAGGCAAGCCGTGGGACCCGGAGAAAAAGCGCCTGGTATGGTGGAACGGCAAGGCCTGGACCGGCACCGACGTGCCCGCCtt                                                                                                                       >  1:59865/1‑133 (MQ=255)
                ccggccAACCGCCGCATCCTTTACAACCGTGCCTCCAGCGACCCGCAAGGCAAGCCGTGGGACCCGGAGAAAAAGCGCCTGGTATGGTGGAACGGCAAGGCCTGGACCGGCACCGACGTGCCGGACTTCAAGGTc                                                                                                                >  2:259947/1‑135 (MQ=255)
                 cggccAACCGCCGCATCCTTTACAACCGTGCCTCCAGCGACCCGCAAGGCAAGCCGTGGGACCCGGAGAAAAAGCGCCTGGTATGGTGGAACGGCAAGGCCTGGAGCGGCCCCGCCGTGCGGGACTTCAAGGTCg                                                                                                               >  2:138481/1‑135 (MQ=255)
                 cggccAACCGCCGCATCCTTTACAACCGTGCCTCCAGCGACCCGCAAGGCAAGCCGTGGGACCCGGAGAAAAAGCGCCTGGTATGGTGGAACGGCAAGGCCTGGACCGGCACCCACGTGCCGGACTTCCACGTCg                                                                                                               >  3:156137/1‑135 (MQ=255)
                  ggccAACCGCCGCATCCTTTACAACCGTGCCTCCAGCGACCCGCAAGGCAAGCCGTGGGACCCGGAGAAAAAGCGCCTGGTATGGTGGAACGGCAAGGCCTGGACCGGCCCCGCCGGGCCGGCCTGCAAGGGCgc                                                                                                              >  5:132594/1‑134 (MQ=255)
                            cgcATCCTTTACAACCGTGCCTCCAGCGACCCGCAAGGCAAGCCGTGGGACCCGGAGAAAAAGCGCCTGGTATGGTGGAACGGCAAGGCCTGGACCGGCCCCGCCGTGCCGGCCGTCAAGGTCGACTCGCCGCCCg                                                                                                   >  8:8538/1‑136 (MQ=255)
                                 ccTTTACAACCGTGCCTCCAGCGACCCGCAAGGCAAGCCGTGGGACCCGGAGAAAAAGCGCCTGGTATGGTGGAACGGCAAGGCCTGGACCCGCACCGCCGTGCCGGACTTCAAGGTCGACTCCCCGCCCGCGGc                                                                                               >  2:289057/1‑135 (MQ=255)
                                                cTCCAGCGACCCGCAAGGCAAGCCGTGGGACCCGGAGAAAAAGCGCCTGGTATGGTGGAACGGCAAGGCCTGGACCGGCCCCGCCGGGCCGGCCGTCAAGGTCGACTCGCCGCCCGCGGCGGGGATGAAcccat                                                                                 >  1:227807/1‑132 (MQ=255)
                                                cTCCAGCGACCCGCAAGGCAAGCCGTGGGACCCGGAGAAAAAGCGCCTGGTATGGTGGAACGGCAAGGCCTGGACCGGCCCCGACGTGGCGGCCTTCAAGGTCGACTCGCCGCCCGCGGCGGGGATGCACCCtt                                                                                 >  1:48751/1‑134 (MQ=255)
                                                cTCCAGCGACCCGCAAGGCAAGCCGTGGGACCCGGAGAAAAAGCGCCTGGTATGGTGGAACGGCAAGGCCTGGACCGGCCCCGACGGGCCGGGCTTCAAGGTCGACGAGCCGCCCGCGGCGGGGATGCAcccgt                                                                                 >  4:355633/1‑132 (MQ=255)
                                                cTCCAGCGACCCGCAAGGCAAGCCGTGGGACCCGGAGAAAAAGCGCCTGGTATGGTGGAACGGCAAGGCCTGGACCGGCACCGCCGTGCCGGACTTCAAGGTCGACTCGCCGCCCGAGGCGGGGATGAACCCtt                                                                                 >  4:76466/1‑134 (MQ=255)
                                                cTCCAGCGACCCGCAAGGCAAGCCGTGGGACCCGGAGAAAAAGCGCCTGGTATGGTGGAACGGCAAGGCCTGGACCGGCACCGACGTGACGGACTTCAAGGTCGCCACGCCCCCCGCGGCCGGGATGCACCCttt                                                                                >  2:322968/1‑135 (MQ=255)
                                                cTCCAGCGACCCGCAAGGCAAGCCGTGGGACCCGGAGAAAAAGCGCCTGGTATGGTGGAACGGCAAGGCCTGGACCGGCACCGACGGGGCGGACTTCAAGGTCGACTCGCCGCCCGAGGCGGGGTTGAAcccat                                                                                 >  3:51261/1‑132 (MQ=255)
                                                                              cccGGAGAAAAAGCGCCTGGTATGGTGGAACGGCAAGGCCTGGACCGGCACCGCCGTGCCGGACTTCAAGGTCGACTCGCCGCCCGAGGCGGGGATGAACCCTTTCATCATGAACCCCGAAGGCGTGGCGCGttt                                                  >  6:204921/1‑135 (MQ=255)
                                                                                                           ccGGCACGCCCGGGACCGCCCCCGCCGTCCCGGACTTCAAGGGCGACTCGCCGCCCGAGGCGGGGATGAACCCTTTCATCATGAACCCCGAAGGCGTGGCGCGTTTCTTCGCCATCGACAAGATGGCCGAAGGcc                     <  3:33587/134‑1 (MQ=255)
                                                                                                           acggcacggcCGGGCCCGGCACCGACGTGCCGGACTTCAAGGTCGACTCGCCGCCCGAGGCGGGGATGAACCCTTTCATCATGAACCCCGAAGGCGTGGCGCGTTTCTTCGCCATCGACAAGATGGCCGAAGGcc                     <  5:108473/135‑1 (MQ=255)
                                                                                                              gCAAGGCCTGGACCGGCCCCCCCGTGCCGGCCTTCAAGGTCGACGCGCCGCCCCCGGCGGGGCTGAACCCGTTCATCATGAACCCCGCAGGCGTGGCGCGGGTCTTAGCCCTCGACAAAATTGCCCCAGGCCCgg                  >  2:354771/1‑134 (MQ=255)
                                                                                                               cAAGGCCTGGAGCGGCACCGACGTGCCGGACTTCAAGTTCGAATCGCCGCCCGAGGCGGGGATGAACCCTTTCATAATGAACCCCGAAGGCGTGGCGCGTTTGTTAGCCCTCGACAAGAAGGCCCACGGCCCGTTc                >  8:333327/1‑136 (MQ=255)
                                                                                                               cAAGGCCGGCAGGGGCACCGCCGTGCGGGACTTCCAGGTCGACTCGCCGCCCGAGGCGGGGATGAACCCTTTCATCATGAACCCCGAAGGCGTGGCGCGTTTCTTCGCCATCGACAAGATGGCCGAAGGCCCGtt                 <  2:138054/135‑1 (MQ=255)
                                                                                                                acggCCGGGCCCGGCACCGACGGGCCGGACTTCAAGGTCGACTCGCCGCCCGAGGCGGGGATGAACCCTTTCATCATGAACCCCGAAGGCGTGGCGCGTTTCTTCGCCATCGACAAGATGGCCGAAGGCCCGTTcc               <  7:37430/134‑1 (MQ=255)
                                                                                                                   gcgtGGACCGGCACCGAGGTGCCGGACTTCAAGGTCGACTCGCCGCCCGAGGCGGGGATGAACCCTTTCATCATGAACCCCGAAGGCGTGGCGCGTTTCTTCGCCATCGACAAGATGGCCGAAGGCCCGTTcccc             <  4:51261/132‑1 (MQ=255)
                                                                                                                   gccggcaCCGGCACCGAGGGGCCGGACTTCAAGGTCGACTCGCCGCCCGAGGCGGGGATGAACCCTTTCATCATGAACCCCGAAGGCGTGGCGCGTTTCTTCGCCATCGACAAGATGGCCGAAGGCCCGTTcccc             <  5:176942/129‑1 (MQ=255)
                                                                                                                   gCCTGGACCGGCACCGACGTGCCGGACTTCAAGGTCGACTCGCCGCCCGAGGCGGGGATGAACCCTTTCATCATGAACCCCGAAGGCGTGGCGCGTTTCTTCGCCATCGACAAGATGGCCGAAGGCCCGTTcccc             <  3:76466/135‑1 (MQ=255)
                                                                                                                   gCCTGGACCGGCACCGACGTGCCGGACTTCAAGGTCGACTCGCCGCCCGAGGCGGGGATGAACCCTTTCATCATGAACCCCGAAGGCGTGGCGCGTTTCTTCGCCATCGACAAGATGGCCGAAGGCCCGTTcccc             <  4:196219/135‑1 (MQ=255)
                                                                                                                   gCCGGGCCCGGCACCGACGTGCCGGACTTCAAGGTCGACTCGCCGCCCGAGGCGGGGATGAACCCTTTCATCATGAACCCCGAAGGCGTGGCGCGTTTCTTCGCCATCGACAAGATGGCCGAAGGCCCGTTcccc             <  6:132594/135‑1 (MQ=255)
                                                                                                                          ccGGCACCGCCGTGCCGGACTTCAAGGTCGACTCGCCGCCCGCGGCGGGGATGAACCCTTTCATCATGAACCCCCCACGCGTGGCGCGTTTCTTAGCCATCGACAAGATTGCCGAAGGCCCGGTCCCCCCGCAc       >  4:177120/1‑134 (MQ=255)
                                                                                                                            ggCACCGCCGTGCCGGACTTCAAGGGCGACTCGCCGCCCGAGGCGGGGATGAACCCTTTCATCATGAACCCCGAAGGCGTGGCGCGTTTCTTCGCCATCGACAAGATGGCCGAAGGCCCGTTCCCCGAGCACTAc    <  1:289057/135‑1 (MQ=255)
                                                                                                                              cACCGAGGGGCGGCACTGCATGGTCGAGGCGCCGCCCGAGGGGGGGATTAACACTTTCATCATGAACCCCGAAGGCGTGGCTCGTTTCTTCGCCATCGACAAGATGGCCGAAGGCCCGTTCCCCGAGCACTACGa  <  4:355480/135‑1 (MQ=255)
                                                                                                                              cACCGACGTGCCGGACTTCAAGGTCGACTCGCCGCCCGAGGCGGGGATGAACCCTTTCATCATGAACCCCGAAGGCGTGGCGCGTTTCTTCGCCATCGACAAGATGGCCGAAGGCCCGTTCCCCGAGCACTACGa  <  4:135067/135‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                 
CTGGGCCTGGGCCTGGCCGGCCAACCGCCGCATCCTTTACAACCGTGCCTCCAGCGACCCGCAAGGCAAGCCGTGGGACCCGGAGAAAAAGCGCCTGGTATGGTGGAACGGCAAGGCCTGGACCGGCACCGACGTGCCGGACTTCAAGGTCGACTCGCCGCCCGAGGCGGGGATGAACCCTTTCATCATGAACCCCGAAGGCGTGGCGCGTTTCTTCGCCATCGACAAGATGGCCGAAGGCCCGTTCCCCGAGCACTACGA  >  ACB122/570317‑570577

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