| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 324,518 | 0 | A | T | 57.1% | ‑3.2 / 18.4 | 21 | D176E (GAT→GAA) | betI | transcriptional regulator BetI |
| Reads supporting (aligned to +/- strand): ref base A (1/8); new base T (12/0); total (13/8) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.42e-05 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.90e-04 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGC‑GGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCAT‑AATCCATCAATCAGCGC‑GGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGG > NZ_CP009273/324392‑324656 | gacagggTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGC‑GGGTTTTATCCAGCGGTTTg < 1:457276/134‑1 (MQ=255)gCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGC‑GGGTTTTATCCAGCGGTTTg < 1:227408/139‑1 (MQ=255) caTGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCCGGGGAGTTGGGC‑GGGGGTTTTTCAGCGGGTTGGCGCTCcgaggccggc > 2:457276/1‑133 (MQ=255) aTGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGCAGGGGAATAACGC‑GGGGTTTTTTCCCCGGGTTTCCGCCCACAGCcgcgcg > 1:332833/1‑139 (MQ=255) ataAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGCAGGGGAATAGCGC‑GGGGGTTTTTTAGCGGGTTTCCGCCCCGAGCCGCGCgccgcc > 1:417790/1‑138 (MQ=255) atgctgctGTTCGGCCATGCGGGACTTCTCGTCCCCTCGGTTAATCGGTGTGTAGATGCTGAGTGTTAAAGTGGCGGGTCGGGAAATTAGGGC‑GGGTTTTATCCAGCGGTTTGCGGCTCGGAGCCGCGCGCAGCCAT‑a < 2:324824/136‑1 (MQ=255) tgctgTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGGATTAGCGC‑GGGGGTTTTTCAGCGGGTTGCCGCCCCCAGACGCGCGCGGCCCT‑AAACCa > 2:88628/1‑139 (MQ=255) ctgTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGCAGGGGAATAGCGC‑GGGGGTTTTTCAGCGGTTTGCCGCTCCGAGCCGCGCGCAGCCAT‑AATCCat > 1:86231/1‑138 (MQ=255) cGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCCGGGGAATATCGC‑GGGGTTTTTTCACCGGGTTTTCGCCCCCAGCCGCGCGCAGCCCT‑AATCCACCAAACAACGC‑GGcc > 1:615661/1‑139 (MQ=255) cACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGCAGGGGATTAACGCGGGGGTTTTTCCACGGGTTTG‑CGCCCCGAGCCGCGCGCCGCCATAAATCCTCCAAT‑AACGC‑GGCCAgggcc > 1:394496/1‑134 (MQ=255) cACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGCAGGGGAGTAGCGG‑GGGGTTTTTTCAGCCGGTTGGCGCCCCGAGCCGCGCGGCGCCAT‑AAAACATCCATCAAAGC‑GGCCCGGCCggc > 2:688758/1‑137 (MQ=255) ccctcGGTTAATCGGTGGTGAGATGCTGAGGGATAAAGTGGCGGGTGAGGGAATTAGCGC‑GGGTTTTATCAAGCGTTTTGCCGCTCAGGGCCGCGCGCAGCCAT‑AATCCATCAATCAGCGC‑GGCCAGGCCGTAGCCCg < 2:624450/137‑1 (MQ=255) aTCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGGAATAGCGC‑GGGGTTTTTTCAGCGGGTTGTCGCTCCGAGCCGCGCGCCGCCCT‑AATCCATCAATCAGCGC‑GGCCCGGCCGTAGCCCGCCTCCTGTGcc > 1:114041/1‑139 (MQ=255) aTGCTGAGTGATAAAGTGGCGGGGGAGGGGATTAAGGG‑GGGGTTTTTTCAGCGGGTTTGCGCCCACAGGCGCGGCCCGCCAC‑AAACCATCAATAAACCCGGGCCCGGCCG‑AGCCCCCCTCCTGTTTCCGTTTGCGCggg > 2:443990/1‑138 (MQ=255) cTGAGTGATAAAGTGGCGGGGCAGGGGATTAACGC‑GGGGGTTTTCCAACGGGTTTCCGCCCAGAGCCGCGCGCCGCCCT‑AAACCCTCCATCAGCGC‑GGCCCGGCCGCAGCCCGCCTCCTGTGCCCGTTTGCGCgggcat > 1:608364/1‑135 (MQ=255) gTGATAAAGTGGCGGGTCCGGGGATTAGCGC‑GGGGTTTTTTCCGCGGGTTGCCGCCCAGAGGCGCGCGCAGCCCT‑AATCCCTCAATCAGCGC‑GGCCAGGCCGGAGCCCCCCTCCTGTTCCTGTTTCCCCGGCAATTcgc > 1:348351/1‑139 (MQ=255) gTGGCGGGTCAGGGAATTAGCGC‑GGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCAT‑AATCCATCAATCAGCGC‑GGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAAc < 1:411108/139‑1 (MQ=255) gTCAGGGAATTAGCGC‑GGGTTTTATCCAGCGTTTTGCCGCTCAGAGCCGCGCGCAGCCAT‑AATCCATCAATCAGCGC‑GGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCa < 1:772065/139‑1 (MQ=255) ggTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCAT‑AATCCATCAATCAGCGC‑GGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAg > 2:242808/1‑139 (MQ=255) tttATCCAGGGGTTTGCCGCTCAGAGCCGCGCGCAGCCAT‑AATCCATCAATCAGCGC‑GGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAgcg < 2:86231/138‑1 (MQ=255) tATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCAT‑AATCCATCAATCAGCGC‑GGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAgcggc < 2:332833/138‑1 (MQ=255) aTCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCAT‑AATCCATCAATCAGCGC‑GGCCAGGCCGTAACCCGCTTCCTGTGCCTGTTAGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCACGcggcgg > 1:389907/1‑139 (MQ=255) | GCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGAATTAGCGC‑GGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCAT‑AATCCATCAATCAGCGC‑GGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGG > NZ_CP009273/324392‑324656 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |