| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 232,635 | T→C | 27.9% | N43D (AAC→GAC) | hemB ← | porphobilinogen synthase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 232,635 | 0 | T | C | 27.9% | 91.3 / 45.0 | 61 | N43D (AAC→GAC) | hemB | porphobilinogen synthase |
| Reads supporting (aligned to +/- strand): ref base T (23/21); new base C (9/8); total (32/29) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01 | |||||||||||
CAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACA > minE/232566‑232703 | cAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTATGTAGTTg > 1:3514560/1‑71 (MQ=255) tttCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCa < 1:4029115/71‑1 (MQ=255) tCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCaa > 1:2745029/1‑70 (MQ=255) gTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATttct < 1:3596263/71‑1 (MQ=255) gATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAAttt > 1:397188/1‑63 (MQ=255) gATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATttcttc < 1:3730568/67‑1 (MQ=255) gATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATttcttc < 1:3725539/67‑1 (MQ=255) gATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTCGTCaa > 1:3821746/1‑60 (MQ=255) tttCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTcaaca < 1:1218100/70‑1 (MQ=255) tctcTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTcaaca > 1:3550298/1‑68 (MQ=255) gAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCgg > 1:2541897/1‑70 (MQ=255) gAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGATCGGc < 1:2994965/71‑1 (MQ=255) aaTGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCa < 1:189627/71‑1 (MQ=255) aaTGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCa < 1:817807/71‑1 (MQ=255) aTGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGct < 1:3364228/70‑2 (MQ=255) aTGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGATCGGCaa < 1:895103/71‑1 (MQ=255) gcgcATCACGCCTGGCATGGCTTGAACGTCTTTGTAGTTGTCAATTTCTTCTTcaacaa < 1:2587487/59‑1 (MQ=255) cgcATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTCGTCAATttct > 1:2621499/1‑48 (MQ=255) cgcATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTcaac > 1:1207659/1‑56 (MQ=255) cgcATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTcaac > 1:3696434/1‑56 (MQ=255) cgcATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGATCGGCAacac < 1:3601780/71‑1 (MQ=255) cgcATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGATCGGCAacac < 1:1300199/71‑1 (MQ=255) aTCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAg > 1:3709354/1‑71 (MQ=255) tCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATCTCTTCTTCAACAAAg < 1:1119852/56‑1 (MQ=255) tCACGCCTGGCATGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGATCGGCAACAcc < 1:2017924/68‑1 (MQ=255) cACGCCTGGCATGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGATCGGCAACAcc < 1:827548/67‑1 (MQ=255) aCGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTcaacaa < 1:3225315/52‑1 (MQ=255) aCGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCgg > 1:3731995/1‑59 (MQ=255) aCGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCgg > 1:3726976/1‑59 (MQ=255) aCGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCaa > 1:460421/1‑62 (MQ=255) aCGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGt > 1:2805862/1‑70 (MQ=255) aCGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTc > 1:1247135/1‑71 (MQ=255) aCGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTc > 1:2625155/1‑71 (MQ=255) aCGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTc > 1:2031183/1‑71 (MQ=255) aCGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATAGGCAACACCAGGTc > 1:2068467/1‑71 (MQ=255) aCGCCTGGCATGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTAAACAAAg > 1:2630362/1‑54 (MQ=39) ccTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGtt < 1:3611799/71‑1 (MQ=255) ccTGGCATGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTCAACaaa > 1:3631373/1‑50 (MQ=255) tGGCATGGCTTCAACGGCTTTGTAGTTGTCAATttct < 1:272338/37‑1 (MQ=255) ggCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTcaaca < 1:1117105/45‑1 (MQ=255) aTGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGc > 1:2572721/1‑70 (MQ=255) aTGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCt < 1:31764/71‑1 (MQ=255) tGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGa < 1:744857/45‑1 (MQ=255) gCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTcaa > 1:1394120/1‑37 (MQ=255) gCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTaa > 1:860752/1‑64 (MQ=255) gCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAg > 1:3352920/1‑71 (MQ=255) cTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTaa > 1:1480448/1‑63 (MQ=255) ttCAACGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAgtg > 1:889999/1‑71 (MQ=255) tCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAaca > 1:846084/1‑49 (MQ=255) tCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAaca > 1:3657059/1‑49 (MQ=255) cAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCaa > 1:252868/1‑46 (MQ=255) cAACGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTaa > 1:1188874/1‑60 (MQ=255) aaCGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGtt < 1:3353381/57‑1 (MQ=255) aCGTCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTa < 1:854963/57‑1 (MQ=255) aCGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACGCCAGGTCGTTa < 1:3444839/57‑1 (MQ=255) ggCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTGAAGGCTAAGTGTTGtct > 1:2967480/1‑71 (MQ=255) gCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGt < 1:930752/52‑1 (MQ=255) gCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAgtgt > 1:1352837/1‑65 (MQ=255) tttGTAGTCGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGtt < 1:203741/51‑1 (MQ=255) gTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCa < 1:359113/70‑1 (MQ=255) gTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACa < 1:988924/71‑1 (MQ=255) ttGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTa < 1:4106304/51‑1 (MQ=255) ttGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTa < 1:3422767/51‑1 (MQ=255) ttGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACa < 1:1454654/70‑1 (MQ=255) tGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGtgttgt > 1:1434409/1‑58 (MQ=255) | CAATTTCGCGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTTGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACA > minE/232566‑232703 |
| 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: ‑ |