Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,001,918 | 0 | A | C | 41.1% | 88.8 / 12.4 | 57 | intergenic (‑129/+192) | fliA/fliC | RNA polymerase, sigma 28 (sigma F) factor/flagellar filament structural protein (flagellin) |
Reads supporting (aligned to +/- strand): ref base A (26/7); new base C (0/23); total (26/31) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.85e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAAT‑CGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGG > NC_000913/2001770‑2002065 | gCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAAT‑CGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAAcc > 2:743952/1‑151 (MQ=255) taGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAAT‑CGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCccgccg > 2:744944/1‑151 (MQ=255) taGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAAT‑CGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCccgcc < 1:918287/150‑1 (MQ=255) atttCACGAAAAACAGCCCGGCTTTATAGGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAAT‑CGGACGATTAGTGGGTAAAATGAGGGTTTATTTGGGGGTTACAGGAAAATTCCAGCCAGAAAAAAACCCCGCCGGTg < 1:537108/145‑1 (MQ=255) aTAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCGGGTTTTTTAAT‑CGGACGATTAGTGGGTGAAAGGGGGGGTTATTTGGGGGTTACAGGTAAATCCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGaa < 2:960706/147‑1 (MQ=255) tAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAAT‑CGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACg > 1:877614/1‑151 (MQ=255) accAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCTTTTCTGCAGGGTTTTTAAT‑CGGACGATTAGTGGGTGAAATGGGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGt < 1:966119/148‑1 (MQ=255) aaCAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAAT‑CGGACGATTAGGGGGTGAAAGGAGGGTTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGtt < 2:820153/151‑1 (MQ=255) aGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGGTTTTAAT‑CGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGGTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCt > 2:414863/1‑151 (MQ=255) gCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAAT‑CGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTg > 2:721375/1‑151 (MQ=255) cccTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAAT‑CGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTg > 1:828276/1‑150 (MQ=255) ccTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAAT‑CGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGAc > 2:852255/1‑151 (MQ=255) cTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAAT‑CGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGAc > 1:185117/1‑150 (MQ=255) gggtttATGTGAGTTAGGGGCAGGATTATCCTTTTCGGGAGGTTGTTGATG‑GGGACGATTGGGGGGAGAAGGGGGGGTTTTTTGGGGGGTAAGGGGAAAATTCGAGGAAGAAAAAACCGCCGCCGGTGGCGGGGAAGCACGTTGCTGACaa < 1:597509/147‑1 (MQ=255) gttttATATGGTTTATCGGAATGATTACCCTTTTCTGCGGGTTTTTAAAT‑GGGACGATTAGTGGGTAAAATGAGGGTTTATTGGGGGGTTACGGGTAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACa < 2:692143/146‑1 (MQ=255) gCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAAT‑CGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACaa > 1:544135/1‑150 (MQ=255) gTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAAT‑CGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAAtt > 2:805081/1‑151 (MQ=255) gTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAAT‑CGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGGTACAGGTAAATTTCAGGCAGAAAAAAACCCCGCCGGGGGCGGGGACGCACGTTGCGGACAAAtt > 1:986251/1‑151 (MQ=255) tatGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAAT‑CGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgcgc > 2:908723/1‑151 (MQ=255) tGATTTATCGGCATGATTATCCTTTTCTGCGGGGTTTTTAAT‑CGGACGTTTAGTGGGGGAAATGAGGGTTTATTTGGGGGTTACAGGAAAATCCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgcgc < 2:820721/149‑1 (MQ=255) tgaTCGGAGTGTTTTTCCTTTTCTGCGGGTTTTTAAAT‑CGACCGATTAGTGGTTTATATGGGGGTTTATTTGGGGGTTAAGGGAAAGTCCCAGGCAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCtt < 1:443948/145‑1 (MQ=255) ggCATGATTATCCTTTTCGGC‑GGGTTTTTAATCCGGACGATGGGGGGGTGAAGTGAGGGTTTATTTGGGGTTTAAAGGTAAATTCCAGGAAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGtt < 1:1000463/148‑1 (MQ=255) ggCATGATTATCCGTTTCGGCAGGTTTTTTAAT‑CGGACGTTAAGGGGGTAAAAGGGGGGGTGTTTTGGGGTTTAGAGGTAAATCCCAGGGAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGt < 1:705179/147‑1 (MQ=255) cATGATTATCCGTTTCTGCAGGGTTTTTAAT‑CGGACGACTAGTGGGTGAAATGAGGGGTTATTTGGGGGATATAGGTAAATTACAGGGAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGAAAAATTGCGCTTTATGTTGCCg > 1:243858/1‑150 (MQ=255) tgccggtgtttttCGGCGGGTTTTTTAAT‑CGGACGATTGTTGGGTAAAATGAGGTGTTTTTTGGGGGTAAGGGGTAAATCCGAATCAGAAAAACACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCg < 1:669426/138‑1 (MQ=255) atttcccGTTTTTGCAGGTTTTTAAAG‑CGGACATTTAGTGGGTAAAATGGGGGGTTGTTGGGGGTTTACGGGAAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGc < 2:539973/139‑1 (MQ=255) gTTTCTGCAGGGTTTTTAAT‑CGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCg > 2:669774/1‑139 (MQ=255) tttcgtcAGGGTTTGAACA‑CGGACGATTAGGGGGTAAAATGAGGGGTTATTTGGGGGTTAGGGGTAATTGCGGGGCAGAAAAACCCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCg < 1:1013994/132‑1 (MQ=255) tgtCTGGAGGTTTTTGAAT‑CGGACGATTAGTGGGTGACATGAGGAGTTATTTGGGGGTTACGGTTAACTTCCCGGAAAAAAAAATCCCCGCCGGTGTCGGGGAAGCACGTTGCTGACAAAGTGCGCTTTATGTTGCCGGATGCGGCGt < 1:32214/146‑1 (MQ=255) tCTGCAGGGTTTTTAAT‑CGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACg > 2:974693/1‑151 (MQ=255) tttAAT‑CGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGc > 2:279135/1‑151 (MQ=255) aT‑CGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAc > 1:472043/1‑151 (MQ=255) ggACGATTGGTGGGTGAAATGGGGGGTTATTGGGGGTTTACAGGTAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaa < 2:1106388/150‑1 (MQ=38) ggACGATTAGGGGGTAAAATGGGGGTTTTTTGGGGGGTTACAGGTAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACa < 2:301058/149‑1 (MQ=38) gACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaa > 2:462432/1‑150 (MQ=255) aaagATTAGTGGGTGAAATGAGGGTTTTTTTGGGGTTTACAGGTAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaa < 1:897544/148‑1 (MQ=38) aCGATTAGTGGGGGAAAGGAGGTGTTTTTTGGGTGTAACAGGTAAATCCCAGACAGAAAAACACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaa < 1:128319/150‑1 (MQ=37) aTTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGGAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaaa > 1:657768/1‑148 (MQ=255) ggggggaaaaGGAGGGGTTATTGGGGGATAAAAGGTAAATCCCGGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCaa < 2:838126/144‑1 (MQ=25) gggggTGAAAGGGGGGTTTTTTTGGGGGTTACAGGTAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGc < 1:479590/147‑1 (MQ=37) gggaaaaaTGGGGGTTTTTTGGGGGGTTACAGGTAAATCCGAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAAt < 2:857187/146‑1 (MQ=25) aaaTGGGGGTTTATTTGGGGTTTAAGGGTAATTTCTAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCa < 1:372578/149‑1 (MQ=37) tGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATaa > 2:837417/1‑150 (MQ=255) gggTTTTTTGGGGGTTACAGGTAATTTCCAGGAAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGc < 1:979010/151‑1 (MQ=255) tATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAtt > 2:736657/1‑151 (MQ=255) aTTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTc > 2:963900/1‑151 (MQ=255) tttGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTc > 2:678134/1‑150 (MQ=255) gggggTTACAGTTAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACt < 2:8681/151‑1 (MQ=255) gggggTTACAGGTAAATTCCAGGAAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACt < 2:277239/151‑1 (MQ=255) ggggTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCACATTCAATAAATTGCAATTccattt > 1:766808/1‑146 (MQ=255) aGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCt > 2:764769/1‑151 (MQ=255) gttAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAgg < 2:401577/145‑1 (MQ=255) ggaaaaTCCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCt < 2:416674/147‑1 (MQ=255) ggTAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCt < 2:644055/150‑1 (MQ=255) aaaTTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATa < 1:705493/151‑1 (MQ=255) aaTTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATaa < 2:758397/151‑1 (MQ=255) ccAAGCAGAAAAAAACCCAGCCGGTGGCGGGGAAGGACGATCCTGGCAAAAGATGCTTTCTGTCGCCGGATGCGGCGTAAACGCCTTATCCGACCTCAGCAGCTTTGCAAGTTGAATCAGTTCCAATTCAACGTGGAGGCCCGAGAagcgc > 1:526451/1‑151 (MQ=255) aaacCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAgg < 1:919766/151‑1 (MQ=255) | GCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAAT‑CGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGG > NC_000913/2001770‑2002065 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |