| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 835,291 | T→C | 26.8% | I3I (ATT→ATC) | 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 | 835,291 | 0 | T | C | 26.8% | 128.8 / 54.9 | 82 | I3I (ATT→ATC) | narH | nitrate reductase 1, beta (Fe‑S) subunit |
| Reads supporting (aligned to +/- strand): ref base T (27/33); new base C (9/13); total (36/46) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 8.05e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01 | |||||||||||
AGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCT > minE/835222‑835358 | aGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAATTc > 1:614527/1‑71 (MQ=255)aGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAATTc < 1:1001081/71‑1 (MQ=255)aGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAATTc > 1:3828491/1‑71 (MQ=255) gATGAAGAACGTTAACTGGCTGGATGATGAAGGTCGCGATCAGGTACAGGAG‑GCGAA‑AAAATGAAAATCCg < 1:1923858/71‑1 (MQ=255) gATGAAGAACGTTAACTGGCTGGATGATGAAGGTCGCGATCAGGTACAGGAG‑GCGAA‑AAAATGAAAATCCg < 1:148799/71‑1 (MQ=255) aagaaCGTTAACTGGCTGGATGATGAAGGTCGCGATCAGGTACAGGAG‑GCGAA‑AAAATGAAAATCCGTTca < 1:1018540/71‑1 (MQ=255) agaaCGTTAACTGGCTGGATGATGAAGGTCGCGATCAGGTACAGGAG‑GCGAA‑AAAATGAAAATCCGTTcac < 1:1445858/71‑1 (MQ=255) gaaCATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAATTCGTTcac < 1:1311226/70‑1 (MQ=255) aCATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAATTc > 1:2275607/1‑62 (MQ=255) gACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAATTCGTTCACAAGTCgg < 1:1714093/70‑1 (MQ=255) gACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAATTCGTTCACAAGTCgg < 1:3132414/70‑1 (MQ=255) cTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAATTCGTTCACAAGTCGGCAt > 1:1855724/1‑71 (MQ=255) cTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAATTCGTTCACAAGTCGGCAt > 1:1350295/1‑71 (MQ=255) tGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAATTc > 1:1875490/1‑54 (MQ=255) tGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAATTCGtt > 1:2995808/1‑57 (MQ=255) tGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAATTCGTTCACAAGTCGGCATg > 1:595237/1‑71 (MQ=255) gTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAATTCGTTCACAAGTCGGCATgg < 1:3760649/70‑1 (MQ=255) gTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAATTCGTTCACAAGTCGGCATgg < 1:2762370/70‑1 (MQ=255) gTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAATTCGTTCACAAGTCGGCATgg < 1:1763620/70‑1 (MQ=255) gTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAATTCGTTCACAAGTCGGCATg < 1:1628791/69‑1 (MQ=255) gTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAATTCGTTCACAAGTCGGCATGGt < 1:3780801/71‑1 (MQ=255) gTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAATTCGTTCACAAGTCGGCATGGt > 1:2223794/1‑71 (MQ=255) gTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAATTCGTTCACAAGTCGGCATGGt > 1:2219106/1‑71 (MQ=255) ttAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAATTCGTTCACAAGTCGGCATGGTg > 1:3238161/1‑71 (MQ=255) tAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAATTCGTTcaca > 1:2932218/1‑57 (MQ=255) tAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAATTCGTTcaca > 1:1305391/1‑57 (MQ=255) aGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAATTCGTTCACAAGTCGGCATGGTGCt > 1:265481/1‑71 (MQ=255) gATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAATTCGTTCACAAGTCGGCa < 1:1642333/63‑1 (MQ=255) gATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAATTCGTTCACAAGTCGGCATGGTGCt < 1:3464655/70‑1 (MQ=255) gCGAAGGCAATGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAATTCGTTCACAAGTCgg < 1:1154401/57‑1 (MQ=255) gCGAAGGCAATGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAATTCGTTCACAAGTCGGCa < 1:1689228/59‑1 (MQ=255) gCGAAGGCAATGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAAt < 1:3189656/70‑1 (MQ=255) cGAAGGCAATGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAAtct < 1:1833539/71‑1 (MQ=255) cGAAGGCAATGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAAtct > 1:3818181/1‑71 (MQ=255) gAAGGCAATGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAATTCGTTCACAAGTCGGCa > 1:713048/1‑57 (MQ=255) aaGGCAATGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAAtctc > 1:1295489/1‑70 (MQ=255) aGGCAATGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAAtct > 1:1676601/1‑68 (MQ=255) aGGCAATGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGa < 1:665169/71‑1 (MQ=255) aGGCAATGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGa > 1:632673/1‑71 (MQ=255) ggCAATGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAATTCGTTCACAAGTCGTCATGGTGCTGAAtc < 1:1098645/66‑1 (MQ=255) gCAATGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAATTCGt < 1:2951322/40‑1 (MQ=255) cAATGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGAt > 1:2548940/1‑69 (MQ=255) cAATGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATaa < 1:2397116/71‑1 (MQ=255) aaTGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAATTCGtt < 1:1026738/39‑1 (MQ=255) aaTGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAATTCGTTcaca > 1:3176505/1‑43 (MQ=255) aTGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAATTCGTTCACAAGTCGGCATgg < 1:260240/53‑1 (MQ=255) tGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAATTCGTTcac > 1:1793006/1‑40 (MQ=255) tGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTg > 1:2816562/1‑71 (MQ=255) tcAGGTACAGGAG‑GCGAA‑AAAATGAAAATCCGTTCACaa < 1:1640459/38‑1 (MQ=9) tcAGGTACAGGAG‑GCGAA‑AAAATGAAAATCCGTTCACaa < 1:2193019/38‑1 (MQ=9) tcAGGTACAGGAG‑GCGAA‑AAAATGAAAATCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTAt < 1:1038941/70‑1 (MQ=255) cAGGTACAGGAG‑GCGAA‑AAAATGAAAATCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTAt > 1:322274/1‑70 (MQ=255) aGGTACAGGAGAGTG‑T‑AAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCAt < 1:3222746/69‑1 (MQ=255) aGGTACAGGAG‑GCGAA‑AAAATGAAAATCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTATCg > 1:580056/1‑71 (MQ=255) ggTACAGGAGAGCG‑T‑AAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCgg < 1:2976711/71‑1 (MQ=255) gTACAGGAG‑GCG‑AAAAAATGAAAATCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTATCGGc > 1:3140096/1‑71 (MQ=255) aCAGGAGAGCG‑T‑AAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCa < 1:3797985/64‑1 (MQ=255) aCAGGAGAGCG‑T‑AAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCa < 1:1988319/64‑1 (MQ=255) cAGGAGAGCG‑T‑AAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTg < 1:2270425/70‑1 (MQ=255) gagaGCG‑T‑AAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGcc > 1:3111585/1‑69 (MQ=255) agagCG‑T‑AAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGc < 1:549269/58‑1 (MQ=255) gagCG‑T‑AAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATaa > 1:681826/1‑53 (MQ=255) cG‑T‑AAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTg < 1:1840511/71‑1 (MQ=255) gaaaaaaTGAAAATCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTg < 1:1483715/68‑1 (MQ=255) aaaaaaTGAAAATCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTGt < 1:3036642/69‑1 (MQ=255) aaaaaTGAAAATCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTGtt < 1:2435410/70‑1 (MQ=255) aaaaaTGAAAATCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTGt < 1:3089650/69‑1 (MQ=255) aaaaTGAAAATCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTg > 1:969075/1‑59 (MQ=255) aaaaTGAAAATCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTGtt > 1:3637829/1‑70 (MQ=255) aaaaTGAAAATCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTGTTc > 1:3290229/1‑71 (MQ=255) aaaTGAAAATCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTGTTcg < 1:2871943/71‑2 (MQ=255) aaaTGAAAATCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTGTTcg > 1:2703013/1‑70 (MQ=255) aaTGAAAATCCGTTCCCAAGTCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTGTTcg < 1:3485491/70‑2 (MQ=255) aTGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGt > 1:1483404/1‑71 (MQ=255) tGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGt > 1:2539524/1‑70 (MQ=255) aaaTTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTAc < 1:1770798/70‑1 (MQ=255) aaaTCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTGTTcg > 1:2777663/1‑64 (MQ=255) aaaTCCGTTCACAAGTCGGCATGGTGCTTAACCTCGACAAATGTATCGGCTGCCATACCTGTTCGGTGAcc > 1:2993079/1‑71 (MQ=255) aaTTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCacac < 1:290989/57‑1 (MQ=255) aaTTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCt < 1:1588623/71‑1 (MQ=255) aTTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATa < 1:1946984/37‑1 (MQ=255) aTTCGTTCACAAGTCGGCATGGTGCGGAATCTCGATAAGTGCATCGGCTGc < 1:997862/51‑1 (MQ=255) | AGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCCACACCTGTTCAGTTACCT > minE/835222‑835358 |
| 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: ‑ |