| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,222,233 | 0 | G | A | 42.9% | 27.9 / 12.2 | 28 | pseudogene (60/75 nt) | ycgH | pseudogene; putative ATP‑binding component of a transport system |
| Reads supporting (aligned to +/- strand): ref base G (9/7); new base A (0/12); total (9/19) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.83e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
AACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCA‑GTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTC > NC_000913/1222095‑1222367 | aaCAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCA‑GTGATGGTGGGCCTGAAATATAAGTTc < 1:261968/151‑1 (MQ=255) aCAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCA‑GTGATGGTGGGCCTGAAATATAAGTTCt < 2:270615/151‑1 (MQ=255) cAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCA‑GTGATGGTGGGCCTGAAATATAAGTTCTg < 2:633532/151‑1 (MQ=255) aTGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCA‑GTGATGGTGGGCCTGAAATATAAGTTCTGATcc < 2:855288/151‑1 (MQ=255) ggCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCA‑GTGATGGTGGGCCTGAAATATAAGTTCTGATCccgcc > 1:196661/1‑151 (MQ=255) ttGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCA‑GTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTc > 1:84640/1‑151 (MQ=255) ttGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGCGATAATGGCTACAATGACACCGCA‑GTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTc > 1:441445/1‑151 (MQ=255) ggCGTAAATGGCACCCTGAACCCAGCGGCCAGCGTCGGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTAAAATGACCCCGCA‑GTGATGGGGGGCATAAAATATAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTc < 1:383188/149‑1 (MQ=255) ggCGAAAAGGGCAACCTGAACCCGGGGGCCGGCGTCGGGGGAGTGGGGGCGGGCAGGTGGGTAATAAGGCATAAAATCCCCCCGAA‑GTAATGGTGGGCTTAAAATATATGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCAt < 2:557712/151‑1 (MQ=255) aaTGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCA‑GTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAcgc > 1:389093/1‑151 (MQ=255) ccggAATCCAGCGGCCAGCGTCTGGGGAAGGTGGGCGTGCGGCTGGGTGATAATGGCTACAATGACCCCGCG‑GTGATGGGGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACat < 1:324751/147‑1 (MQ=255) acccggCGGCAGGCGTTGGGGGAAGGGGGGGGGGCA‑CTGGGGTATAAAGGCCTAAATTACCCCGCATGTGAGGGTGGCCATGAAAAATAAGTCCTGCACCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGaaa < 1:835610/146‑1 (MQ=255) ccgccGGCCAGCGTCTGGGGAAGGGGGGCGGGCAGCGGGGTGATAATGACTACAAGGCCCCCGAA‑GTGTTGGGTGGCTTAAAATATAGGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAata < 2:112165/147‑1 (MQ=255) aGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCA‑GTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAAtatt > 2:251174/1‑151 (MQ=255) gcccggCGTCTGGGGATGGGGGGCGTGCAGGTGGGTAAAAATGCCTACAATGCCCCCGCA‑GTGAGGGGGGGCTTAAAATATAAGTTCGGTTCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACatat < 2:62747/135‑1 (MQ=255) cggCGTGTGGGGAGGGGGGGCGTGCAGGTGGGTGATATTGGCTACAGTGACGCGGGG‑GGGAGGGGGGGCTGGAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAttttt < 2:580759/149‑1 (MQ=255) cAGCGTGGGGGGATTGTGGGGGTGCAGTTGGGTGAAAATGGCAAAAAGCACCCCGCA‑GGGGTGGGGGGCCTAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAttttt < 2:809863/151‑1 (MQ=255) gCGTCTGGGGAATGTGGGGGTGCAGCGGGGTGATAATGGCTACAATGACCCCGCA‑GTGATGGTGGGCCTAAAATATAAGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAtttttt < 1:190818/150‑1 (MQ=255) gggCGTGCAGTGGGGTAAAAATGCCAAAAATGCCCCCGAA‑GTGAGGGGGGGCTGAAAATAAAATTTCTAATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCa < 1:811802/145‑1 (MQ=255) gggTGCAGCGGGGTAATAATGGCTACAATGACCCCGCA‑GTGAGGGTGGGCCTAAAATATAAGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGgt < 1:863657/149‑1 (MQ=255) aGCTGGGTGATAATGGCTACAATGACACCGCA‑GTGATGGTGGGCCTGAAATATAAGTTCTGATCCAGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGtcat > 1:371995/1‑150 (MQ=255) gggTGATAATGGCTAAAATGACACCGCA‑GTGATGGGGGGCATAAAATATAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCaa < 2:841664/149‑1 (MQ=255) ggTGATAATGGCTAAAATGACCCCGCG‑GTGATGGTGGGCTTAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTg < 2:362554/151‑1 (MQ=255) attacccccGCA‑GTGATGGGGGGCCTAAAATAAAGTTTTTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACa < 1:580426/141‑1 (MQ=255) cacCGCA‑GTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTg > 2:490080/1‑150 (MQ=255) acCGCA‑GTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCg > 1:777597/1‑151 (MQ=255) cA‑GTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTc > 1:398720/1‑151 (MQ=255) gTGATGGTGGGCCTAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGt < 1:373229/148‑1 (MQ=255) | AACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCA‑GTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTC > NC_000913/1222095‑1222367 |
| 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: ‑ |