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% | 32.2 / 18.1 | 38 | pseudogene (60/75 nt) | ycgH | pseudogene; putative ATP‑binding component of a transport system |
Reads supporting (aligned to +/- strand): ref base G (14/6); new base A (1/17); total (15/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.01e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01 |
TGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAA‑TATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGC > NC_000913/1222084‑1222373 | tgtgAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTg > 2:32604/1‑150 (MQ=255) gtgAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGa < 2:903616/150‑1 (MQ=255) cGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAA‑TATAAGt > 2:830791/1‑151 (MQ=255) gAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAA‑TATAAGtt < 2:38992/151‑1 (MQ=255) aCAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAA‑TATAAGTTc < 1:267562/150‑1 (MQ=255) cAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAA‑TATAAGTTCTg < 1:679501/151‑1 (MQ=255) aGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAA‑TATAAGTTCTGAt < 1:347702/151‑1 (MQ=255) tGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAA‑TATAAGTTCTGATccc > 1:563086/1‑151 (MQ=255) gggggAAATTAACCTGGGCGTAAGTGGCACCCTGAACCCAGCGGCCAGCGTTTGGGGAAGGTGGGCGTGCAGCTGGGTGATAAGGGCTACAATGACACCGAAGTGATGGTGGGCATGAAA‑TATAAGTTCTGATCCCGCCGTTAGCTaaa < 2:55664/149‑1 (MQ=255) aaaCTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAA‑TATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGcgtcgt > 1:555758/1‑151 (MQ=255) ggCGTAAATGGCACCCTGAACCCAGCGGCCAGCGTCGGGGGAATGTGGGCGTGCAGTGGGGTAATAATGGCTACAATGACACCGCAGTGAGGGTGGGCCTAAAA‑TATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTc < 1:464640/149‑1 (MQ=255) gtgaaaaggacacccTAAACCCGGCGGCCAGCGTTGGGGGAAGGGGGGGGGGCGGGGGGGTGATAGTGGCAACAATCCCCCCGAGGTAAGGGGGGGCATAAAA‑TATAAGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATc < 2:279624/138‑1 (MQ=255) cGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAA‑TATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCg > 1:300449/1‑151 (MQ=255) gTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAA‑TATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCg > 1:337008/1‑150 (MQ=255) gTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAA‑TATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCg > 2:845549/1‑150 (MQ=255) tAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTAATA‑TATATGATCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAc > 1:299552/1‑151 (MQ=255) tGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAA‑TATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCgg > 2:340897/1‑151 (MQ=255) ggCACCCTGACCCCAGGGGCCAGCGTGTGGGGAAGGGGGGCGTGCAGTGGGGTGAAAATGGCTACAATGCCCCCGCAGTGATGGTGGGCCTAAAA‑TAAAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGt < 1:196771/151‑1 (MQ=255) acccTGAACCCAGCGGCCAGCGTCTGGGGAAGGGGGGCGTGCAGCGGGGTGATAATGGCTACAATGACCCCGCAGTGATGGTGGGCCTAAAA‑TATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTAc < 1:844665/148‑1 (MQ=255) gAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAA‑TATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATg > 1:831070/1‑150 (MQ=255) aaaccccggggcaggCGTGTGGGAAAGGGGGACGTGCAGGTGTTTGATATCCGATAAAATGACCCCGCAGGGATGGTGGCTTTAAAA‑TATATTTTCTAATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGa < 2:204823/137‑1 (MQ=255) cccggCGGCCAGCGTGTGGGGAATGTGGGCGGGCAGCTGGGTGATAATGGCTAAAATGACACCGCAGTGAGGGTGGCCCTAAAA‑TATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAat < 1:910781/147‑1 (MQ=255) cggCGGCCAGCGTCGGGGAAATGTGGGGGGGCAGCTGGGTAATAATGTATACAATGCCCCCGTGGTGATGGTGGGCTTAAAA‑TATATTTTCTACCCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAatat < 2:380280/149‑1 (MQ=255) ggCGAGCGTGTGGGAAATGGGGGCGTGCGGCGGGGTGATAATGGTTACAATGCCCCCGAAGTGGGGGTGGCCCTAAAA‑TATAAGTTCTGCTCCCGCCGTTAGCCAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAAtatt < 1:610320/148‑1 (MQ=255) tCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAA‑TATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCg > 2:87024/1‑151 (MQ=255) gggggATTGGGGGGGTGCACGTGGGTAAAAATGTCACAAATGCCACGGAGGTGATGTGGGGCTTAAAA‑TAAAATTCCTGCCCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAtttt < 1:786449/142‑1 (MQ=255) ggggATTGTGGGCGTGCAGTTGGGTGATAAGGGCTAAAATGACCCCGCAGTGATGGTGGGCCTAAAA‑TATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTAATCGACGCGGTACATATGAAATATTAttttt < 1:126086/143‑1 (MQ=255) ggggAAGGGGGGGGGGCAGGGGGGTGATAGTGGCAAAAATGACCCGGAGGTGAGGGTGGGCATAAAA‑TATAAGTGCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGAt < 2:877408/150‑1 (MQ=255) gggAGTGGGGGGGGGCAGGGGGGGTATAAGGGCTACAAGCACCCCGCGGTGATGGTGGGCTAAAAA‑TATATGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAg < 2:246752/151‑1 (MQ=255) tgtgGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAA‑TATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGc > 2:116668/1‑151 (MQ=255) gtgGGCGTGCAGCTGGGTGATAATGGCTACAATGACCCCGCAGTGATGGTGGGCCTAAAA‑TATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCa < 1:957515/151‑1 (MQ=255) gtgGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAA‑TATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGc > 2:977911/1‑150 (MQ=255) gggggggtgaggcggggTAATAAGGGCTACAATCCCCCCGCAGGGATGGGGGGCGTAAAA‑TATAATTTCTAATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCAc < 2:792417/134‑1 (MQ=255) ggggggtgaaggggggTGATAAGGGCTAAATGGACCCCGAGGGGGTGGGGGGCCTAAAA‑TAAAAGTTTTGCTCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCAt < 1:386749/138‑1 (MQ=255) gTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAA‑TATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGgtgt > 2:36294/1‑151 (MQ=255) aGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGGAA‑TATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGGTAGCACGCATGGGGtcat > 2:426390/1‑150 (MQ=255) ggagtgtgggggGGCTAAAATTATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCggg < 1:521365/144‑1 (MQ=255) tggtggGCCTGAAA‑TATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAgc > 2:498976/1‑151 (MQ=255) ggggggCCTAAAA‑TATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAgc < 1:397855/147‑1 (MQ=255) | TGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAA‑TATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGC > NC_000913/1222084‑1222373 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |