Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 646,747 | 0 | T | G | 56.2% | 11.5 / 23.9 | 32 | T255P (ACC→CCC) | citG | 2‑(5''‑triphosphoribosyl)‑3'‑dephosphocoenzyme‑A synthase |
Reads supporting (aligned to +/- strand): ref base T (12/2); new base G (18/0); total (30/2) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.83e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.21e-02 |
CCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGT‑TCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACC > NC_000913/646606‑646896 | ccAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGG‑TTGGATGcc > 1:103852/1‑151 (MQ=255) ttATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGT‑GCGAAGCCCCCCtt > 2:536971/1‑151 (MQ=255) ttATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGG‑TCGAATGCCCCCctt > 2:139117/1‑151 (MQ=255) ttATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGG‑TCGAATGCCCCCctt > 1:530721/1‑151 (MQ=255) cAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGG‑TCGAATACCCCCTTTTTgg > 1:513351/1‑150 (MQ=255) atTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGG‑TTGAATGCCCCCCTTTTTCaataaaattt > 2:408494/1‑142 (MQ=255) tAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTGTGGATAGCCCCCATTTTTCAA‑AATGtttt > 2:42519/1‑145 (MQ=255) tAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGT‑TCGAATGCCCCCTTTTTGCAATAATGTTTGCGcc < 2:530695/150‑1 (MQ=255) tAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGT‑TCGAAAGCCCCCCTTTTTCAATAATGtttttgcc > 1:404561/1‑145 (MQ=255) aaTCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGG‑TCCAATGACCCCCTTTTTCAATAATTTTTTCGcccc > 2:352381/1‑148 (MQ=255) cTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTAGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGG‑TCTCAATCCCCCTTTTTGTACAAATATTTGCGcccc > 2:491639/1‑145 (MQ=255) aCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGG‑TCGAATGCCCCCCTTTTTCAAAAATTTTTTCGCCTCGCGCCGTAGcc > 2:450800/1‑147 (MQ=255) cAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGT‑TTTAAT‑CGCCCTCTTTTCAAAAATATTTTGGCCACGCGCCGtttccca > 2:311794/1‑141 (MQ=255) cAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGG‑TTGAATGCCCCCCTTTTTCCATAATATTTTTGCCTCGCGCTGTAGCCAgcgc > 1:312472/1‑150 (MQ=255) cAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGG‑TCGAATGCCCCCCTTTTTCAATAATGTTTTCGCCCCGCGCCGTAGcc > 1:235049/1‑145 (MQ=255) aGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGG‑TCGAATGCCCCCCTTTTTCAAAAATGGTTGTGGGTCGGGCTTGAGGCAGCGGAggggcc > 2:66011/1‑147 (MQ=255) cAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGG‑TCGAATTCCCCCCTTTTTTAATAATGTTTGTGCCCCCCCCCGTAGCCCACGCAGGcccccccc > 1:313195/1‑148 (MQ=255) aGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGG‑TCGAATGCCCCCCTTTTTTAATAAATTTTTCGCCCCCCGCTGTAGCCAGCGCCGGcccccccccccccccccc > 1:630300/1‑139 (MQ=255) gccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGT‑TCGAATGCCCCCCTTTTTCAATAATGTTTTCGCCCCGCGCTGTAGCCAGCGCCGGCCCCCCCCCCCCCcccgc > 2:161016/1‑142 (MQ=255) agaTTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGG‑TCGAAAGCCCCCCTTTTTCAATAAAGTTTTCGCCCCGCGCTGTCGGCCGCGCCGGCCCCCCCCCCCCCCCCGCGATGCAACGtttggat > 1:390402/1‑143 (MQ=255) tcgttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGT‑TTGAAAGCCCCCCTTTTTGAATAATTTTTTCGCCCCGCGCTGTAGCCAGCGCAGGCCCCCCCCCCCACCCCCCGCTGCAACGTTGGGATCGc > 1:54461/1‑146 (MQ=255) tcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGG‑TTGAATGCCCCCCTTTTTTAATAAATTTTTCGCCCCGCGGCGTTGCCAGCGGAGGGCCCCCCCCCCCCCCCGCGATGCCACGCTGGTATCGCCGGTGAttg > 1:651109/1‑149 (MQ=255) tcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGG‑TTGAAAGCCCCCCCTTTTTCATAAATTTTTTGCCCCGCCCTGGTGGCCGCGCAGGCCCCCCCCGCCAACCCCCGATGCAACCTTGGGTTCGCCGGTGAt > 1:108238/1‑149 (MQ=255) acacTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGG‑TCGCATGCCCCCCTTTTTCAATAATTTTTTCGCCCCCCGCTGTAGCCAGCGCCGGCCCCCCCCCCCCCCCCGCGATGCAACGTTGGGATCGCCGCTGATCGcc > 1:355945/1‑148 (MQ=255) gAACTGCCGGAGATAATCGAGATCGGC‑GGGGGT‑TCCAATTCCCCCCTTTTTCAATAAAATTTTCGCCCCCCGCCGTAGCCAGCGCAGGCCCCCCCCCCCCCCCCGCGCTGCAACGCTTGTATTGCCCGTTATTGCCCTCCGTAgg > 1:710379/1‑145 (MQ=255) aaCTGCCGGAGATAATCGAGATCGGC‑GGGGGT‑TTGAATGGCCCCCTTTTTCAATAATGTTTTCGCCCCCCCCCGTAGCCCGCGCAGGCCCCCCCCCCCCCCCCCCGCTGCCACGGTGGTATCGCCCGTGATCGCCCTCCGTAGGGGCAAgg > 1:30760/1‑151 (MQ=255) agaTCGGC‑GGGGGT‑TCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCa < 1:754018/148‑1 (MQ=255) tCGGC‑GGGGGT‑TCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTCGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGGTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTc > 2:59864/1‑149 (MQ=255) tCGGC‑GGGGGG‑TCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTc > 2:110850/1‑149 (MQ=255) gggT‑TCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTa > 2:538476/1‑150 (MQ=255) gggT‑TCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTa > 2:637006/1‑150 (MQ=255) gT‑TCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAAcc > 1:591721/1‑151 (MQ=255) | CCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGT‑TCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACC > NC_000913/646606‑646896 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |