| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 328,807 | 0 | A | T | 50.0% | 16.6 / 17.8 | 28 | D176E (GAT→GAA) | betI | choline‑inducible betIBA‑betT divergent operon transcriptional repressor |
| Reads supporting (aligned to +/- strand): ref base A (4/10); new base T (14/0); total (18/10) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.53e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
CGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGT‑CAGGGAATTAGCGC‑GGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTG > NC_000913/328661‑328951 | cGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGG‑CGGGGAATGAACAC‑GGGGTTTTTcc > 2:165334/1‑151 (MQ=255) gCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGT‑CAGGGAATTAGCGC‑GGGTTTTATCCAGCGGttt < 1:351046/151‑1 (MQ=255) ccGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGT‑CAGGGAATTAGCGC‑GGGTTTTATCCAGCGGTTTGCCg < 1:213257/151‑1 (MQ=255) tggtggCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGG‑CAGGGAAATAGCGC‑GGGGGTTTTCCCGCGGGTTTTCGCTc > 2:184857/1‑150 (MQ=255) ggtggCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGT‑GTGGGAAGTAGGGG‑GGGTGTTTTTCAGCCGGTTTCCGCTcc > 1:226814/1‑149 (MQ=255) ggAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGT‑GCGGGAAATAAAGC‑GGGGTTTTTTCACCGGTTTTGCGCCCCGAGCcgc > 2:288446/1‑151 (MQ=255) ggAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGG‑GCGGGGATTAACGC‑GGGGGTTTTTCAACGGGTTGCCGCTCCGAGCcgc > 1:93441/1‑151 (MQ=255) accaccATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGG‑CAGGGGAATAGAGG‑GGGGGTTTTTCAGCGGGTTTCCGCCCCCAGCCGCGCgccgcc > 1:624437/1‑150 (MQ=255) caccaTGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGG‑TAGGGAAATAAGGC‑GGGGTTTTTTCCGCGGGTTGGCGCCCCGAGCCGCGCGGAggc > 2:239515/1‑146 (MQ=255) caTGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGG‑CAGGGGAGTTGCGC‑GGGGGTTTTTCAGCGGGTTGTCGCCCCCAGCCGCGCCCCGccc > 2:434861/1‑145 (MQ=255) aTGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGT‑CAGGGAATTAGCGC‑GGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATcc < 1:574564/151‑1 (MQ=255) aTGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGG‑CAGGGAAATAAGGC‑GGGGTTTTTTCCGCGGGTTTTCGCCCCGAGCcgcgccggggcccaa > 2:52999/1‑138 (MQ=255) gTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGT‑CAGGGGAATAGCGCGGGGTTTTTTTCAAGGGTTT‑CCCCCCCGAGCCCCGCGCCGCCCTaa > 1:403015/1‑146 (MQ=255) gTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGG‑GAGGGAAATAGCGG‑GGGGTTTATTTAGCCGTTTGCCGCCCCGACCCgcgcgccggcccaa > 2:440793/1‑138 (MQ=255) gCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGT‑CAGGGAATTAGCGC‑GGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGc < 2:205646/149‑1 (MQ=255) ccATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGT‑CAGGGAATTAGCGC‑GGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGt < 2:350300/151‑1 (MQ=255) ttCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGT‑CAGGGGAATAGGGC‑GGGGGTTTTTCACCGGGTTTCCGCCCCCAGCCGCGCGCAGCCACAAAACATCAAACAACCcggcccggccccaac > 1:401655/1‑146 (MQ=255) cGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGT‑CAGGGAATTAGCGC‑GGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGccc < 1:184862/151‑1 (MQ=255) tcGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGT‑CCGGGGAGTAGCGC‑GGGGTTTTTTCAACGGGTTTTCGCTCCGAGCCCCGCGCGGCCATAATTCATCAAACAGCGCGGGCCGGCCGGAGCCCCCCTCCTTTGCCCGt > 1:390303/1‑149 (MQ=255) cGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGT‑CAGGGAATTAGCGC‑GGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGtt < 1:119710/149‑1 (MQ=255) ttAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGT‑TAGGGGAATAGAGC‑GGGGTTTTTTCAGCGGGTTGTCGCCCACAGCCGCGCGCCGCCCCAAACCCACAAACAACGCCGCCCGGGCGGAGCCCCCCTccccttcc > 1:623230/1‑136 (MQ=255) gTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGTCCGGGGAATTGGGCG‑GGGTTTTTTCCGCGGGTTTGCC‑CCCCCAGGCGCGCGGCGACCAAATCCCTCAATAAGCGCGGCCCCGCCGTCGGCCCCCTCCCGTGCCtttt > 1:283999/1‑139 (MQ=255) tGAGTAATAAAGGGGA‑GGGT‑GAGGGAATTGGCGC‑GGTTTTGATCCGGTGGTTTCCCGCTCAGCGCCGTGGCCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGGTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGaaa < 2:316794/149‑1 (MQ=255) gggT‑CAGGGAATTAGCGC‑GGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCg > 1:226592/1‑150 (MQ=255) cAGGGAATTAGCGC‑GGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGc > 2:503049/1‑151 (MQ=255) gggAATTAGCGC‑GGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCaa < 2:390296/151‑1 (MQ=255) aaTTAGCGC‑GGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAgc < 2:538035/150‑1 (MQ=255) gTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGctgctg > 2:455165/1‑150 (MQ=255) | CGAAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGC‑GGGT‑CAGGGAATTAGCGC‑GGGTTTTATCCAGCGGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAGCGGCGGCTGCTG > NC_000913/328661‑328951 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |