Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,222,242 | 0 | G | T | 44.1% | 49.0 / 10.1 | 34 | pseudogene (69/75 nt) | ycgH | pseudogene; putative ATP‑binding component of a transport system |
Reads supporting (aligned to +/- strand): ref base G (13/6); new base T (0/15); total (13/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.87e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
AAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGG‑TGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAG > NC_000913/1222098‑1222372 | aaGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGG‑TGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTg < 2:705997/150‑1 (MQ=255) cGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGG‑TGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGt > 1:613441/1‑151 (MQ=255) aaCGTGGGGTAAAGGGCAACCTGAACCCAGCGGCCGGCGTGGGGGGAAGGGGGGCGTGCAGCGGGG‑TGATAAGGGCTACAACGCCCCCGAGGTGATGGGGGGCTTAAAATATAATTCCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTa < 1:404927/151‑1 (MQ=255) ttGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGG‑TGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTc > 2:336205/1‑151 (MQ=255) gCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGG‑TGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATc > 2:169887/1‑151 (MQ=255) cGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGG‑TGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCg > 1:829055/1‑151 (MQ=255) gTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGG‑TGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCg > 2:614240/1‑150 (MQ=255) ccTTAACCCAGCGGCCGGGGTCGGGGGAATGGGGGCGTCGGGTTTGG‑TGATAATGCATACAAGGACCCCGTAGTGAGGGGGGGCTGAAAATAAATGTTCTGGACCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACat < 1:682712/150‑1 (MQ=255) gAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGG‑TGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATg > 2:641488/1‑150 (MQ=255) aaTCCAGCGGCCAGCGTCTGGGGAATGTGGGCTCGCAGCTGGG‑TGATAATGGCTACAATGACACCGGAGAGATGGCGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGAGTCGTATTCATCGACGCGGTACatat > 1:548611/1‑148 (MQ=255) cccggCGGCCAGGGTCTGGGAATGGGGGGCGGCCGGCGGGG‑TAAAAATGGCTAAAAGAACCCCGAGGTGAGGGTGGGCCAAAAATATAATTTCTAATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAat < 1:429075/147‑1 (MQ=255) ccagcaccgagctgcTGGGTGATGTGGGCGTG‑AACTGGGATAATAATGCCTAAAATGACGCGCGGGTAATGGTGAGCATAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGaaa < 1:592717/135‑1 (MQ=255) ccAGCGGCCGGCGTGGGGGAAATGTGGGGGGGCAGCTGGG‑TAATAATGGCTACAATGACCCCGCAGTGAGGGGGGGCCTAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGaa < 2:537154/148‑1 (MQ=255) ggggCCGGCGTCGGGGGAATGGGGGGGTGCAGTGGGG‑AGATAATGTCAAAAATGACACCGCAGTAATGGGGGGCATAAAATATAATTTCTTACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGaaa < 1:764951/144‑1 (MQ=255) ggCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGG‑TGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATattatt > 2:76582/1‑151 (MQ=255) ggCCAGCGTCGGGGGAATGTGGGCGGGCAGCTGGG‑TAATAATGGCTAAAATGACACCGGAGTGAGGGTGGGCGTAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAAtatta < 2:222230/149‑1 (MQ=255) ggCCAGCGTCGGGGGAATGGGGGCGTGCAGGGGGG‑TGATAATGCCTAAAATGCCCCCGCAGTGATGGGGGGCCTAAAATATAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATattatt < 2:357753/151‑1 (MQ=255) gcccAGCGTGTGGGGAATGTGGGGGTGCAGCGGGG‑TAATAATGGCTACAATGACACCGCAGTGATGGTGGGCTTAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAAtatt < 1:310664/146‑1 (MQ=255) gTCTGGGGAATGGGGGCGTCCAGGGGGG‑AGAAAATGGCTACAATGCCCCCGAGGTGATGGTGGCCTTAAAAAATAGTTTCTGCCCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGc < 2:173667/150‑1 (MQ=255) tgtgGGCGTGCAGCTGGT‑TGATAATGACAACAAGACCACCGTAGTGGTGGTGGGCTTAAAATATATGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACg < 2:687932/150‑1 (MQ=255) ggggCGTGCAGCTGGG‑TAATAATGGCTACAATGACCCCGGAGTGATGGGGGGCCTAAAATAAAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCAt < 2:502830/150‑1 (MQ=255) gggCGTGCGGCGGGG‑TAATAATGGCTAAAAAGACCCCGGAGGGATGGTGGGCCTAAAATATAATTTTTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCAt < 2:568288/150‑1 (MQ=255) gggCGTGCAGCTGGG‑TAATAATGGCAACAATGACACCGCAGTGATGGTGGGCTTAAAATATAAGTTCTAATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATg < 1:93028/151‑1 (MQ=255) tGGG‑TGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAAt > 1:383372/1‑151 (MQ=255) gg‑TGATAAGGGCTACAATGCCCCCGCAGTGATGGGGGGCCTAAAATAAAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTg < 2:646882/151‑1 (MQ=255) ggCTACAATGACCCCGGAGTGATGGTGGCCCTAAAATATAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGAt < 2:430351/149‑1 (MQ=255) aCAATGACGCGGGAGTGATGTGGCGCCGAAAATATAATTTCTATCCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAg < 1:731683/151‑1 (MQ=255) cggagggggtgtggTCTTGAAATATAATTTATAATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACa < 1:720179/130‑1 (MQ=255) cGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCggg > 1:103048/1‑151 (MQ=255) aGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCt > 1:395750/1‑151 (MQ=255) tGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCttt > 2:776843/1‑151 (MQ=255) gATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCttt > 1:741671/1‑150 (MQ=255) tggtggGCCTGAAATATAAGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAg < 2:385625/150‑1 (MQ=255) ggggggCCAAAAATAAAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCt < 2:127863/142‑1 (MQ=255) | AAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGG‑TGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAG > NC_000913/1222098‑1222372 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |