Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 328,816 | 0 | T | G | 52.0% | 30.1 / 38.0 | 25 | K173N (AAA→AAC) | betI | choline‑inducible betIBA‑betT divergent operon transcriptional repressor |
Reads supporting (aligned to +/- strand): ref base T (7/5); new base G (13/0); total (20/5) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.49e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.51e-04 |
GTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCA‑GGGAAT‑TAGCGC‑GGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACG > NC_000913/328687‑328965 | 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 > 2:210576/1‑148 (MQ=255) ggTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGgctg > 1:509218/1‑148 (MQ=255) gTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGa > 2:166758/1‑151 (MQ=255) gTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACg > 2:605092/1‑151 (MQ=255) | GTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCA‑GGGAAT‑TAGCGC‑GGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTGACCTGCTGTAAACG > NC_000913/328687‑328965 |
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: ‑ |