| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 2,346,508 | G→A | 18.9% | P21S (CCA→TCA) | hemN ← | coproporphyrinogen III oxidase, SAM and NAD(P)H dependent, oxygen‑independent |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 2,346,508 | 0 | G | A | 18.9% | 78.2 / 12.1 | 37 | P21S (CCA→TCA) | hemN | coproporphyrinogen III oxidase, SAM and NAD(P)H dependent, oxygen‑independent |
| Reads supporting (aligned to +/- strand): ref base G (13/17); new base A (3/4); total (16/21) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.30e-01 | |||||||||||
GTAAAAACGCCTGTTCGCCGAAGTCTTCTGAAAACTCCAGCGCGGTCGGGTACGAGGTGTATCGTGGCCCGGAATAGTTATATTTCTGGATCAGGGCCAGATCCCAGTCGATTTGCTGTACAGACATGCTCACTCCT > minE/2346442‑2346578 | gTAAAAACGCCTGTTCGCCGAAGTCTTCTGAAAACTCCAGCGCGGTCGGGTACGAGGTGTATCGTGAcc < 1:2275577/69‑1 (MQ=255) tAAAAACGCCTGTTCGCCGAAGTCTTCTGAAAACTCCAGCGCGGTCGGGTACGAGGTGTATCGTGGCCCg > 1:330162/1‑70 (MQ=255) tAAAAACGCCTGTTCGCCGAAGTCTTCTGAAAACTCCAGCGCGGTCGGGTACGAGGTGTATCGTGACCCg > 1:2644885/1‑70 (MQ=255) aaaaaCGCCTGTTCGCCGAAGTCTTCTGAAAACTCCAGCGCGGTCGGGTACGAGGTGTATCGTGGCCCgg > 1:3812989/1‑70 (MQ=255) aaaaCGCCTGTTCGCCGAAGTCTTCTGAAAACTCCAGCGCGGTCGGGTACGAGGTGTATCGTGGCCCGGaa < 1:3940108/71‑1 (MQ=255) aaaaCGCCTGTTCGCCGAAGTCTTCTGAAAACTCCAGCGCGGTCGGGTACGAGGTGTATCGTGGCCCGGaa < 1:296845/71‑1 (MQ=255) aaaCGCCTGTTCGCCGAAGTCTTCTGAAAACTCCAGCGCGGTCGGGTACGAGGTGTATCGTGGCCCGGAAt < 1:2570890/71‑1 (MQ=255) aaaCGCCTGTTCGCCGAAGTCTTCTGAAAACTCCAGCGCGGTCGGGTACGAGGTGTATCGTGACCCGGaa > 1:100937/1‑70 (MQ=255) aaCGCCTGTTCGCCGAAGTCTTCTGAAAACTCCAGCGCGGTCGGGTACGAGGTGTATCGTGGCCCGGAATa < 1:702721/71‑1 (MQ=255) aCGCCTGTTCGCCGAAGTCTTCTGAAAACTCCAGCGCGGTCGGGTACGAGGTGTATCGTGGCCCGGAATAg > 1:2299423/1‑71 (MQ=255) cTGTTCGCCGAAGTCTTCTGAAAACTCCAGCGCGGTCGGGTACGAGGTGTATCGTGGcc < 1:1621148/59‑1 (MQ=255) cTGTTCGCCGAAGTCTTCTGAAAACTCCAGCGCGGTCGGGTACGAGGTGTATCGTGGcc < 1:3540338/59‑1 (MQ=255) attcttctGAAAACTCCAGCGCGGTCGGGTACGAGGTGTATCGTGGCCCGGAATAGTTATATTTCTGGATc < 1:9107/69‑1 (MQ=255) ttctGAAAACTCCAGCGCGGTCGGGTACGAGGTGTATCGTGGCCCGGAATAGTTATATTTCTGGATCAggg > 1:3435/1‑71 (MQ=255) ctGAAAACTCCAGCGCGGTCGGGTACGAGGTGTATCGTGGCCCGGAATAGTTATATTTCTGGATCAGGGc > 1:2083714/1‑70 (MQ=255) aaaaCTCCAGCGCGGTCGGGTACGAGGTGTATCGTGGCCCGGAATAGTTATAttt > 1:2296408/1‑55 (MQ=255) aaaaCTCCAGCGCGGTCGGGTACGAGGTGTATCGTGGCCCGGAATAGTTATATTTCTGGATCAGGGCCAg < 1:3294088/70‑1 (MQ=255) aaaCTCCAGCGCGGTCGGGTACGAGGTGTATCGTGGCCCGTAATAGTTATATTTCTGGATCAGGGCCAGAt < 1:1879943/71‑1 (MQ=255) cgGTCGGGTACGAGGTGTATCGTGGCCCGGAATAGTTATATTTCTGGATCAGGGCCAg < 1:1827705/58‑1 (MQ=255) cgGTCGGGTACGAGGTGTATCGTGACCCGGAATAGTTATATTTCTGGATCAGGGCCAGATCCCAGTCGAtt < 1:2048083/71‑1 (MQ=255) ggTCGGGTACGAGGTGTATCGTGGCCCGGAATAGTTATATTTCTGGATCAGGGCCa < 1:3479229/56‑1 (MQ=255) ggTCGGGTACGAGGTGTATCGTGACCCGGAATAGTTATATTTCTGGAt < 1:3419677/48‑1 (MQ=255) gggTACGAGGTGTATCGTGGCCCGCAATAGTTATATTTCTGGATCAGGGCCAGATCCCAGTCGATTTGCTg > 1:3718407/1‑71 (MQ=255) gTACGAGGTGTATCGTGGCCCGGAATAGTTATATTTCTGGATCAGGGCCAGAt < 1:232388/53‑1 (MQ=255) gTACGAGGTGTATCGTGGCCCGGAATAGTTATATTTCTGGATCAGGGCCAGATCCCAGTCGATTTGCTGTa > 1:3576520/1‑71 (MQ=255) aCGAGGTGTATCGTGGCCCGGAATAGTTATATTTCTGGATCAGGGCCAGATCCCAGTCGATTTGCTGTa > 1:2018666/1‑69 (MQ=255) aCGAGGTGTATCGTGACCCGGAATAGTTATATTTCTGGATCAGGGCCAGATCCCAGTCGATTTGCTGTa > 1:2293458/1‑69 (MQ=255) cGAGGTGTATCGTGGCCCGGAATAGTTATATTTCTGGATCAGGGCCAGATCCCAGTCGATTTGCTGTACAg < 1:2143671/71‑1 (MQ=255) gAGGTGTATCGTGGCCCGGAATAGTTATATTTCTGGAt > 1:3554573/1‑38 (MQ=255) gAGGTGTATCGTGGCCCGGAATAGTTATATTTCTGGATc < 1:1656563/39‑1 (MQ=255) ggTGTATCGTGGCCCGGAATAGTTATATTTCTGGATCAGGGCCAGATCCCAGTCGATTTGCTGTACAGACa < 1:1674858/71‑1 (MQ=255) ggTGTATCGTGGCCCGGAATAGTTATATTTCTGGATCAGGGCCAGATCCCAGTCGATTTGCTGAACAGACa > 1:1294113/1‑71 (MQ=255) tCGTGGCCCGGAATAGTTATATTTCTGGATCAGGGCCAGATCCCAGTCGATTTGCTGTACAGa < 1:3303952/63‑1 (MQ=255) tCGTGGCCCGGAATAGTTATATTTCTGGATCAGGGCCAGATCCCAGTCGATTTGCTGTACAGACATGCTCa > 1:168681/1‑71 (MQ=255) tCGTGACCCGGAATAGTTATATTTCTGGATCAGGGCCAGATCCCAGTAGATTTGCTGTACAGACATGCTCa < 1:1941336/71‑1 (MQ=255) gTGGCCCGGAATAGTTATATTTCTGGATCAGGGCCAGATCCCAGTCGAttt > 1:1314094/1‑51 (MQ=255) tGGCCCGGAATAGTTATATTTCTGGATCAGGGCCAGATCCCAGTCGATTTGCTGTACa < 1:1295431/58‑1 (MQ=255) gCCCGGAATAGTTATATTTCTGGATCAGGGCCAGATCCCAGTCGATTTGCTGTACAGACATGCTCACTCCt < 1:3237790/71‑1 (MQ=255) | GTAAAAACGCCTGTTCGCCGAAGTCTTCTGAAAACTCCAGCGCGGTCGGGTACGAGGTGTATCGTGGCCCGGAATAGTTATATTTCTGGATCAGGGCCAGATCCCAGTCGATTTGCTGTACAGACATGCTCACTCCT > minE/2346442‑2346578 |
| 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: ‑ |