| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | ACB122 | 1,123,238 | C→A | 100% | intergenic (‑240/+64) | gcvH_1 ← / ← CJ019_01062 | Glycine cleavage system H protein/hypothetical protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | ACB122 | 1,123,238 | 0 | C | A | 91.7% | 14.5 / ‑3.7 | 20 | intergenic (‑240/+64) | gcvH_1/CJ019_01062 | Glycine cleavage system H protein/hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base C (1/1); new base A (17/0); total (19/1) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.05e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.09e-01 | |||||||||||
| Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
| Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
ACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGA > ACB122/1123113‑1123331 | aCGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCTacgac > 3:49404/1‑134 (MQ=255) gTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACACTTGTGCACCGCCGGCGAGGAGGCTGTGTATCGTGTTATCTaagacgca > 2:160718/1‑133 (MQ=255) gAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGGGACGACGTAGCGCAAACGCTACATGTGGTGgctgggctgg > 4:223334/1‑133 (MQ=255) aaTGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTTGCTACGACGTTGCGCAAACCATAAATGTGGTGCCGGGAcatgcgt > 3:356343/1‑129 (MQ=255) aaTGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGGGACGGACTATCGCAAACGAAAAATGTTGTGGCGGGGCTGGCGa > 8:426958/1‑135 (MQ=255) aaTGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCTACGCCGTAGCGCAAAAGATACATGTGGTGGCGGGACTGGGGa > 1:45851/1‑135 (MQ=255) cGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGGGGAGGCTGTGTATCGTTTTAGGTAAGCAGTTGCGCAAACGTTTAAAGTGGTGGCTTGACTGGCGAGGTGCCaa > 3:131230/1‑135 (MQ=255) cGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTATGTTCGCCGTGGCGGAACAGATACCAGGGGTGCCGGGCGTGGCGGGGTGCCaa > 2:24853/1‑135 (MQ=255) cGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCTTAGGCGTTGCGCAAAACCGACAAGTGGTGGGGGGGCGGGCGAGGTGGCaa > 5:199868/1‑135 (MQ=255) cAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCGAAGCCGTTGCCGAAACGCCACACGTGGTGGCGGGGCTGGCGCGGGTCCAAAGGcc > 8:63681/1‑135 (MQ=255) cAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGGTTTAGGGAAGTCGGATCGACAAAGAAACCAGTGGTGCCGGGGCGGGCGAGGTGGCCATGGcc > 8:236532/1‑135 (MQ=255) ggTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCTACGACCTAGCCAAAACGATACAAGTGGTGGCAGGACTGGCGCGGTGGCAAAGGccgt > 2:31311/1‑133 (MQ=255) ggTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGGTTTAGCTACGGCGGTGGGCGAACGCGACATGTGGGGGCGGGGCTGGGGAGGGGGCCATGGccgt > 3:397711/1‑133 (MQ=255) gggCCTCTGAGAGCAAACGGGTTTGAGCATAAAACCCTTGTGCACCGCGGGCGAGGAGGCTGTGTATCGTTTTCGATCCGCCGTAGCTAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTc < 4:49404/133‑1 (MQ=255) ctGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATAAATGTGGTGCCGGGCCTGGCGAGGTGCCAATTGCCATCGGTCAt > 4:331628/1‑135 (MQ=255) ctGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCGTCGACGTAGGGCAAACGCAACATGTGGTGGCTGGACTGGCGTGGTGGCAAAGGCCTTCGGtgct > 8:127044/1‑132 (MQ=255) ctGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCGACGACGTAGCGAAAACGCTACATGTGGTGCCGGGACTGGCGAGGGGGCAATGGCCGTCggctgt > 3:37154/1‑131 (MQ=255) ctGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCGAAGACGTAGCGAAAAAGATAAATTTTGTGCCTGGACTGGCGCGGTTCCAATGGCCTTAGGTc > 2:248199/1‑133 (MQ=255) gAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGGGTCGGCGGAGCGCAAACGATACATGTGGTGGCGGGGCTGGCGGGGTGGCAATGGCCGTCGGTGCTCTCACGTTGAGCAgatg > 4:251963/1‑134 (MQ=255) aCCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTAGCTTCGACGGAGGGAAAACGATAAATGTGGTGGCCGGGCTGGCGGGGGGGCAAATGCCTTCGGTGATCTAACTTTTAGCAGCTTATGCAGCGa > 3:319675/1‑134 (MQ=255) | ACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGA > ACB122/1123113‑1123331 |
| 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: ‑ |