Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,213,525 | 0 | C | T | 57.9% | ‑2.3 / 16.0 | 19 | R11K (AGA→AAA) | iclR | glyoxylate bypass operon transcriptional repressor IclR |
Reads supporting (aligned to +/- strand): ref base C (3/5); new base T (11/0); total (14/5) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.82e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.68e-02 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCT > NZ_CP009273/4213441‑4213609 | tCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTTCcg > 1:160111/1‑90 (MQ=255) cGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCcgc > 1:82515/1‑90 (MQ=255) cTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCCCGTTTTGCggg > 1:176586/1‑90 (MQ=255) ttCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGGTTCGCGGGGATGGGGGCg > 2:560026/1‑90 (MQ=255) cgTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGt < 2:618827/90‑1 (MQ=255) tAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGTTTCGCGGGAATGGGGGGGGCCATGACAGtctc > 2:179446/1‑90 (MQ=255) tGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTc < 1:131968/90‑1 (MQ=255) tGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGTTTTCCGGGGAAGGGGGGGGCCCCGACAGGCTCttttttt > 1:46396/1‑89 (MQ=255) gTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGGGCGACCATGACAGGCTCTTTTTTTTGTATCg > 2:114199/1‑90 (MQ=255) gTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTTCCGCGGTTTGCGGGGGAGGGGGGGGCCCCGAAAGGCTCTTTTTTTTGtgttg > 1:788781/1‑86 (MQ=255) gTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGTTTTGCGGGGATGGGGGCGGCCATGACAGTCTCTTTTTTTTTTATCg > 1:317666/1‑90 (MQ=255) gTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTGCCGCGTTTCGCGGGGATGGGGGCGGCCATGACAGTCTCTTTTTTCTGTTTCGTggga > 2:233639/1‑88 (MQ=255) gTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGGGCGACCATGACAGTCTCTTTTTTTTGTATCGTGGaa > 2:271196/1‑90 (MQ=255) gTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGTTTCGCGGGGATGGGGGCGACCATGACAGTCTCTTTTTTTTTTATCGTggga > 2:22469/1‑88 (MQ=255) tGGTGCGGTGGCAACGGCGGGTTTTTTGCCGCGTTTCGCGGGAATGGGGGCGACCCTGACAGTCTCTTTTTTTTGTATCGTGGGAAacat > 2:233873/1‑90 (MQ=255) tGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGGATGGGGGGGACCATGACAGGCTCTTTTTTTTGTATCGTGGAAAacat > 2:426472/1‑90 (MQ=255) ggCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATcattttcatttt < 2:259792/90‑1 (MQ=255) gCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCAttttt < 2:309177/90‑1 (MQ=255) gTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCt < 1:105856/90‑1 (MQ=255) | TCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCT > NZ_CP009273/4213441‑4213609 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |