Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 646,765 | 0 | G | T | 54.3% | 23.6 / 54.0 | 35 | Q249K (CAA→AAA) | citG | 2‑(5''‑triphosphoribosyl)‑3'‑dephosphocoenzyme‑A synthase |
Reads supporting (aligned to +/- strand): ref base G (14/2); new base T (19/0); total (33/2) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.02e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.95e-06 |
CAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGT‑TCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCCAGC > NC_000913/646616‑646907 | 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:637006/1‑150 (MQ=255) gggT‑TCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTa > 2:538476/1‑150 (MQ=255) gT‑TCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAAcc > 1:591721/1‑151 (MQ=255) gAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTg > 2:175719/1‑151 (MQ=255) aTGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGAt > 2:447334/1‑151 (MQ=255) gCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATcc > 1:350030/1‑151 (MQ=255) ccccTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATcc > 1:194124/1‑149 (MQ=255) cccTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCcagc > 2:10963/1‑151 (MQ=255) cTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCcagc > 2:336429/1‑149 (MQ=255) cTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCcagc > 2:149210/1‑149 (MQ=255) tttttGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCcagc > 1:609148/1‑148 (MQ=255) | CAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGT‑TCGAATGCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACCCCTGATCCAGC > NC_000913/646616‑646907 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |