| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 232,901 | C→T | 24.6% | intergenic (‑140/+30) | hemB ← / ← ykiB | porphobilinogen synthase/hypothetical protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 232,901 | 0 | C | T | 24.6% | 114.5 / 35.9 | 61 | intergenic (‑140/+30) | hemB/ykiB | porphobilinogen synthase/hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base C (34/12); new base T (9/6); total (43/18) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.40e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.92e-01 | |||||||||||
GTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTGGAGTTTGGTTG > minE/232831‑232971 | gTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAAc < 1:1229786/71‑1 (MQ=255)gTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAAc < 1:1522977/71‑1 (MQ=255) tCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGATACAGCAATAATTTta < 1:1848662/71‑1 (MQ=255) gcgATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTtata < 1:2570341/71‑1 (MQ=255) aTGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTc < 1:130233/71‑1 (MQ=255) ttGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAAc < 1:1664713/59‑1 (MQ=255) gaagaGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAACtt > 1:443288/1‑54 (MQ=255) gaagaGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAACtt > 1:458939/1‑54 (MQ=255) gaagaGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAACTTta > 1:1171627/1‑56 (MQ=255) aagaGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGAAAatat > 1:3436557/1‑71 (MQ=255) gagaATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATAttt < 1:4143084/71‑1 (MQ=255) tGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAat > 1:2355444/1‑62 (MQ=255) tGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAAt > 1:1302713/1‑71 (MQ=255) agagagGCCAGGACATCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCtt > 1:666324/1‑70 (MQ=255) agagagGCCAGGACATCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCtt > 1:510138/1‑70 (MQ=255) gagagGCCAGGACATCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGaaaa < 1:66533/56‑1 (MQ=255) agagGCCAGGACATCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGAAAata < 1:2088067/57‑1 (MQ=255) gagGCCAGGACATCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCtttatt > 1:3738782/1‑71 (MQ=255) agGCCAGGACATCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCtttattt < 1:497253/71‑1 (MQ=255) ggACATCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCt > 1:4141553/1‑71 (MQ=255) gACATCTTTTTAATAGAAACAGCAATAATTTTatat > 1:101235/1‑36 (MQ=255) gACATCTTTTTAATAGAAACAGCAATAATTTTatat > 1:342421/1‑36 (MQ=255) aCATCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGt > 1:3218142/1‑71 (MQ=255) aCATCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGt > 1:1973331/1‑71 (MQ=255) cATCTTTTTAATAGAAACAGCAATAACTTTATATTCAATGAAAATATTTTTAATCTTTATTTATAg > 1:2626175/1‑66 (MQ=255) aTCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGAAAata < 1:3857963/45‑1 (MQ=255) aTCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTa < 1:255407/51‑1 (MQ=255) aTCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCt < 1:322255/67‑1 (MQ=255) aTCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTg > 1:2403844/1‑71 (MQ=255) tCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCtt < 1:1311247/55‑1 (MQ=255) tCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGAAAata < 1:22712/44‑1 (MQ=255) cTTTTTAATAGAAACAGCAATAACTTTATATTCACTGAAAATAtt < 1:381912/45‑1 (MQ=255) tttttAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGt > 1:658627/1‑70 (MQ=255) tttttAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGt > 1:2303917/1‑70 (MQ=255) tttttAATAGAAACAGCAATAACTTTATATTCACTGaaaa > 1:1868136/1‑40 (MQ=255) tttttAATAGAAACAGCAATAACTTTATATTCACTGAAAATAtttt > 1:3567385/1‑46 (MQ=255) tttttAATAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGt > 1:240102/1‑70 (MQ=255) ttttAATAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGtta > 1:2233202/1‑71 (MQ=255) ttttAATAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGtta > 1:2604129/1‑71 (MQ=255) tttAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGttat > 1:782560/1‑71 (MQ=255) tttAATAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGttat > 1:3412808/1‑71 (MQ=255) tttAATAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGttat > 1:2041298/1‑71 (MQ=255) tttAATAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGttat > 1:666106/1‑71 (MQ=255) ttAATAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCttta > 1:1858322/1‑52 (MQ=255) aaTAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCtttattt > 1:828676/1‑53 (MQ=255) tAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCtt < 1:1369360/46‑1 (MQ=255) tAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAttt > 1:2226230/1‑71 (MQ=255) aGAAACAGCAATAACTTTATATTCACTGAAAATAttt > 1:459529/1‑37 (MQ=255) gAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCttta > 1:979518/1‑46 (MQ=255) gAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAttttt < 1:3558059/71‑1 (MQ=255) aaCAGCAATAACTTTATATTCACTGAAAATATTTTTAAt > 1:2744409/1‑39 (MQ=255) aaCAGCAATAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTg > 1:4132755/1‑71 (MQ=255) aaCAGCAATAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTg > 1:1333617/1‑71 (MQ=255) cAATAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTGGAGtt > 1:4015896/1‑71 (MQ=255) cAATAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTGGAGtt > 1:501244/1‑71 (MQ=255) aTAATTTTATATTCACTGAAAATATTTTTAATCtttattta > 1:3137668/1‑41 (MQ=255) aTAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTGGAg < 1:3099836/67‑1 (MQ=255) tAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGttat < 1:3619040/54‑1 (MQ=255) tAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTGGAGTTTg < 1:1600777/70‑1 (MQ=255) aaCTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGc > 1:685292/1‑43 (MQ=255) aaCTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTGGAGTTTGGt > 1:1232209/1‑71 (MQ=255) aCTTTATATTCACTGAAAATATTTTTAATCtttattta > 1:1696969/1‑38 (MQ=255) aCTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTGGAGTTTGGtt > 1:2853009/1‑71 (MQ=255) aCTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTGGAGTTTGGtt > 1:2608303/1‑71 (MQ=255) cTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTGGAGTTTGGTTg > 1:130107/1‑71 (MQ=255) | GTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTGGAGTTTGGTTG > minE/232831‑232971 |
| 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: ‑ |