Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 266,349 | 0 | C | A | 44.4% | 99.5 / 31.4 | 55 | V69F (GTT→TTT) | yafY | lipoprotein, inner membrane; degP regulator; CP4‑6 prophage |
Reads supporting (aligned to +/- strand): ref base C (4/26); new base A (24/0); total (29/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.48e-11 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
TCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTAC‑GTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGAT > NC_000913/266204‑266499 | tCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGGAC‑GGAAAAACTGc > 1:347162/1‑150 (MQ=255) cACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTAC‑GTAAAACCTGCAAGGTg < 1:654718/151‑1 (MQ=255) cACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGGAC‑GTAAAAACTGCAAgggg > 1:136460/1‑149 (MQ=255) tGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGGAC‑GTAAAAACCGCCAAGGGGGGCCgg > 2:468891/1‑151 (MQ=255) gTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTAC‑GTAAAACCTGCAAGGTGGGACCGGa < 2:735653/150‑1 (MQ=255) gTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGGGC‑GTAAAAACCGCAAGGTGGGGCCgg > 1:368760/1‑149 (MQ=255) aTTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTAC‑GTAAAAACTGCAAGGGGGGGCCGGATa > 1:532513/1‑150 (MQ=255) tCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGGAC‑GGAAAAACCGCAAAGGGGGGCCgg > 2:777144/1‑145 (MQ=255) ccAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTAC‑GTAAAACCTGCAAGGTGGGACCGGATAagt < 2:150616/150‑1 (MQ=255) ccAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTAC‑GTAAAACCTGCAAGGTGGGACCGGATAagt < 2:945274/150‑1 (MQ=255) cAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGGAC‑GTAAAACCCGCAAGGGGGGGCCggggaa > 1:449654/1‑143 (MQ=255) tttGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTACGGTAAAAGTTGCGGGGGGGGTCCggcgagcg > 2:397377/1‑141 (MQ=255) tGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTAC‑GTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCaa < 1:288156/150‑1 (MQ=255) tGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTAC‑GTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCaa < 1:840724/150‑1 (MQ=255) tATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTAC‑GTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCtgtg < 1:94214/150‑1 (MQ=255) cGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTAC‑GTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCa < 1:257814/151‑1 (MQ=255) cGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGGAC‑GGAAAAACCGCAAGGGGGGGCCGGGTAAGGAGGCAAACTGGGaac > 2:612510/1‑148 (MQ=255) cAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTAC‑GTAAAAACTGCAAGGGGGGGCCGGATAAAGAGTCAACCTGTGACCaat > 1:495315/1‑148 (MQ=255) aaTACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTAC‑GTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCAttct < 2:669612/151‑1 (MQ=255) cATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTAC‑GTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGAcatc < 1:202315/151‑1 (MQ=255) tcatcaGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGGAA‑GTAAAAACCGCCAGGGGGGACCCGGTAAAAAGTCAAACCGTGGCCCATCTTCGCCGAcatca > 1:657127/1‑149 (MQ=255) atcaGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTAC‑GTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCt < 2:619993/151‑1 (MQ=255) tGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGGGA‑GGAAAAACAGCAAGGCGGGGACGGGCAAGGAGTGCACCTGTGCCCATTCTTTACTGTCATAATCATGGGGGCTg > 2:829255/1‑150 (MQ=255) ggACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGGAC‑GTAAAAACCGCAAGGGGGGGCCGGGTAAGTAGTCAACCCGTGACCATTCTTTGCTGACATCATCCTGGGGGGTg > 2:783671/1‑149 (MQ=255) ttttttCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGGAC‑GTAAAAACTTCAAAGTGGGGCCGGACAAAAAGGCAACCCTTGTCCATTCTTTGCTGGTATTATCCCGGTTGCTGa > 1:305332/1‑144 (MQ=255) tttCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGGAA‑AGAAAAACCGCAAGGTGGGGCCGGATAAGGAAACAACCTTTGACCCTTTTTTGCTGTCAATATCCCGGGGGGTGGGCACATc > 1:141953/1‑148 (MQ=255) ttCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTAC‑GTAAAAACTGCAAGGGGGGGCCGGGTAAAAAGGCAACCTGTGACCATTTTTCGCTGACATCATCCTTGGGGCTGAcccc > 1:624432/1‑142 (MQ=255) ccTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGGAC‑GTAAAAACTGCAAGGGGGGGCCGGATAAAGAGGCAAACTGTGACCATTTTTTCCTGACATCATCCTGGTGGCGGGCCACCTCACtta > 1:171347/1‑148 (MQ=255) tGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTAC‑GTAAAACCTGCAAGGTGGGACCTGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGaa < 2:253538/149‑1 (MQ=255) gCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTAC‑GTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGaaa < 2:935844/149‑1 (MQ=255) cATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGGAC‑GTAAAAACTGCAAGGGGGGGCCGGGGAAGTAGTCAACCTGTGACCATTCTTCGCCTACATCATCCTCGTGGGTGACCACATCACCGaaa > 1:562725/1‑148 (MQ=21) ccACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTAC‑GTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTg < 2:436110/150‑1 (MQ=255) aCGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGGAA‑GTAAAAACCGCCAAGGGGGACCGGGAAAGTAGTCCACCTGTGGCCATTCCTCCCTGACCTCATCCTGGTGGGTGACCCCATAACTCAAACGCGGCgggg > 2:191775/1‑151 (MQ=255) cGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTAC‑GTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGt < 2:268869/151‑1 (MQ=255) cGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTAC‑GTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGt < 2:851356/151‑1 (MQ=255) agtagtaACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGGAA‑GTAAAAACGGCAAGGGGGGGCCCGGTAAGAAGTCCAACTGTGACCATTCTTCGCTGACATCCCCCTGGCGGCAGACCACATCACTTAAAAGCGgggggggg > 1:351168/1‑144 (MQ=12) gtagtaACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTAC‑GTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTcc < 2:950856/151‑1 (MQ=255) gtaACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTAC‑GTAAAAACCGCAAGGGGGGACCGGATAAGGAGTCAACCCGTGACCATTCTTCGCTGACATCATCCTGGTGGGTGACCACATCACTGAAACGCGGGGGGGGGAGGTCCttt > 1:581822/1‑150 (MQ=18) gtaACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGGAC‑GGAAAAACCGCAAGGTGGGGCCGGGTAAAGAGTCAAACCGTGACCATTCTTCGCTGACATCATCCTGGGGGCGGACCACATCACTGAAACGCGGGGGGGGGGGGGCCt > 2:172290/1‑149 (MQ=18) aaCTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTAA‑GTAAAAACCGCAAGGTGGGGCCCGCTAAGAAGTCAACCTTTGACCCTTCTTCTCTTACATCATCCCCGTGGGTGGCCACCACACTTAAAAGCGGTGGGGGGGGGGCCttta > 2:588874/1‑148 (MQ=255) aCTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTAC‑GTAAAAACTGCAAAGGGGGACCGGATAAATAATCAACCCGTGAACACTTTTTCCTGACATCATCCTGGGGGGTGACCAAATCACTGAAACGCGGGGGGGGGAGGGCCTTaaa > 1:787370/1‑150 (MQ=255) tGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTAC‑GGAAAAACTGCAAGGGGGGGCCGGATAAATAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGGGGGTGGCCACATCACTGAAACGCGGGGGGGGGGGGGCCTTTAAttt > 2:358009/1‑151 (MQ=18) ccATCCTCGCAGCCATCGGGCGACTGGC‑GGGTAC‑GTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTAc < 1:231285/151‑1 (MQ=255) ccATCCTCGCAGCCATCGGGCGACTGGC‑GGGTAC‑GTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTAc < 1:201791/151‑1 (MQ=255) ccTCGCAGCCATCGGGCGACTGGC‑GGGTAC‑GTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGAt < 2:406526/151‑1 (MQ=255) ccTCGCAGCCATCGGGCGACTGGC‑GGGTAC‑GTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGAt < 2:788675/151‑1 (MQ=255) cTCGCAGCCATCGGGCGACTGGC‑GGGTAC‑GTAAAAACTGCAAGGTGGGGCCCGATAAGTAGTCAACCTGTGACCATTCTTCTCTGACATCATCCTGGGGGCTGACCACATCACTGAAAAGCGGTGGGGGGAGGGCCTTTAATTTACTGAt > 2:622324/1‑150 (MQ=255) aGCCATCGGGCGACTGGCGGGGCAG‑ATAAAACC‑GCAAAGGGGGGCCGGATAAAAAGTTAACCCGTGTCCACTCTTCGCTGACATCCTTCTGGTGGGTGACCACCTCAACGAAACACGGAGGGGGGAGGGGCTTAAAATTAATGATGGCCtt > 2:961345/1‑151 (MQ=255) aGCCATCGGGCGACTGGC‑GGGTAC‑GTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCtt < 2:814675/151‑1 (MQ=255) cGGGCGACTGGC‑GGGTAC‑GTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGtc < 2:722556/151‑1 (MQ=255) ggC‑GGGTAC‑GTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTc < 2:622846/150‑1 (MQ=255) ggC‑GGGTAC‑GTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCt < 1:603455/151‑1 (MQ=255) c‑GGGTAC‑GTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGt > 2:198156/1‑151 (MQ=255) ggTAC‑GTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCa < 2:431756/151‑1 (MQ=255) tAC‑GTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATc > 1:250643/1‑151 (MQ=255) ccTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGAt > 2:826254/1‑151 (MQ=255) | TCCGTCACCTCCTGGTATTCCAGTTTGATGTTATCGGCAATACACAGCGCATTCATCAGTGGCTGGACGGTTTTTTCCTGCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGC‑GGGTAC‑GTAAAACCTGCAAGGTGGGACCGGATAAGTAGTCAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTGGGGTGAGGTCCTTAAATTTACTGATGGCCTTCAGGTCATCACTTCTGTCATCGCATGCGAT > NC_000913/266204‑266499 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |