Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 834,248 | T→G | 14.3% | S902S (TCT→TCG) | narG → | nitrate reductase 1, alpha subunit |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 834,248 | 0 | T | G | 14.3% | 81.5 / 12.2 | 42 | S902S (TCT→TCG) | narG | nitrate reductase 1, alpha subunit |
Reads supporting (aligned to +/- strand): ref base T (26/10); new base G (4/2); total (30/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01 |
AGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATG > minE/834179‑834317 | aGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTa > 1:378712/1‑71 (MQ=255) aGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTa > 1:3951602/1‑71 (MQ=255) ggTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTAt > 1:3084154/1‑71 (MQ=255) aaaaCCGCGCCGAATATTGTGGCGGTGGAGCGCGACTACCCTGCTACGTATGAACGCTTTACCTCGCTCgg < 1:19573/71‑1 (MQ=255) aaaaCCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCgg > 1:2444626/1‑71 (MQ=255) aaaCCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGc > 1:1980575/1‑71 (MQ=255) aaaCCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGc > 1:186387/1‑71 (MQ=255) aaCCGCGCCGAATATTGTGGCGGTGGAGCGCGACTACCCTGCTACGTATGAACGCTTTACCTCGCTCgggc > 1:1391442/1‑69 (MQ=255) aCCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGccc < 1:2561972/71‑1 (MQ=255) cgcgCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCg > 1:1932590/1‑70 (MQ=255) gcgcCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCt > 1:1281191/1‑71 (MQ=255) cacacaTTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGAt < 1:3158329/70‑1 (MQ=255) cacaTTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATg > 1:2937784/1‑69 (MQ=255) aTGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAAt > 1:320999/1‑71 (MQ=255) ggTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACctct < 1:3806373/46‑1 (MQ=255) agagCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATg > 1:603530/1‑71 (MQ=255) agagCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATg > 1:3275776/1‑71 (MQ=255) gagCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATgg < 1:3911881/71‑1 (MQ=255) agCGCGACTACCCTGCTACGTATGAACGCTTTACCTCGCTCGGGCCATTGATGGACAAACTTGGCAATGGc > 1:1367208/1‑71 (MQ=255) cgcgATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCgg > 1:2996914/1‑71 (MQ=255) cgcgACTACCCTGCTACGTATGAACGCTTTACCTCGCTCGGGCCATTGATg > 1:3911127/1‑51 (MQ=255) cgATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGt > 1:558700/1‑70 (MQ=255) gATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTaa < 1:1628483/71‑1 (MQ=255) gATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTa > 1:1968555/1‑70 (MQ=255) ttATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGa > 1:2904507/1‑44 (MQ=255) ttATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGa > 1:2842278/1‑44 (MQ=255) ttATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATg > 1:102027/1‑46 (MQ=255) ttATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATg > 1:132678/1‑46 (MQ=255) ttATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTaa > 1:306251/1‑60 (MQ=255) ttATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTaa > 1:2718844/1‑60 (MQ=255) aTCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGa > 1:75292/1‑42 (MQ=255) gCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGc < 1:3776216/71‑1 (MQ=255) cGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGcc < 1:103031/71‑1 (MQ=255) cTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTaa > 1:1922197/1‑49 (MQ=255) cTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTgg < 1:746364/71‑1 (MQ=255) ttACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTaa > 1:2619149/1‑48 (MQ=255) aaCGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTaaa < 1:3743074/53‑1 (MQ=255) aaCGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTaaa < 1:890923/53‑1 (MQ=255) cGCTTTACCTCGCTCGGGCCATTGATGGACAAACTTGGCAATGGCGGTAAAGGGAttt < 1:1136861/58‑2 (MQ=255) cTTTACCTCTATCGGCCCGCTGATGGAGAATATCGGTAATGGCGGTAAAGGGATTGcc < 1:570014/58‑1 (MQ=255) ttACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGc > 1:3054743/1‑71 (MQ=255) tACCTCGCTCGGGCCATTGATGGACAAACTTGGCAATGGCGGTAAAGGGATTTCGTGGAATACCCaggat > 1:897088/1‑67 (MQ=255) ccTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGGCTGGaa > 1:3735170/1‑58 (MQ=255) ctATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATg > 1:911210/1‑71 (MQ=255) | AGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAGGGATTGCCTGGAACACCCAGAGCGAGATG > minE/834179‑834317 |
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: ‑ |