| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,858,847 | C→T | 18.6% | K245K (AAG→AAA) | cysI ← | sulfite reductase, beta subunit, NAD(P)‑binding |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,858,847 | 0 | C | T | 18.6% | 130.0 / 21.4 | 59 | K245K (AAG→AAA) | cysI | sulfite reductase, beta subunit, NAD(P)‑binding |
| Reads supporting (aligned to +/- strand): ref base C (23/25); new base T (8/3); total (31/28) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.87e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGG > minE/1858781‑1858915 | gTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGc > 1:4558323/1‑71 (MQ=255) tttCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGtttt > 1:4877252/1‑71 (MQ=255) ttgCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGa < 1:2507445/70‑1 (MQ=255) ttgCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGa < 1:60340/70‑1 (MQ=255) ttgCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGa < 1:1665613/70‑1 (MQ=255) tgtTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCAc > 1:382797/1‑71 (MQ=255) ttCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACg < 1:174326/70‑1 (MQ=255) ttCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGa < 1:65929/71‑1 (MQ=255) ttCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACAAGCTTGCCGtttt > 1:4569764/1‑56 (MQ=255) gATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAg < 1:3141286/70‑1 (MQ=255) gATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGt < 1:963666/71‑1 (MQ=255) gATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGt < 1:4121802/71‑1 (MQ=255) gATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGt < 1:2186528/71‑1 (MQ=255) gATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGt < 1:2062524/71‑1 (MQ=255) gATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGt < 1:1304954/71‑1 (MQ=255) aTGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCg < 1:2324579/57‑1 (MQ=255) aTGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCg < 1:5001757/57‑1 (MQ=255) aTGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCg < 1:507387/57‑1 (MQ=255) aTGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCg < 1:1492520/57‑1 (MQ=255) tGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTc < 1:1632057/71‑1 (MQ=255) cccaccgcccaccAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGt < 1:1493628/71‑1 (MQ=255) ccaccgcccaccAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGtt < 1:3715132/71‑1 (MQ=255) gcccaccAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCG‑‑GATCGCCACGAAGTTCATGTCGTTGGCGTg > 1:4621927/1‑70 (MQ=255) caccaACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGtt < 1:4064512/54‑1 (MQ=255) caccaACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGtt < 1:2739174/54‑1 (MQ=255) caccaACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGtt < 1:4186436/54‑1 (MQ=255) caccaACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTg > 1:2119353/1‑69 (MQ=255) ccaACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGc > 1:4226924/1‑37 (MQ=255) ccaACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCa > 1:5020576/1‑54 (MQ=255) ccaACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCa > 1:3137108/1‑54 (MQ=255) ccaACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCa > 1:2957937/1‑54 (MQ=255) ccaACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCa > 1:2921073/1‑54 (MQ=255) ccaACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCa > 1:4070660/1‑54 (MQ=255) ccaACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCa > 1:4661521/1‑54 (MQ=255) ccaACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCa > 1:1522611/1‑54 (MQ=255) caACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGtt < 1:565259/51‑1 (MQ=255) ggTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGata > 1:1275562/1‑71 (MQ=255) ggTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGata > 1:1075576/1‑71 (MQ=255) ggTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGata > 1:135824/1‑71 (MQ=255) ggTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGata > 1:3186400/1‑71 (MQ=255) ggTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGata > 1:605313/1‑71 (MQ=255) ggTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGata > 1:5055894/1‑71 (MQ=255) ggTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGa > 1:3245145/1‑69 (MQ=255) ggTTAAAGCCCACCAGCTTGCCGTTTTCGGAGATCCCCACGAAGTTCATGTc > 1:631277/1‑52 (MQ=39) gTTAAAGCCCACCAGCTTGCCTTTTTCGGCGATCTCCACGAAGTTCATGTCGTTGGCGTGCAGa < 1:511506/64‑1 (MQ=255) gTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGa < 1:3179119/64‑1 (MQ=255) gTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGa < 1:3372639/64‑1 (MQ=255) gTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGa < 1:2891826/64‑1 (MQ=255) tAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATc > 1:3394163/1‑64 (MQ=255) aaaGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGata > 1:1441128/1‑67 (MQ=255) aaaGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGt > 1:2663674/1‑71 (MQ=255) aaaGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGt > 1:2357715/1‑71 (MQ=255) aaGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGt < 1:232850/70‑1 (MQ=255) aGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTc < 1:1152356/71‑1 (MQ=255) accaGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTgg > 1:2798108/1‑71 (MQ=255) caGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCg > 1:1731380/1‑71 (MQ=255) caGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCg > 1:1355055/1‑71 (MQ=255) caGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCg > 1:95249/1‑71 (MQ=255) aGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTATGTGGCgg > 1:3220066/1‑71 (MQ=255) cTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGa < 1:3773389/53‑1 (MQ=255) cTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGa < 1:3059993/53‑1 (MQ=255) | GTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGG > minE/1858781‑1858915 |
| 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: ‑ |