Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,222,233 | 0 | G | A | 43.8% | 53.1 / 20.1 | 48 | pseudogene (60/75 nt) | ycgH | pseudogene; putative ATP‑binding component of a transport system |
Reads supporting (aligned to +/- strand): ref base G (21/6); new base A (0/21); total (21/27) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.66e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01 |
GTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGC‑GTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGC > NC_000913/1222083‑1222373 | gtgtGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGC‑GTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCC‑Tg < 2:911847/151‑1 (MQ=255) gtgtGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGC‑GTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCC‑Tg < 1:685004/151‑1 (MQ=255) aaaGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGC‑GTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCC‑TGAAat > 2:253237/1‑150 (MQ=255) gAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGC‑GTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCC‑TGAAATATAAGtt < 2:59281/151‑1 (MQ=255) aCAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGC‑GTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCC‑TGAAATATAAGTTCt > 2:962912/1‑151 (MQ=255) gATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGC‑GTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATc < 1:158056/151‑1 (MQ=255) aTGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGC‑GTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATcc < 2:199722/151‑1 (MQ=255) tGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGC‑GTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATccc > 1:320629/1‑151 (MQ=255) gaaaTTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGC‑GTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATTccgccgat > 2:748534/3‑147 (MQ=255) tAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATGCAGCGGGCAGCGTCGGGGGACTGGGGGC‑GTGCAGCTGGTAGCTAATGGGTACAATGACCCAGCATTGATGGGGGGCC‑TGAGATATAAGTTCTGATCACGCCGTTa > 2:33884/1‑149 (MQ=255) aaTTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGC‑GTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCt > 2:961189/1‑151 (MQ=255) aaTTTGGGGAAAATTAAACTTGGCGAAAATGGCAACCTGAACCCAGCGGCCAGCGTCTGGGGAATGTGGGC‑GGGCAGCTGGGTGATAATGGCTACAATGCCCCCGCATGGTTGTGGGGCC‑TAAAATATAATTTCTGACCCCGCCGTTAGCt < 2:808672/151‑1 (MQ=255) ggggAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGC‑GTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTaaaaa > 2:41464/1‑150 (MQ=255) ggggAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGC‑GTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTaaaaa > 2:376013/1‑150 (MQ=255) aaTTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGC‑GTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACcgcg > 1:550132/1‑151 (MQ=255) ttAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCGGGGGAATGTGGGC‑GTGCAGCTGGGTGATAATGGCTACAATGACACCGAAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGcgt < 1:715443/150‑1 (MQ=255) aaaCTTGGCGTAAAGGGCACCCTGAATCCAGCGGCCAGCGTCGGGGGAATGTGGGC‑GTGCAGCGGGGTGATAATGGCTACAATGCCCCCGCAGTGATGGTGGGCC‑TAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGcgtcg < 1:578564/150‑1 (MQ=255) ggCGTAAAGGGCACCCGAAACCCGGCGGCCAGCGTCTGGGGAATGTGGGC‑GTGCAGCTGGGTGATAATGGCTACAATGCCCCCGCAGTGATGGGGGGCC‑TGAAATATAATTTCTAACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTc < 1:342204/149‑1 (MQ=255) cGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGC‑GTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCg > 2:815382/1‑151 (MQ=255) cGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGC‑GTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCg > 2:304629/1‑151 (MQ=255) aaggggCACCCGGAATCCAGCGCCCAGCGTCTGGGGAAGGGGGGG‑TTGCAGCGGGTTAATAATGTCTACAACGCCCCCGAAGTGAT‑GTGGGCCATAAAAAATAAGTCCTTATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAcg < 2:500611/147‑1 (MQ=255) aaagggcccccggAACCCGGCGGCATGCGTCGGGGGAGTGGGGGG‑GTCCGGGTGGGTGATATGGGTAACAATGACCCCGCAGTGATGGGGGGCC‑GAAAAAATAGTTGCGGTCCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAc < 1:787416/138‑1 (MQ=255) gactCCAGCGGCAAGTGCGGGGGGATGGGGGGG‑GTGCGGAGGGGTGAAAATGTCAACAATAACCCGGCAGTGTTGGTGGGCT‑TAAAAAAGAATTTCTTTCCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACatat < 1:521607/146‑1 (MQ=255) aacccggCGGCCGGCGTCGGGGGAATGGGGGC‑GTGCGGCGGGGTGAAAATGGCTAAAATGACCCCGCAGTGAGGGTGGGCC‑TAAAATAAAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGaa < 1:649054/145‑1 (MQ=255) aacccAGCGCCCGGCGTCTGGGGAATGGGGGC‑GTGCAGTTGGGTGAAAAGGGCTAAAATGACCCCGCGGTGATGGTGGGCT‑TAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACatat < 1:1341/145‑1 (MQ=255) aaTCCAGCGGCCAGCGTCTGGGGAATGTGGGC‑GTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGa > 2:886716/1‑150 (MQ=255) acccAGCGGCCGGCGTCGGGGGATTGTGGGC‑GTGCAGCTGGGTGATAATGGCTAAAATGCCCCCGCAGTGAGGGGGGGCC‑TAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATg < 1:385031/146‑1 (MQ=255) ccAGCGGCGGGCGTCGGGGGATTGTGGTG‑GGGCAGGGGGGTTATAGGGGTTAAAATGACGCCGTGGTAGTGGGGGGCT‑TAAAATATAAGTTCAGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATg < 2:522276/146‑1 (MQ=255) ccAGCGGCCGGCTTCGGGGGATTGGGGGC‑GTGCAGCTGGGTGATAATGGCAAAAATGCCCCCGCAGTGATGGGGGGCC‑TAAAAAATAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAata < 1:732414/151‑1 (MQ=255) ccAGCGGCCAGCGTCTGGGAAAGGTGGGC‑GTGCAGCTGGGTGATAATGGCTACAATGACCCCGCAGTGATGGTGGCCC‑TAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAata < 1:813075/151‑1 (MQ=255) ggCCAGCGTCTGGGGAAGGGGGGC‑GTGCGGCTGGGTAATAATGCCTACAAGGACCCCGAGGTGATGGTGGGCT‑TAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGa < 1:530181/142‑1 (MQ=255) gCCAGCTTCTGGGGTTTGTGGGC‑TTGGGGCTGGGATATAATGGCAACAATGACCCCGCAGGGGTGGTGGGCC‑GAAAATAAGCGTTTTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATg < 1:691597/140‑1 (MQ=255) ccggCGTCTGGGAAATGTGGGC‑GTGCAGCTGGGTAATAAGGGCTAAAATGACACCGCAGTGATGGTGGGCC‑TAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAtttt < 1:931640/148‑1 (MQ=255) ccAGCGTCGGGGGAATGTGGGC‑GTGCAGCTGGGTGATAATGGCTAAAATGACCCCGCAGTGATGGTGGGCC‑TAAAATATAAGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATattatt < 1:664901/149‑1 (MQ=255) ccggcgtgggagAATGGGGGCTGTGCGGATGGGTAATATGGGCTAAAATACCCCAGCAGGGATGGGGGGCG‑TAAAAAATAAGTTCTGCTCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAttt < 1:908821/139‑1 (MQ=255) cAGCGTCTGGGGAATGGGGGG‑GTGCAGGGGGGTAATAAGGGCTAAAATGACCCCGCAGTGGTGGGGGGCC‑TAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAttttt < 1:226299/151‑1 (MQ=255) gggaAAGGTGGGG‑GTGAGGCTGGGTAATAATGGTTACAATGACCCGGGAGTGGGGGTGGGCA‑AGAAATAGAAGTTCTTACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAtttttt < 1:950425/144‑1 (MQ=255) gC‑GTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATgg > 2:352594/1‑150 (MQ=255) gTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGgtgt > 2:748516/1‑151 (MQ=255) gCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGtc > 2:203011/1‑150 (MQ=255) gggTGATAAGGCCAACATTGACCCCGCAGGGATGGGGGGCT‑TAAAATATAAGTTTGGCTCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAAtt < 1:983308/151‑1 (MQ=255) gggTAATAATGGCTACAATGACCCCGCAGTGATGGGGGGCC‑TAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGtcatca < 1:694093/148‑1 (MQ=255) tACAATGACACCGCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAAc > 1:49420/1‑150 (MQ=255) cAATGACACCGCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGc > 1:996382/1‑151 (MQ=255) acCGCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCg > 1:204532/1‑151 (MQ=255) tGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCttt > 2:232640/1‑151 (MQ=255) gATGGAGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCttt > 1:133669/1‑150 (MQ=255) tggtggGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAgc > 1:929835/1‑151 (MQ=255) gggtggGCC‑TAAAATATATTTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCttt < 1:255178/147‑1 (MQ=255) ggtggGCC‑TAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAg < 1:377438/149‑1 (MQ=255) | GTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGC‑GTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGC > NC_000913/1222083‑1222373 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |