Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 1,602,012:1 | +A | 100% | coding (853/1305 nt) | yfeW → | predicted periplasmic esterase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 1,602,012 | 1 | . | A | 84.6% | 29.6 / 1.9 | 13 | coding (853/1305 nt) | yfeW | predicted periplasmic esterase |
Reads supporting (aligned to +/- strand): ref base . (2/0); new base A (11/0); total (13/0) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.39e-01 | |||||||||||
Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch. |
ATCCATTTCCCGAATATCCGCACCTCCACTCTCTGGGGTCAGGTGCACGATG‑AAAAAGCCTTTTATTCGATGGGCGGCGTTTCCGGG > minE/1601961‑1602047 | aTCCATTTCCCGAATATCCGCACCTCCACTCTCTGGGGTCAGGTGCACGATG‑AAAAAGCCTTTTATTCGa > 1:640958/1‑70 (MQ=255) tCCATTTCCCGAATATCCGCACCTCCACTCTCTGGGGTCAGGTGCACGATG‑AAAAAGCCTTTTATTCGATg > 1:1746134/1‑71 (MQ=255) cGCACCTCCACTCTCTGGGGTCAGGTGCACGATGAAAAAAGCCTTTTATTCGATGGGCGGCGTTTCCggg > 1:1546306/1‑70 (MQ=255) cGCACCTCCACTCTCTGGGGTCAGGTGCACGATGAAAAAAGCCTTTTATTCGATGGGCGGCGTTTCCggg > 1:1644358/1‑70 (MQ=255) cGCACCTCCACTCTCTGGGGTCAGGTGCACGATGAAAAAAGCCTTTTATTCGATGGGCGGCGTTTCCggg > 1:1674577/1‑70 (MQ=255) cGCACCTCCACTCTCTGGGGTCAGGTGCACGATGAAAAAAGCCTTTTATTCGATGGGCGGCGTTTCCggg > 1:2656268/1‑70 (MQ=255) cGCACCTCCACTCTCTGGGGTCAGGTGCACGATGAAAAAAGCCTTTTATTCGATGGGCGGCGTTTCCggg > 1:2780840/1‑70 (MQ=255) cGCACCTCCACTCTCTGGGGTCAGGTGCACGATGAAAAAAGCCTTTTATTCGATGGGCGGCGTTTCCggg > 1:317654/1‑70 (MQ=255) cGCACCTCCACTCTCTGGGGTCAGGTGCACGATGAAAAAAGCCTTTTATTCGATGGGCGGCGTTTCCggg > 1:479529/1‑70 (MQ=255) cGCACCTCCACTCTCTGGGGTCAGGTGCACGATGAAAAAAGCCTTTTATTCGATGGGCGGCGTTTCCggg > 1:49296/1‑70 (MQ=255) cGCACCTCCACTCTCTGGGGTCAGGTGCACGATGAAAAAAGCCTTTTATTCGATGGGCGGCGTTTCCggg > 1:642915/1‑70 (MQ=255) cGCACCTCCACTCTCTGGGGTCAGGTGCACGATGAAAAAAGCCTTTTATTCGATGGGCGGCGTTTCCggg > 1:871978/1‑70 (MQ=255) cGCACCTACACTCTCTGGGGTCAGGTGCACGATGAAAAAAGCCTTTTATTCGATGGGCGGCGTTTCCggg > 1:2456977/1‑70 (MQ=255) | ATCCATTTCCCGAATATCCGCACCTCCACTCTCTGGGGTCAGGTGCACGATG‑AAAAAGCCTTTTATTCGATGGGCGGCGTTTCCGGG > minE/1601961‑1602047 |
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: ‑ |