| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 223,754 | A→C | 50.9% | V447V (GTT→GTG) | pepD ← | aminoacyl‑histidine dipeptidase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 223,754 | 0 | A | C | 50.9% | 6.0 / 63.6 | 59 | V447V (GTT→GTG) | pepD | aminoacyl‑histidine dipeptidase |
| Reads supporting (aligned to +/- strand): ref base A (16/13); new base C (14/16); total (30/29) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.06e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.92e-01 | |||||||||||
ACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTG > minE/223688‑223819 | aCCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAAc < 1:771991/69‑1 (MQ=255)aCCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAAc < 1:530628/69‑1 (MQ=255) tACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCAt > 1:2822349/1‑69 (MQ=255) tACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCAt > 1:2143720/1‑69 (MQ=255) tACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCAt > 1:1357405/1‑69 (MQ=255) tACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCAt > 1:1377902/1‑69 (MQ=255) tACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCAt > 1:1431329/1‑69 (MQ=255) tACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCAt > 1:1534298/1‑69 (MQ=255) ttCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAAc > 1:2626795/1‑60 (MQ=255) ttCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAAc > 1:2805459/1‑60 (MQ=255) ttCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAAc > 1:2870050/1‑60 (MQ=255) ttCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAAc > 1:2909261/1‑60 (MQ=255) ttCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAAc > 1:402019/1‑60 (MQ=255) ttCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAAc > 1:456705/1‑60 (MQ=255) ttCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAAc > 1:505456/1‑60 (MQ=255) ttCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAAc > 1:1338379/1‑60 (MQ=255) ttCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAAc > 1:772333/1‑60 (MQ=255) tCATCCGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGtt < 1:2326173/67‑1 (MQ=255) tCATCCGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGtt < 1:2419012/67‑1 (MQ=255) tCATCCGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGtt < 1:2280958/67‑1 (MQ=255) tCATCCGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGtt < 1:2115368/67‑1 (MQ=255) tCATCCGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGtt < 1:340507/67‑1 (MQ=255) tCATCCGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGtt < 1:448726/67‑1 (MQ=255) tCATCCGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGtt < 1:580694/67‑1 (MQ=255) tCATCCGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGtt < 1:839464/67‑1 (MQ=255) tCATCCGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGtt < 1:1664091/67‑1 (MQ=255) tCATCCGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGtt < 1:1411383/67‑1 (MQ=255) tCATCCGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGtt < 1:1202462/67‑1 (MQ=255) tCATCCGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGtt < 1:1736025/67‑1 (MQ=255) tCATCCGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGtt < 1:1754786/67‑1 (MQ=255) tCATCCGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGtt < 1:1841450/67‑1 (MQ=255) tCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAAc < 1:1876366/41‑1 (MQ=255) ccGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGTTTTTTg > 1:1591842/1‑68 (MQ=255) ccGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGTTTTTTg > 1:4209/1‑68 (MQ=255) ccGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGTTTTTTg > 1:1674299/1‑68 (MQ=255) ccGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGTTTTTTg > 1:58482/1‑68 (MQ=255) ccGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGTTTTTTg > 1:3037908/1‑68 (MQ=255) ccGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGTTTTTTg > 1:1693694/1‑68 (MQ=255) ccGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGTTTTTTg > 1:2692502/1‑68 (MQ=255) ccGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGTTTTTTg > 1:263461/1‑68 (MQ=255) ccGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGTTTTTTg > 1:187153/1‑68 (MQ=255) ccGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGTTTTTTg > 1:2587173/1‑68 (MQ=255) ccGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGTTTTTTg > 1:1892480/1‑68 (MQ=255) ccGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGTTTTTTg > 1:227151/1‑68 (MQ=255) ccGGAGACTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGTTTTTTg > 1:1162406/1‑68 (MQ=255) acGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGTTTTTTg > 1:2784648/2‑68 (MQ=255) cGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGaa < 1:2072798/69‑1 (MQ=255) cGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGaa < 1:1658458/69‑1 (MQ=255) gagTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGtt < 1:520597/58‑1 (MQ=255) gtgtGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGtt < 1:1011388/56‑1 (MQ=255) tgtgGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACca < 1:1905813/69‑1 (MQ=255) tgtgGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACca < 1:2951569/69‑1 (MQ=255) tgGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACcaca > 1:1655019/1‑69 (MQ=255) ggCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACcac < 1:3178154/52‑1 (MQ=255) ggCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACcac < 1:3098397/52‑1 (MQ=255) ggCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACcac < 1:1057455/52‑1 (MQ=255) aGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTg < 1:325792/69‑1 (MQ=255) aGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTg < 1:126315/69‑1 (MQ=255) aGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTg < 1:2260190/69‑1 (MQ=255) | ACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACATTCCAGGCCCGCGTGGATAATCTG > minE/223688‑223819 |
| 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: ‑ |