| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 328,807 | 0 | A | T | 56.0% | 5.2 / 28.0 | 25 | D176E (GAT→GAA) | betI | choline‑inducible betIBA‑betT divergent operon transcriptional repressor |
| Reads supporting (aligned to +/- strand): ref base A (5/6); new base T (14/0); total (19/6) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.61e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.22e-04 | |||||||||||
CTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCA‑GGGAAT‑TAGCGC‑GGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGA > NC_000913/328659‑328952 | ctcGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCA‑GGGAAT‑TAGCGC‑GGGTTTTATc < 1:441316/151‑1 (MQ=255) gTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGCC‑GGGGAT‑TAGCGC‑GGGGGTTTTCCAGCGGGTTGCCGCTCCCAGCcccgcg > 1:295890/1‑150 (MQ=255) taACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCA‑GGGAAT‑TAGCGC‑GGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGc < 1:26259/151‑1 (MQ=255) accaTGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCA‑GGGGAA‑TTGCGC‑GGGGGTTTTTCAGCGGTGTTGCGCCCCGCGCCGCGCgccgccccca > 1:417656/1‑147 (MQ=255) accaTGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCA‑GGGAAT‑TAGCGC‑GGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATaa < 1:34920/151‑1 (MQ=255) ccaTGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGCAGGGGAATATAGCG‑‑GGGTTTTTTCCACCGGTTTTCC‑CCCCCAGCCGCGCgccgcccaaaa > 2:366395/1‑146 (MQ=255) aTGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCA‑GGGAAT‑TAGCGC‑GGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATcc < 2:24938/151‑1 (MQ=255) tGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGGC‑GGGGAT‑TAGCGG‑GGGGTTTTTTTCCCGGGTTGCCGCCCCCAGCCccgcgcgggcccaaatc > 1:183770/1‑139 (MQ=255) atataAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCG‑GGGGAT‑TAGTGG‑GGGGTTTTTTCAACGGGTTGTCGCCCCGAGCCGCGCGCAGCCCTAATccacca > 1:245014/1‑151 (MQ=255) atataAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCA‑GGGGAT‑TAGCGC‑GGGGTTTTTCCAGCGGGTTGCCGCTCAGAGCCGCGCGCCGCCCTAATCCatca > 2:156961/1‑151 (MQ=255) atataAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGGA‑GGGGAA‑TTGGGC‑GGGGGTTTTTCCCCGGGTTGCCGCTCAGAGGCGCGCGCCGCCCAAAACCatc > 2:51444/1‑150 (MQ=255) ctgctgTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCA‑GGGAAT‑TAGCGC‑GGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAgcgc < 1:133340/151‑1 (MQ=255) gTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCA‑GGGGAT‑TAGGGC‑GGGGTTTTTTCAGCGGGTTTTCGCTCCGAGCCGCGCGGAGCCCTAATACATCAATAAACGCGGcc > 2:139511/1‑150 (MQ=255) gTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCA‑GGGGAA‑TAGCGCGGGGGTTTTTCCA‑CGGGTTGTCGCCCCGAGCCCCGCGCAGGCCCAATCCATCAATCAGCgccgcc > 1:452455/1‑150 (MQ=255) gTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGCA‑GGGGAA‑TAGCGG‑GGGGTTTTTTCAGCGGGTTGCCGCCCCGAGCCCCCCGCAGCCCCAATACATCCAACAGCGCGGcc > 1:578822/1‑150 (MQ=255) cGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGCA‑GGGGAA‑TAGCGG‑GGGGGTTTTTCAGCGGGTTGCCGCCCCGAGCCCCGCGCGGCCCCAAACCATCAAACAACGcggcccggccccggcc > 1:274856/1‑144 (MQ=255) tCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGCC‑GGGGAT‑TAACGC‑GGGGTTTTTTCAACGGGTTTTCGCCCCCAGCCGCGCGCGGCCCTAATTCATCAATCAGCGGGGCCCGGCCGGCGcccccctc > 2:148404/1‑144 (MQ=255) gTGGGTAGATGCTGAGTGATAAAGTGGCGGGGGA‑GGGGAA‑TAGCGG‑GGGGTTTTTTCAGGGGGTTGGCGCCCCGAGACGCCCGCCGCCCCAAAACCACCATCAACGCGGCCCGGCCGGAGGCCCCCCCCCCTGCCCGTTTCCGCGGccaat > 1:117180/1‑147 (MQ=255) gggTAGATGCTGAGTGATAAAGTGGCGGGGCA‑GGGGAT‑TAGCGC‑GGGGTTTTTTCCGCGGGTTTGCGCTCCGAGGCGCGGGCGGGCCTAATTCATCAATCAGAGCGGCCCGGCCGGAGCCCGCCTCCCGTGGCCGTTTGCGCGGGAATTc > 1:517610/1‑150 (MQ=255) ggTAGATGCTGAGTGATAAAGTGGCGGGGCA‑GGGAAT‑TAGAGC‑GGGGGTTATCCAGCGGGTTGCCGCTCCGTGCCGCGGGCGGGCCTAATCCCACCAAAAGGGTGGCCCGGGCCCAGCCCCCCTCCCGTGCCTGGTCGCGCCGCGAATcg > 1:290310/1‑150 (MQ=255) ggTCA‑GGGAAT‑TAGCGC‑GGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGCGACTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCCGCTTCGa > 1:303170/1‑150 (MQ=255) gTCA‑GGGAAT‑TAGCGC‑GGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGAc < 2:100300/150‑1 (MQ=255) ggTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGgctg > 1:509218/1‑148 (MQ=255) ggTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGgctg > 2:210576/1‑148 (MQ=255) gTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGa > 2:166758/1‑151 (MQ=255) | CTCGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCA‑GGGAAT‑TAGCGC‑GGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGA > NC_000913/328659‑328952 |
| 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: ‑ |