Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 151,656 | 0 | T | G | 20.0% | 180.5 / 32.2 | 65 | A192A (GCA→GCC) | yadL | putative fimbrial protein YadL |
Reads supporting (aligned to +/- strand): ref base T (0/52); new base G (13/0); total (13/52) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.09e-14 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.84e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTG‑TTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCAC > NC_000913/151594‑151720 | aTGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGc < 1:160452/65‑1 (MQ=255) gCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAg < 1:1053793/65‑1 (MQ=255) cATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGa < 1:1308901/65‑1 (MQ=255) tatataGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGaaa < 1:2267642/65‑1 (MQ=255) atataGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGGCGaaaa > 1:2883616/1‑65 (MQ=255) tataGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAAc < 1:645291/65‑1 (MQ=255) taGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGGGGAAAAAtt > 1:1404375/1‑65 (MQ=255) taGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGCGGGAAAAtt > 1:1764022/1‑65 (MQ=255) aGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGGAGGAAAATTa > 1:2880612/1‑65 (MQ=255) aGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTa < 1:2904923/65‑1 (MQ=255) aGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGGAGGAAAAttt > 1:2684889/1‑64 (MQ=255) aGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGGAGAAAACttt > 1:1384197/1‑64 (MQ=255) aGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGGAGAAAACTTa > 1:1629744/1‑65 (MQ=255) aGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGGAGAAAACTTa > 1:478604/1‑65 (MQ=255) aGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGGAGAAAACTTa > 1:2311160/1‑65 (MQ=255) aGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGGAGAAAAATTa > 1:2273345/1‑65 (MQ=255) aGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGCGGAAAACATa > 1:1181911/1‑65 (MQ=255) aGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGCAGAAGACCTa > 1:1698580/1‑65 (MQ=255) aGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGCAGAAAACTTa > 1:394666/1‑65 (MQ=255) aGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGCAGAAAACTTa > 1:94756/1‑65 (MQ=255) aGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGCAGAAAAAttt > 1:1432704/1‑64 (MQ=255) aGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGCAGAAAAATTa > 1:1036183/1‑65 (MQ=255) aGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGCAGAAAAATTa > 1:2881020/1‑65 (MQ=255) ttCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTAcc < 1:182072/65‑1 (MQ=255) ccTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGGCGAGAACATACCCg > 1:458298/1‑65 (MQ=255) cTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTTCCCGc < 1:2879646/65‑1 (MQ=255) tGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCa < 1:2657848/65‑1 (MQ=255) tGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCa < 1:870974/65‑1 (MQ=255) aaTCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGt < 1:2738245/65‑1 (MQ=255) aaTCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGt < 1:2017915/65‑1 (MQ=255) aaTCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGt < 1:2258272/65‑1 (MQ=255) aaTCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGt < 1:1354351/65‑1 (MQ=255) gTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACAt < 1:353480/65‑1 (MQ=255) aGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTG‑t < 1:588117/65‑1 (MQ=255) aGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTG‑t < 1:1161665/65‑1 (MQ=255) aGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTG‑t < 1:278013/65‑1 (MQ=255) aGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTG‑t < 1:189618/65‑1 (MQ=255) aGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTG‑t < 1:135297/65‑1 (MQ=255) aGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTG‑t < 1:2611155/65‑1 (MQ=255) aGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTAGTGCAGAAAACTTACCCGCAGTGACATTG‑t < 1:1632953/65‑1 (MQ=255) tCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGGGGAAAAATTACCCGCCGGGAAATTG‑ttt > 1:743746/1‑65 (MQ=255) tCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGGAGAAAAATTACCCGCAGGGGCATTG‑ttt > 1:2257415/1‑65 (MQ=255) tCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGGGCAGAAAAATTACCCCCAGGGACATTG‑ttt > 1:853692/1‑65 (MQ=255) ttagttaTAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTG‑TTTACgg < 1:1570355/65‑1 (MQ=255) tagttaTAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTG‑TTTACGGt < 1:2674415/65‑1 (MQ=255) tagttaTAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTG‑TTTACGGt < 1:2677743/65‑1 (MQ=255) tagttaTAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTG‑TTTACGGt < 1:2529636/65‑1 (MQ=255) agttaTAGGTTACTGTAAATGTGGCTGGGGGAGGAAAAATACCCCCCGGGGCATTC‑TTTTCgggt > 1:1694604/1‑63 (MQ=255) ataGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTG‑TTTACGGTTTTAt < 1:2867281/65‑1 (MQ=255) ataGGTTACTGTAAATGTGGCTGGGGGAGAAAAAATACCCCCAGGGGCATTGTTTTAGGGttttt > 1:2349743/1‑64 (MQ=255) ggTTACTGTAAATGTGGCTGGGGGGGAAAAACTACCCCCCGGGGCCTTT‑TTTACGGGtttttttt > 1:1881768/1‑65 (MQ=255) ggTTACTGTAAATGTGGCTGGGGGAGGAAACTTACCCGCCGGGGCATTT‑TTTTTGGTtttttttt > 1:2090483/1‑65 (MQ=255) ggTTACTGTAAATGTGGCTGGGGGAGGAAAATTACCCGCAGCGGCATTG‑TTTTCGGTtttttttt > 1:1876625/1‑65 (MQ=255) ggTTACTGTAAATGTGGCTGGGGGAGAAAAATTATCCGCAGTGACGATG‑TTTACGGCtttttttt > 1:1709500/1‑65 (MQ=255) ggTTACTGTAAATGTGGCTGGGGGAGAAAAATTACCCGCCGGGGCATTT‑TTTTTGGTtttttttt > 1:1956279/1‑65 (MQ=255) ggTTACTGTAAATGTGGCTGGGGGAGAAAAATTACCCCCGGTGTCATTT‑TTTTTGGGtttttttt > 1:1385516/1‑65 (MQ=255) ggTTACTGTAAATGTGGCTGGGGGAGAAAAATTACCCCCCGGGGCATTG‑TTTTCGGGtttttttt > 1:808212/1‑65 (MQ=255) ggTTACTGTAAATGTGGCTGGGGGAGAAAAATTAACCCCCGTGGCATTG‑TTTTCGGGtttttttt > 1:1060275/1‑65 (MQ=255) ggTTACTGTAAATGTGGCTGGGGCAGGAAACATTCCCCCCGTGACAATG‑TTTTCGGGTTTCtttt > 1:950365/1‑65 (MQ=255) ggTTACTGTAAATGTGGCTGGGGCAGAAAAATTAACCCCAATGCCACTG‑TTTACGGTTTTAtttt > 1:78952/1‑65 (MQ=255) aCTGTAAATGTGGCTGGGGGGGAAAACTTACCCGCAGGGACACTT‑TTTTCGGGTttttttttttt > 1:934185/1‑64 (MQ=255) aCTGTAAATGTGGCTGGGGGGGAAAAATTACCCGCAGTGGCATTG‑TTTTTGGttttttttttttt > 1:2303432/1‑64 (MQ=255) tGTAAATGTGGCTGGGGGCGAAAAATTACCCGCCGCGACATTT‑TTTTTGGtttttttttttttta > 1:2608919/1‑62 (MQ=255) tGTAAATGTGGCTGGGGCCGAAAAATTACCCGCCGGGACAATGCTTTAGGGTTTTTTTTTTTTcc > 1:1341050/1‑65 (MQ=255) gTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTG‑TTTACGGTTTTATTTTTTTCCAc < 1:656322/65‑1 (MQ=255) gTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTG‑TTTACGGTTTTATTTTTTTCCAc < 1:2898430/65‑1 (MQ=255) aaaTGTGGCTTGTGCAGAAAACTTACCCGCAGTGACATTG‑TTTACGGTTTTATTTTTTTCCACGa < 1:207934/65‑1 (MQ=255) aaaTGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTG‑TTTACGGTTTTATTTTTTTCCACGa < 1:1403673/65‑1 (MQ=255) aaaTGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTG‑TTTACGGTTTTATTTTTTTCCACGa < 1:1277807/65‑1 (MQ=255) aaaTGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTG‑TTTACGGTTTTATTTTTTTCCACGa < 1:2336555/65‑1 (MQ=255) aaaTGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTG‑TTTACGGTTTTATTTTTTTCCACGa < 1:2353854/65‑1 (MQ=255) aaaTGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTG‑TCTACGGTTTTATTTTTTTCCACGa < 1:1152761/65‑1 (MQ=255) aaTGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTG‑TTTACGGTTTTATTTTTTTCCACGAc < 1:593742/65‑1 (MQ=255) aTGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTG‑TTTACGGTTTTATTTTTTTCCACGAcc < 1:1669419/65‑1 (MQ=255) tgGCTGGTGCAGAAAACTTACCCGCAGTGACATTG‑TTTACGGTTTTATTTTTTTCCACGACCAgg < 1:1442333/65‑1 (MQ=255) tgGCTGGTGCAGAAAACTTACCCGCAGTGACATTG‑TTTACGGTTTTATTTTTTTCCACGACCAgg < 1:1891688/65‑1 (MQ=255) tgGCTGGTGCAGAAAACTTACCCGCAGTGACATTG‑TTTACGGTTTTATTTTTTTCCACGACCAgg < 1:846417/65‑1 (MQ=255) tgGCTGGTGCAGAAAACTTACCCGCAGTGACATTG‑TTTACGGTTTTATTTTTTTCCACGACCAgg < 1:2495524/65‑1 (MQ=255) tgGCTGGTGCAGAAAACTTACCCGCAGTGACATTG‑TTTACGGTTTTATTTTTTTCCACGACCAgg < 1:2516976/65‑1 (MQ=255) tgGCTGGTGCAGAAAACTTACCCGCAGTGACATTG‑TTTACGGTTTTATTTTTTTCCACGACCAgg < 1:2172477/65‑1 (MQ=255) tgGCTGGTGCAGAAAACTTACCCGCAGTGACATTG‑TTTACGGTTTTATTTTTTTCCACGACCAgg < 1:1109092/65‑1 (MQ=255) tgGCTGGTGCAGAAAACTTACCCGCAGTGACATTG‑TTTACGGTTTTATTTTTTTCCACGACCAgg < 1:399933/65‑1 (MQ=255) tgGCTGGTGCAGAAAACTTACCCGCAGTGACATTG‑TTTACGGTTTTATTTTTTTCCACGACCAgg < 1:1219812/65‑1 (MQ=255) tgGCTGGTGCAGAAAACTTACCCGCAGTGACATTG‑TTTACGGTTTTATTTTTTTCCACGACCAgg < 1:2054612/65‑1 (MQ=255) tgGCTGGTGCAGAAAACTTACCCGCAGTGACATTG‑TTTACGGTTTTATTTTTTTCCACGACCAgg < 1:2932237/65‑1 (MQ=255) tgGCTGGTGCAGAAAACTTACCCGCAGTGACATTG‑TTTACGGTTTTATTTTTTTCCACGACCAgg < 1:1226025/65‑1 (MQ=255) tgGCTGGTGCAGAAAACTTACCCGCAGTGACATTG‑TTTACGGTTTTATTTTTTTCCACGACCAgg < 1:1019299/65‑1 (MQ=255) tgGCTGGTGCAGAAAACTTACCCGCAGTGACATTG‑TTTACGGTTTTATTTTTTTCCACGACCAgg < 1:1249561/65‑1 (MQ=255) tgGCTGGTGCAGAAAACTTACCCGCAGTGACATTG‑TTTACGGTTTTATTTTTTTCCACGACCAgg < 1:131265/65‑1 (MQ=255) tgGCTGGTGCAGAAAACTTACCCGCAGTGACATTG‑TTTACGGTTTTATTTTTTTCCACGACCAgg < 1:1920542/65‑1 (MQ=255) tgGCTGGGTCAGAAAACTTACCCGCAGTGACATTG‑TTTACGGTTTTATTTTTTTCCACGACCAgg < 1:715278/65‑1 (MQ=255) ggTGCAGAAAACTTACCCGCAGTGACATTG‑TTTACGGTTTTATTTTTTTCCACGACCAGGCGTGc < 1:2910120/65‑1 (MQ=255) ggTGCAGAAAACTTACCCGCAGTGACATTG‑TTTACGGTTTTATTTTTTTCCACGACCAGGCGTGc < 1:1089248/65‑1 (MQ=255) tGCAGAAAACTTACCCGCAGTGACATTG‑TTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCAc < 1:1806454/65‑1 (MQ=255) | ATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAACTTACCCGCAGTGACATTG‑TTTACGGTTTTATTTTTTTCCACGACCAGGCGTGCAC > NC_000913/151594‑151720 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |