| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | W3110S.gb | 4,292,897 | C→T | L21L (CTG→TTG) | nrfB → | nitrite reductase, formate‑dependent, penta‑heme cytochrome c |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | W3110S.gb | 4,292,897 | 0 | C | T | 94.3% | 100.6 / ‑1.0 | 35 | L21L (CTG→TTG) | nrfB | nitrite reductase, formate‑dependent, penta‑heme cytochrome c |
| Reads supporting (aligned to +/- strand): ref base C (2/0); new base T (14/19); total (16/19) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.02e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.83e-01 | |||||||||||
TGAGCGTATTACGTTCGTTATTAACTGCCGGGGTGCTGGCGTCGGGCCTGTTGTGGAGCCTGAACGGAATTACCGCTACCCCTGCGGCGCAGGCATCTGATGATCGCTACGAAGTTACCCAG > W3110S.gb/4292838‑4292959 | tGAGCGTATTACGTTCGTTATTAACTGCCGGGGTGCTGGCGTCGGGCCTGTTGTGGAGCTTGAACGGAAt > 1:2571586/1‑70 (MQ=255)tGAGCGTATTACGTTCGTTATTAACTGCCGGGGTGCTGGCGTCGGGCCTGTTGTGGAGCTTGAACGGAAt > 1:2494375/1‑70 (MQ=255)tGAGCGTATTACGTTCGTTATTAACTGCCGGGGTGCTGGCGTCGGGCCTGTTGTGGAGCTTGAACGGAAt > 1:2044056/1‑70 (MQ=255) aGCGTATTACGTTCGTTATTAACTGCCGGGGTGCTGGCGTCGGGCCTGTTGTGGAGCTTGAACGGAATTAc > 1:867946/1‑71 (MQ=255) aGCGTATTACGTTCGTTATTAACTGCCGGGGTGCTGGCGTCGGGCCTGTTGTGGAGCTTGAACGGAATTAc < 1:3065766/71‑1 (MQ=255) gCGTATTACGTTCGTTATTAACTGCCGGGGTGCTGGCGTCGGGCCTGTTGTGGAGCTTGAACGGAATTAcc < 1:3157949/71‑1 (MQ=255) aCGTTCGTTATTAACTGCCGGGGTGCTGGCGTCGGGCCTGTTGTGGAGCCTGAACGGAATTACCGCTAccc > 1:1949622/1‑71 (MQ=255) aCGTTCGTTATTAACTGCCGGGGTGCTGGCGTCGGGCCTGTTGTGGAGCCTGAACGGAATTACCGCTAccc > 1:3443559/1‑71 (MQ=255) gttATTAACTGCCGGGGTGCTGGCGTCGGGCCTGTTGTGGAGCTTGAACGGAATTACCGCTACCCCTgcgg < 1:1731521/71‑1 (MQ=255) gttATTAACTGCCGGGGTGCTGGCGTCGGGCCTGTTGTGGAGCTTGAACGGAATTACCGCTACCCCTgcgg < 1:1799015/71‑1 (MQ=255) gttATTAACTGCCGGGGTGCTGGCGTCGGGCCTGTTGTGGAGCTTGAACGGAATTACCGCTACCCCTgcgg < 1:1903948/71‑1 (MQ=255) gttATTAACTGCCGGGGTGCTGGCGTCGGGCCTGTTGTGGAGCTTGAACGGAATTACCGCTACCCCTgcgg < 1:1301494/71‑1 (MQ=255) cGGGGTGCTGGCGTCGGGCCTGTTGTGGAGCTTGAACGGAATTACCGCTACCCCTGCGGCGCAGGCATCtg < 1:1528607/71‑1 (MQ=255) cGGGGTGCTGGCGTCGGGCCTGTTGTGGAGCTTGAACGGAATTACCGCTACCCCTGCGGCGCAGGCATCtg < 1:1780319/71‑1 (MQ=255) cGGGGTGCTGGCGTCGGGCCTGTTGTGGAGCTTGAACGGAATTACCGCTACCCCTGCGGCGCAGGCATCtg < 1:3905422/71‑1 (MQ=255) cGGGGTGCTGGCGTCGGGCCTGTTGTGGAGCTTGAACGGAATTACCGCTACCCCTGCGGCGCAGGCATCtg < 1:1523734/71‑1 (MQ=255) cGGGGTGCTGGCGTCGGGCCTGTTGTGGAGCTTGAACGGAATTACCGCTACCCCTGCGGCGCAGGCATCtg < 1:2980971/71‑1 (MQ=255) cGGGGTGCTGGCGTCGGGCCTGTTGTGGAGCTTGAACGGAATTACCGCTACCCCTGCGGCGCAGGCATCtg < 1:2837603/71‑1 (MQ=255) cGGGGTGCTGGCGTCGGGCCTGTTGTGGAGCTTGAACGGAATTACCGCTACCCCTGCGGAGCAGGCATCtg < 1:2871511/71‑1 (MQ=255) cGGGGTGCGGGCGTCGGGCCTGTTGTGGAGCTTGAACGGAATTACCGCTACCCCTGCGGCGCAGGCATCtg < 1:1040011/71‑1 (MQ=255) ggggTGCTGGCGTCGGGCCTGTTGTGGAGCTTGAACGGAATTACCGCTACCCCTGCGGCGCAGGCATCtga > 1:4047675/1‑71 (MQ=255) ggggTGCTGGCGTCGGGCCTGTTGTGGAGCTTGAACGGAATTACCGCTACCCCTGCGGCGCAGGCATCtg > 1:769256/1‑70 (MQ=255) ggggTGCTGGCGTCGGGCCTGTTGTGGAGCTTGAACGGAATTACCGCTACCCCTGCGGCCCAGGCATCtg < 1:1252661/70‑1 (MQ=255) gCTGGCGTCGGGCCTGTTGTGGAGCTTGAACGGAATTACCGc > 1:4076867/1‑42 (MQ=255) gCTGGCGTCGGGCCTGTTGTGGAGCTTGAACGGAATTACCGCTACCCCTGCGGCGCAGGCATCTgatgat > 1:1690720/1‑70 (MQ=255) gCTGGCGTCGGGCCTGTTGTGGAGCTTGAACGGAATTACCGCTACCCCTGCGGCGCAGGCATCTGATGATc > 1:2571561/1‑71 (MQ=255) gCTGGCGTCGGGCCTGTTGTGGAGCTTGAACGGAATTACCGCTACCCCTGCGGCGCAGGCATCTGATGATc > 1:3563303/1‑71 (MQ=255) gCTGGCGTCGGGCCTGTTGTGGAGCTTGAACGGAATTACCGCTACCCCTGCGGCGCAGGCATCTGATGATc > 1:1783528/1‑71 (MQ=255) gCGTCGGGCCTGTTGTGGAGCTTGAACGGAATTACCGCTACCCCTGCGGCGCAGGCAt < 1:750909/58‑1 (MQ=255) gggCCTGTTGTGGAGCTTGAACGGAATTACCGCTACCCCTGCGGCGCAGGCATCTGATGATCGCTACGaa < 1:3463407/70‑1 (MQ=255) gggCCTGTTGTGGAGCTTGAACGGAATTACCGCTACCCCTGCGGCGCAGGCATCTGATGATCGCTACGaa < 1:51804/70‑1 (MQ=255) gggCCTGTTGTGGAGCTTGAACGGAATTACCGCTACCCCTGCGGCGCAGGCATCTGATGATCGCTACGaa < 1:731258/70‑1 (MQ=255) gCCTGTTGTGGAGCTTGAACGGAATTACCGCTACCCCTGCGGCGCAGGCATCTGATGATCGCTACGAAGt > 1:3658261/1‑70 (MQ=255) gtgGAGCTTGAACGGAATTACCGCTACCCCTGCGGCGCAGGCATCTGATGATCGCTACGAAGTTACCcag > 1:1487973/1‑70 (MQ=255) gtgGAGCTTGAACGGAATTACCGCTACCCCTGCGGCGCAGGCATCTGATGATCGCTACGAAGTTACCcag > 1:2163334/1‑70 (MQ=255) | TGAGCGTATTACGTTCGTTATTAACTGCCGGGGTGCTGGCGTCGGGCCTGTTGTGGAGCCTGAACGGAATTACCGCTACCCCTGCGGCGCAGGCATCTGATGATCGCTACGAAGTTACCCAG > W3110S.gb/4292838‑4292959 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |