Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 807,177 | 0 | T | G | 59.1% | 22.5 / 26.2 | 22 | V67G (GTG→GGG) | bioC | malonyl‑ACP O‑methyltransferase BioC |
Reads supporting (aligned to +/- strand): ref base T (2/7); new base G (13/0); total (15/7) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.11e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.40e-05 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCA‑GAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGA > NZ_CP009273/807091‑807261 | tCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGAc < 1:288476/90‑1 (MQ=255) ccACAGCGTAAATACCCCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACg < 1:339969/90‑1 (MQ=255) cccccAGGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGAtctctc < 1:75571/88‑1 (MQ=255) cACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGGCACGCGCGGGGGGCGGCCCTAGATCTCTccccg > 2:29286/1‑90 (MQ=255) aCGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGGACGTCACGCGCAGGGGACGGCCCTAGATCTCTcgccgc > 1:339395/1‑90 (MQ=255) aCGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGGACGGCAGGCGCCGGGGGCGGCCTTTGATCTCTcccccc > 2:353368/1‑88 (MQ=255) aCGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTcgccgc > 1:143475/1‑90 (MQ=255) aCGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGGGACGGCCTTAGATCTCTcgccgc > 1:3682/1‑90 (MQ=255) aCGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGGGACGGCCCTAGATCTCTcgccgc > 1:606952/1‑90 (MQ=255) aCGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAAAGTCACGCGCAGGGGACGGCCCTAAAACTTTcgcccc > 1:289254/1‑88 (MQ=255) ggACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGGACGTCACGCGCCGGGGACGGCCTTTGATCTCTCGCCGCCCATGCt > 2:187005/1‑90 (MQ=255) ggACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGGGGCGGCCTTAGATCTCTCGCCGCCCATGCt > 1:160361/1‑90 (MQ=255) ctggctggATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGAc > 1:161034/1‑46 (MQ=255) ctggctggATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGAc < 2:161034/46‑1 (MQ=255) tggctggATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGGGACGGCCTTAGATCTCTCGCCGCCAAAGCTTGTTCAGGCACGCCA‑Gaa > 1:412400/1‑90 (MQ=255) tggctggATGAGCCGCCACTGGCGGGAAAGTCACGCGCAGGGGACGGCCTTTGATCTCTCCCCGCCCATGCTTGTTTAGGCACGCCA‑Gaa > 2:15165/1‑90 (MQ=255) ctggATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCA‑GAAGGa < 2:143475/90‑1 (MQ=255) ggATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGGGACGGCCCTAGATCTCTCCCCGCCCATGCTTGTTCAGGCACGCCC‑GAAGGATg > 2:81973/1‑90 (MQ=255) gAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCA‑GAAGGATGCCGc < 1:81973/90‑1 (MQ=255) aGCCGCCACTGGCGGGGAAGTCACGCGCAGGGGACCGGCCTAGATCTCTCCCCCCCCAAGCTTTTTCACGCACCCCACGAAGG‑TGCCGcc > 2:310378/1‑89 (MQ=255) gAACGTCACGCGCAGGGGACGGCCCTAGATCTCTCGCCGCCCAAGCTTGTTCAGGCGCGCCA‑GAAGGAGGCCGCAGACCATTATCTgggg > 2:527548/1‑88 (MQ=255) gcAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCA‑GAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGa < 2:160361/90‑1 (MQ=255) | TCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCA‑GAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGA > NZ_CP009273/807091‑807261 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |