| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,222,233 | 0 | G | A | 47.2% | 50.7 / 29.4 | 53 | pseudogene (60/75 nt) | ycgH | pseudogene; putative ATP‑binding component of a transport system |
| Reads supporting (aligned to +/- strand): ref base G (24/4); new base A (0/25); total (24/29) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.71e-11 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAT‑ATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTG > NC_000913/1222083‑1222382 | gtgtGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTg < 1:1352108/151‑1 (MQ=255) gAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGaaa < 1:779283/150‑1 (MQ=255) aaaGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAat > 1:186049/1‑150 (MQ=255) aaaGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAat > 2:525330/1‑150 (MQ=255) gAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAT‑ATAAGtt > 2:85818/1‑151 (MQ=255) aCAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAT‑ATAAGTTCt < 2:1044011/151‑1 (MQ=255) gAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAT‑ATAAGTTCTGATCCCGCCGtt > 2:1157696/1‑151 (MQ=255) gAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAT‑ATAAGTTCTGATCCCGCCGtt > 1:1296074/1‑151 (MQ=255) aGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAT‑ATAAGTTCTGATCCCGCCGtt > 1:1204916/1‑150 (MQ=255) gggggAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAT‑ATAAGTTCTGATCCCGCCGTTAGCTaaaa > 2:171176/1‑150 (MQ=255) gggAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAT‑ATAAGTTCTGATCCCGCCGTTAGCTAAAAAAc > 2:1097443/1‑151 (MQ=255) aaaTTAATCTGGGGAAAAAGGACAACTTAAACCCAGGGGCCGGCTTGTGGGGAGTGTGGGCGTGCGGCTGGGTAAAAGTGAAAAAAATGACCCCGAAGTGGTGGTGGTCCTAAAAA‑TAAATTTTTGATCCCGCCGTTAGCTAAAAAAcc < 1:388534/149‑1 (MQ=255) aaaTTAACCTGGGCTTAAAGGGCAACATGAACCCAGCGGCCAGGGTGAGGGGATGGGGGGCGGTAAGCGGGGAAATATGGTCTACAATGCCCCCCGAGTGATGGTGGCCAAAAAAT‑AGATGTTCTGGTCCCGCCGTTAGCTAAAAAACcg < 1:291987/150‑1 (MQ=255) aaaCTGGGCGTAAAGGGCACCCTAAACCCAGGGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCGGGGTGATAATGGCTACAATGCCCCCGAGGTGATGGTGGGCCTAAAAT‑ATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGcgtcg < 2:653903/150‑1 (MQ=255) aCTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAT‑ATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTAt > 2:1359004/1‑151 (MQ=255) ggCTTAAATGGCAACCTGAACCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCGGGTGGGTGATAAGGGCTAAAATGACCCCGCAGTGTGGGTGGGCCAAAAAT‑ATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCa < 2:920698/150‑1 (MQ=255) ggCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAT‑ATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCAt > 1:408712/1‑151 (MQ=255) ggCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAT‑ATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCAt > 1:28067/1‑151 (MQ=255) ggCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAT‑ATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCAt > 2:337377/1‑151 (MQ=255) gTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAT‑ATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCg > 2:767854/1‑150 (MQ=255) aaaggacacccTGAATCCAGCGGCCGGCGTCGGGGGATGGGGGGCGGGCAGCTGGGTAATAATGGCTAAAAGGACCCCGCGGTGATGGGGGGCCTAAAAT‑ATAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAc < 2:616872/141‑1 (MQ=255) aaaGGGCACCCTGAACCCAGCGGCCGGGGTCGGGGGATTGGGGGGGTGCAGTGGGGTAATAAGGGCTACAACGCCCCCGAGGTGATGGTGGGCCTAAAAT‑ATATGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAc < 1:566584/150‑1 (MQ=255) tGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAT‑ATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCgg > 2:975639/1‑151 (MQ=255) ggccacccTAAACCCAGCGGCCGGCGTCGGGGGAATGTGGGGGTGCGGGTGGGTTATAAGGGCTAAAATGACCCCGCGGTGTTGGTGGGCCAAAAAT‑AAAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCgg < 1:857503/145‑1 (MQ=255) aaCCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAT‑ATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTAc > 2:371566/1‑150 (MQ=255) cctaaacccgggggCCGGCGTGTGGGGATTGGGGGCGTGCGGTTGGGTAATAATGCCTACAATGCCCCCGAAGTGAGGGTGGCCCTAAAAT‑ATAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACa < 2:1309967/137‑1 (MQ=255) ccaaaacccggCGGCCAGCGTCTGGGGATGGTGGGCGGGCAGGTGGGTAATAAGGGCAAAAATGACCCCGCGGTGATGGTGGCACAAAAAT‑ATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACata < 2:169487/141‑1 (MQ=255) aTCCGGCGGGCAGCTTTGGGGGATTGTGGGGGGGCTGGTGGGTGATAGTGGCTAAAATACCCCCGGGGTGATGGTGGGCTAAAAAT‑ATAAGTTCTGTTCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATg < 1:1254041/148‑1 (MQ=255) ggCCGGCGTTTGGGGAAGGTGGGCGTGGAGCGGGGTGAAAATGGCTAAAATGACCCCGAGGTGAGGGTGGGCTAAAAAT‑AAAAGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAAtatt < 2:1266996/148‑1 (MQ=255) ggCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAT‑ATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATattatt > 2:755904/1‑151 (MQ=255) cAGCGTCTGGGGAAGGTGGGCGTGCAGCTGGGTGATAATGGCTAAAATACCCCCGAAGTGATGGGGGGCCTAAAAT‑ATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAttt < 1:664697/149‑1 (MQ=255) tCTGGGGATGGGGGGGGTGCAGTTGGGTAATAATGGCTAAAATGACCCCGAAGTGATGGTGGGCTTAAAAT‑ATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTg < 1:727928/148‑1 (MQ=255) ggggAATGTGGGGGTGCAGGTGGGTAAAAAGGTCTAAAATCACCCCGAAGTGAGGGTGGGCTAAAAAT‑AAAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATa < 1:781328/151‑1 (MQ=255) gggaAATGTGGGGGTGCGGGTGGGTAAAAAGGGATACAATGACCCCGAGGTGATGG‑GGGGATAAAATAATAAGTTCTGACTCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATa < 1:965674/151‑1 (MQ=255) gggaAATGTGGGCGTCCGGTGGGGTAAAAATGGCTACAATGACCCCGAGGTGTTGG‑GGTGCTGAAATAATAAGTTCGCATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAAtattat < 2:860811/139‑1 (MQ=255) ggAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAT‑ATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGc > 2:702472/1‑151 (MQ=255) aaTGGGGGCGGGCAGCTGGGTGATAATGGCTACAATGACCCCGAGGTGAGGGTGGGCTTAAAAT‑ATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCa < 1:1385560/150‑1 (MQ=255) gtgGGCGTGAAGCTGGGTGATAAGGACAACAATCCCCCCGTGGGGGTGGGGGGCCAAAAAT‑ATAATTTTTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGc < 1:1125790/150‑1 (MQ=255) gggCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAT‑ATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATg > 1:627929/1‑151 (MQ=255) aGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAT‑ATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGtcatc > 1:577528/1‑151 (MQ=255) gggTGAAAATGGCTACAAGCCCCCCGAAGGGAGGGGGGGCTGAAAAA‑AAAATTTCTGTTCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAAtt < 2:712698/151‑1 (MQ=255) ggTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAT‑ATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTg > 2:601349/1‑151 (MQ=255) aaTGGCAACAATGACCCCGCAGGGAGGGTGGGCTTAAAAT‑AAAAGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGc < 2:123742/146‑1 (MQ=255) aaaaTGACCCCGCAGTGATGGGGGCCCAAAAAT‑ATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAg < 2:1271717/149‑1 (MQ=255) aCAATGCCACCGTAGTGAGGGTGGCCCTAAAAT‑AAAGGTTCTAATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACa < 2:1169419/150‑1 (MQ=255) aaTGACCCCGTAGTGATGGGGGGCCTAAAAT‑ATATGTTCTGAGCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACa < 1:923063/148‑1 (MQ=255) aaTGACACCGCAGTGATGGTGGGCCTGAAAT‑ATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGc > 1:350513/1‑150 (MQ=255) ccccGCGGTGGTGGGGGGCGTAAAAT‑ATAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTg < 1:975353/148‑1 (MQ=255) tGATGGTGGGCCTGAAAT‑ATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCttt > 1:161701/1‑151 (MQ=255) ggtggGCTTAAAAT‑ATAGGTTCGGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAgcg < 2:747710/151‑1 (MQ=255) ggggggCCTAAAAT‑ATAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAgc < 2:443377/147‑1 (MQ=255) ggCCTGAAAT‑ATAATTTCGGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCaa < 1:635619/150‑1 (MQ=255) ccTGAAAT‑ATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAAcc > 1:80033/1‑151 (MQ=255) tGAAAT‑ATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTg > 1:929439/1‑151 (MQ=255) gAAAT‑ATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTg > 1:94977/1‑150 (MQ=255) | GTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAT‑ATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTG > NC_000913/1222083‑1222382 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |