Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 3,591,290 A→G V398A (GTC→GCC)  yihF ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,591,2900AG89.5% 152.7 / 6.9 57V398A (GTC→GCC) yihFconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (6/0);  new base G (0/51);  total (6/51)
Fisher's exact test for biased strand distribution p-value = 2.76e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

TCCGCATCCATTCCTTCAGATAAATTAAGCTGTTTTGCGGTTTCTGTGACGACATTAAGCGGTAACTTCACATCAAA  >  W3110S.gb/3591242‑3591318
                                                |                            
tCCGCATCCATCCCTTCATATAAATTAAGCTGTTTTGCGGTTTCTGTGGCGACATTAAGc                   <  1:2046320/60‑1 (MQ=255)
tCCGCATCCATCCCTTCAGATAAATTAAGCTGTTTTGCGGTTTCTGTGGCGACATTAAGc                   <  1:33617/60‑1 (MQ=255)
tCCGCATCCATCCCTTCAGATAAATTAAGCTGTTTTGCGGTTTCTGTGGCGACATTAAGc                   <  1:1099250/60‑1 (MQ=255)
tCCGCATCCATCCCTTCAGATAAATTAAGCTGTTTTGCGGTTTCTGTGGCGACATTAAGc                   <  1:225504/60‑1 (MQ=255)
tCCGCATCCATCCCTTCAGATAAATTAAGCTGTTTTGCGGTTTCTGTGGCGACATTAAGc                   <  1:2255147/60‑1 (MQ=255)
tCCGCATCCATCCCTTCAGATAAATTAAGCTGTTTTGCGGTTTCTGTGGCGACATTAAGc                   <  1:2264264/60‑1 (MQ=255)
tCCGCATCCATCCCTTCAGATAAATTAAGCTGTTTTGCGGTTTCTGTGGCGACATTAAGc                   <  1:234673/60‑1 (MQ=255)
tCCGCATCCATCCCTTCAGATAAATTAAGCTGTTTTGCGGTTTCTGTGGCGACATTAAGc                   <  1:2527599/60‑1 (MQ=255)
tCCGCATCCATCCCTTCAGATAAATTAAGCTGTTTTGCGGTTTCTGTGGCGACATTAAGc                   <  1:2543096/60‑1 (MQ=255)
tCCGCATCCATCCCTTCAGATAAATTAAGCTGTTTTGCGGTTTCTGTGGCGACATTAAGc                   <  1:2614248/60‑1 (MQ=255)
tCCGCATCCATCCCTTCAGATAAATTAAGCTGTTTTGCGGTTTCTGTGGCGACATTAAGc                   <  1:2726542/60‑1 (MQ=255)
tCCGCATCCATCCCTTCAGATAAATTAAGCTGTTTTGCGGTTTCTGTGGCGACATTAAGc                   <  1:2745641/60‑1 (MQ=255)
tCCGCATCCATCCCTTCAGATAAATTAAGCTGTTTTGCGGTTTCTGTGGCGACATTAAGc                   <  1:275287/60‑1 (MQ=255)
tCCGCATCCATCCCTTCAGATAAATTAAGCTGTTTTGCGGTTTCTGTGGCGACATTAAGc                   <  1:2863446/60‑1 (MQ=255)
tCCGCATCCATCCCTTCAGATAAATTAAGCTGTTTTGCGGTTTCTGTGGCGACATTAAGc                   <  1:2217127/60‑1 (MQ=255)
tCCGCATCCATCCCTTCAGATAAATTAAGCTGTTTTGCGGTTTCTGTGGCGACATTAAGc                   <  1:363553/60‑1 (MQ=255)
tCCGCATCCATCCCTTCAGATAAATTAAGCTGTTTTGCGGTTTCTGTGGCGACATTAAGc                   <  1:420163/60‑1 (MQ=255)
tCCGCATCCATCCCTTCAGATAAATTAAGCTGTTTTGCGGTTTCTGTGGCGACATTAAGc                   <  1:586677/60‑1 (MQ=255)
tCCGCATCCATCCCTTCAGATAAATTAAGCTGTTTTGCGGTTTCTGTGGCGACATTAAGc                   <  1:631926/60‑1 (MQ=255)
tCCGCATCCATCCCTTCAGATAAATTAAGCTGTTTTGCGGTTTCTGTGGCGACATTAAGc                   <  1:802080/60‑1 (MQ=255)
tCCGCATCCATCCCTTCAGATAAATTAAGCTGTTTTGCGGTTTCTGTGGCGACATTAAGc                   <  1:810696/60‑1 (MQ=255)
tCCGCATCCATCCCTTCAGATAAATTAAGCTGTTTTGCGGTTTCTGTGGCGACATTAAGc                   <  1:85751/60‑1 (MQ=255)
tCCGCATCCATCCCTTCAGATAAATTAAGCTGTTTTGCGGTTTCTGTGGCGACATTAAGc                   <  1:863787/60‑1 (MQ=255)
tCCGCATCCATCCCTTCAGATAAATTAAGCTGTTTTGCGGTTTCTGTGGCGACATTAAGc                   <  1:896596/60‑1 (MQ=255)
tCCGCATCCATCCCTTCAGATAAATTAAGCTGTTTTGCGGTTTCTGTGGCGACATTAAGc                   <  1:925167/60‑1 (MQ=255)
tCCGCATCCATCCCTTCAGATAAATTAAGCTGTTTTGCGGTTTCTGTGGCGACATTAAGc                   <  1:950981/60‑1 (MQ=255)
tCCGCATCCATCCCTTCAGATAAATTAAGCTGTTTTGCGGTTTCTGTGGCGACATTAAGc                   <  1:981868/60‑1 (MQ=255)
tCCGCATCCATCCCTTCAGATAAATTAAGCTGTTTTGCGGTTTCTGTGGCGACATTAAGc                   <  1:998248/60‑1 (MQ=255)
tCCGCATCCATCCCTTCAGATAAATTAAGCTGTTTTGCGGTTTCTGTGGCGACATTAAGc                   <  1:218754/60‑1 (MQ=255)
tCCGCATCCATCCCTTCAGATAAATTAAGCTGTTTTGCGGTTTCTGTGGCGACATTAAGc                   <  1:1158370/60‑1 (MQ=255)
tCCGCATCCATCCCTTCAGATAAATTAAGCTGTTTTGCGGTTTCTGTGGCGACATTAAGc                   <  1:1163785/60‑1 (MQ=255)
tCCGCATCCATCCCTTCAGATAAATTAAGCTGTTTTGCGGTTTCTGTGGCGACATTAAGc                   <  1:1226942/60‑1 (MQ=255)
tCCGCATCCATCCCTTCAGATAAATTAAGCTGTTTTGCGGTTTCTGTGGCGACATTAAGc                   <  1:1364080/60‑1 (MQ=255)
tCCGCATCCATCCCTTCAGATAAATTAAGCTGTTTTGCGGTTTCTGTGGCGACATTAAGc                   <  1:1575011/60‑1 (MQ=255)
tCCGCATCCATCCCTTCAGATAAATTAAGCTGTTTTGCGGTTTCTGTGGCGACATTAAGc                   <  1:1590872/60‑1 (MQ=255)
tCCGCATCCATCCCTTCAGATAAATTAAGCTGTTTTGCGGTTTCTGTGGCGACATTAAGc                   <  1:1592763/60‑1 (MQ=255)
tCCGCATCCATCCCTTCAGATAAATTAAGCTGTTTTGCGGTTTCTGTGGCGACATTAAGc                   <  1:1678055/60‑1 (MQ=255)
tCCGCATCCATCCCTTCAGATAAATTAAGCTGTTTTGCGGTTTCTGTGGCGACATTAAGc                   <  1:1701953/60‑1 (MQ=255)
tCCGCATCCATCCCTTCAGATAAATTAAGCTGTTTTGCGGTTTCTGTGGCGACATTAAGc                   <  1:1707083/60‑1 (MQ=255)
tCCGCATCCATCCCTTCAGATAAATTAAGCTGTTTTGCGGTTTCTGTGGCGACATTAAGc                   <  1:1864088/60‑1 (MQ=255)
tCCGCATCCATCCCTTCAGATAAATTAAGCTGTTTTGCGGTTTCTGTGGCGACATTAAGc                   <  1:1894295/60‑1 (MQ=255)
tCCGCATCCATCCCTTCAGATAAATTAAGCTGTTTTGCGGTTTCTGTGGCGACATTAAGc                   <  1:2138464/60‑1 (MQ=255)
tCCGCATCCATCCCTTCAGATAAATTAAGCTGTTTTGCGGTTTCTGTGGCGACATTAAGc                   <  1:2097444/60‑1 (MQ=255)
tCCGCATCCATCCCTTCAGATAAATTAAGCTGTTTTGCGGTTTCTGTGGCGACATTAAGc                   <  1:2090032/60‑1 (MQ=255)
tCCGCATCCATCCCTTCAGATAAATTAAGCTGTTTTGCGGTTTCTGTGGCGACATTAAGc                   <  1:2045647/60‑1 (MQ=255)
tCCGCATCCATCCCTTCAGATAAATTAAGCTGTTTTGCGGTTTCTGTGGCGACATTAAGc                   <  1:1999865/60‑1 (MQ=255)
tCCGCATCCATCCCTTCAGATAAATTAAGCTGTTTTGCGGTTTCTGTGGCGACATTAAGc                   <  1:1968530/60‑1 (MQ=255)
tCCGCATCCATCCCTTCAGATAAATTAAGCTGTTTTGCGGTTTCTGTGGCGACATTAAGc                   <  1:1922506/60‑1 (MQ=255)
tCCGCATCCATCCCTTCAGATAAATTAAGCTGTTTTGCGGTTTCTGTGGCGACATTAAGc                   <  1:1889560/60‑1 (MQ=255)
tCCGCATCCATCCCTTCAGATAAATTAAGCTGTTCTGCGGTTTCTGTGGCGACATTAAGc                   <  1:751929/60‑1 (MQ=255)
               tCAGATAAATTAAGCTGTTTTGCGGTTTCTGTGACGACATTAAGCGGTAACTTCACATCaaa  >  1:147591/1‑62 (MQ=255)
               tCAGATAAATTAAGCTGTTTTGCGGTTTCTGTGACGACATTAAGCGGTAACTTCACATCaaa  >  1:1357367/1‑62 (MQ=255)
               tCAGATAAATTAAGCTGTTTTGCGGTTTCTGTGACGACATTAAGCGGTAACTTCACATCaaa  >  1:2113203/1‑62 (MQ=255)
               tCAGATAAATTAAGCTGTTTTGCGGTTTCTGTGACGACATTAAGCGGTAACTTCACATCaaa  >  1:1201666/1‑62 (MQ=255)
               tCAGATAAATTAAGCTGTTTTGCGGTTTCTGTGACGACATTAAGCGGTAACTTCACATCaaa  >  1:983011/1‑62 (MQ=255)
               tCAGATAAATTAAGCTGTTTTGCGGTTTCTGTGACGACATTAAGCGGTAACTTCACATCaaa  >  1:2202696/1‑62 (MQ=255)
                  gATAAATTAAGCTGTTTTGCGGTTTCTGTGGCGACATTAAGc                   <  1:375075/42‑1 (MQ=255)
                                                |                            
TCCGCATCCATTCCTTCAGATAAATTAAGCTGTTTTGCGGTTTCTGTGACGACATTAAGCGGTAACTTCACATCAAA  >  W3110S.gb/3591242‑3591318

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: