| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 836,323 | T→C | 13.9% | Y347Y (TAT→TAC) | narH → | nitrate reductase 1, beta (Fe‑S) subunit |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 836,323 | 0 | T | C | 13.9% | 90.5 / 7.7 | 36 | Y347Y (TAT→TAC) | narH | nitrate reductase 1, beta (Fe‑S) subunit |
| Reads supporting (aligned to +/- strand): ref base T (21/10); new base C (2/3); total (23/13) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.28e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
CACAGCAGTCGCCGGTTTATAAAATGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGAATATCGCACACTGCCGATGGTCTGGTACGTGCCGCCTCTGTCTCCGATTCAGTCTGCAGCAGACGCGGGT > minE/836256‑836389 | cacaGCAGTCGCCGGTTTATAAAATGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGAATATCg > 1:1272532/1‑70 (MQ=255) agTCGCCGGTTTATAAAATGGCAATGGAATGGAAACTGGCGCTGACGCTGCATCCGGAATATCGCACACTg < 1:3267119/71‑1 (MQ=255) tCGCCGGTTTATAAAATGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGAATATCGCACACTGcc < 1:2447661/71‑1 (MQ=255) gCCGGTTTATAAAATGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGAATATc > 1:3258717/1‑59 (MQ=255) gCCGGTTTATAAAATGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGAATATCGc > 1:629750/1‑61 (MQ=255) gTTTATAAAATGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGAATATCGca < 1:3130053/58‑1 (MQ=255) ttATAAAATGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGAATATCGca > 1:3518329/1‑56 (MQ=255) ttATAAAATGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGAATATCGCACACTGCCGATGGTCt > 1:1026838/1‑71 (MQ=255) ttATAAAATGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGAATATCGCACACTGCCGATGGTCt > 1:2688942/1‑71 (MQ=255) aaaTGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGAATATc > 1:3387021/1‑48 (MQ=255) aaaTGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGAATATCGCACACTGCCGATGGTCTGGTa < 1:2336277/70‑1 (MQ=255) aaTGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGAATATCg < 1:3333269/48‑1 (MQ=255) aaTGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGAATATCGCACACTGCCGATGGTCTGGTAc < 1:413084/70‑1 (MQ=255) gCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGAATATCGca > 1:1122780/1‑46 (MQ=255) gCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGAATATCGCACACTGCCGATGGTCTGGTACGTgc > 1:423653/1‑70 (MQ=255) gCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGAATATCGCACACTGCCGATGGTCTGGTACGTgc > 1:321505/1‑70 (MQ=255) atggaatggaAACTGGCGCTGCCGCTGCATCCGGAATATCGCACACTGCCGATGGTCTGGTACGTgccgcc < 1:3299052/71‑1 (MQ=255) gaatggaaACTGGCGCTGCCGCTGCATCCGGAATATc > 1:2523542/1‑37 (MQ=255) gaatggaaACTGGCGCTGCCGCTGCATCCGGAATATCGCACACTGCCGCTGGTCTGGTACGTGCCGCctct > 1:312905/1‑71 (MQ=255) gaatggaaACTGGCGCTGCCGCTGCATCCGGAATATCGCACACTGCCGATGGTCTGGTACGTgcc > 1:3808693/1‑65 (MQ=255) gaatggaaACTGGCGCTGCCGCTGCATCCGGAATATCGCACACTGCCGATGGTCTGGTACGTg > 1:2426313/1‑63 (MQ=255) gaatggaaACTGGCGCTGCCGCTGCATCCGGAATATCGCACACTGCCGATGGTCTGGTACGTGCCGCctct > 1:1394244/1‑71 (MQ=255) tggaaAGTGGCGCTGCCGCTGCATCCGGAATATCGCACACTGCCGATGGTCTGGTACGTGCCGCct < 1:3051941/66‑1 (MQ=255) aaaCTGGCGCTGCCGCTGCATCCGGAATATCGCACACTGcc > 1:2818521/1‑41 (MQ=255) aCTGGCACTACCGTTGCACCCTGAATACCGCACCCTGCCGATGGTCTGGTACGTTCCTCCGCTGt > 1:1250512/1‑65 (MQ=255) tGGCGCTGCCGCTGCATCCGGAATATCGCACACTGCCGATGGTCTGGTACGTg > 1:1917899/1‑53 (MQ=255) tGGCGCTGCCGCTGCATCCGGAATATCGCACACTGCCGATGGTCTGGTACGTGCCGCCTCTGTCTcc < 1:131169/67‑1 (MQ=255) ctaccgttgcACCCTGAATACCGCACCCTGCCGATGGTCTGGTACGTTCCTCCGCTGTCACCGATTCAg < 1:119849/66‑1 (MQ=255) ctaccgttgcACCCTGAATACCGCACCCTGCCGATGGTCTGGTACGTTCCTCCGCTGTCACCGATTCAGTc > 1:3446171/4‑71 (MQ=255) gccgctgcATCCGGAATATCGCACACTGCCGATGGTCTGGTAc > 1:94659/1‑43 (MQ=255) ccgttgcACCCTGAATACCGCACCCTGCCGATGGTCTGGTACGTTCCTCCGCTGTCACCGATTCAGTcct < 1:3321228/70‑3 (MQ=255) ttgcACCCTGAATACCGCACCCTGCCGATGGTCTGGTa < 1:819627/37‑1 (MQ=19) aTCCGGAATATCGCACACTGCCGATGGTCTGGTACGTGCCGCCTCTGt > 1:3156697/1‑48 (MQ=255) cGGAATATCGCACACTGCCGATGGTCTGGTACGTGCCGCCTCTGTCTCCGATTCAGTCTGCAGCAGa > 1:2144561/1‑67 (MQ=255) cGGAATATCGCACACTGCCGATGGTCTGGTACGTGCCGCCTCTGTCTCCGATTCAGTCTGCAGCAGAcgcg > 1:2219229/1‑71 (MQ=255) aaTATCGCACACTGCCGATGGTCTGGTACGTGCCGCCTCTGTCTCCGATTCAGTCTGCAGCAGACGCGGGt < 1:2874373/71‑1 (MQ=255) tatCGCACACTGCCGATGGTCTGGTACGTGCCGCCTCTGTCTCCGATTCAGTCTgcagca > 1:11336/1‑60 (MQ=255) atCGCACACTGCCGATGGTCTGGTACGTGCCGCCTCTGTCTCCGATTCAGTCTGCAGCAGAcgc < 1:1494690/64‑1 (MQ=255) tCGCACACTGCCGATGGTCTGGTACGTGCCGCCTCTGTCTCCGa > 1:2768524/1‑44 (MQ=255) tCGCACACTGCCGATGGTCTGGTACGTGCCGCCTCTGTCTCCGATTCAGTCTgcag > 1:2857568/1‑56 (MQ=255) | CACAGCAGTCGCCGGTTTATAAAATGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGAATATCGCACACTGCCGATGGTCTGGTACGTGCCGCCTCTGTCTCCGATTCAGTCTGCAGCAGACGCGGGT > minE/836256‑836389 |
| 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: ‑ |