| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 324,502 | 0 | A | G | 52.6% | 10.2 / 25.4 | 19 | S182P (TCC→CCC) | betI | transcriptional regulator BetI |
| Reads supporting (aligned to +/- strand): ref base A (0/9); new base G (10/0); total (10/9) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.08e-05 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.65e-04 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTC > NZ_CP009273/324378‑324616 | tGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGCCAGGGATTAGTGGGGGgtt > 1:355210/1‑139 (MQ=255) caTGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTTGGTAGATGCTGAGTGATAAAGTGGCGGGTTAGGGGATTAGGGGGGGGTTTATACTGCGGGTTGTTCCTCTGAGCcgcgc > 2:450333/1‑139 (MQ=255) caTGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGCAGGGGAATATGGCGGGGTTTTTTCAGCGGGTTGTCGCGCAGAGAcgggg > 1:162547/1‑136 (MQ=255) aTGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGAGGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCcgcgcg < 1:472917/139‑1 (MQ=255) aaaGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATa < 2:355102/139‑1 (MQ=255) aaGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCCGGGGATTAGGGCGGGGTTTTTTCAGCGGGTTTTCGCTCAGAGACGCGCGCCGCCCTaa > 1:416259/1‑139 (MQ=255) aaGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATa < 1:412024/138‑1 (MQ=255) aaGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATa < 1:422258/138‑1 (MQ=255) tgctgTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCCGGGGATTTGCGGGGGGTTTTTTCAGCGGGTTGGCGCCCCGAGCCCCGCGCCGCCCTCATCCa > 2:46026/1‑139 (MQ=255) tgctgTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGGATTAGCGCGGGGTTTTTTCAGCGGGTTTCCGCTCCGAGCCGCGCGCAGCCCCAATCCa > 2:189131/1‑139 (MQ=255) ctgTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCGGGGGGTTAGCGCGGGGGTTTTTCAGCGGGTTGCCGCTCAGAGCCGCGCGCAGCCACAATccct > 2:327995/1‑136 (MQ=255) ctgTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCGGGGGATTAGGGGGGGGTTTTTTCCGAGGGTTTCCGGCCCGGGCCGCGCGCAGCCCCAAACCat > 2:412024/1‑138 (MQ=255) ctgTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGGATTATCGCGGGGTTTTTTCCGCGGGTTTTCCCTCCGAGCCGCGCGCCGCCCTAAttcct > 1:292542/1‑134 (MQ=255) cGGGACATCACGTCTCCTCGGTTAAGCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGcc < 2:322366/139‑1 (MQ=255) gACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGGGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAGCGGTTTCCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAgg < 1:132845/139‑1 (MQ=255) ctcGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCg < 1:189131/136‑1 (MQ=255) ggTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTc < 1:467454/137‑1 (MQ=255) gATGCTGAGTGATAAAGTGGCGGGGCGGGGGAATAGTGCGGGGTTTTTTTCCCGGGTTGTCGCTCCGAGCCGCGCGCCGCCATAATCCATCCATCAGCGCGGCCCGGGCGGGACCCGCCTCCTGTGCCTGTTCGCGCgg > 2:392419/1‑139 (MQ=255) aGTGATAAAGTGGCGGGTCAGGGATTTAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTc < 2:416259/138‑1 (MQ=255) | TGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGCGGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTC > NZ_CP009273/324378‑324616 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |