Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | KT2440 | 1,666,173 | 0 | A | C | 52.9% | 5.8 / 13.7 | 19 | D54A (GAC→GCC) | xerD | Tyrosine recombinase XerD |
Reads supporting (aligned to +/- strand): ref base A (3/5); new base C (9/0); total (14/5) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 9.05e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.66e-02 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TGGAAAAAGGCCTGTCCGACAACACCCGCGTGTCCTATCGCAGCGACCTGGCCCTGTTCAATGGCTGGTTGCAGGAGCATTCGGTATCGCTGCCGGACGCCGGTCGCGACCTGATCCTCGATCACTTGGCCTGGCGCCTGGATCAGGGCTACAAGCCACGTTCCACGGCGCGTTTCCTCTCCGGCCTGCGCGGCTTCTTCCGTTACCTGCTGCGGGAAAAGCTGGTGGC > KT2440/1666065‑1666293 | tGGAAAAAGGCCTGTCCGACAACACCCGCGTGTCCTATCGCAGCGACCTGGCCCTGTTCAATGGCTGGTTGCAGGAGCATTCGGTTTCGCTGCCGGACGCCGGTCGCGACCTGATCCTCGATCACTTGGCCTGgc < 7:100820/135‑1 (MQ=255) aaaGGCCTGTCCGACAACACCCGCGTGTCCTATCGCAGCGACCTGGCCCTGTTCAATGGCTGGTTGCAGGAGCATTCGGTATCGCTGCCGGACGCCGGTCGCGCCCTGCGCCGCGCTCACGTGGGCTGGCGcccg > 1:98476/1‑133 (MQ=255) aaaGGCCTGTCCGACAACACCCGCGTGTCCTATCGCAGCGACCTGGCCCTGTTCAATGGCTGGTTGCAGGAGCATTCGGTATCGCTGCCGGACGCCGGTCGCGACCAGATCCTCGATCACTTGGCCTGGCGCCt > 4:96791/1‑134 (MQ=255) gACAACACCCGCGTGTCGTTTCGCAGCGACCTGGCCCTGATCAATGGCTGGTTGCAGGAGCATTCGGTATCGCTGCCGGACGCCGGTCGCGACCTGATCCTCGATCACTTGGCCTGGCGCCTGGATCAGGGCTAc < 6:23492/135‑1 (MQ=255) cgTGTCCTATCGCAGCGACCTGGCCCTGTTCAATGGCTGGTTGCAGGAGCATTCGGTATCGCTGCCGGACGCCGGTCGCGACCTGATCCTCGATCACTTGGCCTGGCGCCTGGATCAGGGCTACAAGCCACGTTc < 6:89040/135‑1 (MQ=255) ccTGGCCCTGTTCAATGGCTGGTTGCAGGAGCATTCGGTATCGCTGCCGGACGCCGGTCGCGGCCAGATCCTCGATCACTTTGCCTGGCGCCTGGATCCGGGCTACAAGCCACGTTCCACGGCGCGGTTCCtctc > 3:39738/1‑135 (MQ=255) tGGCCCTGTTCAATGGCTGGTTGCAGGAGCATTCGGTATCGCTGCCGGACGCCGGTCGCGGCCTGCTGCTCGATCACTTGGCCTGGCGCCGGGATCAGGGCTACAAGCCACGTTCCACGGCGCGGTTGGTCTCCg > 5:60116/1‑135 (MQ=255) tGGCCCTGTTCAATGGCTGGTTGCAGGAGCATTCGGTATCGCTGCCGGACGCCGGTCGCGCCCTGCTCCTCGATGACTTGGCCTGGCGCCTGGATCAGGGCGACAAGCCACGTTCCACGGCGCGTTTCCTCTACg > 7:143804/1‑135 (MQ=255) cTGTTCAATGGCTGGTTGCAGGAGCATTCGGTATCGCTGCCGGACGCCGGTCGCGCCCTGATCCGCGCTCACTTTGCCTGGCGCCGGGGTCAGGGCGACAAGCCCCGTTCCACGGCGCGTTTACTCTCCGGccgg > 1:125806/1‑133 (MQ=255) aaTGGCTGGTTGCAGGAGCATTCGGTATCGCTGCCGGACGCCGGTCGCGCCCAGATACTCGCTCCAATTGCCTGGCGCCTGGCTCAGGGGGACAACCCCAGTTCCACGGCGCGGTTTCACAGCGGCCAGCCCGGc > 6:147604/1‑135 (MQ=255) cTGGTTGCAGGAGCATTCGGTATCGCTGCCGGACGCCGGTCGCGCCCTGCTCCTCGCTCACTTGGCCTGGCGCCTGGATCAGGGCTACCACCCCCCTTCCACGGCGCGGTTGCTCTCCGGCCTGCGCGGcttcgt > 2:3720/1‑133 (MQ=255) gTTGCAGGAGCATTCGGTATCGCTGCCGGACGCCGGTCGCGCCCTGATCCTCGATCACTTGGCCTGGCGCCTGGATCAGGGCTACAAGCCACGTTCCACGGCGCGTTTCCTCTCCGGCCTGCGCGGCTTCTTCCg > 8:30568/1‑135 (MQ=255) tGCAGGAGCATTCGGTATCGCTGCCGGACGCCGGTCGCGCCCTGATCCTCGATCACTTGGCCTGGCGCCTGGATCAGGGCTACAAGCCCCGTTCCACGGCGCGTTTCCTCTCCGGCCTGCGCGGCTTCTTCCGtt > 4:159668/1‑135 (MQ=255) tGCAGGAGCATTCGGTATCGCTGCCGGACGCCGGTCGCGCCCTGATCCTCGATCACTTGGCCTGGCGCCTGGATCAGGGCTACAAGCCACGTTCCACGGCGCGGTTCCGCTCCGGCCTGCGCGGCGTCTTCCGtt > 5:206290/1‑135 (MQ=255) tGCAGGAGCATTCGGTATCGCTGCCGGACGCCGGTCGCGACCTGATcc > 4:177011/1‑48 (MQ=255) tGCAGGAGCATTCGGTATCGCTGCCGGACGCCGGTCGCGACCTGATcc < 3:177011/48‑1 (MQ=255) tATCGCTGCCGGACGCCGGTCGCGCCCTGCTGCTCGATCACTTGGCCTGGCGCCTGGATCAGGGCTACAAGCCCCGGTGCACGGCGCGGTTCCTCGCCGGCCTGCGCCGCTTATTGCGTTTCCTGCTGCgggc > 8:10017/1‑132 (MQ=255) ccGGACGCCGGTCGCGACCTGATCCTCGATCACTTGGCCTGGCGCCTGGATCAGGGCTACAAGCCACGTTCCACGGCGCGTTTCCTCTCCGGCCTGCGCGGCTTCTTCCGTTACCTGCTGCGGGAAAAGCtggtg > 7:19640/1‑135 (MQ=255) ggACGCGGGTCGCGACCTGATCCTCGATCACTTGGCCTGGCGCATGGATGAGGGCTACAAGCGAGGGTACACGGCGAGTTTCGTCTCAGGCCTGCGCGGCTTCTTCCGTTACCTGCTGCGTGAAATGGTGGTGGc < 7:156210/135‑1 (MQ=255) | TGGAAAAAGGCCTGTCCGACAACACCCGCGTGTCCTATCGCAGCGACCTGGCCCTGTTCAATGGCTGGTTGCAGGAGCATTCGGTATCGCTGCCGGACGCCGGTCGCGACCTGATCCTCGATCACTTGGCCTGGCGCCTGGATCAGGGCTACAAGCCACGTTCCACGGCGCGTTTCCTCTCCGGCCTGCGCGGCTTCTTCCGTTACCTGCTGCGGGAAAAGCTGGTGGC > KT2440/1666065‑1666293 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |