Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | KT2440 | 5,914,046 | 0 | A | C | 52.4% | 9.1 / 19.9 | 22 | K266T (AAG→ACG) | CJ019_05537 | hypothetical protein |
Reads supporting (aligned to +/- strand): ref base A (7/3); new base C (0/11); total (8/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.03e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.23e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GAACCACGCCGCGAACAGGGCCAGCTGCGCTACCGTATCGATGGCCTGCTGCACACCGTCTACGCCTTCCCCGCCGGGGTCACCCTGGCCCTCGTCAGCCGCCTGAAACACCTGGG‑CCGCATGGACGTTGCAGAGAAGC‑GCCGGCCGCAGGACGGCCGGCTGCAAAGCCGCCTACCCGGTGGTGGCGAAGTGGAGCTGCGGCTTTCGACCCTGCCCACCCCATTTGGCGAGAAACTGGTGTTGCGCCTGTTCGACCCGCACCAGTTGC > KT2440/5913910‑5914177 | gAACCACGCCGCGAACAGGGCCAGCTGCGCTACCGTATCGATGGCCTGCTGCACACCGTCTACGCCTTCCCCGCCGGGGTCACCCTGGCCCTCGTCAGCCGCCTGAAACACCTGGG‑CCGCATGGACGTTGCAGAGaa > 6:396087/1‑137 (MQ=255) cccgcGAACAGGGCCAGCTGCGCTACCGTATCGATGGCCTGCTGCACACCGTCTACGCCTTCCCCGCCGGGGTCACCCTGGCCCTCGGCAGCCGCCTGAAACACCTGGG‑CCGCATGGACGTTGCAGAGAAGC‑gccgg < 8:310326/136‑1 (MQ=255) cgcgAACAGGGCCAGCTGCGCTACCGTATCGATGGCCTGCTGCACACCGTCTACGCCTTCCCCGCCGGGGTCACCCTGGCCCTCGTCCGCCGCCTGAAACACCTGGG‑CCGCATGGACGTTGCAGAGAAGC‑ggcggcc > 4:19582/1‑137 (MQ=255) cgcgAACAGGGCCAGCTGCGCTACCGTATCGATGGCCTGCTGAACACCGTCTACTCCTTCCCCGCCGGGGTCACCCTGGCCCTCGTCAGCCGCCTGAAACACCTGGG‑CCGCATGGACGTTGCAGAGAAGc‑gcccgcc > 8:428321/1‑137 (MQ=255) gcgAACAGGGCCAGCTGCGCTACCGTATCGATGGCCTGCTGCACACCGTCTACGCCTTCCCCGCCGGGGTCACCCTGGCCCTCGTCAGCCGCCTGAAACACCTGGG‑CCGCATGGACGTTGCAGAGAAGC‑ggccgccc > 6:166200/1‑136 (MQ=255) aCAGGGCCAGCTGCGCTACCGTATCGATGGCCTGCTGCACACAGTCTCCGCCTTCCCCGCCGGGGTCACCCTGGCCCTCGTCAGCAGCCGGAAACACCTGGG‑CCGCATGGACGTAGCAGACAAGC‑ggcggcgcacgg > 6:78217/1‑131 (MQ=255) ccGTATCGATGGCCTGCTGCACACCGTCTACGCCTTCACCGCCGGGGTCACCCTGGCCCTCGTCAGCCGCCTCAAACACCTGGG‑CCGCATTGACGTTGCAGAGAAGCGGCGGGGCGCAGG‑CGGCCGCCTGccacccc > 5:250292/1‑131 (MQ=255) gggcgccacaCGGTCACCGCCTGCCCCGCGGGGGCCCCCCTGCCCCTCGGCCCCGCCCTCAAACCCCGGGC‑CCGCATTGACGTTGCAGTGACGG‑GCGGGCCGCTGGCGGGCCGGCTGCAAAGCCGCCTACCCggtgg < 8:425943/131‑1 (MQ=255) tgctgcACACCGTCTACGCCTTCCCCGCCGGGGTCACCCTGGCCCTCGTCAGCCGCCTGAAACACCTGGG‑CCGCATGGACGTTGCAGAGATGC‑GCCGGCCGCCGGACGGCCCGCTGCAAAGCCGCCGACCCcgtggg > 5:207653/1‑136 (MQ=255) gcggcACCCCGTCTACGCCGGCCCCGCCGGGGCCACCCTGCCCCGCGCCCCCCCCCTCAAACCCCTGGG‑CCGCTGGGACGTTGCACACACCC‑GCCGGCCCCCGGCCGCCCGGCTGCAAAGCCGCCTACCCGgtggtg < 6:192562/134‑1 (MQ=255) gcACACCGTCTACGCCTTCCCCGCCGGGGTCACCCTGGCCCTCGTCAGCCGCCTGAAACACCTGGG‑CCGCATGGACGTTGCAGAGAAGC‑gcc > 5:90834/1‑92 (MQ=255) gcACACCGTCTACGCCTTCCCCGCCGGGGTCACCCTGGCCCTCGTCAGCCGCCTGAAACACCTGGG‑CCGCATGGACGTTGCAGAGAAGC‑gcc < 6:90834/92‑1 (MQ=255) tACGCCTCCCCCGCGGGGGACCCCGTGCCCCCCGCCCCCCGCCTCAACCCCGGGGG‑CCGCATGGCGGTTCCACACCCGC‑GCCGGCCGCAGCCCGGCCGGCTGCAAAGCCGCCTACCCGGTGGTGGCGAAGTGGAGCt < 4:195059/137‑1 (MQ=255) accGCCTGCCGCGCGGGGGCCCCCGGGGCGCTGGGCAGCCGCCT‑TAACCCGTGGGCCCGCTTGAAGTTTGCATACACGC‑GCGGGCCGCAGTACGGCCGGCTGCAAAGCCGCCTACCCGGTGGTGGCGAAGTGGAGCt < 8:37824/135‑1 (MQ=255) aCGCCTTCCCCGCCGGGGTCACCCTGGCCCTCGTCAGCCGCCTGAAACACCTGGG‑CCGCATGGACGTTGCAGAGAAGC‑gcc < 2:321885/81‑1 (MQ=255) aCGCCTTCCCCGCCGGGGTCACCCTGGCCCTCGTCAGCCGCCTGAAACACCTGGG‑CCGCATGGACGTTGCAGAGAAGC‑gcc > 1:321885/1‑81 (MQ=255) ccccccCGGGGCCACCCGGCCCCTCGCCCCCCGCCTCAAACACGGGGC‑CCGCAGGGCGGTGCCAGAGACGC‑GCGGGCCGCCGGCCGGCCGGCTGCAAAGCCGCCTACCCGGTGGTGGCGAAGTGGAGCTGCGGCttt < 5:79226/137‑1 (MQ=255) ccgccgGGGGCACCCTGCCCCTCGACCGCCCCCTCAAACACGTGGG‑CCGCAGGGCGGTTGCAGACACGC‑GCCGGCCGCAGGACGGCCGGCTGCAAAGCCGCCTACCCGGTGGTGGCGAAGTGGAGCTGCGGCTTTCg < 6:360688/137‑1 (MQ=255) ccTCGTCCGCCGCCTCAACCACGTGGG‑CCGCATGTCGGTTGCAGTGACGC‑GCCGGCGGCAGGACGGCCGGCTGCAAAGCCGCCTACCCGGTGGTGGCGAAGTGGAGCTGCGGCTTTCGACCCTGCCCACCCCATTTg < 6:249322/137‑1 (MQ=255) cgcccccTCAAACCCCTGGG‑CCGCAGGGCCGTTCCCGAGCCGC‑GCCGGCCGCAGGACGGCCGGCTGCAAAGCCGCCTACCCGGTGGTGGCGAAGTGGAGCTGCGGCTTTCGACCCTGCCCACCCCATTTGGCGAGaa < 7:421888/136‑1 (MQ=255) aaCCCCGGGG‑CGGCTTGTACGTGGCAGACACGC‑CCGGGCGGCCGTCCGGCCGGCTGCAAAGCCGCCTACCCGGTGGTGGCGAAGTGGAGCTGCGGCTTTCGACCCTGCCCACCCCATTTGGCGAGAAACTGGTGTTg < 8:152983/137‑1 (MQ=255) ccGCAGGGCGGTTGCCGACACGC‑GCCGGCCGCAGGCCGCCCGGCTGCAAAGCCGCCTACCCGGTGGTGGCGAAGTGGAGCTGCGGCTTTCGACCCTGCCCACCCCATTTGGCGAGAAACTGGTGTTGCGCCTGTTCg < 8:4065/137‑1 (MQ=255) cgagaCGC‑GCGGGCCGCTGGCCGCCCGGCTGCAAAGCCGCCTACCCGGTGGTGGCGAAGTGGAGCTGCGGCTTTCGACCCTGCCCACCCCATTTGGCGAGAAACTGGTGTTGCGCCTGTTCGACCCGCACCAGTTGc < 4:410705/136‑1 (MQ=255) | GAACCACGCCGCGAACAGGGCCAGCTGCGCTACCGTATCGATGGCCTGCTGCACACCGTCTACGCCTTCCCCGCCGGGGTCACCCTGGCCCTCGTCAGCCGCCTGAAACACCTGGG‑CCGCATGGACGTTGCAGAGAAGC‑GCCGGCCGCAGGACGGCCGGCTGCAAAGCCGCCTACCCGGTGGTGGCGAAGTGGAGCTGCGGCTTTCGACCCTGCCCACCCCATTTGGCGAGAAACTGGTGTTGCGCCTGTTCGACCCGCACCAGTTGC > KT2440/5913910‑5914177 |
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: ‑ |