Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 52,880 T→G A179A (GCA→GCC pdxA ← 4‑hydroxy‑L‑threonine phosphate dehydrogenase, NAD‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb52,8800TG87.5% 143.4 / 7.9 48A179A (GCA→GCCpdxA4‑hydroxy‑L‑threonine phosphate dehydrogenase, NAD‑dependent
Reads supporting (aligned to +/- strand):  ref base T (6/0);  new base G (0/42);  total (6/42)
Fisher's exact test for biased strand distribution p-value = 8.15e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.47e-01

ATGCAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGTGATAGCGTCTGCGATATCGCGCAGCGGTAAATGCGTCGTTGCCAGC  >  W3110S.gb/52850‑52934
                              |                                                      
aTGCAAAATAGCAATCAGTTCGTGCAAAAGGGCAgg                                                   <  1:2492339/36‑1 (MQ=37)
aTGCAAAATAGCAATCACTTCGTGCAAAAGGGCAgg                                                   <  1:2768591/36‑1 (MQ=255)
aTGCAAAATAGCAATCACTTCGTGCAAAAGGGCAgg                                                   <  1:227049/36‑1 (MQ=255)
aTGCAAAATAGCAATCACTTCGTGCAAAAGGGCAgg                                                   <  1:1008765/36‑1 (MQ=255)
aTGCAAAATAGCAATCACTTCGTGCAAAAGGGCAgg                                                   <  1:2317189/36‑1 (MQ=255)
aTGCAAAATAGCAATCACTTCGTGCAAAAGGGCAgg                                                   <  1:2394850/36‑1 (MQ=255)
aTGCAAAATAGCAATCACTTCGTGCAAAAGGGCAgg                                                   <  1:2491863/36‑1 (MQ=255)
aTGCAAAATAGCAATCACTTCGTGCAAAAGGGCAgg                                                   <  1:2571336/36‑1 (MQ=255)
aTGCAAAATAGCAATCACTTCGTGCAAAAGGGCAgg                                                   <  1:2613852/36‑1 (MQ=255)
aTGCAAAATAGCAATCACTTCGTGCAAAAGGGCAgg                                                   <  1:2652131/36‑1 (MQ=255)
aTGCAAAATAGCAATCACTTCGTGCAAAAGGGCAgg                                                   <  1:2673100/36‑1 (MQ=255)
aTGCAAAATAGCAATCACTTCGTGCAAAAGGGCAgg                                                   <  1:227619/36‑1 (MQ=255)
aTGCAAAATAGCAATCACTTCGTGCAAAAGGGCAgg                                                   <  1:282560/36‑1 (MQ=255)
aTGCAAAATAGCAATCACTTCGTGCAAAAGGGCAgg                                                   <  1:2899291/36‑1 (MQ=255)
aTGCAAAATAGCAATCACTTCGTGCAAAAGGGCAgg                                                   <  1:377071/36‑1 (MQ=255)
aTGCAAAATAGCAATCACTTCGTGCAAAAGGGCAgg                                                   <  1:381936/36‑1 (MQ=255)
aTGCAAAATAGCAATCACTTCGTGCAAAAGGGCAgg                                                   <  1:421127/36‑1 (MQ=255)
aTGCAAAATAGCAATCACTTCGTGCAAAAGGGCAgg                                                   <  1:668531/36‑1 (MQ=255)
aTGCAAAATAGCAATCACTTCGTGCAAAAGGGCAgg                                                   <  1:67198/36‑1 (MQ=255)
aTGCAAAATAGCAATCACTTCGTGCAAAAGGGCAgg                                                   <  1:718613/36‑1 (MQ=255)
aTGCAAAATAGCAATCACTTCGTGCAAAAGGGCAgg                                                   <  1:734445/36‑1 (MQ=255)
aTGCAAAATAGCAATCACTTCGTGCAAAAGGGCAgg                                                   <  1:968280/36‑1 (MQ=255)
aTGCAAAATAGCAATCACTTCGTGCAAAAGGGCAgg                                                   <  1:1630042/36‑1 (MQ=255)
aTGCAAAATAGCAATCACTTCGTGCAAAAGGGCAgg                                                   <  1:1212983/36‑1 (MQ=255)
aTGCAAAATAGCAATCACTTCGTGCAAAAGGGCAgg                                                   <  1:1300350/36‑1 (MQ=255)
aTGCAAAATAGCAATCACTTCGTGCAAAAGGGCAgg                                                   <  1:1203472/36‑1 (MQ=255)
aTGCAAAATAGCAATCACTTCGTGCAAAAGGGCAgg                                                   <  1:1345401/36‑1 (MQ=255)
aTGCAAAATAGCAATCACTTCGTGCAAAAGGGCAgg                                                   <  1:1391202/36‑1 (MQ=255)
aTGCAAAATAGCAATCACTTCGTGCAAAAGGGCAgg                                                   <  1:1430498/36‑1 (MQ=255)
aTGCAAAATAGCAATCACTTCGTGCAAAAGGGCAgg                                                   <  1:1506289/36‑1 (MQ=255)
aTGCAAAATAGCAATCACTTCGTGCAAAAGGGCAgg                                                   <  1:1605801/36‑1 (MQ=255)
aTGCAAAATAGCAATCACTTCGTGCAAAAGGGCAgg                                                   <  1:2221033/36‑1 (MQ=255)
aTGCAAAATAGCAATCACTTCGTGCAAAAGGGCAgg                                                   <  1:1608515/36‑1 (MQ=255)
aTGCAAAATAGCAATCACTTCGTGCAAAAGGGCAgg                                                   <  1:1754285/36‑1 (MQ=255)
aTGCAAAATAGCAATCACTTCGTGCAAAAGGGCAgg                                                   <  1:1796028/36‑1 (MQ=255)
aTGCAAAATAGCAATCACTTCGTGCAAAAGGGCAgg                                                   <  1:1966744/36‑1 (MQ=255)
aTGCAAAATAGCAATCACTTCGTGCAAAAGGGCAgg                                                   <  1:2012471/36‑1 (MQ=255)
aTGCAAAATAGCAATCACTTCGTGCAAAAGGGCAgg                                                   <  1:1118077/36‑1 (MQ=255)
aTGCAAAATAGCAATCACTTCGTGCAAAAGGGCAgg                                                   <  1:2027279/36‑1 (MQ=255)
aTGCAAAATAGCAATCACTTCGTGCAAAAGGGCAgg                                                   <  1:2184316/36‑1 (MQ=255)
aTGCAAAATAGCAATCACTTCGTGCAAAAGGGCAgg                                                   <  1:1056992/36‑1 (MQ=255)
aTGCAAAATAGCAATCACTTCGAGCAAAAGGGCAgg                                                   <  1:1707075/36‑1 (MQ=38)
                       gCAAAAGTGCAGGGGTGATAGCGTCTGCGATATCGCGCAGCGGTAAATGCGTCGTTGCCAgc  >  1:273648/1‑62 (MQ=255)
                       gCAAAAGTGCAGGGGTGATAGCGTCTGCGATATCGCGCAGCGGTAAATGCGTCGTTGCCAgc  >  1:724840/1‑62 (MQ=255)
                       gCAAAAGTGCAGGGGTGATAGCGTCTGCGATATCGCGCAGCGGTAAATGCGTCGTTGCCAgc  >  1:1302027/1‑62 (MQ=255)
                       gCAAAAGTGCAGGGGTGATAGCGTCTGCGATATCGCGCAGCGGTAAATGCGTCGTTGCCAgc  >  1:2536287/1‑62 (MQ=255)
                       gCAAAAGTGCAGGGGTGATAGCGTCTGCGATATCGCGCAGCGGTAAATGCGTCGTTGCCAgc  >  1:2185087/1‑62 (MQ=255)
                       gCAAAAGTGCAGGGGTGATAGCGTCTGCGATATCGCGCAGCGGTAAATGCGTCGTTGCCAg   >  1:235218/1‑61 (MQ=255)
                              |                                                      
ATGCAAAATAGCAATCACTTCGTGCAAAAGTGCAGGGGTGATAGCGTCTGCGATATCGCGCAGCGGTAAATGCGTCGTTGCCAGC  >  W3110S.gb/52850‑52934

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