| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,222,233 | 0 | G | A | 46.2% | 34.3 / 19.6 | 39 | pseudogene (60/75 nt) | ycgH | pseudogene; putative ATP‑binding component of a transport system |
| Reads supporting (aligned to +/- strand): ref base G (16/5); new base A (0/18); total (16/23) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 9.98e-07 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
AAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGG‑TGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAG > NC_000913/1222088‑1222372 | aaaGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGG‑TGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAat > 1:173652/1‑150 (MQ=255) aaGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGG‑TGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAatat > 2:805515/1‑151 (MQ=255) aaGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGG‑TGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAatat < 2:450411/151‑1 (MQ=255) aGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGG‑TGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAatat > 2:253709/1‑150 (MQ=255) gTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGG‑TGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATaa < 2:217521/151‑1 (MQ=255) 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) | AAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGG‑TGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAG > NC_000913/1222088‑1222372 |
| 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: ‑ |