Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | KT2440 | 4,143,814 | 0 | A | C | 57.9% | 3.4 / 14.1 | 19 | intergenic (+168/‑229) | lrp_3/CJ019_03903 | Leucine‑responsive regulatory protein/hypothetical protein |
Reads supporting (aligned to +/- strand): ref base A (3/5); new base C (11/0); total (14/5) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.82e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.68e-02 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTACCCCAAGGCCTTCAGCCAGCTGGCACTCTCCTTGCTACCGTCACCCGCTTAAGCTTTTGTCAGCCTTAGTTTTTTTAACGCTTATGGGTACATATTTACTAATTTCT > KT2440/4143688‑4143938 | ccTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAACCGCGCACATCGCACCCACCCAGCCCACCGTGGTCCGGTTGCCCACCCccccgc > 2:168374/1‑137 (MQ=38) tGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACAACGCACCAACCCAGCCCACCGTGAACCGGTTGCCACCCCCCCCGccc > 7:42855/1‑137 (MQ=255) aGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCCTGATCCGGTTGCCCCCCCCGCccccccgcga > 2:58903/1‑133 (MQ=255) aGGAGCGGGTTTACACGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCACCCCAGCCCCCCGGGCTCCGGGTGCCCCCCCCGCCGCCCCGCGCCCCCCACGCCGTcagccac > 6:8058/1‑136 (MQ=37) ggTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTACCCCAAGGCCTTCAGCCAGCTGGCAc < 5:8058/137‑1 (MQ=255) agGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTACCCCAAGGCCTTCAGCCAGCTGGCACTCTCCTTGCTACCGt < 5:281954/137‑1 (MQ=255) ggccCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGCTCCGGTTGCCCCCCCCGCCGCCCCGCTACCCCAAGGCCGTCCGCCCGCTGGGCCGCTCCGTGCTACCGTCcccc > 3:81450/1‑137 (MQ=255) aCCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCCGTTGCCACCCCCCCCGCCCCTATAACCCAACGCCGTAAGCCAGCAGGCACGCTCAGTTCTAACGTGACCCCCTTAAGcat > 5:62078/1‑135 (MQ=255) ccTTTCATCTGCCCTGAAAACGCGCAGATCGCACCAACCCAGCCCACCTTGATCGGGTTGCGACCCCCGCCGCCCCTCTACGCCAAGGCCTTCAGCGAGCTGGCACTCTCCTTGCTACCGTCTCCCGCTTAAGCttt < 5:180409/137‑1 (MQ=255) tGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTACCCCAACGCCTTCATCCAGCTGGCACTCTCCTTGCTACCGTCACCCGCTTAAGCTTTTGTCAGCCTTAGtt > 8:241102/1‑137 (MQ=255) aaaCGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTACCCCAAGGCCTTCAGCCAGCTGGCACTCTCCTTGCTACCGACACCCGCTTAAGCTTTTGTCAGCCTTAGttttt < 3:87352/137‑1 (MQ=255) aaCGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCCCCCCCCCCGCCCCTCTACCCCAAGGCCTTCAGCCAGCTGGCACTCTCCTTGCTACCGCCACCCGCATAAGCGTTTGTCAGCCTTAGtttttt > 4:87352/1‑137 (MQ=255) caTCGCACCAACCCAGCCCACCGTGATCCGGTTGCCCCCCCCGCCGCCCCTCTACCCCAAGGCCGTCCGCCAGCTGGCCAGCTCCTTGCTACCGTCACCCCCGTAAGCTTTTTTCAGGCGTAGTTTTTTTAACGCtt > 2:350179/1‑137 (MQ=255) gCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTACCCCAAGGCCTTCAGCCAGCTGGCACTCTCCTTGCTACCGTCACCCGCTTAAGCTTTTGTCAGCCTTAGTTTTTTTAACGCTTATgg < 6:55835/137‑1 (MQ=255) aCCAACCCAGCCCACCGTGATCCGGTTGCCCCCCCCGCCGCCCCTCTACCCCAAGGCCATCAGCCCGCTGGCACTCTCCGTGCGACCGTCACCCGCGTAAGCTTTTGTCAGCCGTAGTTTTTTTAACGCTTATgggg > 2:7579/1‑136 (MQ=255) aGCCCACCGTGATCCGGTTGCCCCCCCCCCCGCCCCTCTACCCCAACGCCTTCAGCCAGCTGGCACTCTCCTTGCTACCGTCACCCCCGTAAGCTTTTTTCAGCCTTAGTTTTTTTTAAGCTTATGGGTACATAttt > 5:55835/1‑137 (MQ=255) aGCCCACCGTGATCCGGTTGCCCCCCCCCCCGCCCCTCAACCCCAACGCCATAAGCCACCTGGCACTCTCCATGCTACCGTACACCGCGTAAGCTTTTTTAAGACTTAGTTTTTTTAACGCTTATTGGTACATtttt > 6:110739/1‑137 (MQ=255) aCCGTGATCCGGTTGCCCCCCCCCCCCCCCCGCGACCCCCAGGCCTTCAGCCACCTGGCCCAATCCTTGCTACCGTCACCCGCTTAAGCAGTTGTGCGCCTTAGTTTTTTTTAAGCTTAAGGGTAAATAATtaataa > 4:3062/1‑137 (MQ=255) gATCCGGTTGCCCCCCCCCCCGCCCCTCTACCCCAAGGCCTTAAGCCAGCTGGCACTCTACTTTCTACCGTCACCCGCGTAAGCTTTTTTAAGCCTTAGTTTTTTTAAAGCTTATGGGTAAATATTTAATAAtttat > 2:211210/1‑135 (MQ=255) | CCTGCAAAGGTGACGTAGCCCGTCGCAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCTTTCATCTGCCCTGAAAACGCGCACATCGCACCAACCCAGCCCACCGTGATCCGGTTGCCACCCCCGCCGCCCCTCTACCCCAAGGCCTTCAGCCAGCTGGCACTCTCCTTGCTACCGTCACCCGCTTAAGCTTTTGTCAGCCTTAGTTTTTTTAACGCTTATGGGTACATATTTACTAATTTCT > KT2440/4143688‑4143938 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |