| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | KT2440 | 381,781 | 0 | A | C | 58.3% | 7.0 / 18.8 | 25 | D201A (GAC→GCC) | soxB | Sarcosine oxidase subunit beta |
| Reads supporting (aligned to +/- strand): ref base A (1/9); new base C (14/0); total (16/9) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.65e-06 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.29e-03 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTC > KT2440/381668‑381901 | gaagaaCACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGCCCTGATCCAGCAAACCGAAGt > 4:75719/1‑135 (MQ=255)gaagaaCACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCCCCGCCCGGGGCGGGGCCCGGATCCAGCAAAccccag > 2:230293/1‑130 (MQ=255)gaagaaCACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCGTCGCCCGCGCTGCCGCCGCCCCGGGCGGCGCCCTGCTCCAGCAAACCGaac > 4:189336/1‑133 (MQ=255)gaagaaCACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTCGTGGCGTTGCCCGTCACGAAGCCGAGGCCTGGGGCCGCGCCCGCGCTGACCACGCCCTGGGCCTCGCCCTGATCAAGCAAACCgcaat > 2:220333/1‑131 (MQ=255) cTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGGTTCGCCCGCGCAGGCCCCCCCCAGGGGCTCGCCCGGGTCCAGCAAACCCAAGGGAACGGCGTccc > 7:273920/1‑134 (MQ=255) cTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACCCCCTGGGCGTCGCCCAGATGCAGCAAACCGCAATGATAGGCTTACGCAAGGaaaaa > 7:227782/1‑134 (MQ=255) tCGGCGCCCCCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCACCAAACCGAAGTGATCGGCTTCCGCAAGGAAAAc < 2:183067/134‑1 (MQ=255) ggCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGGGGCCTGGGGCTTCGCCCGCGCGGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGgc < 6:135107/135‑1 (MQ=255) cgcCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGCCGCCCTGGGCGTCGCCCTGATCCAGCAAACCGCAGTGCTCGGCTTCCGCCAGGAAAACGgcgc > 2:214115/1‑135 (MQ=255) gcCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGCCGCCCTGGGCGGCGCCCCGCTACACGCAACCGCAGTGATCGGCGTCCGCCAGGAAAACGgcgcc > 5:121370/1‑134 (MQ=255) ggggCGGTGCCCGTCACGATGCGGGGGCCTGGGGCTGCGCCGGCGCGGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACca < 1:130719/133‑1 (MQ=255) tggCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCAGCCGCCCCCCGGGGCGGCGCCCCGCTCCAGCAAACCGAAGTGCTCGGCGTCCGCAATGAAAACGGCGCGGTGATCGGTGTGGCAAccc > 4:87025/1‑133 (MQ=255) cGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTc < 4:110899/135‑1 (MQ=255) tGCCGTGGCCGGGGGCTGCGCCCGCGCTGCCGACGCCGTGGGCGTCGACCTGATCCAGCAAACCGAAGTGAGCGGCGTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGgcgc < 7:186343/135‑1 (MQ=255) ccGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGGGGCGCCCGGCTCCAGCAAACCGAAGTGCTCGGCGTCCGCAAGGAAAACCGCGCGGGGATCGGTGTGGGAACCAACAAAGGCGTCATCGgcgc > 1:26148/1‑133 (MQ=255) cGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGCCCGGATCCCGCAAACCGAAGTGCTCCGCTTCCGCAAGGAAAACGGCGCGGGCATCGGTGTGGAAACCAACAAAGGCTTCAACGGCGccca > 7:157229/1‑133 (MQ=255) ccTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGCCCTGATCCAGCCAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCCcgcgcg > 8:167205/1‑135 (MQ=255) ccTGGGGCTTCGCCCGCGCTGCCGACGCCCGGGGCGTCGGCCTGATCCAGCAAACCGCACTGATCGGCTTCCGCCAGGAAAACGGCGCGGTCATCGGTGTGGGAACCAACAAAGGGTTCATCGGCGCCCAGcgcg > 7:12840/1‑135 (MQ=255) gggggCTTCGCCGGGGCGGCCGACGCCCTGGGCGGCGACCTGATCCACCAAACCGAAGGGATCGGCTGCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTc < 7:21363/134‑1 (MQ=255) ccGCGCTGCCGCCGCCCTGGGCCTCGCCCAGATCCACCAAACCGAAGTGATCGGCTTCCGCCAGGCAAACGGCGCCGGCATCCGTGTGGAAACCCACAAAGGCTTCATCGGCGCCAACCGCCTAGGCGTGGtgcc > 8:21363/1‑132 (MQ=255) gcgcTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCAccg < 8:12840/135‑1 (MQ=255) gcTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGACGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCAccgcc < 1:62489/135‑1 (MQ=255) tGCCGACGCCCTGGGCGTCGCCCTGCTCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCCCCGCCgg > 3:187788/1‑135 (MQ=255) tGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCGTCCGCAAGGCAAACGGCGCGGGCATCGGTGTGGAGACCAACAAAGGCTTCATCGGCGGCAAGCGCGTCGGCGTGGTCACCGCCgg > 4:8661/1‑135 (MQ=255) cGCCCTGGGCGGCGACCGGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGGCGGCGTGGTCACCGCCGGTAACTc < 8:250025/135‑1 (MQ=255) | GAAGAACACCCGCTACCCGATCCTCGGCGCCACCGTACAGCGCCGTGGTGGCGTTGCCCGTCACGATGCCGTGGCCTGGGGCTTCGCCCGCGCTGCCGACGCCCTGGGCGTCGACCTGATCCAGCAAACCGAAGTGATCGGCTTCCGCAAGGAAAACGGCGCGGTCATCGGTGTGGAAACCAACAAAGGCTTCATCGGCGCCAAGCGCGTCGGCGTGGTCACCGCCGGTAACTC > KT2440/381668‑381901 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |