| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,222,233 | 0 | G | A | 43.8% | 36.7 / 13.7 | 32 | pseudogene (60/75 nt) | ycgH | pseudogene; putative ATP‑binding component of a transport system |
| Reads supporting (aligned to +/- strand): ref base G (13/5); new base A (0/14); total (13/19) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.94e-05 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACA‑CCGCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGC > NC_000913/1222085‑1222375 | gtgAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACA‑CCGCAGTGATGGTGGGCC‑TGaa < 2:2209/151‑1 (MQ=255)gtgAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACA‑CCGCAGTGATGGTGGGCC‑TGaa < 1:364923/151‑1 (MQ=255) aaGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACA‑CCGCAGTGATGGTGGGCC‑TGAAatat < 2:75359/151‑1 (MQ=255) gAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACA‑CCGCAGTGATGGTGGGCC‑TGAAATATAAGtt < 2:381900/151‑1 (MQ=255) aCAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACA‑CCGCAGTGATGGTGGGCC‑TGAAATATAAGTTCt < 1:792481/151‑1 (MQ=255) gATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACA‑CCGCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATc > 2:60051/1‑151 (MQ=255) gAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACA‑CCGCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGtt > 2:129446/1‑151 (MQ=255) ggggaAAATAAAACTTGGCGTAAAGGGCACCCTAAACCCAGCGGCCGGCGTCGGGGGAATGTGGGCGGGCAGCTGGGTGATAATGCCTACAATGACC‑CCGCAGTGATGGTGGGCC‑TAAAATATAAGTTCTGATCCCGCCGTTAGCTaaa < 1:333977/149‑1 (MQ=255) ggAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACA‑CCGCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAAcc > 1:766256/1‑151 (MQ=255) tAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACA‑CCGCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGcgtcg > 1:735925/1‑151 (MQ=255) aaaCTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACA‑CCGCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGcgtcgt > 1:281369/1‑151 (MQ=255) aaTGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACA‑CCGCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAcgc > 2:619099/1‑151 (MQ=255) ggCACCCTGATCCCGGCGGCAAGCGTCGGGGGAATGTGGGCGTGCAGTGGGGTGATAAGGGCTACAAGGACA‑CCGGAGTGTGGGTGGGCC‑AAAAATATAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAcgcg < 1:670739/149‑1 (MQ=255) cTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACA‑CCGCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACatat > 1:317791/1‑151 (MQ=255) tGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACA‑CCGCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATg > 1:93560/1‑151 (MQ=255) aaacccAGCGGCCAGCGTCGGGGAAATGTGGGCGTCGGGCGGGGTGATAATGCCTACAATAACC‑CCGAAGGGTGGGGGGCCC‑TAAAATATAAGTTCTAACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACata < 1:135419/144‑1 (MQ=255) aaTCCGGCGGCCAGCGTCGGGGGAATGTGGGCGTGCAGGTGGGTGATAAGGGCTAAAATAACC‑CGGCAGTGAGGGTGGCCT‑TAAAATATAAGTTTTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACatat < 1:231981/148‑1 (MQ=255) aaTCCAGCGGCCAGCGTCTGGGGAGTGTGGGCGTGCGGCTGGGTAAT‑AGGGGTAAAAAGGCACCCGGGGTGATGGGGGGCA‑TAAAATATAAGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACatat < 2:429529/148‑1 (MQ=255) aaTCCAGCGGCAAGTGTGGGGGGAATGGGGGCTTGCGGTTGGGTAATAATGGCTAAAACGACA‑CCGCTGTGGTGGGGGGTC‑TAAAATATAATTTCTGTCCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGa < 1:560823/150‑1 (MQ=255) ggCCTGCGTGTGGGGATTGGGGGCTGGCAGCGGGTTAATATGGGCTAAAAGCACA‑CCGCAGTGATGGTGGGCT‑TAAAATAAAAGTTGTATTCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAAtattat < 1:390329/150‑1 (MQ=255) ggCCGGCGTGGGGGGAATGGGGGCGTGCAGTTGGTTAATATGGGCAACAATGACC‑CCGAAGTGAGGGGGGCCA‑TAAAATATAAGTTTTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAAtattat < 2:185956/150‑1 (MQ=255) ggggATTGGGGGCGTGCGGTTGGTTAATAAGGGCTACAAGGCCC‑CCGCAGTGTTGGGGGCCA‑TAAAATATAATTTCTAATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAttttt < 1:328141/143‑1 (MQ=255) ggAAGGGGGGCGTGCAGTGGGGTGATAATGGCTAAAATGACC‑CCGCAGTGAGGGTGGCCC‑AAAAATAAAAGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGc < 1:690169/151‑1 (MQ=255) gggggCGTGCAGTGGGGTGATAATGGTTAAAATGACC‑CCGCAGTGATGGTGGGCC‑TAAAATATAAGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACg < 1:196926/147‑1 (MQ=255) gggCGTGCAGTGGGGGAATAATGGCAACAATGACC‑CCGAAGTGGTGTTGGGCG‑TAAAATATAAGTTTTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCAt < 2:318590/150‑1 (MQ=255) ccGCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCCCGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCgg > 2:454913/1‑151 (MQ=255) cGCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCggg > 2:37368/1‑151 (MQ=255) gTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCtt > 1:227543/1‑151 (MQ=255) gggaggtggGCCTTAAAATAAAAGTTCGGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCt < 2:388295/146‑1 (MQ=255) ggggggCC‑AAAAATATAGTTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAgcg < 2:183110/148‑1 (MQ=255) tggGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAgcgc > 2:403512/1‑150 (MQ=255) tggGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAgcgc > 1:30351/1‑150 (MQ=255) | GTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACA‑CCGCAGTGATGGTGGGCC‑TGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGC > NC_000913/1222085‑1222375 |
| 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: ‑ |