| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | W3110S.gb | 542,820 | C→T | 28.9% | A151T (GCC→ACC) | ylbA ← | conserved hypothetical protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | W3110S.gb | 542,820 | 0 | C | T | 28.9% | 56.4 / 26.5 | 38 | A151T (GCC→ACC) | ylbA | conserved hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base C (14/13); new base T (4/7); total (18/20) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.85e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01 | |||||||||||
GGCAGAAAATCCAGCAGAATAACATCGTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATAAA > W3110S.gb/542756‑542881 | ggCAGAAAATCCAGCAGAATAACATCGTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGTATTGc < 1:105704/70‑1 (MQ=255) aGAAAATCCAGCAGAATAACATCGTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGa < 1:1035107/71‑1 (MQ=255) aGAAAATCCAGCAGAATAACATCGTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGa < 1:2491147/71‑1 (MQ=255) tCCAGCAGAATAACATCGTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGTATTGCCAGaa < 1:957885/66‑1 (MQ=255) gcagAATAACATCGTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCAc < 1:3573489/71‑1 (MQ=255) gAATAACATCGTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGcc < 1:3088299/66‑1 (MQ=255) aaTAACATCGTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGCATTGcc < 1:2873390/54‑1 (MQ=255) tAACATCGTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCACGGTGCa < 1:1369724/71‑1 (MQ=255) gTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGTATTGGCAGAAACCAGCCACGGTGCATAGCCtt > 1:3599243/1‑71 (MQ=255) gTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGTATTGCCAGAAACCAGCCACGGTGCATAGCCt > 1:3247147/1‑70 (MQ=255) gTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGTATTGCCAGAAACCAGCCACGGTGCATAGCCt > 1:2638671/1‑70 (MQ=255) gTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGCa > 1:109631/1‑40 (MQ=255) gTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGCa > 1:2041737/1‑40 (MQ=255) gTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCACGGTGCATAGCCt < 1:2165368/70‑1 (MQ=255) tCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCAc < 1:240939/57‑1 (MQ=255) tCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCAc < 1:2680881/57‑1 (MQ=255) ccATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCACGGTGCATAGCCt < 1:142913/68‑1 (MQ=255) cATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTa > 1:3261781/1‑71 (MQ=255) ccTTCATAATGAATGCGTTCCAGTTCGCTGGTATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTACCg > 1:3684328/1‑70 (MQ=255) ccTTCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTACCg > 1:2497230/1‑70 (MQ=255) ccTTCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTACCg > 1:3739280/1‑70 (MQ=255) ttCATAATGAATGCGTTCCAGTTCGCTGGTATTGCCAGAAACCAGCCAc < 1:79395/49‑1 (MQ=255) ttCATAATGAATGCGTTCCAGTTCGCTGGTATTGCAAg < 1:3433047/38‑1 (MQ=25) cATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTACCGGAAc > 1:2870950/1‑70 (MQ=255) cATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTACCGGAACa > 1:947230/1‑71 (MQ=255) aTAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTACCGGAACAt < 1:2512782/71‑1 (MQ=255) atgaatgCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTACCGGAACATa > 1:3308099/1‑69 (MQ=255) gTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTACCGGAACATAGCGGCGCTta < 1:1654511/71‑1 (MQ=255) ttCCAGTTCGCTGGTATTGCCAGAAACCAGCCACGGTGCATAGCCTTCt < 1:2376682/49‑1 (MQ=255) ttCCAGTTCGCTGGTATTGCCAGAAACCAGCCACGGTGCATAGCCTTCt < 1:2180259/49‑1 (MQ=255) ttCCAGTTCGCTGGCATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTACCGGAACATAGCGGCGCTta > 1:2757560/1‑70 (MQ=255) ccAGTTCGCTGGCATTGCCAGAAACCAGCCACGTTGCATAGCCTTCTACCGGAACATAGCGGCGCTtata > 1:1956680/1‑70 (MQ=255) gTTCGCTGGCATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATaaa < 1:2787966/71‑1 (MQ=255) ttCGCTGGTATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTACCGGAACATAGCGGCGCt < 1:3693086/62‑1 (MQ=255) cGCTGGCATTGCCAGAAACCAGCCACGGTGCATAGCCtt > 1:2588730/1‑39 (MQ=255) cGCTGGCATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTACCGGAACATAgcg > 1:2685938/1‑55 (MQ=255) cGCTGGCATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTACCGGAACATAgcg > 1:2534886/1‑55 (MQ=255) cGCTGGCATTGCCAGAAACCAGCCACGGTGCATAGACTTCTACCGGAAc > 1:1864708/1‑49 (MQ=255) | GGCAGAAAATCCAGCAGAATAACATCGTCCATGCCTTCATAATGAATGCGTTCCAGTTCGCTGGCATTGCCAGAAACCAGCCACGGTGCATAGCCTTCTACCGGAACATAGCGGCGCTTATATAAA > W3110S.gb/542756‑542881 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |