Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,594,167 | 0 | C | T | 57.9% | 3.5 / 19.9 | 19 | S64F (TCT→TTT) | bglJ | DNA‑binding transcriptional activator BglJ |
Reads supporting (aligned to +/- strand): ref base C (0/8); new base T (11/0); total (11/8) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.32e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-05 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGG > NZ_CP009273/4594089‑4594240 | tCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGaa < 1:146179/90‑1 (MQ=255) tataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTCCATGGGAAGGGAgcgc > 2:368515/1‑90 (MQ=255) tataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTTCCCTGAGAAGGGAgcgc > 2:53548/1‑90 (MQ=255) taCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTTCCATGAGAAGGGAGCGCag > 2:202559/1‑90 (MQ=255) aCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGGGAGCGCaga > 2:29616/1‑90 (MQ=255) aCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCaga < 1:53548/90‑1 (MQ=255) cAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCagaga < 1:20936/90‑1 (MQ=255) aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTTAGAAGTGAGCGCagagag > 2:95684/1‑90 (MQ=255) aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGTGAGCGCagagag > 2:108363/1‑90 (MQ=255) aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCagggag > 2:99707/1‑90 (MQ=255) aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTTCCCTGAGAAGTGAGCGCagaggg > 2:244869/1‑88 (MQ=255) gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGGGAGCGCAGAGAgg > 1:345669/1‑90 (MQ=255) cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGa < 1:95684/90‑1 (MQ=255) tCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTGCCATGAGAAGGGAGCGCAGAGAGGGGTTATTTTGCCTGACTGAAATGGCg > 1:248693/1‑90 (MQ=255) gcgGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGGGAGCGCAGAGAGGGGTTATCTTGCCTGACTGAACTGGCGAt > 2:383384/1‑90 (MQ=255) ctcGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGtt < 1:108363/90‑1 (MQ=255) gcggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTAcc < 1:244869/90‑1 (MQ=255) cggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTAccc < 1:368515/90‑1 (MQ=255) gTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCgg < 1:202559/90‑1 (MQ=255) | TCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGG > NZ_CP009273/4594089‑4594240 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |