| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | W3110S.gb | 1,068,568 | T→A | 12.2% | G22G (GGT→GGA) | ymdF → | conserved hypothetical protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | W3110S.gb | 1,068,568 | 0 | T | A | 12.2% | 106.8 / 8.0 | 41 | G22G (GGT→GGA) | ymdF | conserved hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base T (16/20); new base A (3/2); total (19/22) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.49e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01 | |||||||||||
AATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAG > W3110S.gb/1068502‑1068629 | aaTGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTaaa > 1:2498663/1‑70 (MQ=255)aaTGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTaaa > 1:3341909/1‑70 (MQ=255)aaTGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTaaa > 1:1030075/1‑70 (MQ=255) aTGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTaaaaa < 1:1036504/71‑1 (MQ=255) aTGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTaaaaa < 1:3050033/71‑1 (MQ=255) aTGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTaaaaa < 1:3659338/71‑1 (MQ=255) aTGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTaaaa > 1:1381107/1‑70 (MQ=255) aTGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTaaaa > 1:2992961/1‑70 (MQ=255) aTGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTaaaa < 1:3596807/70‑1 (MQ=255) aTGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTaaaa < 1:2014579/70‑1 (MQ=255) aTGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTaaaa > 1:1861953/1‑70 (MQ=255) aTGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTaaaa < 1:1726422/70‑1 (MQ=255) aTGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTaaaa < 1:1487587/70‑1 (MQ=255) aTGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACAGCGAAAGAGCATCAGAAGCAGGTaaaa > 1:2179232/1‑70 (MQ=255) gCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAgg > 1:3098994/1‑71 (MQ=255) gCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAAAAAAAgg > 1:1131980/1‑71 (MQ=255) gCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAAAAAAAgg > 1:3640565/1‑71 (MQ=255) gCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAAAAAAAgg > 1:1526603/1‑71 (MQ=255) aaaCCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTg < 1:1801674/71‑1 (MQ=255) aaaCCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTg < 1:2335236/71‑1 (MQ=255) aaaCCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTg < 1:320458/71‑1 (MQ=255) aaaCCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTg < 1:2950151/71‑1 (MQ=255) aaaCCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAAAAAAAGGTg < 1:1119202/71‑1 (MQ=255) aaaCCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGAAAAAAAGGTg < 1:2989194/71‑1 (MQ=255) aCCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGa < 1:1398417/71‑1 (MQ=255) tCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGacag > 1:1114587/1‑70 (MQ=255) tacGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAgg > 1:3382852/3‑53 (MQ=255) tCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAgg > 1:3318549/1‑53 (MQ=255) tCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAgg > 1:926850/1‑53 (MQ=255) tCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCggggg > 1:1946820/1‑71 (MQ=255) cGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGt < 1:3063523/70‑1 (MQ=255) aTTTTGCAGAAGACCGCGAAAGAGCATCAGAATCAGGTaa < 1:3587095/40‑1 (MQ=255) aTTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTaa < 1:1706404/40‑1 (MQ=255) tGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAg < 1:2600997/70‑1 (MQ=255) agaagaCCGCGAAAGAGCATCAGAAGCAGCTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATg < 1:3019230/71‑1 (MQ=255) cgcgAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGc < 1:2452195/70‑1 (MQ=255) tCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGa < 1:2818361/71‑1 (MQ=255) aGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAg > 1:2355498/1‑71 (MQ=255) aGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAg > 1:779405/1‑71 (MQ=255) aGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAg > 1:1128373/1‑71 (MQ=255) gTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAg < 1:2378598/38‑1 (MQ=255) | AATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAG > W3110S.gb/1068502‑1068629 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |