Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 184,245 | 0 | A | T | 31.6% | 21.0 / 10.4 | 19 | intergenic (‑150/+12) | yaeH/yaeI | UPF0325 family protein/phosphodiesterase with model substrate bis‑pNPP |
Reads supporting (aligned to +/- strand): ref base A (9/4); new base T (0/6); total (9/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.08e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TCAGTAAATTGTGCTTATTTTAGCATTTGGCCTGGCCCGCCCGGCTCGACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATCTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCAAAGGCATTTAATCCGGCGACGTAACGTTTATCTT > NC_000913/184115‑184387 | tCAGTAAATTGTGCTTATTTTAGCATTTGGCCTGGCCCGCCCGGCTCGACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATCTacac > 8:36889/1‑148 (MQ=255) aTTTTAGCATTTGGCCTGGCCCGCCCGGCTCGACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTTTTTTTTTTATTGTTTTGTTTAAATTCTAGATCTACACCAGTTCCAGCATCGt < 2:132259/147‑1 (MQ=255) ttAGCATTTGGCCTGGCCCGCCCGGCTCGACCGTGATCGCAGAAATCCTGTTTCTTCCAGACTTCTATCATAGCGTAAAAATCAAAAATTTTTTTTTTTTATTGTTTTGTTTAAATTCTAGATCTACACCAGTTCCAGCATCGTTACtt < 8:59295/149‑1 (MQ=255) aTTTGGCCTGGCCCGCCCGGCTCGACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTTATTTATTGTGGTGATTAAAATCACGATCTACACCAGTCCCAGCATCGTTACTTCCggg < 2:88217/149‑1 (MQ=255) ccgcccgGCTCGACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATCTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGa < 1:76922/149‑1 (MQ=255) cAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATCTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCcca > 5:48029/1‑149 (MQ=255) aGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATCTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCccac < 7:36889/149‑1 (MQ=255) ttttCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATCTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGGCCACGCCACGGGt > 1:30810/1‑149 (MQ=255) ttcttcCAGACTTCTATCATGGCGTAAAAATCAAAAAATTATTTTTTTTATTGTTTTGATTAAATTCTAGATCTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGGCCACGCCACGGGtt < 2:30810/148‑1 (MQ=255) tcCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATCTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTtgtg > 6:54087/1‑147 (MQ=255) cAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATCTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTAAGACGCAAACCATACAAAATGACAAAGCCAAGGGTtgag > 7:56657/1‑143 (MQ=255) tGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATCTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGcc < 5:54087/149‑1 (MQ=255) aaaaTCAAAAAATTTTTTTTTTTATTGTTTTGTTTAAATTCTAGATCTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGAt < 4:120282/130‑1 (MQ=255) aaaaTCAAAAAATTTTTTTTTTTATTGTTTTGTTTAAATTCTAGATCTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGAt < 4:122284/130‑1 (MQ=255) aaaaTCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCCCGATCTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGAt > 3:120282/1‑130 (MQ=255) aaaaTCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATCTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCCCGCCACGGGTGGTGTAGAt > 3:122284/1‑130 (MQ=255) ggTGATAAAATTCACGCTCTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCAAAGGCATTTAATCCGGCGACGTAACGTTTATCtt > 5:126123/1‑149 (MQ=255) gATAAAATTCACGATCTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCAAAGGCAtt < 5:77875/120‑1 (MQ=255) gATAAAATTCACGATCTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCAAAGGCAtt > 6:77875/1‑120 (MQ=255) | TCAGTAAATTGTGCTTATTTTAGCATTTGGCCTGGCCCGCCCGGCTCGACTGTGATCGCAGAAAGCCTTTTTCTTCCAGACTTCTATCATGGCGTAAAAATCAAAAAATTACCTGCTTTATTCTGGTGATAAAATTCACGATCTACACCAGTTCCAGCATCGTTACTTCCGGGCGGCAGTTCAGACGCAAACCATACAAACTGCCCACGCCACGGGTTGTGTAGATGTGTCTTTCGCCAAAGGCATTTAATCCGGCGACGTAACGTTTATCTT > NC_000913/184115‑184387 |
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: ‑ |