Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 570,448 A→C 42.2% 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,4480AC42.2% 9.9 / 5.1 34D599A (GAC→GCC) fdoG_2Formate dehydrogenase‑O major subunit
Reads supporting (aligned to +/- strand):  ref base A (9/9);  new base C (13/2);  total (22/12)
Fisher's exact test for biased strand distribution p-value = 3.42e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGGCCTGGGCCTGGCCGGCCAACCGCCGCATCCTTTACAACCGTGCCTCCAGCGACCCGCAAGGCAAGCCGTGGGACCCGGAGAAAAAGCGCCTGGTATGGTGGAACGGCAAGGCCTGGACCGGCACCGACGTGCCGGACTTCAAGGTCGACTCGCCGCCCGAGGCGGGGATGAACCCTTTCATCATGAACCCCGAAGGCGTGGCGCGTTTCTTCGCCATCGACAAGATGGCCGAAGGCCCGTTCCCCGAGCACTACGAGCCCT  >  ACB122/570318‑570582
                                                                                                                                  |                                                                                                                                      
tgggcctgggccTGGCCGGCCAACCGCCGCATCCTTTACAACCGTGCCTCCAGCGACCCGCAAGGCAAGCCGTGGGACCCGGAGAAAAAGCGCCTGGTATGGTGGAACGGCAAGGCCTGGACCGGCGccgccgag                                                                                                                                    >  5:118385/1‑133 (MQ=255)
     ctgggcctggCCGGCCAACCGCCGCATCCTTTACAACCGTGCCTCCAGCGCCCCGCAAGGCAAGCCGTGGGACCCGGAGAAAAAGCGCCTGGGATGGTGGAACGGCAAGGCCTGGACCGGCACCGACGTGGCggcc                                                                                                                              >  4:292109/1‑134 (MQ=255)
     ctgggcctggCCGGCCAACCGCCGCATCCTTTACAACCGTGCCTCCAGCGACCCGCAAGGCAAGCCGTGGGACCCGGAGAAAAAGCGCCTGGTATGGTGGAACGGCAAGGCCTGGAGCGGCACCGCCGGGCCGGa                                                                                                                               >  5:369477/1‑135 (MQ=255)
     ctgggcctggCCGGCCAACCGCCGCATCCTTTACAACCGTGCCTCCAGCGACCCGCAAGGCAAGCCGTGGGACCCGGAGAAAAAGCGCCTGGTATGGTGGAACGGCAAGGCCTGGACCGGCCCCgccgggccggcc                                                                                                                              >  7:392697/1‑134 (MQ=255)
     ctgggcctggCCGGCCAACCGCCGCATCCTTTACAACCGTGCCTCCAGCGACCCGCAAGGCAAGCCGTGGGACCCGGAGAAAAAGCGCCTGGTATGGTGGAACGGCAAGGCCTGGACCGGCCCCgccgggccggc                                                                                                                               >  1:381564/1‑134 (MQ=255)
     ctgggcctggCCGGCCAACCGCCGCATCCTTTACAACCGTGCCTCCAGCGACCCGCAAGGCAAGCCGTGGGACCCGGAGAAAAAGCGCCTGGTATGGTGGAACGGCAAGGCCTGGACCGGCACCGCCGTGGCggc                                                                                                                               >  1:162771/1‑134 (MQ=255)
     ctgggcctggCCGGCCAACCGCCGCATCCTTTACAACCGTGCCTCCAGCGACCCGCAAGGCAAGCCGTGGGACCCGGAGAAAAAGCGCCTGGTATGGTGGAACGGCAAGGCCTGGACCGGCACCGACGTGCCgggc                                                                                                                              >  3:125586/1‑134 (MQ=255)
          cctggCCGGCCAACCGCCGCATCCTTTACAACCGTGCCTCCAGCGACCCGCAAGGCAAGCCGTGGGACCCGGAGAAAAAGCGCCTGGTATGGTGGAACGGCAAGGCCTGGACCGGCCCCgccgggccgggcgt                                                                                                                            >  1:70498/1‑129 (MQ=255)
          cctggCCGGCCAACCGCCGCATCCTTTACAACCGTGCCTCCAGCGACCCGCAAGGCAAGCCGTGGGACCCGGAGAAAAAGCGCCTGGTATGGTGGAACGGCAAGGCCTGGACCGGCCCCGACGTGCCGGAcgt                                                                                                                            >  6:379891/1‑131 (MQ=255)
                 ggccAACCGCCGCATCCTTTACAACCGTGCCTCCAGCGACCCGCAAGGCAAGCCGTGGGACCCGGAGAAAAAGCGCCTGGTATGGTGGAACGGCAAGGCCTGGACCGGCACCGCCGGGGCGGACGTCAAGGTCGAc                                                                                                                  >  4:219920/1‑136 (MQ=255)
                     aaCCGCCGCATCCTTTACAACCGTGCCTCCAGCGACCCGCAAGGCAAGCCGTGGGACCCGGAGAAAAAGCGCCTGGTATGGTGGAACGGCAAGGCCTGGACCGGCCCCCCCGTGCCGGCCTGCAAGGTCGACTcgc                                                                                                              >  8:247949/1‑136 (MQ=255)
                            ccacccTTTACAACCGGGCCTCCTCCCACCCCCAAGGCACGCCGGGGGCCCCGCATAAAAAGCGCCTGGTATGGTGGAACGGCAAGGCCTGGACCGGCACCGACGTGCCGGACTTCAAGGTCGACTCGCCGCCCGa                                                                                                       <  8:397060/132‑1 (MQ=255)
                                       aaCCGTGCCTCCAGCGACCCGCAAGGCAAGCCGTGGGACCCGGAGAAAAAGCGCCTGGTATGGTGGAACGGCAAGGCCTGGAGCCGCCCCCCCGGGCCGGCCGTCAAGGTGGACTCCCCCCCCCCCGCGGGGGTg                                                                                             >  1:360251/1‑135 (MQ=255)
                                       aaCCGTGCCTCCAGCGACCCGCAAGGCAAGCCGTGGGACCCGGAGAAAAAGCGCCTGGTATGGTGGAACGGCAAGGCCTGGACCGGCCCCGCCGTGCCGGACTTCAAGGTCGACTCGCCGCCCGCGGCGGGGATg                                                                                             >  6:278401/1‑135 (MQ=255)
                                       aaCCGTGCCTCCAGCGACCCGCAAGGCAAGCCGTGGGACCCGGAGAAAAAGCGCCTGGTATGGTGGAACGGCAAGGCCTGGACCGGCCCCGCCGGGCCGGCCGTCAAGGTCGGCGCGCCGCCCGCGGCGGGGCTg                                                                                             >  5:129266/1‑135 (MQ=255)
                                          cGTGCCCCCCCCCACCCCAAACGCACGCGGGGGGACCCGGATAACAAGCGCCTGGTATGGTGGACGGGCAAGGCCGGGACCGGCACCGTCGTGGCGGACTTCAAGGTCGACTCGCCGCCCGAGGCGGGGATGAAcc                                                                                         <  4:355781/136‑1 (MQ=255)
                                            tGCCTCCAGCGACCCGCAAGGCAAGCCGTGGGACCCGGAGAAAAAGCGCCTGGTATGGTGGAACGGCAAGGCCTGGACCGGCCACGAACTGGCGGCCGTAAAGGTCGACTCGCCGCCCGACGCGGGGAAGAAccc                                                                                        >  6:349219/1‑135 (MQ=255)
                                               cTCCAGCGACCCGCAAGGCAAGCCGTGGGACCCGGAGAAAAAGCGCCTGGTATGGTGGAACGGCAAGGCCTGGACCGGCCCCGCCGGGCCGGGCGTCAAGGTCGCCTCGCCGCCCGAGGCGGGGATGAAcccga                                                                                      >  1:346093/1‑132 (MQ=255)
                                               cTCCAGCGACCCGCAAGGCAAGCCGTGGGACCCGGAGAAAAAGCGCCTGGTATGGTGGAACGGCAAGGCCTGGACCGGCACCGCCGTGCCGGACTTCAAGGTCGACACGCCGCCCGAGGCGGGGCAGAAcccagt                                                                                     >  4:56915/1‑132 (MQ=255)
                                               cTCCAGCGACCCGCAAGGCAAGCCGTGGGACCCGGAGAAAAAGCGCCTGGTATGGTGGAACGGCAAGGCCTGGACCGGCACCGCCGTGCCGCACTTCAAGGTCGACTCGCCGCCCGCGGCGGGGATGAACCCttt                                                                                     >  7:220172/1‑135 (MQ=255)
                                               cTCCAGCGACCCGCAAGGCAAGCCGTGGGACCCGGAGAAAAAGCGCCTGGTATGGTGGAACGGCAAGGCCTGGACCGGCACCGCCGGGCCGGACTTCAAGGTCGCCTCGCCGCCCGAGGCGGGGATGAACCCttt                                                                                     >  7:305221/1‑135 (MQ=255)
                                               cTCCAGCGACCCGCAAGGCAAGCCGTGGGACCCGGAGAAAAAGCGCCTGGTATGGTGGAACGGCAAGGCCTGGACCGGCACCGACGTGCCGGACTTCAAGGTCGACTCGCCGCCCGAGGCGGGGATGAACCCtt                                                                                      >  2:97843/1‑134 (MQ=255)
                                               cTCCAGCGACCCGCAAGGCAAGCCGTGGGACCCGGAGAAAAAGCGCCTGGTATGGTGGAACGGCAAGGCCTGGACCGGCACCCCCGGGCCGGACGTCAAGGTCGACTCGCCGCCCGCGGCGGGGATGAACCCttt                                                                                     >  3:184250/1‑135 (MQ=255)
                                               cTCCAGCGACCCGCAAGGCAAGCCGTGGGACCCGGAGAAAAAGCGCCTGGTATGGTGGAACGGCAAGGCCTGGACCGGCACCCACGTGCCGGGCCTAAAAGTCGACTCCCCCCCCCACGCGGGGGGGAACCCgtt                                                                                     >  7:357677/1‑135 (MQ=255)
                                                 ccAGCGCCCCGCAAGGCAAGCCGTGGGACCCGGAGAAAAAGCGCCTGGTATGGTGGAACGGCAAGGCCTGGACCGGCACCGCCGTGCCGGCCTTCAAGGTCGACTCGCCGCCCGAGGCGGGGATGCACCCTTtca                                                                                   >  6:199354/1‑135 (MQ=255)
                                                   cccgccccGCACGCCACGCGGGGGGACCCGTATAAACCGCGCCTGTTATGGTGGAACGGCAAGGCCGGGACCGGCACCGACGTGCCGGACTTCAAGGTCGACTCGCCGCCCGAGGCGGGGATGAACCCTTtcatc                                                                                 <  1:35705/130‑1 (MQ=255)
                                                           gcaaggcaagCCGTGGGACCCGGAGAAAAAGCGCCTGGTATGGTGGAACGGCAAGGCCTGGACCGGCCCCGACGTGCCGGCCGTCAAGGTCGCCTCGCCCCCCGCGGCGGGGCTGGACCCGTTGATCATGAAcccc                                                                        >  3:40415/1‑136 (MQ=255)
                                                                  acgCCGTGGGACCCGCACAAAACGCGCCTGGTATGGTGGAACGGCAAGGCCGGGCCCGGCACCGACGGGCCGGACTTCAAGGTCGACTCGCCGCCCGAGGCGGGGATGAACCCTTTCATCATGAACCCCGAAGGCg                                                                 <  7:267328/134‑1 (MQ=255)
                                                                       gTGGGACCCGGAGAAAAAGCGCCTGGTATGGTGGAACGGCAAGGCCTGGACCGGCACCCACGTGCCGGACTTCAACGTAGACTCGCCGCCCGAGGCGGGGATGAACCCGTTCATCATGAACCCCGAACGCGTGgcg                                                            >  4:253217/1‑136 (MQ=255)
                                                                                          cgcCGGGTTTGGGGCAACGGCACGGCCTGGACCGGCACCGCCGGGCCGGACTTCAAGGTCGACTCGCCGCCCGAGGCGGGGATGAACCCTTTCATCATGAACCCCGAAGGCGTGGCGCGTTTCTTCGCCATCGACa                                         <  3:79348/136‑1 (MQ=255)
                                                                                                            ggCCAGGCCGGGCCCGGCACCGACGGGCCGGACTTCAAGGTCGACTCGCCGCCCGAGGCGGGGATGAACCCTTTCATCATGAACCCCGAAGGCGTGGCGCGTTTCTTCGCCATCGACAAGATGGCCGAAGGCCCGt                       <  8:305221/136‑1 (MQ=255)
                                                                                                            ggCACGGCGTGCACCGGCACCGACGGCCGGGACTTCAAGGTCGACGCGCCGCCCGAGGCGGGGATGAACCCTTTCATCATGAACCCCGAAGGCGTGGCGCGTTTCTTCGCCATCGACAAGATGGCCGAAGGCCCGt                       <  3:369306/136‑1 (MQ=255)
                                                                                                                aGGCCTGGACCGGCACCGACGTGCCGGGCTTCAAGGTCGACTCGCCGCCCGCGGCGGGGATGAACCCTATAATCATGAACCCCGAAGGCGTGGGGCGTTTGTTAGCCATCCAAAACATGGCCCAACGCCCGGTccc                   >  8:35574/1‑136 (MQ=255)
                                                                                                                 gggcTGGACCGGCACCGAGGGGCCGGACTTCAAGGTCGACTCGCCGCCCGAGGCGGGGATGAACCCTTTCATCATGAACCCCGAAGGCGTGGCGCGTTTCTTCGCCATCGACAAGATGGCCGAAGGCCCGTTcccc                  <  7:13204/133‑1 (MQ=255)
                                                                                                                        aCCGGCACCGACGTGCCGGACTTCAAGGTCGACTCGCCGCCCGAGGCGGGGATGAACCCTTTCATCATGAACCCCGAAGGCGTGGCGCGTTTCTTCGCCATCGACAAGATGGCCGAAGGCCCGTTCCCCGAGCACt           <  4:15801/136‑1 (MQ=255)
                                                                                                                         ccGGCACCGCCGGGCCGGACTTCAAGGTCGACTCGCCGCCCGAGGCGGGGATGAACCCTTTCATCATGAACCCCGAAGGCGTGGCGCGTTTCTTCGCCATCGACAAGATGGCCGAAGGCCCGTTCCCCGAGCACt           <  2:360251/135‑1 (MQ=255)
                                                                                                                         ccGGCACCGACGGGCCGGACTGCAAGGTCGACTCGCCGCCCGAGGCGGGGATGAACCCTTTCATCATGAACCCCGAAGGCGTGGCGCGTTTCTTCGCCATCGACAAGATGGCCGAAGGCCCGTTCCCCGAGCACt           <  2:100341/135‑1 (MQ=255)
                                                                                                                            gCACCGACGTGCCGGACTTCAAGGTCGACTCGCCGCCCGAGGCGGGGATGAACCCTTTCATCATGAACCCCGAAGGCGTGGCGCGTTTCTTCGCCATCGACAAGATGGCCGAAGGCCCGTTCCCCGAGCACTACGa       <  4:185719/136‑1 (MQ=255)
                                                                                                                                  aCGTGCCGGACTTCAAGGTCGACTCGCCGCCCGAGGCGGGGATGAACCCTTTCATCATGAACCCCGAAGGCGTGGCGCGTTTCTTCGCCATCGACAAGATGGCCGAAGGCCCGTTCCCCGAGCACTACGAGCCCt  <  1:97843/135‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                      
TGGGCCTGGGCCTGGCCGGCCAACCGCCGCATCCTTTACAACCGTGCCTCCAGCGACCCGCAAGGCAAGCCGTGGGACCCGGAGAAAAAGCGCCTGGTATGGTGGAACGGCAAGGCCTGGACCGGCACCGACGTGCCGGACTTCAAGGTCGACTCGCCGCCCGAGGCGGGGATGAACCCTTTCATCATGAACCCCGAAGGCGTGGCGCGTTTCTTCGCCATCGACAAGATGGCCGAAGGCCCGTTCCCCGAGCACTACGAGCCCT  >  ACB122/570318‑570582

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: