Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,222,233 | 0 | G | A | 41.7% | 77.2 / 25.2 | 60 | pseudogene (60/75 nt) | ycgH | pseudogene; putative ATP‑binding component of a transport system |
Reads supporting (aligned to +/- strand): ref base G (28/7); new base A (0/25); total (28/32) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.75e-11 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTG > NC_000913/1222085‑1222382 | gtgAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑aa < 1:1156393/151‑1 (MQ=255) aaGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAatat < 2:1076070/151‑1 (MQ=255) gAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGtt > 1:702349/1‑151 (MQ=255) gAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGtt > 2:302317/1‑151 (MQ=255) gAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGtt < 1:934625/151‑1 (MQ=255) aaGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGa < 1:1224208/151‑1 (MQ=255) aGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGAt > 2:119175/1‑151 (MQ=255) ggCGCGTAATTTGGGGGAAATTAAACCTGGCGTAAATGGCACCCTGAATCCAGCGGCCAGCGTCTGGGGAATGTTAGCGCGCAGCTGGGTGATAATGGCCACAATGACACCGCAGCGATGGTGGGCCTG‑ACACATAAGTTCCGATcc > 2:498121/1‑147 (MQ=255) cGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATCCCGCCGt > 2:74943/1‑151 (MQ=255) gAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATCCCGCCGtt > 1:778334/1‑151 (MQ=255) gAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATCCCGCCGtt > 2:1034022/1‑151 (MQ=255) aGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATCCCGCCGtt > 1:488104/1‑150 (MQ=255) aaTTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATCCCGCCGTTAGCt > 1:1159521/1‑151 (MQ=255) ttGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATCCCGCCGTTAGCTaa > 1:305277/1‑150 (MQ=255) ttGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATCCCGCCGTTAGCTaa > 2:1313623/1‑150 (MQ=255) ttGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATCCCGCCGTTAGCTaa > 2:998566/1‑150 (MQ=255) aaaTTAAACTGGGGGAAAATGGCAACCTGAACCCAGCGGCCGGCGTCGGGGGAATGTGGGGTTGCAGCTGGGTGATAATGGCTAAAATGACCCCGAAGTGATGGTGGGCTTA‑AAATATAAGTTCTGACCCCGCCGTTAGCTAAAAAACcg < 2:697339/150‑1 (MQ=255) aaTTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACcgcg > 1:496168/1‑151 (MQ=255) tAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGcgtcg > 1:518490/1‑151 (MQ=255) aCTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTAt > 1:333066/1‑151 (MQ=255) cTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTAtt > 1:544359/1‑151 (MQ=255) ggCGAAAAGGGCACCCTGAACCCAGCGGCCGGCGTGGGGGGAAGGGGGGGGTGCAGTGGGGTGATAATGGCTAAAATGACCCCGCAGTGATGGGGGGCCTA‑AAATATAAGTTTTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTa < 1:759172/146‑1 (MQ=255) ggCGAAAAGGCCACCCTAAACCCAGCGGCCAGGGTGGGGGGAAGGGGGGCGGGCAGGGGGGTAATAAGGGCTAAAATGACCCCGCGGTGAGGGGGGGCTTA‑AAATAAAAGTTTTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCAt < 1:556239/151‑1 (MQ=255) gCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATc > 2:1232023/1‑151 (MQ=255) aaggccacccTGAACCCAGCGGCCGGCGTCGGGGGAATGGGGGGGTGAGGCGGGGTGATAATGGCTAAAATGACCCCGCGGTGAGGGGGGCCCTA‑AAATATAAGTTCTGCTCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAcgc < 1:917309/143‑1 (MQ=255) aaTGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAcg > 2:1206302/1‑150 (MQ=255) cAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTAc > 2:851977/1‑151 (MQ=255) acccaaaaTACAGGGGCCAGCGTGTGGGGATTGGGGGGGTGCAGGTGGTTAAAAATGTCTACAACGACCCCCCAGTGATGGTGGGCCTA‑AAATAAATGTTCTCATTCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTa < 2:211240/143‑1 (MQ=255) aaCCTGAACCCAGGGCCGGGCGTGGGGGGAATGGGGGCGTGGGGCTGGGTGATAAGGGCTAAAATGACCCCGCGGTGAGGGTGGGCCTA‑AAATATATGTTCTAATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTa < 2:39209/149‑1 (MQ=255) aaCCTGAACCCAGCGGCCAGGGTCTGGGAATTGTGGGGGTGCAGCTGGTTGATAAGGGCTACAAGGACCCCGCGGTGGTGGGGGGCCTA‑AAATATAGTTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTa < 1:613136/149‑1 (MQ=255) ccTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACat > 1:928169/1‑150 (MQ=255) ccTAAACCCAGCGCCCGGCGTGTGGGGAATGGGGGCGTGCAGCTGGGTGATAATGGCTACAATGCCCCCGCAGTGATGGGGGGCTTA‑AAATAAAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACata < 1:345268/151‑1 (MQ=255) ctaaatcaggggcccAGCCTCGGGGGTTTGTGGGGGTGCAGTGGGGGTAAAATGGGCAAAATGCCCCCGCATGGTGGGGGGGCCTA‑AAATAAAAGTTCATACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACata < 2:940798/137‑1 (MQ=255) gAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATg > 2:934152/1‑150 (MQ=255) aaccaaGCGGCCGGCTTCGGGGGAATGGGGGCGGGCAGCGGGGTAATAATGCTAACAACAACCCGGCAGTGATGGGGGGCCGA‑AAATAAATTTTCTAACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGaa < 2:378391/146‑1 (MQ=255) tCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAat > 2:867751/1‑151 (MQ=255) gCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGGTGGGTGATAAGGGCTACAATGACCCCGCAGTGATGGTGGGCCTA‑AAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAAtatta < 1:99888/151‑1 (MQ=255) ggCCAGCGTCTGGGGATGGGGGGTGGGCGGTTGGGTGAAAATGGCTAAAAGGACCCGGCAGTGGTGGGGGGCTGA‑AAATATAATTTTTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATTGACGCGGTACATATGAAAtatta < 1:342567/149‑1 (MQ=255) ggCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATattatt > 1:226249/1‑151 (MQ=255) gCCAGCGTCTGGGGAATGGGGGCGTGCAGCGGGGTGATAATGGCTAAAATGACCCCGCAGTGAGGGTGGGCTTA‑AAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAAtattat < 1:912598/149‑1 (MQ=255) ccAGCTTCTGGGGATGGTGGGCGTGCAGCGGGGTAATAATGGTTACAATCCCCCCCCAGTGAGGGGGGGCCTA‑AAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAttt < 2:408732/150‑1 (MQ=255) cAGCGTGTGGGTAATGGGGGCGTGCAGCGGGGTGATAATGGCTACAAGGACCCCTCGGTGATGGTGGGCCTA‑AAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAtttt < 2:226245/150‑1 (MQ=255) gCGTCTGGGGAATGTGGGCGTGCAGCTGGGTAATAATGGCTACAATGACCCCGCAGTGATGGGGGGCCTA‑AAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTg < 1:1322474/151‑1 (MQ=255) tgtGGGGAGGGTGGGAGTGAGGGGGGGTAATATTGGCAAAATTAACACCCGAGGGATGGTGGGACTA‑AAAAATAGTTATTTACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAGTATTATTTTTTg < 1:1213144/146‑1 (MQ=255) ggggAATGGGGGGGTGCGGCGGGTTAATAATGTCTACAATGCCCCCGTAGGGATGG‑GGGCCTGAAAATATAAGTTCTGTCCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATa < 2:837794/151‑1 (MQ=255) ggggGTGCAGCGGGGTTATAATGCCTAAAATGACCCCGGAGTTGTGGTGGCTATG‑AAATATAATTTTTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCAc < 2:197919/146‑1 (MQ=255) gggcgggcAGGGGGGTAAAAATGGCAACAAAGCCCCCGAGGTGAGGGGGGGCAAA‑AATTATAAGTTCGGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATg < 2:111538/151‑1 (MQ=255) gTGGAGCGGGGTGATAATGGTAAAAATGACCCCGGGGTGTGGGTGGGCGTA‑AAAAATAGTTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGgt < 2:517840/149‑1 (MQ=255) ggCTGGGTGATAAGGGCTAAAAGGCCCCCGCGGTGAGGGGGGGCATA‑AAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGtca < 1:74943/148‑1 (MQ=255) aaaaagccTACATTGCCCCCGCAGTGTGGGGGGGCTTA‑AAAAATAGTTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGt < 1:1018052/144‑1 (MQ=255) tAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTg > 1:272596/1‑151 (MQ=255) aaTGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTg > 1:747197/1‑150 (MQ=255) ggCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATc > 1:1032391/1‑150 (MQ=255) aCAATTACCCCCAAGTTATGGTGGGCCAA‑AAATATAATTTTTAATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCa < 1:596674/147‑1 (MQ=255) cAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGc > 1:1338223/1‑151 (MQ=255) gagggggggCTTA‑AAATATAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGt < 1:451503/140‑1 (MQ=255) ggtggtCCTG‑AAATATATGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCtt < 2:1111522/146‑1 (MQ=255) ggtggGCATA‑AAAAATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAg < 1:302323/149‑1 (MQ=255) ggggggCTTA‑AAATATAAGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAgc < 1:119177/147‑1 (MQ=255) tG‑AAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTg > 2:1334169/1‑151 (MQ=255) | GTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTG‑AAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTG > NC_000913/1222085‑1222382 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |