Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,055,679 | 0 | A | C | 55.2% | 24.9 / 38.9 | 29 | T491P (ACC→CCC) | scpA | methylmalonyl‑CoA mutase |
Reads supporting (aligned to +/- strand): ref base A (11/2); new base C (0/16); total (11/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.25e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.92e-06 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACAGCATAACGAAAACCTGCTGGCTGCCGCTGTTAAT > NZ_CP009273/3055591‑3055768 | ggATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCAc > 1:42185/1‑90 (MQ=255) gATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCAcc > 2:465647/1‑90 (MQ=255) aTCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCAccc > 2:67138/1‑90 (MQ=255) aTCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCAccc < 1:512504/90‑1 (MQ=255) aaaCCGATGTACTTGAGATCGACAACGGGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCCCCCGTGATGATgccg < 2:376893/90‑1 (MQ=255) aCTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATGCCGCCGTAAccgc < 2:330044/90‑1 (MQ=255) aCTTGAGATCGACAACGTGATGGTGCGTAACGAGAAAATTGCTTCGCTGGAACGCATTCGCGCCCCCCGTGATGATGCCGCCGTAAccgc < 2:763600/90‑1 (MQ=255) tCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGa > 1:916438/1‑90 (MQ=255) ccGACAACGGGAGGGTGCGAAACGAGAAAATTGCTTCGCTGGAACGCTTTCGCCCCCCCCGTGATGATGCCGCCGTAACCGCCGCGTTGa < 1:231984/89‑1 (MQ=255) gACACCGGGAGGGTGCGAAACGCGAAATTTCCTCCGCGGGACCGCCTCCCCCCCCCCGGTGATGATGCCGCCGTAACCGCCGCGTTGAAc < 1:606494/90‑1 (MQ=255) gACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGAAc > 1:113129/1‑90 (MQ=255) gACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGAAc > 1:861732/1‑90 (MQ=255) aCGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGccc > 2:356613/1‑90 (MQ=255) agggTGCGAAACGAGAAAATTGCTTCGCGGGAACGCATCCGCCCCCCCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACt < 1:290834/88‑1 (MQ=255) gCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCAcgc > 1:504294/1‑90 (MQ=255) gCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCAcgc > 1:631520/1‑90 (MQ=255) gCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCAcgc > 1:545177/1‑90 (MQ=255) tAACGAGAAAATTCCTCCGCGGGACCGCATCCGCCCCCCCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCAcgccgc < 1:262368/90‑1 (MQ=255) accaGCAAATTTCTTCGCCGGACCCCATTCGCGCCCCCCGCGATGATCCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGcac < 2:505944/87‑1 (MQ=255) aCGAGCAATTTGCTCCGCGGGACCGCTTCCGCCCCCCCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGcac < 1:895162/90‑1 (MQ=255) aCGAGAAATTTGCTCCGCTGGAACGCATCCGCGCCCCCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGcac < 2:706514/90‑1 (MQ=255) aCGAGAAAATTGCTTCGCGGGACCCCATTCGCCCCCCCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGcac < 2:870995/90‑1 (MQ=255) aCGAGAAAATTGCTCCGCTGGACCGCTTCCGCCCCCCCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGcac < 1:843723/90‑1 (MQ=255) aCGAGAAAATTGCTCCGCGGGAACGCATCCGCCCCCCCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGcac < 1:330730/90‑1 (MQ=255) gAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACAg > 1:702484/1‑90 (MQ=255) aaaaTTGCTTCGCGGGACCGCATTCGCCCCCCCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACAGCAt < 2:653024/89‑1 (MQ=255) cgcgGGACCGCTTCCGCCCCCCCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACAGCATAACGAAAAcc < 1:32105/90‑1 (MQ=255) ttCGCCCCCCCCGTGAGGAGGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACAGCATAACGAAAACCTGCTGGCTGCCg < 1:356613/90‑1 (MQ=255) tccgccCCCCCGGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACAGCATAACGAAAACCTGCTGGCTGCCg < 1:250316/88‑1 (MQ=255) aCCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGcac > 1:185720/1‑55 (MQ=255) aCCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGcac < 2:185720/55‑1 (MQ=255) aCCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACAGCATAACGAAAACCTGCTGGCTGCCGCTGTTAAt > 2:881790/1‑90 (MQ=255) | GGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACAGCATAACGAAAACCTGCTGGCTGCCGCTGTTAAT > NZ_CP009273/3055591‑3055768 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |