Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 834,218 | T→C | 23.4% | D892D (GAT→GAC) | narG → | nitrate reductase 1, alpha subunit |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 834,218 | 0 | T | C | 23.4% | 110.6 / 31.4 | 64 | D892D (GAT→GAC) | narG | nitrate reductase 1, alpha subunit |
Reads supporting (aligned to +/- strand): ref base T (26/23); new base C (6/9); total (32/32) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.56e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
TGGAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAG > minE/834150‑834288 | tGGAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTa > 1:2027666/1‑71 (MQ=255) gcaaaGGCGAATGCGATCTGATTCCAGGCAAAACCGCGCCGAATATTGTGGCGGTGGAGCGCGACTAccct > 1:2952907/3‑70 (MQ=255) aaaaGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCgg > 1:220715/1‑71 (MQ=255) aaGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTAt < 1:5422594/65‑1 (MQ=255) ggCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGa > 1:1320910/1‑70 (MQ=255) gCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACt < 1:2612768/71‑1 (MQ=255) gCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACt < 1:2308008/71‑1 (MQ=255) aGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTAc < 1:1073281/71‑1 (MQ=255) aGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTAc < 1:4876532/71‑1 (MQ=255) gTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTa > 1:5092015/1‑69 (MQ=255) gTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTAc > 1:1805810/1‑70 (MQ=255) atcTGATTCCAGGCAAAACCGCGCCGAATATTGTGGCGGTGGAGCGCGACTACCCTGCTACGTATGAAc < 1:1171900/67‑1 (MQ=255) atcTGATTCCAGGCAAAACCGCGCCGAATATTGTGGCGGTGGAGCGCGACTACCCTGCTACGTATGAAc < 1:1531459/67‑1 (MQ=255) atcTGATTCCAGGCAAAACCGCGCCGAATATTGTGGCGGTGGAGCGCGACTACCCTGCTACGTATGAAc < 1:1759713/67‑1 (MQ=255) tGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATcc < 1:230778/52‑1 (MQ=255) tGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAAc < 1:3983531/66‑1 (MQ=255) aTCCCAGGTAAAACCGCGCCACACATTATGGTCGTATAGCGCGATTATCCGGCGACTTACGAACGCtt < 1:1439717/68‑1 (MQ=255) aTCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCtt < 1:1762731/68‑1 (MQ=255) aTCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTAc > 1:1664951/1‑71 (MQ=255) aTCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTAc > 1:3137129/1‑71 (MQ=255) cAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACctct < 1:2961830/71‑1 (MQ=255) cAGGCAAAACCGCGCCGAATATTGTGGCGGTGGAGCGCGACTacc > 1:3652708/1‑43 (MQ=255) aGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTa > 1:4131261/1‑71 (MQ=255) aGGCAAAACCGCGCCGAATATTGTGGCGGTGGAGCGCGACTACCCTGCTACGTATGAACGc < 1:2274748/61‑1 (MQ=255) aGGCAAAACCGCGCCGAATATTGTGGCGGTGGAGCGCGACTACCCTGCTACGTATGAACGCTTTACctcgc < 1:1793489/71‑3 (MQ=255) ggTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGaa > 1:4401786/1‑57 (MQ=255) gTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCttt < 1:3205612/62‑1 (MQ=255) aaaaCCTCGCCGAATATTGTGGCGGTGGAGCGCGACTACCCTGCTACGTATGAACGCTTTACCTCGCTCgg > 1:3373877/1‑71 (MQ=255) aaaaCCGCGCCGAATATTGTGGCGGTGGAGCGCGACTACCCTGCTACGTATGAACGCTTTACCTCGCTCgg > 1:4059875/1‑71 (MQ=255) aaaaCCGCGCCGAATATTGTGGCGGTGGAGCGCGACTACCCTGCTACGTATGAACGCTTTACCTCGCTCgg > 1:3264816/1‑71 (MQ=255) aaaaCCGCGCCGAATATTGTGGCGGTGGAGCGCGACTACCCTGCTACGTATGAACGCTTTACCTCGCTCgg < 1:537943/71‑1 (MQ=255) aaaaCCGCGCCGAATATTGTGGCGGTGGAGCGCGACTACCCTGCTACGTATGAACGCTTTACCTCGCTCgg > 1:2028724/1‑71 (MQ=255) gcgcCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTAc > 1:2475642/1‑56 (MQ=255) gcgcCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTAc > 1:2682170/1‑56 (MQ=255) gcgcCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCt < 1:151129/71‑1 (MQ=255) cacacaTTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTAc > 1:1043999/1‑52 (MQ=255) cacaTTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACctc < 1:5323215/53‑1 (MQ=255) acaTTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACctct < 1:38355/53‑1 (MQ=255) acaTTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATg > 1:2989700/1‑68 (MQ=255) caTTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTa > 1:3775983/1‑53 (MQ=255) aTTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGaa > 1:4049950/1‑39 (MQ=255) aTTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGaa > 1:2117296/1‑39 (MQ=255) aTTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGaa > 1:2688210/1‑39 (MQ=255) aTTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGa > 1:1524131/1‑38 (MQ=255) aTTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAAc < 1:4409047/40‑1 (MQ=255) aTTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTAc > 1:1337985/1‑47 (MQ=255) aTTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTAc > 1:4700031/1‑47 (MQ=255) aTTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTAc > 1:4276509/1‑47 (MQ=255) ttATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGaaa < 1:276779/71‑1 (MQ=255) tGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTAt > 1:1417611/1‑49 (MQ=255) gTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCgg < 1:512139/71‑1 (MQ=255) gTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCgg < 1:2863443/71‑1 (MQ=255) tCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGGGAAAATCGGt < 1:918393/71‑1 (MQ=255) tCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGt < 1:4618116/71‑1 (MQ=255) tCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGt < 1:2213556/71‑1 (MQ=255) gTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTaa > 1:2423121/1‑71 (MQ=255) agagCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATg > 1:2292926/1‑71 (MQ=255) agagCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATg > 1:1718139/1‑71 (MQ=255) gagCGCGACTACCCTGCTACGTATGAACGCTTTACCTCGCTCGGGCCATTGATGGACAAACTTGGCAATgg < 1:2770130/71‑1 (MQ=255) gagCGCGACTACCCTGCTACGTATGAACGCTTTACCTCGCTCGGGCCATTGATGGACAAACTTGGCAATgg < 1:2331647/71‑1 (MQ=255) agCGCGACTACCCTGCTACGTATGAACGCTTTACCTCGCTCGGGCCATTGATGGACAAACTTGGCAATGGc < 1:2732854/71‑1 (MQ=255) cgcgATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAAt < 1:2180099/59‑1 (MQ=255) cgATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGt > 1:2191912/1‑70 (MQ=255) aTTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCgg < 1:3378964/58‑1 (MQ=255) ttATCCGGCGACTTACGAACGCTTTACCTCTATCGGc > 1:3829941/1‑37 (MQ=255) ttATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTaa > 1:5053776/1‑60 (MQ=255) ttATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGt > 1:2175296/1‑67 (MQ=255) ttATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAg > 1:4055910/1‑71 (MQ=255) | TGGAAAAAAGGCGAGTGCGACCTGATCCCAGGTAAAACCGCGCCACACATTATGGTCGTAGAGCGCGATTATCCGGCGACTTACGAACGCTTTACCTCTATCGGCCCGCTGATGGAGAAAATCGGTAATGGCGGTAAAG > minE/834150‑834288 |
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: ‑ |