Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,001,911 | 0 | G | A | 42.9% | 56.0 / 33.4 | 56 | intergenic (‑122/+199) | fliA/fliC | RNA polymerase, sigma 28 (sigma F) factor/flagellar filament structural protein (flagellin) |
Reads supporting (aligned to +/- strand): ref base G (28/4); new base A (0/24); total (28/28) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 9.40e-12 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
CCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGC > NC_000913/2001764‑2002059 | ccTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGa < 2:762344/149‑1 (MQ=255) cTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGaaa > 1:34526/1‑150 (MQ=255) cTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGaaa > 2:830911/1‑150 (MQ=255) cGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAAccc > 2:487270/1‑151 (MQ=255) ggTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAAcccc < 1:811763/151‑1 (MQ=255) gTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCGGGCAGAAAAAAAccc > 2:900109/1‑149 (MQ=255) tataGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCccgc > 2:270889/1‑151 (MQ=255) tataGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCccgc > 1:460159/1‑151 (MQ=255) aaTTCACAATAAACAGCCCTGCTTTATATGAGTTATCGGAATGATTATCCGTTTCTGCAGGTTTTTTAACCGGACGATTAGTGGGTAAAATGGGGGTTTTTTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAACCCCCGCCGGTg < 2:1107545/147‑1 (MQ=255) aTAAACAGCCCTGCGTTATATGAGTTATGGGCATGATTTTCCTTTTCTGCAGGTTTTTTAATCGGACGATTGGGGGGTAAAAGGAGGGGTTTTTGGGGGG‑TTACAGGTAAATT‑CCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGc < 1:889027/149‑1 (MQ=255) tAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACg > 2:385599/1‑151 (MQ=255) aGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCt > 2:441296/1‑151 (MQ=255) ccTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGAc > 1:817519/1‑151 (MQ=255) gTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAAtt > 2:707643/1‑151 (MQ=255) ttATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTg > 1:84102/1‑151 (MQ=255) ttATCGGCATGTTTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTAAAATGGGGGGTTATTTGGGGG‑TTACAGGTAAATC‑CCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCtt < 1:539447/147‑1 (MQ=255) atgggtagatttATCCGTTTCTGCGGGTTTTTAAATGAGACGATTAGTGGGTAAAATGAGGGTTTATTGGGGGG‑TTACAGGTAAATT‑CCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGt < 1:616687/140‑1 (MQ=255) aTCGGCATGATTATCCTTTTCTGCAGGTTTTTTAATCGGACGATTAGGGGGTGAAAGGGGGGGTTATTGGGGTT‑TTACGGGTAATTT‑CCGGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGtt < 1:457863/151‑1 (MQ=255) gCATGATTATCCTTTTCTGCAGGTTTTTTAATCGGACGTTTAGTGGGTGAAATGAGGGTTTTTTTGGGGG‑TTAAAGGTAAATC‑CCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCg < 2:476766/151‑1 (MQ=255) cATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCgg > 1:624840/1‑151 (MQ=255) ttATCCGTTTCTGCGGGTTTTTTAATCGGACGATTAGTGGGTAAAAGGGGGGGTTATTTGGGGG‑TTACAGGTAAATC‑CCAGGAAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGAt < 1:353553/148‑1 (MQ=255) ccGTTTCTGCAGGGTTTTTAATCGCACGATTAGTGGGTAAAATGAGGGGTTTTTTGGGGG‑TTAAAGGAAAATT‑CCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcgg < 2:846595/148‑1 (MQ=255) gTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTa > 2:152301/1‑150 (MQ=255) tCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAAc > 2:1103922/1‑150 (MQ=255) gCAGGTTTTTTAATCGGACGATTAGTGGGTGAAAGGAGGGTTTTTTTGGGGG‑TTACAGGTAAATT‑CCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAAc < 2:220883/147‑1 (MQ=255) gggTTTTTAATGGGACGATTGGGGGGTAAAAGGAGGGGTTTTTTGGGGG‑TTACAGGAAAATT‑CCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACg < 2:1077854/145‑1 (MQ=255) tttttAATGGGAAGATGGGGGGGT‑AAAAGGGGGGTTTTTTGGGGGTTTACAGGAAATTC‑CCGGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATcc < 2:181752/150‑1 (MQ=255) tttttAATCGGTCTTTTAGGGGGTAAAATGGGGGGTTTTTTGGGGG‑TAACAGTAAAATT‑CAAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATcc < 1:695974/150‑1 (MQ=255) tttttAATCGGACGATTAGGGGGTAAAAGGAGGGTTTATTGGGGGT‑TTACGGGTAAATC‑CCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATc < 1:57471/149‑1 (MQ=255) tttttAAACGGAAGATTAGGGGGGTAAATGGGGGGTTTTTTGGGGT‑TTACAGTAAAATC‑CCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATcc < 1:596347/150‑1 (MQ=255) ttAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGcc > 1:271769/1‑151 (MQ=255) aaTCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCt > 2:970861/1‑150 (MQ=255) cGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACa > 1:67368/1‑150 (MQ=255) aTTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtgt > 2:43639/1‑151 (MQ=255) ttAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtgtg > 2:151399/1‑151 (MQ=255) gTGGGTGAAAGGGGGGGTTTTTTGGGGG‑TTAAAGAAAAATT‑CCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaa < 2:1046427/143‑1 (MQ=37) gTGGGTAAAAGGAGGGTTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGc < 1:562181/149‑1 (MQ=255) gggaaaaaTGAGGGGTTTTTTGGGGT‑TAGCAGGAAAATT‑CCAGAAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAAt < 1:644904/146‑1 (MQ=25) gggTGAAATGGGGGTTTTTTTGGGGG‑TT‑CAGGTTAATTCCCAGGAAAAAAAACCCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCaaa < 1:529551/150‑1 (MQ=25) gggTAAAATGAGGGTTTTTTTGGGGG‑TTACAGGAAAATT‑CCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCaa < 2:666540/149‑1 (MQ=38) aTGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATa > 2:974410/1‑150 (MQ=255) gAGGGTTTTTTTGGGGG‑TTACAGGTAAATC‑CCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATaa < 2:436311/149‑1 (MQ=37) ttttgtGGGGG‑TTACAGGAAAATC‑CCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCa < 1:502278/144‑1 (MQ=37) aTTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTc > 1:460642/1‑151 (MQ=255) tttGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTc > 2:491753/1‑150 (MQ=255) ttGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCa > 2:732853/1‑150 (MQ=255) tGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAAc > 2:913743/1‑151 (MQ=255) ggggt‑TTACAGGAAAATT‑CCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAAc < 1:326967/150‑1 (MQ=255) tACAGGAAAATT‑CCAGGCAGAAAGAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTTACAAATTGCGCTTTATGTTGCCGGGTGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGGAAATTCAATAAATTGCAATTGAACTTGTa > 2:153840/1‑149 (MQ=38) ggTAAATT‑CCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTg > 2:1071495/1‑151 (MQ=255) aaaTC‑CCGGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGAt < 1:835646/150‑1 (MQ=255) aaaTC‑CCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAgg < 1:965984/144‑1 (MQ=255) aaTC‑CCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATaa < 1:270923/151‑1 (MQ=255) tt‑CCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagc > 1:340119/1‑151 (MQ=255) ggCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgca < 1:441358/149‑1 (MQ=255) aGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcagcgc > 1:82615/1‑150 (MQ=255) | CCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGC > NC_000913/2001764‑2002059 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |