Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,222,233 | 0 | G | A | 44.7% | 42.9 / 24.8 | 47 | pseudogene (60/75 nt) | ycgH | pseudogene; putative ATP‑binding component of a transport system |
Reads supporting (aligned to +/- strand): ref base G (23/3); new base A (0/21); total (23/24) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.81e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01 |
GTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAAC > NC_000913/1222083‑1222379 | gtgtGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACATGGCAGTGATGGTGGGCCTg < 1:592968/151‑1 (MQ=255) gAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAat > 2:466058/1‑151 (MQ=255) aaaGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGAAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAat > 2:856408/1‑150 (MQ=255) aaGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAatat > 1:410591/1‑151 (MQ=255) cAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGGGGGCCTGAAATATAAGTTCTg > 2:142601/1‑151 (MQ=255) aGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGAAAGCGTGTGGGGAATGTGGGGGGGCAGGTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTg > 2:30638/1‑149 (MQ=255) gATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATc > 2:710299/1‑151 (MQ=255) gATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATc < 1:191764/151‑1 (MQ=255) aGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGtt > 1:580487/1‑150 (MQ=255) tGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTaaa > 2:316963/1‑150 (MQ=255) aaaCTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCTTGTGGGGAATGGGGGCGTGCAGCGGGGTGATAATGGCTAAAATGCCACCGGAGGAATGGTGGGCGTGAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGcgtcg < 2:472141/150‑1 (MQ=255) aCTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTAt > 1:576457/1‑151 (MQ=255) ggggaaaaTGGCAACCTGATCCCGGCGCCTGCCGTCTGGGGAATGTGGGCGTGAAGTGGGGTGATAATGGCTAAAATCACCCCGGGTTGATGGTGGGCTTAAAATTTAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTc < 2:748685/144‑1 (MQ=255) ggCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGAGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGAGTCGTATTCAt > 1:322533/1‑151 (MQ=255) acccggaccccAGGGGCCGGCGTCGGGGGAAGGGGGGGGTCCGGCGGGGAGAAAGGGGCTAAATGGACCCCGGAGTGGGGGGGGCCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAcgcg < 2:37402/137‑1 (MQ=255) acccggAACCCAGCGGCCGGCGTGGGGGGAATGGGGGCGTGCAGGGGGGTGATAATGGCTACAATGACACCGCAGTGATGGGGGGCTAAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTa < 2:795285/144‑1 (MQ=255) ccTAAACCCGGCGGCCGGCGTCGGGGGATTGGGGGCGTGCGGCTGGGTAATAATGGCTAAAATGCCCCCGCGGTGAGGGGGGGCCTAAAATATAGTTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACat < 1:147241/150‑1 (MQ=255) gAACCCGGCGCCCAGCTTGTGGGAATTGTGGGCTTGCTGCTGGTTTATATTGCCTCCAATCCCCCCGCAGTGATGGTGGCCTTAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATg < 2:102412/150‑1 (MQ=255) ggCCGGCGTCGGGGGATTGTGGGCGTGCAGTGGGGTAATAATGGCAACAATGACCCCGCAGTGATGGGGGGCCTAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAatat < 2:599789/147‑1 (MQ=255) ccggggtggggggAGTGGGGGGGTGCGGGTGGGTAATATTGCCAACAATAACCCCGAGGTGAGGGTGGGCTTAAAATATAATTTCTTTCCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATattatt < 2:661392/140‑1 (MQ=255) tgTCGGGGGAATGTGGGGGTGCGGTGGGGTAAAAATGGTAAAAAAGACACAGAAGTGTTGGGGGGCCTAAAATAAAATTTCTGATTCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAtttt < 2:364954/146‑1 (MQ=255) ggggAATGTGGGCGTGCAGCGGGGTAATAATGGCTACAATGACACCGCAGTGATGGTGGGCATAAAATATAATTTCTGTTCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCg < 1:133982/148‑1 (MQ=255) ggggAATGGGGGCGTGCAGTTGGGTAAAAATGGCTACAAGGACCCCGCAGTGATGGTGGGCCTAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATa < 2:361017/151‑1 (MQ=255) ggggAATGGGGGCGTGCAGCTGGGTAATAATGGCTACAATGACCCCGCGGTGATGGTGGGCCAAAAATATAATTTTTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCg < 1:38285/148‑1 (MQ=255) ggAAGGGGGGCTGGCGGCTGGGTAATAATGGCTCCAACAACCCCGAAGTGTTGGGGGGCTTAAAAAATAATTTTTTCTCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGc < 1:751326/151‑1 (MQ=255) atggggggCGTGCATATGGGTAATAATGGCTAAAAGGACCCCGAAGTGATGGTGGGATAAAAATAAAAGTTCTAATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGc < 2:323236/144‑1 (MQ=255) ggggggGTGCAGCTGGGTAAAAAGGGCTACAATCCCCCCTCGGTGATGTTGGCCTTAAAATATAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCAc < 2:307838/146‑1 (MQ=255) ggggGTGCGGTTGGGTAATAATGGCTACAATGACCCCGCAGTGAGGGTGGGCTTAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATg < 2:844577/151‑1 (MQ=255) ggCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATgg > 2:491591/1‑151 (MQ=255) tGCAGCTGGGTAATAAGGTCTACAACGCCCCCGCGGTGATGTTGGCCTTAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGgtg < 1:425944/149‑1 (MQ=255) ggttGATAATGGCTAAAATGCCCCCGCGGTGAGGGTGGGCCTAAAATATAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAAt < 1:78960/147‑1 (MQ=255) gggTGATAAGGGCTAAAATGACCCCGCAGTGAGGGGGGGCCTAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAAtt < 2:69135/151‑1 (MQ=255) ggTGATAATGCCTAAAATGCCACCGCAGTGAGGGGGGGCTTAAAATATAATTTCTGCCCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCaa < 2:149679/148‑1 (MQ=255) gATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCg > 1:483328/1‑151 (MQ=255) tAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTg > 2:24029/1‑151 (MQ=255) aaTGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTg > 2:68649/1‑150 (MQ=255) aTGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGAt > 1:453451/1‑151 (MQ=255) gCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCa > 1:864902/1‑150 (MQ=255) aCAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAAc > 2:71866/1‑149 (MQ=255) tGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCAtt > 1:57593/1‑151 (MQ=255) acacCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTg > 1:326963/1‑151 (MQ=255) acCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCg > 1:844191/1‑151 (MQ=255) cggaaGTGATGGTGGCCCAAAAATATAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCAttt < 1:764066/143‑1 (MQ=255) ccGCGGGGATGGGGGGCTTAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCg < 1:36590/150‑1 (MQ=255) cGCAGGGAGGGTGGGCATAAAATAAAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCggg < 2:630437/151‑1 (MQ=255) cAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTc > 1:72703/1‑151 (MQ=255) cAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTc > 2:357119/1‑151 (MQ=255) gCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAAc > 2:511131/1‑151 (MQ=255) | GTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAAC > NC_000913/1222083‑1222379 |
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: ‑ |