| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,001,918 | 0 | A | C | 45.8% | 62.0 / 14.5 | 49 | 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 (22/4); new base C (0/22); total (22/27) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.67e-10 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
CAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGAT‑TAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCA > NC_000913/2001768‑2002067 | cAGCGGTATAGAGTGCATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGGCTGATTACCCGTTTCTGCAGGGTTTTTAATCGGACGAT‑GAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCATGCAGaaaaaaa < 2:343875/151‑1 (MQ=255) aGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGAT‑TAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAAc < 1:34121/151‑1 (MQ=255) cGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGAT‑TAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAAccc > 1:788180/1‑151 (MQ=255) agTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGAT‑TAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTg > 2:753126/1‑151 (MQ=255) aaTTCACGATAAACAGCCCTGCGTTATATGAGTGATCGGCATTCTTATCCGGTTCTGCAGGGTTTTTAATCGGACGAT‑TAGGGGGTGAAATGAGGGGGTATTTGGTGG‑TTACAGGAAAATTCCAGGCAGAAAAAAACACCGCCGGTGGCgg > 2:382772/1‑151 (MQ=255) ttCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCGGGTTTTTTAATCGGACGAT‑TGGGGGGTAAAAGGAGGGGTTTTTTGGGGG‑TTACAGGTAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCggg < 2:946965/150‑1 (MQ=255) tCAAGATAAACAGACCTGCGATATATGAGTTATCGGCATGATTATCCGTTTCTAAAGGGTTTTTAATTGGACGAG‑GATTGGGTGAAAAGAGGGGTTATTTGGGGT‑GAACAGGCAAATTCAAGGCAGAAAAAAACCCCGACGGTGgggggg > 2:689489/1‑150 (MQ=255) aTAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGAT‑TAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGc > 2:712010/1‑149 (MQ=255) aaaCAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGAT‑TAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGt > 1:869945/1‑151 (MQ=255) aaCGGCCCGGCGTTATAGGAGTTATGGGAATGATTATCCGTTTCTGCAGGTTTTTAAATCGGAGAAT‑TGGTGGGTGAAATGAGGGTTTATTTGGGGG‑TTACGGTTAAATTCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGtt < 2:474478/151‑1 (MQ=255) aGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGAT‑TAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCt > 2:346847/1‑151 (MQ=255) cTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGAT‑TAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCAAGTTGCTGa > 2:792667/1‑149 (MQ=255) gctttATATGATTTAGCGGCATGATTATCCGTTGCTGGAGGGTTTTTAATCGGACGAT‑GAGTGGGTGAAATGAGGGTTTATTTGGGGG‑TTACAGATAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACaa < 2:584740/147‑1 (MQ=255) ttATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGAT‑TAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTg > 1:443868/1‑151 (MQ=255) atGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGAT‑TAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCt > 1:1006555/1‑151 (MQ=255) gTTATCGGTATGATTATCCTTTTTGGCGGGGTTTTTAAGCGAACGAT‑GGGTGGGTGAAGGGGGGGGTTATTTGGGGT‑TGACAGGTAAATTCGAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTa < 2:555378/150‑1 (MQ=255) ttATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGAT‑TAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGAAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATg > 2:873645/1‑151 (MQ=255) ggCATGATTATCCGTTTCTGCGGGGTTTTTAATCGGACGAT‑TAGGGGGTGAAATGGGGGGTTATTGGGGGG‑TTACAGGTAAATCCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGc < 1:843970/150‑1 (MQ=255) ttctctttttCTGTAGGTTTTTTAAGCGGACGAT‑TGGTGGGTGAAAGGAGGGGTTTTTTGGGGT‑TTACATGAAAATTCCGGAAAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGAt < 2:992363/141‑1 (MQ=255) acccGTTTCTGAGGGGTTTTTAATCGGACGAT‑TGGTGGGTGAAATGGGGGTTTTTTTGGGGG‑TTACAGGTAAATCCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcg < 2:21335/147‑1 (MQ=255) aTCCGTTTCTGCAGGGTTTTTAATCGGACGAT‑TAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcgg > 2:451062/1‑150 (MQ=255) tttCTGAGGTGTTTTTAATGGGACGAT‑TAGGGGGTGAAGGGAGGGGTTTTTTGGGGG‑TGAGGGGAAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGc < 2:66793/136‑1 (MQ=255) ttgttgtgggTTTTTAAACGGAAGGTGTGGTGGTAGAAAGGAGGAGTTTTTTGGGGT‑TTACAGAAAACTTCCGGGAAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATGGCGCTTTATGTTTCCGGATGCGGCGt < 2:1016068/141‑1 (MQ=255) ttcttcAGGGTTTTTAATCGAATGTT‑TGGTGGGAGAAAGTAGGGTTTATTTGGGGG‑TTAGAGAAAATGTCCGGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTTTGTTGCCGGATGCGGCGTaaa < 2:982283/150‑1 (MQ=255) ggttttttAAGCGGACGAT‑TAGTGGTGGAAAGGAGGGGTTTTTTGGGGG‑TTAAAGTTAAATTCCAGAAAAAAAAACCCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCtt < 2:1020892/146‑1 (MQ=255) ttttAATCGGAAGAT‑TGGGGGGTAAAAGGAGGGTTTTTTTGGGGG‑TAACAGTAAAATTCCGGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCg < 2:232643/150‑1 (MQ=37) tAATCGGACGAT‑TAGTGGGTGCAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCCGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGGTGCCGGATGCGGCGTCAACGCCTTATCCGGCCt > 2:848846/1‑151 (MQ=255) gaagAT‑TGGTGGGTAAAATG‑GGGGTTTTTTGGGGGTTTAAAGGTAATTACCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaa < 1:607682/148‑1 (MQ=25) t‑tAGGGGGTAAAATGGGGGTTTTTTTGGGGG‑TTACAGGTAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAt < 1:309996/148‑1 (MQ=37) tGGGTAAAATGAGGGTTTTTTTGGGGG‑TGACGGGAAAATCCAAGACAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCa < 2:75568/149‑1 (MQ=37) gggggtaaaaggggggTTTTTTGGGGGTTTACAGTAAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtgtg < 1:346996/136‑1 (MQ=25) gtaaaaTGAGTGTTTGTTTGGGGG‑TTACGGGTAAGTTGCAGGAAAAAAAAACCCCCCCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCaa < 2:538893/144‑1 (MQ=25) gTGAAATGGGGGGTTATTTGGGGG‑TTACAGGTAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTc < 1:715124/151‑1 (MQ=255) aaaTGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCa > 2:594259/1‑149 (MQ=255) aaaTGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAAt > 1:430359/1‑151 (MQ=255) ctGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGGGCAAATTCAATa > 2:198077/2‑150 (MQ=255) gggtTTTTTTGGGGT‑TTACAGGTAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTg < 1:1065676/151‑1 (MQ=37) ggggTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTg > 1:459219/1‑151 (MQ=255) gggTTTTTTGGGGG‑TTCCGGGAATTTCCCAGAAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTgg < 2:787915/151‑2 (MQ=25) ggTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGc > 1:1009342/1‑150 (MQ=255) ttGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCa > 1:278419/1‑150 (MQ=255) gggg‑TAACGGATAAATCCCAGGAAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACtt < 2:487250/151‑1 (MQ=37) ttACAGTAAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTa < 2:434052/150‑1 (MQ=255) cAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGcc > 1:817424/1‑151 (MQ=255) agggaaaaGTCCAAGAAGAAAAAATCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAg < 2:736300/143‑1 (MQ=37) agttAAATGCCGGGAAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGcc < 2:649209/147‑1 (MQ=255) ggTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTg > 2:404873/1‑151 (MQ=255) aattaccaggagaaaaaaacCCCCGCCAGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCAGCGTAAACGCCTTAGCCGGCCTACAAAAATGTGCAAATTCAATAAATGGCAATTCAACTTGTAGGCCTGAt < 2:889726/137‑1 (MQ=255) aaaTTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGAt < 1:386544/150‑1 (MQ=255) tcccAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAag < 1:975195/148‑1 (MQ=255) aGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgc > 2:597602/1‑149 (MQ=255) gaaaaaaaaacCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcagcg < 2:1043862/147‑1 (MQ=255) aaaaaaaCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATc > 2:831348/1‑151 (MQ=255) aaaaaaCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCa < 1:866148/151‑1 (MQ=255) aaaaCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAgg > 1:262203/1‑151 (MQ=255) aaaCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGc > 2:193044/1‑151 (MQ=255) aCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCa > 1:529093/1‑150 (MQ=255) | CAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGAT‑TAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCA > NC_000913/2001768‑2002067 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |