Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 650,809 | 0 | T | G | 32.6% | 82.9 / 21.1 | 49 | T350P (ACC→CCC) | citC | citrate lyase synthetase |
Reads supporting (aligned to +/- strand): ref base T (11/22); new base G (16/0); total (27/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.29e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAGCGGATAT > NC_000913/650680‑650947 | ctATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGa > 1:889145/1‑140 (MQ=255) ttGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTCGTTTCTTTTCATTTTTCACCTGATTTTTATGCCGGGGACTTTTGACGAGCTGCcgc > 2:243030/1‑140 (MQ=255) ggTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCcg > 1:420562/1‑134 (MQ=255) ggTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCcg < 2:420562/134‑1 (MQ=255) gTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTg < 1:866150/140‑1 (MQ=255) tCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTgg < 1:714137/140‑1 (MQ=255) cGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGc > 2:556016/1‑140 (MQ=255) tatCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCggg < 1:497842/140‑1 (MQ=255) tGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGc > 2:240167/1‑140 (MQ=255) tGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGGTTTTTTTTATTTTTTACCTGTTTTTCATGCGGGGGGCTTTTTACGAGCTGCCGCGTCCTGGCgggggggg > 2:356831/1‑135 (MQ=255) tGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGGTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCgggggggg > 2:178689/1‑135 (MQ=255) tGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGGTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCgggggggg > 1:156977/1‑135 (MQ=255) tATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGCTCaa > 2:1172900/1‑140 (MQ=255) aGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAgg < 1:147360/140‑1 (MQ=255) ggggCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTTACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGGCCTGGCGGGGGGGCTCCAGCAgggtt > 2:147360/1‑137 (MQ=255) gggCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCt < 1:622113/140‑1 (MQ=255) gggCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCt < 1:1043638/140‑1 (MQ=255) aTGCGTATCATCACATCCCCAGACTCAAGGGTGCTTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCaaa < 2:1421550/140‑1 (MQ=255) tGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTTACGAGCTGCCGCGTCCTGGCGGGGGGGGTCAAGCAGGGTCTTCaaaa > 2:1142672/1‑139 (MQ=255) catcaCATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGGTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGCTCAAACAGGGTCTTCAAATAAAGCAg > 2:799096/1‑140 (MQ=255) atcaCATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGc < 1:556016/140‑1 (MQ=255) caTCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGa < 2:889145/140‑1 (MQ=255) caTCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGGGCTCAAACAGGTTCTGCAAATAATGCAGCGTGa > 1:220351/1‑140 (MQ=255) aTCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGGGCTCAAGCAGGGTCTGCAAATAATGCAGCGTGAc > 1:204678/1‑140 (MQ=255) ccccAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCCCCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGCTCAAGCAGGTTTTGCCAATAAAGCAGCGTGCCtt > 2:220540/1‑139 (MQ=255) aGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAgc < 1:1172900/140‑1 (MQ=255) gggTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTTACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGTGCTCAAGCAGGGTTTGCAAATAATGCAGCGGGACTGCAGGGGCCAgcg > 2:735032/1‑140 (MQ=255) tGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTTACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGGTCAAACAGGGTCTGCAAAAAAAGCAGCGTGACTGCAGGGACCAgcggcg > 1:883473/1‑140 (MQ=255) gCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGCTCAAGCAGGGTCTGCAAATAAAGCAGCGTGACTGCAGGGGCCAGCGgcgc > 2:1345971/1‑140 (MQ=255) gCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCCGGCGGGGGGGGTCAAGCAGGGTTTTCAAAAAATGCAGCGGGGCTGCAGGGGCCCGCGgggc > 2:1053828/1‑140 (MQ=255) ccTGCAACGGCGGGCTGGTTTTTTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGGGCTCAAACAGGTTTTGCCAATAAAACAGCGGGGCTGCAGGGGCCAGCGGcgcg > 2:254269/1‑140 (MQ=255) ccTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGCTCAAGCAGGGTTTGCAAAAAATGCAGGGTGGCTGCAGGGGCCAACGGggcg > 2:1374667/1‑140 (MQ=255) cTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGa < 1:1142672/140‑1 (MQ=255) cTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGa < 1:240167/140‑1 (MQ=255) aCGGCGGGCTGGGTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGCTCAAGCAGGGTTTTCAAATAAAGCAGCGGGACTGCCGGGGCCAGCGGCGCGGTAGcc > 2:632293/1‑140 (MQ=255) cggcggGCCGGTTTATTTTCATTTTTCACCTGTTTTTCATGCGGGGGGCTTTTGACGAGCTGCCGCGTCCTGGCGGGGGGGGTCAAGCAGGGTTTGCAAAAAAAGCAGCGGGGCTGCAGGGGCCAGCGGCGCGAAAGccc > 2:1249476/1‑139 (MQ=255) tttCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGc < 1:73936/140‑1 (MQ=255) tttttCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCcagca < 1:1053828/140‑1 (MQ=255) tttttCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCcagca < 1:1345971/140‑1 (MQ=255) tCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTg > 2:1335393/1‑140 (MQ=255) ccTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACg < 1:254269/140‑1 (MQ=255) cTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGt < 1:1374667/140‑1 (MQ=255) tttCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTAcc > 2:1230914/1‑140 (MQ=255) tttCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTAcc > 1:151593/1‑140 (MQ=255) tttCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTAcc > 2:76845/1‑140 (MQ=255) tttCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGCGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTAcc > 2:1159510/1‑140 (MQ=255) cATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGa < 1:220540/140‑1 (MQ=255) tGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAg > 1:657529/1‑140 (MQ=255) gCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAGc < 2:220351/140‑1 (MQ=255) gggggTTTTTTAACGAGCTGCGGCGCCCTGGGGGGGGTCTTCAAGCAGGTTCTGCAAATAAAGCGGCGGAACTGCGGGGGCACGCGGCGCGACAACCTAGAGAGCGTTTTTCGCAAACAGTTGTCGTACCCGGGAATCgg > 2:904804/1‑140 (MQ=255) gggggTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAg < 2:1193264/137‑1 (MQ=255) gggggTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAg > 1:1193264/1‑137 (MQ=255) gggTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAGCGGat < 1:632293/140‑1 (MQ=255) gTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAGCGGatat > 2:425690/1‑140 (MQ=255) | CTATTGATTTGCAGGTCGATATCCTGCGTATCGAGTGGGGCGATGCGTATCATCACATCCCCAGACTCAAGGGTGCCTGCAACGGCGGGCTGGTTTATTTTCATTTTTCACCTGTTTCTCATGCGGGGGTCTTTTGACGAGCTGCCGCGTCCTGGCGGGAGTGCTCAAGCAGGTTCTGCAAATAATGCAGCGTGACTGCAGGGACCAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAGCGGATAT > NC_000913/650680‑650947 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |