| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 328,807 | 0 | A | T | 54.5% | 10.1 / 41.9 | 33 | D176E (GAT→GAA) | betI | choline‑inducible betIBA‑betT divergent operon transcriptional repressor |
| Reads supporting (aligned to +/- strand): ref base A (7/8); new base T (18/0); total (25/8) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.63e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.29e-04 | |||||||||||
AAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCA‑GGGAATTAGCGC‑GGGTTTTATCCAGC‑GGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAG > NC_000913/328663‑328939 | aaGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCA‑GGGAATTAGCGC‑GGGTTTTATCCAGc < 2:715104/151‑1 (MQ=255) cGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCA‑GGGGATTAACGC‑GGGGTTTTTTCAGCGGGTTTGcg > 2:133580/1‑149 (MQ=255) gCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCAGGGGAATAGGGGG‑GGGTTTTTTCAACG‑GGTTTGcgc > 2:449841/1‑148 (MQ=255) ggTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCA‑GGGGAATAGCGC‑GGGGTTTTTTCAGC‑GGGTTGCCGCCCAGAGCcgcgcg > 2:271497/1‑151 (MQ=255) ggTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGCA‑GGGGATTAGCGC‑GGGGGTTTTTCAGC‑GGGTTGTCGCCCCGAGAcgcgcg > 1:663515/1‑151 (MQ=255) accaccATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGAGGGTAGATGCTGAGTGATAAAGTGGCGGCGAC‑GGGAAGTATGGC‑GGGGTTTTTTCCGC‑GGGTTGGCGCTCagagaccggcgcgcct > 2:581473/1‑138 (MQ=255) aTGTATATAAAGCTGTTGTTCTGCCATTCGGGACATCACGTCTCCACGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGTTCA‑GTGAATTAGCGC‑GGGTTTTATCCTGC‑GGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATcc < 1:816559/151‑1 (MQ=255) tGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGCC‑GGGAAATAGGGC‑GGGGTTTTTTCAAC‑GGTTTGCCGCCCAGAGCCGCGCGCAGccccaaacc > 2:305777/1‑143 (MQ=255) tGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGCA‑GGGGAATAGAGC‑GGGGTTTTTTTCGC‑GGGTTGTCGCGCCCAGCCGCGGGCAGccccaatac > 1:217228/1‑143 (MQ=255) gTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCC‑GGGGATTAGCGC‑GGGGTTTTTTCTGC‑GGGTTTCCGCCCCGAACCGCGCGCCGCCCCAATccc > 2:490339/1‑149 (MQ=255) gTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCA‑GGGAATTAGCGC‑GGGTTTTATCCAGC‑GGttt < 1:409069/119‑1 (MQ=255) atataAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTTA‑GGGAAATAGCGC‑GGGGTTTTTTCAGC‑GGGTTGCCGCCCCGAGGCGCGCGCAGCCATAATCCatca > 2:818428/1‑151 (MQ=255) atataAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCA‑GGGGAATATCGCGGGGTTTTTTCCA‑C‑GGGTTTTCGCCCCCAGCCGCGCGCCGCCCCAATCCctca > 2:598447/1‑151 (MQ=255) taAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCA‑GGGAATTAGCGC‑GGGTTTTATCCAGC‑GGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCatcaatc < 2:7201/151‑1 (MQ=255) gTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGCA‑GGGGATTAACGG‑GGGGGTTTTTCCGC‑GGGTTTCCGCCCCGAGGCGCGCGCAGCCCTAATACCTCAAAAAACGCGGcc > 1:187987/1‑150 (MQ=255) gCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGCA‑GGGGATTAGCGC‑GGGGTTTTTTCCGC‑GGGTTGCCGCCCCGAGCCGCGCGCCGCCCTAATCCATCAAACAACGCGGCCAGGcc > 2:496640/1‑150 (MQ=255) aTTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCC‑GGGGATTAGCGC‑GGGGGTTTTTCAGC‑GGGGTTTCGCCCCAAGACGCCCGCCGCCCCAATACATCAATCACCGGGGCCccgccgct > 1:29237/1‑148 (MQ=255) cGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGCC‑GGGGAATAAGGC‑GGGGTTTTTTCAGC‑GGGTTTTCGCCCAGAGGCGCGCGCGGCCCAAATCCCTCAATCAACGCGGGCCCGGCGCAGCCCCCCTCCcgtg > 2:138476/1‑151 (MQ=255) ctcGGTTAATCGGTGGGTGGATGCTGAGTGATAAAGTGGGGGGTCA‑GGGAATTAGCGC‑GGGTTTTATCCAGC‑GGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGtt < 1:857654/151‑1 (MQ=255) ctcGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCA‑GGGAATTAGCGC‑GGGTTTTATCCAGC‑GGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGtt < 2:337004/151‑1 (MQ=255) gTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGGCC‑GGGGATTAGCGG‑GGGGTTTTTTCAGC‑GGGTTGGCGCTCCGAGGCGCGCGCCGCCCCAAATCATCAATCAGCGCGGCCCGGCCCGAGCCCCCCTCCCCTGCCCGTTtgcg > 1:126749/1‑151 (MQ=255) ttAATCGGTGGGAGGATGCTGAGTGATAAAGTGGCGGGTCG‑GGGAATTAGCGC‑GGGTTTTATCCAGC‑GGTGTCCGGCGCAGAGCCGCCCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTcgcg < 1:457032/150‑1 (MQ=255) gggTAGATGCTGAGTGATAAAGTGGCGGGTCA‑GGGAAATAGCGC‑GGGGTTTTTCCAGC‑CGGTTGCCCCCCAGAGCCCCCCGCACCCATAATCCCTCCAACACCGCGGCCCGGCCCGAGCCCCCCTCCTGTGCCCGtttccccggccaatt > 2:743201/1‑137 (MQ=255) gggTAGATGCTGAGTGATAAAGTGGCGGGGCC‑GGGAATTAGCGC‑GGGGGTTATCCAGC‑GGGTTGTCGCTCAGAGCCGCGCGCAGCCACAATTCATCCAACAACGCGGGCCGGGCGCAGCCCCCCTCCTCTGCCCGGTTGCGCGGGAAATcg > 1:41195/1‑151 (MQ=255) ggTAGATGCTGAGTGATAAAGTGGCGGGGCA‑GGGGAATAGCGC‑GGGGTTTTTTCCGC‑GGGTTTTCGCTCCGAGCCGCGCGCCGCCCTAATCCATCAATCAGCGCGGCCCCGCCCTAGCCCCCCTCCTGTGCCCGTTTGCGCGGcgaattg > 2:275078/1‑144 (MQ=255) aTGCTGAGTGATAAAGTGGCGGGGCA‑GGGGAATAGGGC‑GGGGTTTTTCCAGA‑GGTTTGCCGCTCAGAGGCGCGCGCAGCCCTAAACCATCAATCAGCGCGGCCAGGGCGGAGCCCGCTTCCCGTGCCCCGTTGCGCGGGAAATCGcgacga > 2:189281/1‑151 (MQ=255) cTGAGGGAGAAGGAGGCGGCACA‑GGGAATTAGCGC‑GGGTTTTAACCAGG‑GGTTCGCCGCGCAGCGCCGAGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAAc < 1:268807/151‑1 (MQ=255) gTGGCGGGTCA‑GGGAATTAGCGC‑GGGTTTTATCCAGC‑GGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGa > 2:359275/1‑151 (MQ=255) gggAATTAGCGC‑GGGTTTTATCCAGC‑GGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATTAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCaa > 2:308249/1‑151 (MQ=255) gAATTAGCGC‑GGGTTTTATCCAGC‑GGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAg > 2:45220/1‑150 (MQ=255) gAATTAGCGC‑GGGTTTTATCCAGC‑GGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAg > 2:177623/1‑150 (MQ=255) gAATTAGCGC‑GGGTTTTATCCAGC‑GGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAg > 1:576834/1‑150 (MQ=255) gAATTAGCGC‑GGGTTTTATCCAGC‑GGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAg > 2:842532/1‑150 (MQ=255) | AAGGTGCGACCGCTGGTGGCGGAGGTATAACCACCATGTATATAAAGCTGCTGTTCTGCCATTCGGGACATCACGTCTCCTCGGTTAATCGGTGGGTAGATGCTGAGTGATAAAGTGGCGGGTCA‑GGGAATTAGCGC‑GGGTTTTATCCAGC‑GGTTTGCCGCTCAGAGCCGCGCGCAGCCATAATCCATCAATCAGCGCGGCCAGGCCGTAGCCCGCTTCCTGTGCCTGTTCGCGCGGCAATTCGCGACGAAACTCGCTCACCAGATTCGACAGCAAG > NC_000913/328663‑328939 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |