Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,594,161 | 0 | C | T | 60.0% | ‑5.4 / 12.8 | 20 | S62F (TCT→TTT) | bglJ | DNA‑binding transcriptional activator BglJ |
Reads supporting (aligned to +/- strand): ref base C (4/4); new base T (12/0); total (16/4) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.44e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.35e-02 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGC > NZ_CP009273/4594088‑4594232 | cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATgaga > 1:57344/1‑90 (MQ=255) cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCCTgaga > 1:150722/1‑90 (MQ=255) acagtACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGaa < 1:42759/86‑1 (MQ=255) ccATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCCTGAGAAg > 2:192466/1‑90 (MQ=255) aTACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTg > 1:108648/1‑90 (MQ=255) tataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGGGAgcgc > 2:123388/1‑90 (MQ=255) tataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGGAGTGAgcgc > 1:34120/1‑90 (MQ=255) tataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAgcgc < 1:176529/90‑1 (MQ=255) tataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGCCATTTTTTTTTTTTCTGCCATGAGAAGTGAgcgc > 2:16452/1‑90 (MQ=255) ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTTCCCTGAGAAGGGAGGGCa > 2:264779/1‑90 (MQ=255) aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCTTTTTTTTTTTTTTTGCCATGAGAAGGGAGGGGagagag > 1:364426/1‑90 (MQ=255) aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGTGGGCGCagctgt > 2:42759/1‑86 (MQ=255) aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGGGAGCGCagagag > 1:351228/1‑90 (MQ=255) gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGGGAGGGCAGAgggg > 2:360625/1‑90 (MQ=255) cAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTTCTTTTTTTTTTTTCTCGCCTGAGGAGTTTGCGGAGAGAGGGGTTATCTTGccct > 1:251323/1‑88 (MQ=255) cATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACt < 2:108648/90‑1 (MQ=255) gcgGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGGGAAGTGGGGGCAGAGGGGGGTTATTTTTCCTGACTGAACTGGCggt > 2:172452/1‑88 (MQ=255) ttGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGGGAGCGCCGAGAGGGGTTATCTTGCCTGACTGAACTGGCGATTAAGttttccc > 1:39603/1‑86 (MQ=255) gcggcgGTCATTTTTTCTTTTTCTGCCATGAGAAGTGa < 2:138040/38‑1 (MQ=255) gcggcgGTCATTTTTTCTTTTTCTGCCATGAGAAGTGa > 1:138040/1‑38 (MQ=255) | CTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGC > NZ_CP009273/4594088‑4594232 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |