Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913_3_pae_tpiA | 1,202,957 | 0 | G | . | 25.0% | 36.5 / 12.0 | 16 | intergenic (‑24/‑67) | cohE/croE | e14 prophage; repressor protein phage e14/e14 prophage; putative DNA‑binding transcriptional regulator |
Reads supporting (aligned to +/- strand): ref base G (6/6); new base . (2/2); total (8/8) | |||||||||||
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.87e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTT > NC_000913_3_pae_tpiA/1202828‑1203083 | tCTGGGTTGGCTCGCTATCAATCAGAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGa < 4:110645/150‑1 (MQ=255) cATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGACTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGc > 1:141477/1‑150 (MQ=255) cTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTa > 1:147052/1‑129 (MQ=255) cTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTa < 2:147052/129‑1 (MQ=255) aTGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGaa > 3:100567/1‑60 (MQ=255) aTGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGaa < 4:100567/60‑1 (MQ=255) gATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTg > 4:32019/1‑150 (MQ=255) aaCGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGaa < 1:183837/123‑1 (MQ=255) aaCGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGaa > 2:183837/1‑123 (MQ=255) aaCGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGtt < 3:32019/150‑1 (MQ=255) aaaCCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACAcc > 1:135100/1‑65 (MQ=255) aaaCCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACAcc < 2:135100/65‑1 (MQ=255) ttGTACC‑TTAACGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAg < 1:21139/63‑1 (MQ=38) ttGTACC‑TTAACGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAg > 2:21139/1‑63 (MQ=38) ttGTACC‑TTAACGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCaaa > 1:40934/1‑104 (MQ=255) ttGTACC‑TTAACGGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCaaa < 2:40934/104‑1 (MQ=255) | TCTGGGTTGGCTCGCGATCAATCATAGTGGCAAAGGAAGAATTCCCGCCAACACCATCTCTCAGTTTTCTGGCGTTAGACCGCCGGATGTCATGGATTGTTTTCATAACGAAATTAAAACCCTTGTACCGTTAA‑GGTACAAGTATCTTGAAGGTTCATTTCAATCATGTAATATGTACACCGGAGGTACATATTGTATGAAAGCGTATTGGGACTCTTTAACCAAAGAACAGCAGGGCGAGTTGGCCGGAAAAGTT > NC_000913_3_pae_tpiA/1202828‑1203083 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |