Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | ACB122 | 3,440,367 | 0 | T | G | 68.0% | ‑4.7 / 14.5 | 26 | W26G (TGG→GGG) | CJ019_03266 | hypothetical protein |
Reads supporting (aligned to +/- strand): ref base T (4/4); new base G (17/0); total (22/4) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.53e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.20e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAGCGCTCGATCAGCTGGC > ACB122/3440237‑3440503 | gCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGgtgt > 5:248898/1‑137 (MQ=255) gCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCAGGgtgt > 6:449134/1‑137 (MQ=255) ggaggaTCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGCAGGCCaa > 3:378758/1‑137 (MQ=255) cAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCCGGCACCCGGGGTGTGGCAGGCCACCGCCTATaa > 8:8440/1‑137 (MQ=255) cAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCCCGGGCGGCGCCCCGGCAGCCGGGGGGTGGCAGGCCAACGCCTATaa > 4:125807/1‑137 (MQ=255) aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGGCGCCCGCCCAGGCGGCGCCCCGGCCGCCGGGGGGTGGCCGGCCCCAGCCGATAAg > 2:123152/1‑137 (MQ=255) aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCAGCCTGGGGGTGGCAGGCCCACCCCGATaaa > 2:123148/1‑136 (MQ=255) aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCCCGGGCGGCGCCCGGGCCGCCGGGGTGTGGGAGGCCAACCCCGATaaa > 3:242415/1‑136 (MQ=255) aaGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCACCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGCAGGGCACCGCCGATAAg > 5:232086/1‑137 (MQ=255) tACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCAGCCGGGGTGTGGCAGGCCAACGCCTATAATGCGCAGCTc > 6:72931/1‑137 (MQ=255) aCAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGCTGGCCAACGCCGAAAAGGCGCAGCTCg > 5:30583/1‑137 (MQ=255) aaGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCCCGCTGGCGGCGCCCGCGCCGCCGGGGTGTGGCACGGCAACGCCGATAATGCGCAgcgcgcg > 1:329454/1‑136 (MQ=255) cccgcgcgcTCGTGGTGCGGCAGTTCCGGCTTGGGCCCCCCGGGGGGGCGCGGGGGCGCTGCTGCGGGCCCCCCGGGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTg < 8:194619/132‑1 (MQ=255) gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCGGGCGGCGCCCGGGCCGCCGGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGc > 8:136680/1‑137 (MQ=255) gcTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGACGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCACCGCCTATAAGGCGCAGCTCGCTGATCAGGc > 1:127766/1‑137 (MQ=255) aGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCCGCCGGGGGGTGGGAGGCCCACGCCGATAAGGCGCCGCTCGCTGATCAGGCCGCTg > 6:398540/1‑137 (MQ=255) tgctgcAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCGGGGTGTGGGCGGCCAACGCCGATAAGGCGCAGCGCGGGGATCAGGCCGCTGAcg > 4:191546/1‑136 (MQ=255) gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGTGTGGCCGGCCAACCCCGATAAGGCGCCGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGtt > 8:108020/1‑137 (MQ=255) gcTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGGGTGGCCGGCCCACCCCGATAAGGCGCCGCTCGCTGCTCAGGCCGCTGACGTCGGCAAGCAGtt > 6:225777/1‑137 (MQ=255) cccGCCTGGTCGCCGTGGCGCTGCTGCTGGGCGCCCGCCCTGGCGGGGCCCCGGCCCCCGGGGGGTGGGGCGCCCACGCCTATAACGCGCAGCACGCCGCTCAGGCCGCGGTCATCGGCCAGCACTTGGGGGGGaaa > 4:349608/1‑136 (MQ=255) gCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCCGCGCTGGCGGCGCCCGGGCCGCCGGGGGGTGGGCGGCCCACGCCGATAAGGCGCCGCTCGCTGCTCAGGCCGCTGACTTCGGCAAGCAGTTGGGGGGGAAGGAt > 6:182620/1‑137 (MQ=255) ggcgcgGTGGCGCTGCTGCGGGCCCCCCGGGCGGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATcgcgc < 8:232619/135‑1 (MQ=255) gTGGCGCTGCTGCTGGCCGCCCGCGCGGGCGGCGCCCGGGCCGCCTGGGGGTGGGCGGCCCACGCCGCTAAGGCGCCGCGCGGGGGTCAGGCCGAGGGCTTGGGCGAGCGGTTGGGGGGGAAGGGTcgcgcgggcgg > 2:50988/1‑132 (MQ=255) gccgcggctgGCCGGGGGGGCCGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTgg < 1:92086/131‑1 (MQ=255) ggggCCGCCCGCGGTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGga < 6:373163/135‑1 (MQ=255) gggCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGAAGCAGCGCTCGAt > 5:91857/1‑137 (MQ=255) tGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAGCGCTCGATCAgctggc < 7:136680/137‑1 (MQ=255) | GCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAGCGCTCGATCAGCTGGC > ACB122/3440237‑3440503 |
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: ‑ |