| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 836,317 | G→T | 15.3% | P345P (CCG→CCT) | 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,317 | 0 | G | T | 15.3% | 114.9 / 6.7 | 39 | P345P (CCG→CCT) | narH | nitrate reductase 1, beta (Fe‑S) subunit |
| Reads supporting (aligned to +/- strand): ref base G (23/10); new base T (3/3); total (26/13) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.80e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.70e-01 | |||||||||||
GAAGCCGCACAGCAGTCGCCGGTTTATAAAATGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGAATATCGCACACTGCCGATGGTCTGGTACGTGCCGCCTCTGTCTCCGATTCAGTCTGCAGCAGACGCG > minE/836249‑836386 | gAAGCCGCACAGCAGTCGCCGGTTTATAAAATGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCgg < 1:581930/70‑1 (MQ=255)gAAGCCGCACAGCAGTCGCCGGTTTATAAAATGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCgg < 1:3419046/70‑1 (MQ=255) aGCCGCACAGCAGTCGCCGGTTTATAAAATGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGAat < 1:3188699/71‑1 (MQ=255) cacaGCAGTCGCCGGTTTATAAAATGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGAATATCg > 1:1272532/1‑70 (MQ=255) gcagTCGCCGGTTTATAAAATGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCg > 1:1722871/1‑58 (MQ=255) agTCGCCGGTTTATAAAATGGCAATGGAATGGAAACTGGCGCTGACGCTGCATCCGGAATATCGCACACTg < 1:3267119/71‑1 (MQ=255) tCGCCGGTTTATAAAATGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGAATATCGCACACTGcc < 1:2447661/71‑1 (MQ=255) cGCCAGTCTACAAAATGGCGATGGACTGGAAACTGGCACTACCGTTGCACCCTGAatac > 1:47812/1‑58 (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:2688942/1‑71 (MQ=255) ttATAAAATGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGAATATCGCACACTGCCGATGGTCt > 1:1026838/1‑71 (MQ=255) ttATAAAATGGCAATGGAATGGAAACTGGCGCTGACGCTGCATCCGGa > 1:3479611/1‑48 (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:321505/1‑70 (MQ=255) gCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGAATATCGCACACTGCCGATGGTCTGGTACGTgc > 1:423653/1‑70 (MQ=255) atggaatggaAACTGGCGCTGCCGCTGCATCCGGAATATCGCACACTGCCGATGGTCTGGTACGTgccgcc < 1:3299052/71‑1 (MQ=255) tggaatggaaACTGGCGCTGCCGCTGCATCCGGaa > 1:573187/1‑35 (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) | GAAGCCGCACAGCAGTCGCCGGTTTATAAAATGGCAATGGAATGGAAACTGGCGCTGCCGCTGCATCCGGAATATCGCACACTGCCGATGGTCTGGTACGTGCCGCCTCTGTCTCCGATTCAGTCTGCAGCAGACGCG > minE/836249‑836386 |
| 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: ‑ |