Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,213,521 | 0 | T | G | 56.2% | 16.7 / 20.4 | 16 | K12N (AAA→AAC) | iclR | glyoxylate bypass operon transcriptional repressor IclR |
Reads supporting (aligned to +/- strand): ref base T (0/7); new base G (9/0); total (9/7) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 8.74e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.80e-04 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTT > NZ_CP009273/4213434‑4213599 | aTTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGgtt > 2:256639/1‑90 (MQ=255) tGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGttttt > 2:5476/1‑89 (MQ=255) ggATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGtttttt > 1:283929/1‑88 (MQ=255) cGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTTCcgg > 1:475330/1‑89 (MQ=255) cACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTcgcg < 1:5476/90‑1 (MQ=255) cTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGGTTTGCggg > 2:253437/1‑90 (MQ=255) aTTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTg < 2:283929/90‑1 (MQ=255) aTTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTg < 2:128790/90‑1 (MQ=255) aGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGAcc < 1:277297/90‑1 (MQ=255) cACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGGTTCGCGGGGATGGGGGCGGCCATGa > 1:4779/1‑90 (MQ=255) aCGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGGTTTGCGGGGATGGGGGGGGCCATGaa > 1:416420/1‑89 (MQ=255) tAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGTTTCGCGGGAATGGGGGCGACCATGACAGtctc > 1:47148/1‑90 (MQ=255) cTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCCCGGTTTGCGGGGATGGGGGGGGCCCTGACAGGCTCtttttt > 2:127341/1‑90 (MQ=255) tGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTc < 2:75956/90‑1 (MQ=255) tGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTc < 1:270443/90‑1 (MQ=255) gTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGGTTCGCGGGAATGGGGGCGACCATGACAGGCTCTTTTTTTTTTATCg > 1:75956/1‑90 (MQ=255) aGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTGCCGCGTTTCGCGGGAATGGGGGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGa > 1:150851/1‑90 (MQ=255) aGTCGCTGGTGCGGTGGCAACGGCGGGGGTTGTTCCGCGTTTTCTCGGATTGGGGGCGGCCGTGACCGACTCTTTTTTTTTTTTCGttgg > 2:612534/1‑87 (MQ=255) gCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCAttttt < 1:392293/90‑1 (MQ=255) | ATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTT > NZ_CP009273/4213434‑4213599 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |