| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 223,627 | T→A | 17.2% | intergenic (+47/+10) | prfH → / ← pepD | predicted peptide chain release factor/aminoacyl‑histidine dipeptidase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 223,627 | 0 | T | A | 17.2% | 141.7 / 23.1 | 58 | intergenic (+47/+10) | prfH/pepD | predicted peptide chain release factor/aminoacyl‑histidine dipeptidase |
| Reads supporting (aligned to +/- strand): ref base T (29/19); new base A (3/7); total (32/26) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 9.46e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01 | |||||||||||
ATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTA > minE/223560‑223692 | aTGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATc < 1:1202729/71‑1 (MQ=255) tGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCa < 1:1489167/71‑1 (MQ=255) gCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCaa < 1:2573073/70‑1 (MQ=255) gCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCaaa < 1:2709807/71‑1 (MQ=255) ttATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAaataat > 1:3923440/1‑71 (MQ=255) ttATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAaataat > 1:3568544/1‑71 (MQ=255) ttATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAaataat > 1:3694087/1‑71 (MQ=255) tATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAaataat < 1:3179785/70‑1 (MQ=255) tATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAaataa > 1:3387457/1‑69 (MQ=255) gAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACt < 1:4025775/70‑1 (MQ=255) aGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTc < 1:300274/70‑1 (MQ=255) aGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCg > 1:1564543/1‑71 (MQ=255) aGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCg > 1:608368/1‑71 (MQ=255) aTAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTc < 1:1087631/67‑1 (MQ=255) tAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCg > 1:3269546/1‑70 (MQ=255) tAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCg > 1:3712249/1‑70 (MQ=255) tAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCg > 1:580438/1‑70 (MQ=255) aaTGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCg > 1:2988826/1‑69 (MQ=255) aTGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGa < 1:2102423/70‑1 (MQ=255) tGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAAt > 1:184724/1‑71 (MQ=255) tGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCaaa > 1:291289/1‑53 (MQ=255) gAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAaat < 1:124119/53‑1 (MQ=255) gAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAAtt > 1:3296545/1‑71 (MQ=255) aaGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAAt > 1:1052346/1‑69 (MQ=255) aGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAAtt > 1:422456/1‑55 (MQ=255) aGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAAtatc > 1:63064/1‑71 (MQ=255) aGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAAtatc > 1:4022758/1‑71 (MQ=255) aTGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAAtttcttt < 1:3713623/71‑1 (MQ=255) tGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAAtttcnt > 1:521767/1‑69 (MQ=255) tGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAAtttcttt > 1:463418/1‑70 (MQ=255) gAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGc < 1:1526529/57‑1 (MQ=255) aaaCTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAAtttc > 1:659573/1‑65 (MQ=255) cTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTcag > 1:1664980/1‑68 (MQ=255) ggAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTcagc < 1:2454943/64‑1 (MQ=255) aaGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAAtttct > 1:566485/1‑56 (MQ=255) aGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtcagt > 1:1734534/1‑69 (MQ=255) gCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtcag < 1:3686845/66‑1 (MQ=255) gCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTcagc < 1:182240/71‑1 (MQ=255) gCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTcagc < 1:3928558/71‑1 (MQ=255) cGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtcagtcag > 1:3994632/1‑69 (MQ=255) cGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtcagtcag > 1:2923733/1‑69 (MQ=255) cGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtcagtcag > 1:1956766/1‑69 (MQ=255) gATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTcagca < 1:2957405/58‑1 (MQ=255) ttAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGtt < 1:3829692/59‑1 (MQ=255) ttAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGtt < 1:2133958/59‑1 (MQ=255) aCGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGtt < 1:2351997/55‑1 (MQ=255) cGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAg > 1:3822945/1‑71 (MQ=255) tCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGt < 1:3899407/67‑1 (MQ=255) tCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGt < 1:3894474/67‑1 (MQ=255) tCCGGCAGCAAATCAAATAATTACTTCGCCGGAAtttctttc < 1:1162140/42‑1 (MQ=255) tCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAAt > 1:2984144/1‑70 (MQ=255) ccGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAAt > 1:2108801/1‑69 (MQ=255) cGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTcagcag > 1:2544756/1‑57 (MQ=255) cGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTa > 1:2463241/1‑66 (MQ=255) ggCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGAcc > 1:559146/1‑71 (MQ=255) ggCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGAcc < 1:3423583/71‑1 (MQ=255) gcaTATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTa < 1:4143088/69‑1 (MQ=255) caTATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTc < 1:3729584/55‑1 (MQ=255) | ATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTA > minE/223560‑223692 |
| 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: ‑ |