Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,594,161 | 0 | C | T | 56.5% | 3.0 / 26.2 | 23 | S62F (TCT→TTT) | bglJ | DNA‑binding transcriptional activator BglJ |
Reads supporting (aligned to +/- strand): ref base C (1/9); new base T (13/0); total (14/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.22e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.06e-04 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGA‑TCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCC > NZ_CP009273/4594087‑4594238 | gCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGA‑TCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATgag > 1:203221/1‑90 (MQ=255) tataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGA‑TCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCCTGAGAAGTGGgcgc > 1:326019/1‑90 (MQ=255) tataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGA‑TCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAgcgc < 1:208790/90‑1 (MQ=255) tataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGA‑TCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAgcgc < 2:121042/90‑1 (MQ=255) ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATTCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTTCCCTTAGAAGTGAgcgc > 2:333313/1‑90 (MQ=255) ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGA‑TCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCCTGAGAAGTGAGCGCa > 1:95772/1‑90 (MQ=255) ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGA‑TCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTTCCCTGTGAAGAGAGCGCa > 2:251441/1‑90 (MQ=255) aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGA‑TCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGGGAGCGCagagag > 1:318006/1‑90 (MQ=255) aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGA‑TCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCagaggg > 1:201678/1‑88 (MQ=255) aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGA‑TCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGTGAGCGCagagag > 2:180463/1‑90 (MQ=255) aGTCAGGAGTCATTTCAGGATGCCATGTCGCGGA‑TCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTGCCATGAGAAGTGAGCGCagagag > 1:263753/1‑90 (MQ=255) gTCAGGAGTCATTTCAGGATGCCATGTCGCGGA‑TCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGGGAGGGCAGAGAgg > 1:328234/1‑90 (MQ=255) tCAGGAGTCATTTCAGGATGCCATGTCGCGGA‑TCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCCTGAGAAGTGAGGGCAGAGAggg > 2:153200/1‑90 (MQ=255) agggAGTCATTTCAGGATGCCATGTCGCGGA‑TCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGa < 1:364133/88‑1 (MQ=255) tttCAGGATGCCATGTCGCGGA‑TCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGc < 2:35930/90‑1 (MQ=255) ttCAGGATGCCATGTCGCGGA‑TCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGcc < 2:37049/90‑1 (MQ=255) tCAGGATGCCATGTCGCGGA‑TCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCct < 2:73250/90‑1 (MQ=255) tCGCGGA‑TCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTTCCATGAGAAGGGAGCGCAGAGGGGGGTTATCTTGCCTGACTGAACTGGCg > 1:127793/1‑90 (MQ=255) gcgGA‑TCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGCCATGAGAAGGGAGCGCAGAGAGGGGTTATTTTGCCTGACTGAACTGGCGAt > 2:388554/1‑90 (MQ=255) tcGTTTGCGGCGGTCATTTTTTTTTTTTTTTCCATGAGAAGGGAGCGCAGAGAGGGGTTATCTTGCCTGACTGAACTGGCGGTTAAGttt > 1:391323/1‑90 (MQ=255) cggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTAccc < 2:201678/90‑1 (MQ=255) cggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTAccc < 1:333313/90‑1 (MQ=255) cggTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTAccc < 1:251441/90‑1 (MQ=255) | GCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGA‑TCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCC > NZ_CP009273/4594087‑4594238 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |