Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 223,739 | A→G | 40.0% | T452T (ACT→ACC) | pepD ← | aminoacyl‑histidine dipeptidase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 223,739 | 0 | A | G | 40.0% | 28.8 / 51.7 | 75 | T452T (ACT→ACC) | pepD | aminoacyl‑histidine dipeptidase |
Reads supporting (aligned to +/- strand): ref base A (38/7); new base G (14/16); total (52/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.95e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.84e-01 |
ATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACA > minE/223685‑223795 | aTGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGTCCCGATAGa > 1:362569/1‑69 (MQ=255) aTGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGc > 1:3232653/1‑61 (MQ=255) aTGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGa > 1:1850641/1‑69 (MQ=255) aTGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGa > 1:2605085/1‑69 (MQ=255) aTGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGa > 1:2146257/1‑69 (MQ=255) aTGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGa > 1:2500862/1‑69 (MQ=255) aTGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGa > 1:2507784/1‑69 (MQ=255) aTGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGa > 1:2628711/1‑69 (MQ=255) aTGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGa > 1:2777490/1‑69 (MQ=255) aTGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGa > 1:2899595/1‑69 (MQ=255) aTGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGa > 1:3164543/1‑69 (MQ=255) aTGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGa > 1:3317888/1‑69 (MQ=255) aTGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGa > 1:3433735/1‑69 (MQ=255) aTGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGa > 1:53245/1‑69 (MQ=255) aTGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGa > 1:561076/1‑69 (MQ=255) aTGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGa > 1:637656/1‑69 (MQ=255) aTGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGa > 1:710974/1‑69 (MQ=255) aTGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGa > 1:926195/1‑69 (MQ=255) aTGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGa > 1:1085409/1‑69 (MQ=255) aTGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGa > 1:1146301/1‑69 (MQ=255) aTGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGa > 1:1454564/1‑69 (MQ=255) aCCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAAc < 1:771991/69‑1 (MQ=255) aCCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAAc < 1:530628/69‑1 (MQ=255) tACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCAt > 1:2822349/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) tACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCAt > 1:2143720/1‑69 (MQ=255) ttCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGa > 1:2671643/1‑53 (MQ=255) ttCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAAc > 1:1338379/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:2909261/1‑60 (MQ=255) ttCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAAc > 1:2870050/1‑60 (MQ=255) ttCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAAc > 1:2805459/1‑60 (MQ=255) ttCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAAc > 1:2626795/1‑60 (MQ=255) ttCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAAc > 1:772333/1‑60 (MQ=255) tCATCCGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGtt < 1:2115368/67‑1 (MQ=255) tCATCCGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGtt < 1:1754786/67‑1 (MQ=255) tCATCCGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGtt < 1:340507/67‑1 (MQ=255) tCATCCGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGtt < 1:1411383/67‑1 (MQ=255) tCATCCGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGtt < 1:448726/67‑1 (MQ=255) tCATCCGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGtt < 1:839464/67‑1 (MQ=255) tCATCCGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGtt < 1:2280958/67‑1 (MQ=255) tCATCCGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGtt < 1:2326173/67‑1 (MQ=255) tCATCCGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGtt < 1:2419012/67‑1 (MQ=255) tCATCCGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGtt < 1:1664091/67‑1 (MQ=255) tCATCCGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGtt < 1:1736025/67‑1 (MQ=255) tCATCCGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGtt < 1:1202462/67‑1 (MQ=255) tCATCCGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGtt < 1:1841450/67‑1 (MQ=255) tCATCCGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGtt < 1:580694/67‑1 (MQ=255) tCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAAc < 1:1876366/41‑1 (MQ=255) ccGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGTTTTTTg > 1:1693694/1‑68 (MQ=255) ccGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGTTTTTTg > 1:187153/1‑68 (MQ=255) ccGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGTTTTTTg > 1:3037908/1‑68 (MQ=255) ccGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGTTTTTTg > 1:4209/1‑68 (MQ=255) ccGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGTTTTTTg > 1:58482/1‑68 (MQ=255) ccGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGTTTTTTg > 1:1892480/1‑68 (MQ=255) ccGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGTTTTTTg > 1:1591842/1‑68 (MQ=255) ccGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGTTTTTTg > 1:2692502/1‑68 (MQ=255) ccGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGTTTTTTg > 1:263461/1‑68 (MQ=255) ccGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGTTTTTTg > 1:2587173/1‑68 (MQ=255) ccGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGTTTTTTg > 1:227151/1‑68 (MQ=255) ccGGAGAGTGTGGACCGGTGATGGTTGGCCCGATAGACACCATGTCCATTTCCGGATACGGTTTTTTg > 1:1674299/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) | ATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGGTGATAGTTGGCCCGATAGAAACCATGTCCATTTCCGGATACGGTTTTTTGAACAGACCACA > minE/223685‑223795 |
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: ‑ |