Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | CP009273 | 1,193,873 | G→T | 100% | intergenic (‑180/‑278) | ymfE ← / → lit | e14 prophage; putative inner membrane protein/e14 prophage; cell death peptidase, inhibitor of T4 late gene expression |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | CP009273 | 1,193,873 | 0 | G | T | 82.2% | 25.3 / 5.1 | 17 | intergenic (‑180/‑278) | ymfE/lit | e14 prophage; putative inner membrane protein/e14 prophage; cell death peptidase, inhibitor of T4 late gene expression |
Reads supporting (aligned to +/- strand): ref base G (1/2); new base T (7/7); total (8/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.04e-01 | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. |
TTTAATTAAGTGTCATTATATACTTTTCAGTAGGTTAGTTACAATTTTTTGTATCTGTTCAGGCTGACCTAGCTTCGCTGACAGACAATATTGTGATCAGTAGCACGTATCGAGGAGGAGTAGCGCTACAAATTTGACGCTGGGTGAGAATCTGAAATTGATAGAAATGAAATAATGAAATAATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATGGTGTCTTACGACATGCTTACCTTA > CP009273/1193728‑1194009 | tttAATTAAGTGTCATTATATACTTTTCAGTAGGTTAGTTACAATTTTTTGTATCTGTTCAGGCTGACCTAGCTTCGCTGACAGACAATATTGTGATCAGTAGCACGTATCGAGGAGGAGTAGCGCTACAAATTTGACGCTGGGTGAGaa < 2:4201640/150‑1 (MQ=255) aCTTTTCAGTAGGTTAGTTACAATTTTTTGTATCTGTTCAGGCTGACCTAGCTTCGCTGACAGACAATATTGTGATCAGTAGCACGTATCGAGGAGGAGTAGCGCTACAAATTTGACGCTGGGTTAGAATCTGAAATTGATAgaaatgaa < 2:4412489/150‑1 (MQ=255) ttttCCGTAGGTTAGTTACAATTTTTTGTATCTGTTCAGGCTGAACTAGCTACGCTGACAGACAATATTGTGATCAGTAGCACGTATCGAGGAGGAGTAGCGCTACAAATTTGACGCTGGGTTAGAATCTGAAATAGATAgaaatgaaat > 2:2196069/1‑150 (MQ=255) ttttCAGGAGGTTAGTTACAATTTTTTGTATCTGTTCAGGCTGACCTAGCTTCGCTGACAGACAATATTGTGATCAGTAGCACGTATCGAGGAGGAGTAGCGCTACAAATTTGACGCTGGGTTAGAATCTGAAATTGATAgaaatgaaat > 2:2695481/1‑150 (MQ=255) ttagttaCAATTTTTTGTATCTGTTCAGGCTGACCTAGCTTCGCTGACAGACAATATTGTGATCAGTAGCACGTATCGAGGAGGAGTAGCGCTACAAATTTGACGCTGGGTTAGAATCTGAAATTGATAGAaatgaaataatgaaataat < 2:1213413/150‑1 (MQ=255) agttaCAATTTTTTGTATCTGTTCAGGCTGACCTAGCTTCGCTGACAGACAATATTGTGATCAGTAGCACGTATCGAGGAGGAGTAGCGCTACAAATTTGACGCTGGGTTAGAATCTGAAATTGATAGAaatgaaataatgaaataatga > 1:983640/1‑150 (MQ=255) agttaCAATTTTTTGTATCTGTTCAGGCTGACCTAGCTTCGCTGACAGACAATATTGTGATCAGTAGCACGTATCGAGGAGGAGTAGCGCTACAAATTTGACGCTGGGTTAGAATCTGAAATTGATAGAaatgaaataatgaaataatga > 1:983631/1‑150 (MQ=255) ttttGTATCTGTTCAGGCTGACCTAGCTTCGCTGACAGACAATATTGTGATCAGTAGCACGTATCGAGGAGGAGTAGCGCTACAAATTTGACGCTGGGTTAGAATCTGAAATTGATAGAAATGAAATAATGAAATAatgaaatgatgaaa < 1:1453698/150‑1 (MQ=255) ttttGTATCTGTTCAGGCTGACCTAGCTTCGCTGACAGACAATATTGTGATCAGTAGCACGTATCGAGGAGGAGTAGCGCTACAAATTTGACGCTGGGTTAGAATCTGAAATTGATAGAAATGAAATAATGAAATAatgaaatgatgaaa < 1:238760/150‑1 (MQ=255) gTATCTGTTCAGGCTGACCTAGCTTCGCTGACAGACAATATTGTGATCAGTAGCACGTATCGAGGAGGAGTAGCGCTACAAATTTGACGCTGGGTTAGAATCTGAAATTGATAGAAATGAAATAATGAAATAatgaaatgatgaaatgat > 1:1618270/1‑150 (MQ=255) gTTCAGGCTGACCTAGCTTCGCTGACAGACAATATTGTGATCAGTAGCACGTATCGAGGAGGAGTAGCGCTACAAATTTGACGCTGGGTTAGAATCTGAAATTGATAGAAATGAAATAATGAAATAatgaaatgatgaaatgatgaaatg > 1:1409338/1‑150 (MQ=255) gTTCAGGCTGACCTAGCTTCGCTGACAGACAATATTGTGATCAGTAGCACGTATCGAGGAGGAGTAGAGCTACAAATTTGACGCTGGGTTAGAATCTGAAATTGATAGAAATGAAATAATGAAATAatgaaatgatgaaatgatgaaatg > 1:1409402/1‑150 (MQ=255) tatTGTGATCAGTAGCACGTATCGAGGAGGAGTAGCGCTACAAATTTGACGCTGGGTGAGAATCTGAAATTGATAGAAATGAAATAATGAAATAATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACaga < 1:1227501/150‑1 (MQ=255) ggaggagTAGCGCTACAAATTTGACGCTGGGTTAGAATCTGAAATTGATAGAAATGAAATAATGAAATAATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATGGTGTc < 1:3006981/150‑1 (MQ=255) aCAAATTTGACGCTGGGTGAGAATCTGAAATTGATAGAAATGAAATAATGAAATAATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATGGTGTCTTACGACATGCTTa > 1:4330112/1‑150 (MQ=255) tttGACGCTGGGTTAGAATCTGAAATTGATAGAAATGAAATAATGAAATAATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATGGTGTCTTACGACATGCTTACCTTa < 1:3882376/150‑1 (MQ=255) tttGACGCTGGGTTAGAATCTGAAATTGATAGAAATGAAATAATGAAATAATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATGGTGTCTTACGACATGCTTACCTTa < 1:2845175/150‑1 (MQ=255) | TTTAATTAAGTGTCATTATATACTTTTCAGTAGGTTAGTTACAATTTTTTGTATCTGTTCAGGCTGACCTAGCTTCGCTGACAGACAATATTGTGATCAGTAGCACGTATCGAGGAGGAGTAGCGCTACAAATTTGACGCTGGGTGAGAATCTGAAATTGATAGAAATGAAATAATGAAATAATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGAAATGATGACAGAGTGTCCAGTGGGCACGGATGGTGTCTTACGACATGCTTACCTTA > CP009273/1193728‑1194009 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |