Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,213,521 | 0 | T | G | 55.6% | 30.3 / 37.1 | 27 | K12N (AAA→AAC) | iclR | glyoxylate bypass operon transcriptional repressor IclR |
Reads supporting (aligned to +/- strand): ref base T (0/12); new base G (15/0); total (15/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.75e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.62e-06 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGT > NZ_CP009273/4213448‑4213604 | tCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTcg < 2:246078/90‑1 (MQ=255) aTTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTg < 1:78778/90‑1 (MQ=255) ttCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTTCCGCGGTTTGCGGGAATggggggg > 1:87120/1‑86 (MQ=255) gCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGTTTCGCGGGGATGGGGGCGGCCAt > 1:852998/1‑90 (MQ=255) cACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCCCGTTTCGCGGGAATGGGGGCGAccccga > 1:253436/1‑86 (MQ=255) cgcgTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTTCCGCGTTTTGCGGGGATGTGGGCGACCCTGTCa > 1:177953/1‑90 (MQ=255) cgTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGt < 1:317944/90‑1 (MQ=255) cgTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGt < 1:95822/90‑1 (MQ=255) ttAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGtct < 2:542126/90‑1 (MQ=255) aaGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCtt < 1:579643/90‑1 (MQ=255) aCTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGTTTTGCGGGGATGGGGGGGGCCACGACAGGCTCttttt > 1:538011/1‑90 (MQ=255) tGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGGTTTGCGGGGATGGGGGGGGCCATGACAGTCTCttttttt > 2:634382/1‑89 (MQ=255) tGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGGTTTGCGGGGAAGGGGGGGGCCACGACAGGCTCttttttt > 2:91955/1‑89 (MQ=255) tGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGGTTCGCGGGGAAGGGGGGGGCCCTGACAGTCTCttttttt > 2:833223/1‑89 (MQ=255) gAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTTCCGCGTTTTGCGGGGGTGGGGGGGGCCAAGACAGTCTCtttttttt > 2:138852/1‑88 (MQ=255) gAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTTCCGCGGTTTGCGGGGAAGGGGGGGGCCCTGAAAGGCTCtttttttt > 1:577197/1‑88 (MQ=255) gTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGTTTTGCGGGGATGGGGGGGGCCACGACAGACTCTTTTTTTTTTAttg > 1:43188/1‑88 (MQ=255) gTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGTTTCGCGGGAATGGGGGGGGCCATGACAGGCTCTTTTTTTTGTATCg > 1:854825/1‑90 (MQ=255) aGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGTTTTGCGGGGATGGGGGGGACCCTGAGAGGCTCTTTTTTTTTTATCGTggg > 1:688341/1‑89 (MQ=255) gTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGaa < 1:91955/90‑1 (MQ=255) gCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAAtcattttc < 2:43188/90‑1 (MQ=255) gCGGTGGCAACGGCGGGGTTTTTTCCCCGGTTCCCGGGGAAGGGGGGGGCCATGGCCGGCTCTTTTTTTTGTTTTGTGGGAAtacttttt > 2:257371/1‑89 (MQ=255) ggCAACGGCGGGGTTTTTGCCGCGTTTCGCGGGAATGGGGGCGACCATGACAGTCTCTTTTTTTTGTATCGTGGGAAAcatttttatttt > 1:749135/1‑90 (MQ=255) gCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCAttttt < 2:177953/90‑1 (MQ=255) gCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCAttttt < 2:570494/90‑1 (MQ=255) ggcggGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGt < 1:404147/90‑1 (MQ=255) ggcggGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGt < 2:749135/90‑1 (MQ=255) | TCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGT > NZ_CP009273/4213448‑4213604 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |