Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | KT2440 | 2,094,808 | 0 | A | C | 51.7% | ‑3.5 / 10.2 | 30 | L186R (CTG→CGG) | CJ019_01995 | hypothetical protein |
Reads supporting (aligned to +/- strand): ref base A (9/5); new base C (12/3); total (21/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.27e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATGCGGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGC‑GCGCC‑TGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAG > KT2440/2094688‑2094936 | aTGCGGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGCAGGCCCGGC‑GCGCC‑GGGGCg > 7:49827/1‑135 (MQ=255) ggCGGTCCAGTCCCGGGGCATGGGCCCGGGCCAACAGTTTCTCAAGATACGTCCGCTGCTCTTCATACCCCTTCGGCTCCGGCTGCTTGCTCAAATACTGCACGCCGCGCTCGACGGCCCGGC‑GCgcc‑gggccgg < 8:49827/134‑1 (MQ=255) gTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCCGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGGCCCGC‑GCCCC‑Ggggcgggcgg > 6:36718/1‑133 (MQ=255) gTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCCGC‑GCGcc‑gggccgggggg > 8:257514/1‑132 (MQ=255) aGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCGGC‑GCGCC‑cgggcgggcggggcg > 7:114798/1‑132 (MQ=255) cATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGCAGGCCCCGC‑GCGCC‑GGGCCGGGCGGGGCGTTTTGGAggg > 2:156974/1‑135 (MQ=255) cATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCGGC‑GCGCCGGGGCCGGGGGGGCGGTTTTgaggg > 3:164835/1‑131 (MQ=255) cATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTAGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGGCACGCGGGGCC‑GGGCCGGGGGGGCGGTTTTGCTgg > 2:131407/1‑135 (MQ=255) gCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCCGC‑GCGCC‑TGGGCGGGGGGGCCGTTTTGCAGGGGCGGGTGGCCg > 4:56357/1‑135 (MQ=255) tcAAGATCCGGCAGGTGCTCTCCATCCCCCTGCGGGCCCGGCTGCTGGCGCAAATCCTCCCCGCCGGGCGCGCCGGCCAGGC‑GCGCG‑GGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTgaa < 4:113981/134‑1 (MQ=255) cAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCGGC‑GCGCC‑TGGCCGGGCGGGCCTTTTTGCAGGGCCGGGTGGGCGAGGGTGTTgaaga > 6:96824/1‑135 (MQ=255) gTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCGGC‑GCGCC‑TGGCCGGGCTGGCCGTTTTTCAGGGCCGGTTGGCCGCGGTTGTTGAAGAACTCtata > 3:34343/1‑135 (MQ=255) tCAGCGGCTCTTCAACCCCCTCCGGCTCCGGCGGCTGGCGCAAATCCTGCACGCCGCGGTCGACGGCCAGGC‑GCGCC‑GGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTatat < 4:78343/135‑1 (MQ=255) ctctTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGC‑GCCCC‑GGGCCGGGGTGGGCGGTTTTCAGGGGCGGTTTGGCGAGGGTTTTGAAGAACTCTAAATgatttct > 3:26591/1‑129 (MQ=255) tctTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGC‑GCGCC‑GGGGCCGGCGGGCCGTTTTTGAGGGCCGGTTTGGCGCGGTTGTTGTAGAAAACAATATtaagcagc > 7:249644/1‑135 (MQ=255) acacCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGC‑GCGCC‑TGGCCGGGCAGGCCTGTTTGCAGGGGCGGTTGGCCGCGGGTGTTGAAGACCTCTATATGGAGCAGCACCa > 2:196788/1‑132 (MQ=255) acacCTTCGGCTCCGGCTGCTTGCGCAAATAATGCACGCCGCGCTCGACGGCCCGGC‑GCGAC‑TGGGCGCGCTGGACGTTTTGCAGGGCCGTTGGGCGGTGGGTGTTGAAGTAATGTCTATTCAGCAGCACCAGGa > 8:216284/1‑135 (MQ=255) ggCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCGGC‑GCGCC‑TGGCCGGGCGGGCCGTTTTGGAGGGCCGGTTGGGCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCCGGATATGgggg > 3:189123/1‑135 (MQ=255) cTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCCGC‑GCGCC‑TGGCCGGGCTGGCCTTTTTGCACGGCCGGTTGGCCGCGGTTGTTGAAGAACGCTATATGCAGCAGCACCAGGAGATGGCGGag > 1:47851/1‑134 (MQ=255) tgctTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGC‑GCGCC‑TGGGCGGGCTGGCCGTTTTTGAGGGGCGGGTGGGCGAGGTTGTTGAAGAAATCTATATTGAGCACCCCCAGGAGTTTGCGGAttttcatg > 8:18551/1‑130 (MQ=255) ctggcgcAAATACTCCCCGCCCGGCTCGACGCCCAGGC‑GCGCC‑GGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGAc < 8:66369/132‑1 (MQ=255) cgcAAATACTGCACGCCGCGCTCGAAGGCCAGGC‑GCGCC‑TGGCCGGGCTGGCCTTTTTGCAGGGCCGGGTGGCCGGGGTTGTTGAAGAAATCTATATGCCCCCGCACCAGGATATGGCGGATTTCCTTGccccca > 6:289665/1‑133 (MQ=255) gccgcGCTCGAAGGCCAGGC‑GCGCC‑TGGCCGGGCTGGCCGTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGTAGAACTCTATATGCAGCAGCCCCAGGATATGGCGGATTTCCTTGACCCAATGGTTGGCCTCg > 2:235917/1‑133 (MQ=255) gccgcGCTCGAAGGCCAGGC‑GCGCC‑TGGCCGGGCTGGCCGGTTTGCCGGGGCTGTTGGGCGAGGTTGTTGAAGAACTCTATATTCAGCAGCCCCAGGATATGGCGGCTTTCCTTGACCCAATGCTTGGCcgcg > 2:232145/1‑133 (MQ=255) gccgcGCTCGAAGGCCAGGC‑GCGCC‑GGGCCGGGCTGGCCGTTTTTCAGGGCCGGGTGGGCGAGGTTGTTGAAGAACTCTATATGCAGCAGCCCCAGGATATGGCGGGTTTCCTTGACCCAATGCTTGGCCTCGAc > 1:19453/1‑135 (MQ=255) cgcgcgcgACGGCCCGGC‑GCGCC‑GGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGAccg < 7:23226/135‑3 (MQ=255) gcGAGGGCCCGGC‑GCGCG‑GGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGc < 6:95901/134‑1 (MQ=255) gccggCCGGGC‑GCGCG‑GGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCaa < 3:222618/132‑1 (MQ=255) acggCCAGGC‑GCGCC‑GGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCaag < 4:183210/133‑1 (MQ=255) acgcccAGGC‑GCGCC‑GGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCaag < 2:246526/131‑1 (MQ=255) | ATGCGGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGC‑GCGCC‑TGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAG > KT2440/2094688‑2094936 |
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: ‑ |