Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,222,233 | 0 | G | A | 56.5% | 4.7 / 26.6 | 23 | pseudogene (60/75 nt) | ycgH | pseudogene; putative ATP‑binding component of a transport system |
Reads supporting (aligned to +/- strand): ref base G (9/1); new base A (0/13); total (9/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.22e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.06e-04 |
GTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACC > NC_000913/1222083‑1222380 | gtgtGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTg < 2:220263/151‑1 (MQ=255) gtgAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGa > 2:162763/1‑150 (MQ=255) gAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGtt > 1:478905/1‑151 (MQ=255) gTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGtt > 2:197203/1‑149 (MQ=255) tAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGcgtcg > 2:35025/1‑151 (MQ=255) ttGGCTAAAATGGCAACCTGAATCCGGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCGGGGTGATAATGGCTACAATGACCCCGCAGTGATGGTGGGCCTAAAATATAAGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTAt < 2:38863/149‑1 (MQ=255) ggCGAAAATGCCACCCTGAACCCGGCGCCCGGCGTCGGGGGAAGGGGGGGGGGCAGCGGGGTGATAATGGCAACAATGACCCCGAAGTGATGGTGGGCGTAAAATATAAGTTCTAACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTc < 1:146063/149‑1 (MQ=255) aacccAGCGGCCAGCTTCTGGGGAATGGGGGCGTGCAGCGGGGTGAAAATGGCTACAATGACCCCGCAGTGTTGGTGGGCCTAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGa < 2:89544/147‑1 (MQ=255) aacccAGCGGCCAGCGTCGGGGGAATGGGGGCGTGCAGCTGGGTGATAATGGCTAAAATGACCCCGCAGTGATGGGGGGCCTAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATg < 2:237649/146‑1 (MQ=255) cGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAAtattat > 2:151357/1‑151 (MQ=255) ggCCAGCGTCGGGGAAATGGGGGGGTGCAGCTGGGTAATAATGCTTAAAATGACACCGAAGTGAGGGGGGGCTTAAAATAAAAGTTCTTACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATattatt < 1:492580/151‑1 (MQ=255) ccAGCGTCTGGGGAATGGGGGCGTGCAGCTGGGTGATAATGGCAACAATGACCCCGCGGTGATGGGGGGCCTGAAATATAGGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAtttt < 2:492666/151‑1 (MQ=255) gtgTCTGGGAAGTGTGGGCGTGCAGTGGGGTAATAATGACTACAATGACCCCGCAGGGAGGGGGGGCATAAAATAAAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAtttttt < 1:504577/148‑1 (MQ=255) gCGTCGGGGGAATGGGGGCGTGAAGCGGGGTGATAATGGCAAAAATGACCCCGGGGTGAGGGGGGGCCTAAAATATAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAtttttt < 2:111909/150‑1 (MQ=255) tCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCg > 1:90770/1‑151 (MQ=255) ggAATGTGGGCGTACAGCTGGGTGATAATGGCTAAAATGACACCGCAGGGATGGGGGGCCTAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGc < 2:120768/151‑1 (MQ=255) gtgGGCGTGCGGCTGGGTGATAAGGCCAACAATGACGCCGCAGTGATGGTGGGCTTAAAAAAAAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACg < 1:529287/149‑1 (MQ=255) gggCGTGCAGCTGGGTAATAATGGCTACAATGCCCCCGGTGTGATGGGGGGCAAAAAATAAATTTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCAt < 2:80758/150‑1 (MQ=255) gggCGTGCAGCGGGGTGATAATGGCAAAAATGACCCCGCGGTGATGGTGGGCTTAAAATATAATTTTTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCAt < 2:341087/150‑1 (MQ=255) acCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCg > 1:503798/1‑151 (MQ=255) gCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGt > 1:366560/1‑151 (MQ=255) ggtggGCCGAAAATATATTTTTTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAg < 2:438703/149‑1 (MQ=255) ggggggCCTAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTa < 2:479025/145‑1 (MQ=255) ccTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAAcc > 1:421792/1‑151 (MQ=255) | GTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACC > NC_000913/1222083‑1222380 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |