| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | ACB122 | 3,440,363 | 0 | A | C | 57.7% | 4.3 / 18.4 | 26 | A24A (GCA→GCC) | CJ019_03266 | hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base A (7/4); new base C (15/0); total (22/4) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.21e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.86e-03 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAGCGCTCGAT > ACB122/3440237‑3440495 | 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) | GCAACGGTCTGGAGGATCGCAAGGCCCTCTACAAGCGAGCGCTCGAGGTGCTGCAGTGAGCGCTTGGGCCACCCGCCTGGTCGCCGTGGCGCTGCTGCTGGCCGCCAGCGCTGGCGGCGCCCGGGCAGCCTGGGTGTGGCAGGCCAACGCCTATAAGGCGCAGCTCGCTGATCAGGCCGATGACTTCGGCAAGCAGTTGGCGGAGAAGGATCGCGCGCACGGTCTGGAGCGGGAGAAGGCTGCAGCAGCAGCGCTCGAT > ACB122/3440237‑3440495 |
| 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: ‑ |