| 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.4% | 55.7 / 30.8 | 57 | pseudogene (60/75 nt) | ycgH | pseudogene; putative ATP‑binding component of a transport system |
| Reads supporting (aligned to +/- strand): ref base G (25/5); new base A (0/27); total (25/32) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.32e-11 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTG > NC_000913/1222087‑1222382 | gAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCC‑TGAAat < 1:682457/151‑1 (MQ=255) gTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCC‑TGAAATATaa < 1:306782/151‑1 (MQ=255) cGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCC‑TGAAATATAAGt > 2:546567/1‑151 (MQ=255) aaGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCC‑TGAAATATAAGTTCTg > 1:479518/1‑150 (MQ=255) aTGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATcc > 2:657838/1‑151 (MQ=255) ggggCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATccc < 2:229027/150‑1 (MQ=255) gCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCccgccg > 1:140034/1‑151 (MQ=255) cGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGt > 2:798606/1‑151 (MQ=255) tAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGc > 1:985944/1‑151 (MQ=255) tAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGc > 1:199434/1‑151 (MQ=255) aTTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCt > 1:475582/1‑150 (MQ=255) ggaaaaTTAAACTGGGCGAAAATGGCAACCTGAATCCGGCGCCCCGCGTCGGGGAAATGGGGGGGGGCAGCGGGGGAAAAATGCCAACAATGACCCCGCAG‑GAAGGGGGGCCATAAAATATAGGTTCTGACCCCGCCGTTAGCTAAAAAAc < 2:536869/148‑1 (MQ=255) aaaTTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACcgc > 1:743412/1‑151 (MQ=255) aaTTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACcgcg > 2:593455/1‑151 (MQ=255) ttAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGcgtc > 2:695306/1‑151 (MQ=255) taaaaCTGGGGGTAAATGGCACCCTGAATCCGGCGGCCAGCGTCGGGGGAATGGGGGCGTGCGGCTGGGTGATAAGGCCAAAAATGACCCCGCAGTGTGGGTGGCCT‑TAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGcgt < 1:144321/148‑1 (MQ=255) aaaTGGCACCCTGAACCCGGCGGCCAGCGTCGGGGAAAGGGGGGGGGGCGGCGGGTTGATAATGGCTACAATCACCCCGAAGGGAGGGTGGCCT‑TAAAATATAGTTTTTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAcg < 2:89507/151‑1 (MQ=255) aaaTGGCAACCTGAATCCAGCGCCCGGCGTCGGGGGAATGTGGGCGTGCAGCTGGGTAAAAATGGCTACAATGCCCCCGCAGTGATGGGGGGCC‑TAAAATATAAGTTCTGCCCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAc < 1:317162/150‑1 (MQ=255) agtGGCACCCAGAACCGCGCGGCCGTCGTCGGGGGAATGGGGGCGTGCGGCGGGGGAATAGGGGCTAAAAGACCACCGCAATGAGGGGGGGCA‑TAAAATATACTTTTGGTCCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATc < 2:948262/144‑1 (MQ=255) aagggcaccttaaacccAGCGCCCGGGGTCGGGGAAATGGGGGGGTGCAGGTGGGTGAAATGGGCAACAATGCCCCCGGGGTGTTGGGGGGCC‑AAAAATAAAATTTCTACTCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAcgc < 1:317916/136‑1 (MQ=255) aaTGGCACCCTAAACCCGGCGGCCGGCGTGGGGGGAATGGGGGCGTGCGGGTGGGTGAAAATGCCAACAATGCCCCCGCAGTGAGGGTGGGCC‑TAAAAAATAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAcgc < 1:269349/151‑1 (MQ=255) acccTAAATCCAGCGGCCAGCGTCGGGGGAATGTGGGGGTGCAGCTGGGTGATAATGCCAAAAATGCCCCCGCGGGGGTGGTGGCCC‑TAAAATAAAAGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACa < 1:572340/149‑1 (MQ=255) acccTAAACCCGGCGGCGGGCGTTGGGGAAATGGGGGCGTGCAGCGGGGTAATAATGCCTAAAATGCCCCCGCAGTGAGGGGGGCCT‑TAAAATATAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTa < 2:140938/147‑1 (MQ=255) aaCCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTAc > 2:379981/1‑150 (MQ=255) cccTGAACCCAGCGGCCAGCGTCGGGGGAATGGGGGCGTGCAGGTGGGGGATAAGGGCTACAATAACCCCGCAGTGATGGTGGGCC‑TGAAATATAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACat < 2:100537/150‑1 (MQ=255) ccTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACat > 1:211590/1‑150 (MQ=255) ccTGAACCCAGCGGCCGGCTTTGGGGAAATGGGGGCGTGCAGCGGGGTGATAATGGCTAAAATGACACCGCAGTGATGGGGGGCT‑TAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACata < 1:704513/151‑1 (MQ=255) tGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATg > 2:130396/1‑151 (MQ=255) aGCGGCCAGGGTCGGGGGATTGGGGGCGTGCAGCTGGGTGATATTGGTTAAAATACCCCCGCATTGTTGGTGGCCT‑TAAAATATAAGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATg < 1:1040338/144‑1 (MQ=255) ggCCCGCGTGGGGGGAGGGGGGGGGGGCGGTTGGGTGATAATGGCTAAAATGCCCCCGATGTGGGGGTGGGCC‑TAAAATATAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAAtattat < 2:493866/150‑1 (MQ=255) ggCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATattatt > 2:1015503/1‑151 (MQ=255) cAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAttttt > 1:872726/1‑151 (MQ=255) cAGCGTAGGGAGACTGGGGGGTTGCAGAGGGGTTAAAATGGCAAAAAAGACCCCGCAATGTTGGTGGTCC‑AAAAATATAATTGATAATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAtttt < 1:488698/150‑1 (MQ=255) aGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAtttttt > 1:375517/1‑151 (MQ=255) gggTCTGGGAAAGGGGGGCGTGAAGGTGGGTGATAATGGCTAAAATGACCCCGCAGTGAGGGTGGGCC‑TAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAtttttt < 1:586654/148‑1 (MQ=255) gCGTCGGGGGAGTGTGGGCGTGCAGCTGGGTAATAATGGCTACAATGCCCCCGCAGTGATGGGGGGCT‑TAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAtttttt < 2:21438/150‑1 (MQ=255) ggggAATGGGGGGGTGCGGCTGGGTAATAATGGCTACAATACCCCCGCAGTGATGGTGGGCC‑TAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCg < 2:494378/148‑1 (MQ=255) gggaAATGTGGGCGTGCAGTTGGGTGATAATGGCTAAAATGACACCGCAGTGATGGTGGGCT‑TAAAATATAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGa < 1:453594/149‑1 (MQ=255) ggggCGTGCAGGGGGGTAAAAATGGCTACAATGCCCCCGCAGTGATGGGGGGCC‑TAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCAt < 1:1092967/150‑1 (MQ=255) cGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGgtg > 1:12081/1‑151 (MQ=255) aGCGGGGTGAAAATGCCTAAAATAACCCCGTGGGGATGGTGGGCC‑TAAAATATATTTTTTGTTCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGtca < 2:325792/149‑1 (MQ=255) aTGGCTACAATGACACCGCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGAt > 2:1057135/1‑151 (MQ=255) ggCTACAAGCACCCCGCAGTGTGGGTGGGCC‑AGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCa < 1:798725/151‑1 (MQ=255) cAATGACACCGCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGc > 1:983771/1‑151 (MQ=255) aaTGACACCGCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGc > 2:725552/1‑150 (MQ=255) ccgcggtgtgtgggtggGCT‑AAAAAAATAATTTCGGCTCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGc < 2:850537/139‑1 (MQ=255) acCGCAGTGATGGTGGCCT‑AAAAAAATAATTTGTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCAttt < 2:978621/148‑1 (MQ=255) ccGCGGGGATGGTGGGCA‑TAAAATATATGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCgg < 2:65501/151‑1 (MQ=255) ccGCAGTGATGGTGGGCC‑TAAAATATAAGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCgg < 1:279247/151‑1 (MQ=255) ccGCAGTGAGGGTGGGCA‑TAAAATATAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTg < 2:1028954/148‑1 (MQ=255) cAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTc > 1:788935/1‑151 (MQ=255) ggtggCCA‑TAAAATAAAAGTCCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTc < 1:643901/144‑1 (MQ=255) gggCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCa > 2:809124/1‑150 (MQ=255) gggCA‑TAAAATATAGGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCaa < 2:140007/151‑1 (MQ=255) ggCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCaa > 2:330569/1‑150 (MQ=255) ggCC‑TAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCaaa < 1:237721/151‑1 (MQ=255) cc‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAAcc > 2:1053726/1‑151 (MQ=255) gAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTg > 1:1137667/1‑150 (MQ=255) | GAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTG > NC_000913/1222087‑1222382 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |