| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 835,275 | Δ1 bp | 18.7% | coding (3733/3744 nt) | narG → | nitrate reductase 1, alpha subunit |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 835,275 | 0 | A | . | 18.7% | 191.4 / 51.2 | 80 | coding (3733/3744 nt) | narG | nitrate reductase 1, alpha subunit |
| Reads supporting (aligned to +/- strand): ref base A (30/35); new base . (6/9); total (36/44) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.77e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01 | |||||||||||
TTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCC > minE/835207‑835340 | ttgttgTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGgaga > 1:2702740/1‑69 (MQ=255)ttgttgTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGgaga > 1:780615/1‑69 (MQ=255) tgttgtAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCg > 1:541158/1‑71 (MQ=255) gttgtAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCg > 1:3855325/1‑70 (MQ=255) tgtAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑T‑a > 1:3829155/1‑70 (MQ=255) gtAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑T‑aa < 1:2938315/70‑1 (MQ=255) gtAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑T‑aa > 1:328648/1‑70 (MQ=255) gtAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑T‑aa > 1:3759474/1‑70 (MQ=255) tAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑T‑aaa < 1:3301694/70‑1 (MQ=255) aGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑T‑aaaa < 1:2017284/70‑1 (MQ=255) aGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑T‑AAAAt > 1:3537989/1‑71 (MQ=255) gCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑T‑AAAATGa > 1:2445047/1‑70 (MQ=255) cgcaaGATGAAGAACGTTAACTGGCTGGATGATGAAGGTCGCGATCAGGTACAGGAG‑GCGAA‑AAAATGaaa < 1:1657235/68‑1 (MQ=255) cgcaaGATGAAGAACGTTAACTGGCTGGATGATGAAGGTCGCGATCAGGTACAGGAG‑GCGAA‑AAAATGaaa < 1:2209642/68‑1 (MQ=255) tAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑T‑AAAATGa < 1:3772670/67‑1 (MQ=255) tAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAAt < 1:751456/71‑1 (MQ=255) tAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGAG‑T‑AAAATGAAAAt < 1:3862461/71‑1 (MQ=255) caaGATGAAGAACGTTAACTGGCTGGATGATGAAGGTCGCGATCAGGTACAGGAG‑GCGAA‑AAAATGAAAAt > 1:2576504/2‑71 (MQ=255) caaGATGAAGAACGTTAACTGGCTGGATGATGAAGGTCGCGATCAGGTACAGGAG‑GCGAA‑AAAATGAAAAt > 1:3320715/2‑71 (MQ=255) aGATGAAGAACGTTAACTGGCTGGATGATGAAGGTCGCGATCAGGTACAGGAG‑GCGAA‑AAAATGAAAAtcc > 1:2944020/1‑69 (MQ=255) aGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAATTc < 1:1001081/71‑1 (MQ=255) aGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAATTc > 1:3828491/1‑71 (MQ=255) aGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAATTc > 1:614527/1‑71 (MQ=255) gATGAAGAACGTTAACTGGCTGGATGATGAAGGTCGCGATCAGGTACAGGAG‑GCGAA‑AAAATGAAAATCCg < 1:148799/71‑1 (MQ=255) gATGAAGAACGTTAACTGGCTGGATGATGAAGGTCGCGATCAGGTACAGGAG‑GCGAA‑AAAATGAAAATCCg < 1:1923858/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) aTTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑T‑AAAATGaaaa < 1:3940656/57‑1 (MQ=255) gACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAATTCGTTCACAAGTCgg < 1:3132414/70‑1 (MQ=255) gACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAATTCGTTCACAAGTCgg < 1:1714093/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:1763620/70‑1 (MQ=255) gTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAATTCGTTCACAAGTCGGCATgg < 1:2762370/70‑1 (MQ=255) gTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAATTCGTTCACAAGTCGGCATgg < 1:3760649/70‑1 (MQ=255) gTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAATTCGTTCACAAGTCGGCATg < 1:1628791/69‑1 (MQ=255) gTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAATTCGTTCACAAGTCGGCATGGt > 1:2223794/1‑71 (MQ=255) gTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAATTCGTTCACAAGTCGGCATGGt > 1:2219106/1‑71 (MQ=255) gTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAATTCGTTCACAAGTCGGCATGGt < 1:3780801/71‑1 (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‑AAAATGaaaa < 1:484874/46‑1 (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:3818181/1‑71 (MQ=255) cGAAGGCAATGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAAtct < 1:1833539/71‑1 (MQ=255) gAAGGCAATGACCAGGTACAGGAGAGCG‑T‑AAAATGa < 1:123019/36‑1 (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:632673/1‑71 (MQ=255) aGGCAATGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGa < 1:665169/71‑1 (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) | TTGTTGTAGTGCGTAAGATGAAGAACATTGACTGGTTAGATGGCGAAGGCAATGACCAGGTACAGGAGAGCG‑T‑AAAATGAAAATTCGTTCACAAGTCGGCATGGTGCTGAATCTCGATAAGTGCATCGGCTGCC > minE/835207‑835340 |
| 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: ‑ |