Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 223,627 | T→A | 51.9% | 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 | 51.9% | ‑2.6 / 37.3 | 27 | intergenic (+47/+10) | prfH/pepD | predicted peptide chain release factor/aminoacyl‑histidine dipeptidase |
Reads supporting (aligned to +/- strand): ref base T (7/6); new base A (6/8); total (13/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.06e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01 |
TGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACG > minE/223561‑223694 | tGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATc > 1:533248/1‑70 (MQ=255) ttATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAaataat > 1:2515293/1‑71 (MQ=255) ttATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAaataat > 1:2070993/1‑71 (MQ=255) tCGAAGGATAATGAAGGAAGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTAc > 1:3207112/1‑71 (MQ=255) aaGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAAttt > 1:2837535/1‑71 (MQ=255) gATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAAtttctt < 1:1569891/71‑1 (MQ=255) gAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAAtttcttt < 1:599184/69‑1 (MQ=255) aaaCTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTcag < 1:880056/71‑1 (MQ=255) aCTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTcagca < 1:1731403/71‑1 (MQ=255) cTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTcagcag > 1:170806/1‑71 (MQ=255) ccGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGt < 1:1481138/69‑1 (MQ=255) aTTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAAtttctt < 1:1763287/51‑1 (MQ=255) ttAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAgtg > 1:3270125/1‑71 (MQ=255) tAAACGCCATCCGGCAGCATATGAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAgtg < 1:3448158/70‑1 (MQ=255) tAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAAtttcttt > 1:1564573/1‑50 (MQ=255) tAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtcagt > 1:342379/1‑62 (MQ=255) tAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGaa > 1:1234902/1‑43 (MQ=255) tAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAAtttctttc > 1:188391/1‑51 (MQ=255) aaCGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTc < 1:1196873/70‑1 (MQ=255) cATCCGGCAGCAAATCAAATAATTACTTCGCCGGAAtttctt < 1:1595159/42‑1 (MQ=255) aTCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCa > 1:3441972/1‑66 (MQ=255) tCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTtcagtcag < 1:2413048/56‑1 (MQ=255) gcagcaAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCt < 1:2855291/71‑1 (MQ=255) gcaTATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACg > 1:3441905/1‑71 (MQ=255) gcaAATCAAATAATTACTTCGCCGGAATTTCTTTcagc < 1:135424/38‑1 (MQ=255) gcaAATCAAATAATTACTTCGCCGGAATTTCTTTcagc < 1:961532/38‑1 (MQ=255) caTATCAAATAATTACTTCGCGGGAATTTCTTTcag < 1:3392643/36‑1 (MQ=255) | TGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACG > minE/223561‑223694 |
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: ‑ |