| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,222,233 | 0 | G | A | 51.6% | 9.9 / 39.2 | 31 | pseudogene (60/75 nt) | ycgH | pseudogene; putative ATP‑binding component of a transport system |
| Reads supporting (aligned to +/- strand): ref base G (13/2); new base A (0/16); total (13/18) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.09e-07 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.48e-05 | |||||||||||
TCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCT‑GAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTGA > NC_000913/1222092‑1222383 | tCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCT‑GAAATATAAg < 2:393702/151‑1 (MQ=255) cAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGGGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCT‑GAAATATAAGtt > 2:305513/1‑148 (MQ=255) aGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCT‑GAAATATAAGTTCTGAt > 1:114396/1‑151 (MQ=255) tgggttAAATGGCACCCTAAACCCCGGGCCCAGGGTGTGGGAAAGGGGGGGGTGAGGCGGGGTAATAATGCCAACAACCACACCGAAGTGATGGGGGGCTT‑GAAATATAAGTTTTAAACCCGCCGTTAGCTAAAAAACCGCGTCGTATTc < 1:901271/145‑1 (MQ=255) ggCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCT‑GAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCAt > 1:725349/1‑151 (MQ=255) aaTGGCACCCGAAAGCCAGGGCCCGGCGTCGGGGGAATGTGGGGGTGCAGTTGGGTGATAATGGCTAAAAGCCCCCCGCAGTGGTGGTGGGCCT‑AAAATAAAAGTTGTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAcgc < 2:617102/151‑1 (MQ=255) aacccTGAAGCCAGGGCCCCGGGTGGGGGGAATGGGGGGGTGCAGTGGGGTGATAATGCCTACAATGACCCCGAAGTGATGGTGGGCCT‑AAAATAAAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTAc < 2:441851/148‑1 (MQ=255) aacccTGAAACCGGGGCCCGGCGTCGGGGGAGGGGGGGGGGGAAGTGGGGTAATATTGGCAAAAAGGACCCCGAAGGGGGGGTGGCCCT‑AAAATAAAAGTGCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTAc < 2:257086/148‑1 (MQ=255) aCCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCT‑GAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACat > 1:576499/1‑151 (MQ=255) taaacccAGGGGCCGGCGTTGGGGGAAGGGGGGCGTGCGGCTGGGTAAAAAGGCCTACAATCCCCCCGAATTGAGGGGGGGCCT‑AAAATAAAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATg < 2:702880/146‑1 (MQ=255) gaaccgggCGGCCAGCGGCGGGGGAAGGGGGGGGTGAAGTGGGGTAAAAAGGGCTAAAAGGACACCGCAGTGATGGGGGGCTT‑AAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATg < 1:733198/143‑1 (MQ=255) gAACCCGGCGGCCGGCGTGTGGGAAATGGGGGGGTGCGGGGGGGTAATAATGGCTAAAATCACACCGCGGGGATGGTGGCCTT‑AAAATATAAGTTCTTACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATg < 1:770743/150‑1 (MQ=255) aTCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCT‑GAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGaaa > 2:999494/1‑151 (MQ=255) cccAGGGGCCAGCGTCGGGGAATTGGGGGCGTGCGGTGGGGAAATATGGTCAAAAACGACCCCGCAGGGAGGGGGGGCCA‑AAAATAAAAGTTCTGAACCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAat < 2:1029273/150‑1 (MQ=255) ccAGCGGCCGGCTTCTGGGAATTGGGGGGGTGCAGCTGGTGAATAAGGGCTAGAAGGACCCCGTGGTGAGGGTGGGCAT‑AAAAAAAAGTTTTTAATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAata < 1:753061/151‑1 (MQ=255) ggCCGGCGTTGGGGAAATGGGGGCGTGCAGGGGGGTGAAAATGGCAAAAAGGCCCCCGCGGTGAGGGGGGGCCT‑AAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAAtattat < 2:195810/150‑1 (MQ=255) cGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCT‑GAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTg > 2:32257/1‑150 (MQ=255) tgtgGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCT‑GAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACg > 2:314751/1‑150 (MQ=255) gggtcaaGCTGGGTGAAAATGCCAACAAAGACCCGTCAGGGAGGGGGGGCAT‑AAAATAAAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATgg < 1:599585/144‑1 (MQ=255) cgggcgggggggTGAAAATGGCAAAAATAACCCCGCAGTGAGGGGGGGCCT‑AAAATTAAAGTTCTTATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGgtg < 2:34758/142‑1 (MQ=255) ggggggTTAAAATGCCTACAATGACACCGCAGTGATGGGGGGCCT‑AAAATATAAGTTTTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGtcatc < 2:207042/147‑1 (MQ=255) ggggggTAATTATGGCTAAAAGGACCCCGAGGTGATGG‑GGGCTTAAAAATAAATGTTCAGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGtcatc < 1:463249/147‑1 (MQ=255) cggggTGATAATGGCTAAAAGGCCACCGAGGTGATGGTGGGCCT‑AAAATAAAATTTCGGTTCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGtcatca < 1:347406/148‑1 (MQ=255) gtaatatTGTCTACAATGCCCCCACAGTGATGGTGGGCTT‑AAAATATAGGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTg < 1:594429/147‑1 (MQ=255) tGGCTACAATGACACCGCAGTGATGGTGGGCCT‑GAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGAt > 1:187554/1‑150 (MQ=255) ggCTACAATGACACCGCAGTGATGGTGGGCCT‑GAAATATAAGGTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATG‑AATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATaaa > 2:628995/1‑151 (MQ=255) cTAAAATGCCACGGGAGTGAGGGGGGCCCT‑AAAATATAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAAc < 2:443702/151‑1 (MQ=255) cAATGACACCGCAGTGATGGTGGGCCT‑GAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAg > 2:247462/1‑150 (MQ=255) cAATGACACCGCAGTGATGGTGGGCCT‑GAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAg > 1:248777/1‑150 (MQ=255) cAGTGATGGTGGGCCT‑GAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGt > 1:499835/1‑150 (MQ=255) gTGATGGTGGGCCT‑GAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCtt > 1:907458/1‑151 (MQ=255) gAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTGa > 2:628074/1‑151 (MQ=255) | TCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCT‑GAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTGA > NC_000913/1222092‑1222383 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |