Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,001,911 | 0 | G | A | 40.0% | 46.1 / 19.1 | 41 | 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 (21/3); new base A (0/16); total (21/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.54e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
ACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATT‑CCA‑GG‑CAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGC > NC_000913/2001763‑2002059 | acCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATT‑CCA‑GG‑CAGa > 2:658576/1‑150 (MQ=255) ccTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATT‑CCA‑GG‑CAGaaa < 1:302448/151‑1 (MQ=255) ataGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATT‑CCA‑GG‑CAGAAAAAAACCccgcc > 2:15791/1‑151 (MQ=255) ataGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATT‑CCA‑GG‑CAGAAAAAAACCccgcc > 1:145093/1‑151 (MQ=255) agTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATT‑CCA‑GG‑CAGAAAAAAACCCCGCCGGTg > 2:346284/1‑151 (MQ=255) gAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATT‑CCA‑GG‑CAGAAAAAAACCCCGCCGGTGGCg > 2:401547/1‑151 (MQ=255) aaaaaCGGCCCTGCTTTATATGAGTTAGCGGCATGTTTATCCGTTTCTGCGGGTTTTTTAAGCGGACGATTAGTGGGTAAAAGGAGGGTTTATTGGGGGGTTAAAGGTAAATT‑CCA‑GG‑CAAAAAAAACCCCCGCCGGTGGCGGGGAAGCa < 2:477292/148‑1 (MQ=255) aaacgccctggctttATATGATTTATCGGGATGTTTATCCGTTTCTGCAGGTTTTTAAATCTGACGTTTAGTGGGTAAAAGGAGGGGTTATTGGGGGGTTAAAGG‑AAATT‑CCC‑GGACAAAAAAAACCCCCGCCGGTGGCGGGGAAGCAc < 1:516297/136‑1 (MQ=255) aaaCAGCCCTGCGTTATATGAGTTACCGACATAATTATCCGTTTCTGGAGGTTTTTTAAGCGGACGATTGGTGGGTGAAAGGAGGGGTTTTTGGGGGGTTAAAGGTAAATC‑CCA‑GC‑CAGAAAAAAACCCCGCCGGTGGCGGGGAAGc < 1:166032/147‑1 (MQ=255) aaCAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATT‑CCA‑GG‑CAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGtt > 2:188961/1‑151 (MQ=255) cGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATT‑CCA‑GG‑CAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAAt > 1:586158/1‑151 (MQ=255) tatGAGTTTCCGGATTGATTATCCGTTTCTGCGGGTTTTTAAATCGGACGATTAGTGGGTGAAAGGGGGGTTTATTGGGGGGTTACAGGTAAATT‑CCA‑GG‑CAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgc < 1:376965/149‑1 (MQ=255) ggCATGTTTCTCCTTTTCCGGAGGTTTTTTAAGCGAACGATTAGTGGGTAAAAGGGGGGGTTTTTTGGGGGTTAC‑GTTAAATTCCCA‑GT‑AAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTg < 1:73961/149‑1 (MQ=255) tttCTGCGGGGTTTTTAATCGGACAATTGGTGGGTAAAATGGGGGGTTATTTGGGGGTTACAGGAAAATC‑CCA‑GG‑CAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcggcg < 1:17880/147‑1 (MQ=255) tttCTGCAGGGTTTTTAATCGGACGATTGGGGGGTAAAATGGGGGTTTTTTTGGGGTTTACAGGTAAATT‑CCA‑GG‑CAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTa < 1:222736/149‑1 (MQ=255) tttCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATT‑CCA‑GG‑CAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTaa > 2:533730/1‑150 (MQ=255) tttCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATT‑CCA‑GG‑CAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTaa > 2:84982/1‑150 (MQ=255) gcggggTTTTTAAGCGGAGGATTAGTGGTTGAAATGGGGGTTTTTTGGGGGGTTAAAGGAAATTC‑CCA‑GA‑AAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCt < 1:29070/148‑1 (MQ=255) aGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATT‑CCA‑GG‑CAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCtt > 1:206263/1‑150 (MQ=255) ggttttttAATCGGACGATTAGTGGGTAAAATGGGGGTTTATTGGGGGGTTACAGGTAAATT‑CCA‑GG‑CAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCtt < 2:262905/146‑1 (MQ=255) taagggggTGAAAGGAGGGGTTATTTGGGGGTAAAAGGAAAATT‑CCA‑GG‑CACAAAAACACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACa < 1:363295/138‑1 (MQ=38) tAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATT‑CCA‑GG‑CAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGc > 1:460968/1‑151 (MQ=255) aGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATT‑CCA‑GG‑CAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGc > 1:60281/1‑150 (MQ=255) gTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATT‑CCA‑GG‑CAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCa > 2:370589/1‑150 (MQ=255) gTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATT‑CCA‑GG‑CAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAAtt > 1:686972/1‑150 (MQ=255) gAAATGGGGTTTTTTTTGGGGTTTACAGTTAAATT‑CCA‑GA‑CAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtgtg < 1:710263/142‑1 (MQ=37) aTGAGGGGTTATTTGGGGGTTACAGGTAAATT‑CCA‑GG‑CAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATa > 2:344171/1‑150 (MQ=255) gAGGGGTTTTTGGGGGGTTACAGTTAAATC‑CAA‑GA‑CAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAt < 1:401455/151‑1 (MQ=37) ggggtttaaagaaaaTT‑CCCGGG‑CAAAAAAACCCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAt < 1:360697/133‑1 (MQ=25) ggggTTACAGGTAAATC‑CCA‑GG‑AAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACtt < 2:481829/151‑1 (MQ=255) ggttaaaaGGAAAATC‑CCA‑GG‑AAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTg < 1:553801/144‑1 (MQ=25) gtttACAGGAAAATT‑CCG‑GG‑CAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGt < 2:320318/149‑1 (MQ=255) ttACAGGTAAATT‑CCA‑GG‑CAGAAAAAAACCCCGCCGGAGGCGGGGAAGCACGTTGCTGACAAATTGAGCTTTATGTTGCCGGAAGCGGCGCAAACGCCTTATCCGGCCTGCAAAAATGTGCAAATTGAATAAATTGCAATTCAAc > 1:520384/1‑145 (MQ=255) cAGGTAAATT‑CCA‑GG‑CAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGcc > 2:129710/1‑151 (MQ=255) aGGTAAATT‑CCA‑GG‑CAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCt > 2:599846/1‑151 (MQ=255) ggaaaaTC‑CCA‑GG‑AAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTg < 2:111897/148‑1 (MQ=255) ggTAAATT‑CAA‑GC‑AAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCt < 1:658489/150‑1 (MQ=255) gTAAATT‑CCA‑GG‑CAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTg > 1:580641/1‑150 (MQ=255) aaaTC‑CCA‑GG‑CAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATa < 2:180189/151‑1 (MQ=255) aaTT‑CCA‑GG‑CAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATa > 2:259527/1‑150 (MQ=255) aGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcagcgc > 2:594506/1‑150 (MQ=255) | ACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATT‑CCA‑GG‑CAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGC > NC_000913/2001763‑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: ‑ |