Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,068,173 | 0 | T | G | 23.1% | 109.6 / 10.0 | 53 | G31G (GGT→GGG) | ymdF | conserved protein YmdF |
Reads supporting (aligned to +/- strand): ref base T (18/22); new base G (12/0); total (31/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.25e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGTTAATTCTATGCAATATGATGTCTATACCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACA‑GCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCG > NC_000913/1068034‑1068312 | cgTTAATTCTATGCAATATGATGTCTATACCCAGACGGAGGTCAGTAATGTCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACA‑GCGGGGGt < 2:1031648/140‑1 (MQ=255) tGCAATATGATGTCTATACCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTGCGGGAATTTTGTAGAAGAGCGCGGAAAAGCAACCATAGCAAGTAAAAAAAGAGGGCAACACAGGCGGGGAAAAtttaacaa > 2:513928/1‑135 (MQ=255) gatgTCTATACCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACA‑GCGGGGGTAATTTCAAGAATGACCCGc > 2:844959/1‑140 (MQ=255) gatgTCTATACCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGGAAAAAAGGAGGACAGCACA‑GCGGGGGGAAATTCAAGAATGACCCGc > 2:1132655/1‑140 (MQ=255) tgTCTATACCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACA‑GCGGGGGGAATTTCAAGAATGACCCGCAg > 2:953121/1‑140 (MQ=255) tgTCTATACCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGATCCCGCCAATTTTTTGAAGACCACGCAAGAGCATCCGAAGGCGGGAAAAAAGGAGGACAGGACA‑AGGGGGGGAATTTTAAGGATGACCCGccg > 2:990821/1‑138 (MQ=255) tCTATACCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACA‑GCGGGGGTAATTTCAAGAATGACCCGCAgcg < 1:228750/140‑1 (MQ=255) tCTATACCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACA‑GCGGGGGTAATTTCAAGAATGACCCGCAgcg < 2:496705/140‑1 (MQ=255) tCTATACCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTACGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACA‑GCGGGGGTAATTTCAAGAATGACCCGCAgcg < 2:1157834/140‑1 (MQ=255) cTATACCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACA‑GCGGGGGTAATTTCAAGAATGACCCGCAgcgc < 2:136316/140‑1 (MQ=255) cTATACCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACA‑GCGGGGGTAATTTCAAGAATGACCCGCAgcgc < 2:730867/140‑1 (MQ=255) ataCCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAACTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACA‑GCGGGGGGAATTTCAAGAATGACCCGCAgcgcgc > 2:190088/1‑140 (MQ=255) aCCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACA‑GCGGGGGGAATTTCAAGAATGACCCGCAGCGCGCAt > 2:607126/1‑140 (MQ=255) ccAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACA‑GCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCt > 2:314019/1‑140 (MQ=255) gACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACA‑GCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGaa < 2:885755/140‑1 (MQ=255) gACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACA‑GCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGaa < 1:587689/140‑1 (MQ=255) gACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACA‑GCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGaa > 1:308228/1‑140 (MQ=255) gACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACA‑GCGGGGGGAATTTCAAGAATGACCCGCAGCGCGCATCTGaa > 1:467134/1‑140 (MQ=255) gACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAAAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACA‑GCGGGGGGAATTTCAAGAATGACCCGCAGCGCGCATCTGaa > 1:687316/1‑140 (MQ=255) gACGGAGGTCAGTAATGGCAAACCATCGAGGCGGGTCCGGCAATTTTGCAGAAGACCGCGAAAGAACATCAGAAGCACGGAAAAAAGGTGGACAGCACA‑GCGGGGGCAATTTTAAGAAAAACCCCCAGCGCGCCTCTGaa > 2:671929/1‑140 (MQ=255) cGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGGAAAAAAAGTGGACAGCACA‑GCGGGGGGAATTTCAAGAATGACCCGCAGCGGGCATCTGAAGc > 2:408863/1‑140 (MQ=255) ggAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACA‑GCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCa < 1:65400/140‑1 (MQ=255) tCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACA‑GCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTaa < 2:1052862/140‑1 (MQ=255) gTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGGAAAAAAGGAGGACAGCACA‑GCGGGGGGAATTTCAAGAATGACCCGCCGCGCGCATCTGAAGCAGGGaaaaa > 1:1041821/1‑140 (MQ=255) gCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAAAGCATCAGAAGGAGGGAAAAAAGGGGGACAGCACA‑GCGGGGGGAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGGAAAAAAggtgg > 1:632303/1‑140 (MQ=255) gCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGAACGCGAAAGAGCATCAGAAGCAGGGAAAAAAGGTGGACAGCACA‑GCGGGGGGAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGGAAAAAAggggg > 1:1151274/1‑140 (MQ=255) aTCGAGGCGGTTCCGGCAATTTTGCAGAAAACCGCGAAAAAGCATCAGAAGCCGGGAAAAAAGGGGGACAGCACA‑GCGGGGGGAATTTCAAAAATGACCCGCAGCGCGCATCTGAAGCAGGGAAAAAAGGGGGGAaaaga > 2:292113/1‑139 (MQ=255) gAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACA‑GCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGt < 2:864742/140‑1 (MQ=255) aGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGGGGACAGCACA‑GCGGGGGTAATTTCAAGAAAGACCCGCAGCGCGCATCTGAAGCAGGGAAAAAAGGTGGGAAGAGCAGTc > 2:955329/1‑140 (MQ=255) cGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACA‑GCGGGGGGAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACg > 1:427209/1‑140 (MQ=255) ttCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACA‑GCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCa < 2:308228/140‑1 (MQ=255) tCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGGGGACAGCACA‑GCGGGGGGAATTTCAAGAATGACCCGCAGCGCGCATCTTAAGGAGGGAAAAAAGGGGGGAAGAGCAGGCACGGCaa > 2:158082/1‑140 (MQ=255) ttttGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAACACA‑GCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACa > 2:577642/1‑140 (MQ=255) cAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACA‑GCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAg < 2:170724/140‑1 (MQ=255) cAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACA‑GCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAg < 1:1132655/140‑1 (MQ=255) cAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACA‑GCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAg < 1:314019/140‑1 (MQ=255) gaagaCCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACA‑GCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGcc < 1:440298/140‑1 (MQ=255) gaCCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACA‑GCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCggg < 1:844959/140‑1 (MQ=255) cgcgAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACA‑GCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTa > 2:647062/1‑140 (MQ=255) cgcgAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACA‑GCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTa > 2:513943/1‑140 (MQ=255) gagCATCAGAAGCAGGTAAAAAAGGTGGACAGCACA‑GCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATg > 2:517558/1‑140 (MQ=255) gCATCAGAAGCAGGTAAAAAAGGTGGACAGCACA‑GCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGAc < 1:190088/140‑1 (MQ=255) aTCAGAAGCAGGTAAAAAAGGTGGACAGCACA‑GCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGa < 2:687316/140‑1 (MQ=255) aaGCAGGTAAAAAAGGTGGACAGCACA‑GCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCa < 1:607126/140‑1 (MQ=255) aaGCAGGTAAAAAAGGTGGACAACACA‑GCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCa < 1:577642/140‑1 (MQ=255) gCAGGTAAAAAAGGTGGACAGCACA‑GCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCAtt > 2:205955/1‑140 (MQ=255) gCAGGTAAAAAAGGTGGACAGCACA‑GCGGGGGAAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGCGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCAtt > 2:436399/1‑140 (MQ=255) aGGTAAAAAAGGTGGACAGCACA‑GCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCAACAACTAGCCGGGCTAATCAATGACGAATGCAtttt > 1:1045894/1‑140 (MQ=255) ggTAAAAAAGGTGGACAGCACA‑GCGGGGGTAATTTCAAGAATGCCCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCAttttt < 1:990821/140‑1 (MQ=255) ggTAAAAAAGGTGGACAGCACA‑GCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCAttttt > 2:849398/1‑140 (MQ=255) gTAAAAAAGGTGGACAGCACA‑GCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTg > 1:17228/1‑140 (MQ=255) gTAAAAAAGGTGGACAGCACA‑GCGGGGGTAATTTCAAGAATGACACGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTg > 1:736834/1‑140 (MQ=255) ggTGGACAGCACA‑GCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGc > 1:269148/1‑140 (MQ=255) tGGACAGCACA‑GCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTc > 2:648315/1‑140 (MQ=255) caca‑gcGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCaaaa > 2:607229/1‑140 (MQ=255) caca‑gcGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCaaaa > 1:131009/1‑140 (MQ=255) gggggTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGAAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTAAAAAAGCCAt < 2:269148/140‑1 (MQ=255) gTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCAcc < 1:75964/140‑1 (MQ=255) gTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCAcc > 2:219479/1‑140 (MQ=255) tAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCAccg < 1:648315/140‑1 (MQ=255) | CGTTAATTCTATGCAATATGATGTCTATACCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACA‑GCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCG > NC_000913/1068034‑1068312 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |