Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | CP009273 | 2,434,063 | T→C | 12.8% | Q341R (CAG→CGG) | fabB ← | 3‑oxoacyl‑[acyl‑carrier‑protein] synthase I |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | CP009273 | 2,434,063 | 0 | T | C | 12.8% | 61.5 / 3.7 | 33 | Q341R (CAG→CGG) | fabB | 3‑oxoacyl‑[acyl‑carrier‑protein] synthase I |
Reads supporting (aligned to +/- strand): ref base T (9/17); new base C (2/2); total (13/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.11e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.58e-01 |
GACATAACGGTGGTCAGTTCGCGATCGGTCGTTTCGGTCACGATGTTCAGACCCGCAGCCTGCTCGTCCAGCTCTTCAATGTTGATGCTCGGGGCGATAAAGCCGTGTTCCAGCATCAGCAGAGAGTAGATAGCTTCCTGTACGCCAGCAGCGCCCAGAGAGTGACCGGTCATGGCTTTGGTTGCAGAAATCGCCGGGCTCTTATCGCCGAACACTTCACGGATAGCTGCCAGCTCTTTCACGTCGCCAACCGGAGTCGAAGTACCGTGGGAGTTCAGGTAA > CP009273/2433925‑2434206 | gACATAACGGTGGGCAGTTCGCGATCGGTCCTTTCGGTCACGATGTTCAGACCCGCAGCCTGCTCGTCCAGCTCTTCAATGTTGATGCTCGGGGCGATAAAGCCGTGTTCCAGCATCAGCAGAGAGTAGATAGCTTCCTGTACg < 2:118352/144‑1 (MQ=255) aaCGGTGGTCAGTTCGCGATCGGTCGTTTCGGTCACGATGTTCAGACCCGCAGCCTGCTCGTCCAGCTCTTCAATGTTGATGCTCGGGGCGATAAAGCCGTGTTCCAGCATCAGCAGAGAGTAGATAGCTTCCCGTACGCcagc > 2:270782/1‑144 (MQ=255) aaCGGTGGTCAGTTCGCGATCGGTCGTTTCGGTCACGATGTTCAGACCCGCAGCCTGCTAGTCCAGATCTTCAATGTTGATGCTCGGGGCGATAAAGCCGTGTTCCAGCATAAGCAGAGAGAAGATAGCTTCCCGTACGCctgc > 2:35055/1‑144 (MQ=255) tCAGTTCGCGCTCGGTCGTTTCGGTCACGATGTTCAGACCCGCAGCCTGCTCGTCCAGCTCTTCAATGTTGATGCTCGGGGCGATATAGCCGTGTTCCAGCATCCGCAGAGAGTAGATAGCTTCCTGTACGCCCGCAGCGCCCa < 2:156738/144‑1 (MQ=255) cgcgATCGGTCGTTTCGGTCACGATGTTCAGACCCGCAGCCTGCTCGTCCAGCTCTTCAATGTTGATGCTCGGGGCGATAAAGCCGTGTTCCAGCATCAGCAGAGAGTAGATAGCTTCCTGTACGCCAGCAGCGCCCAGAGAGTGACCg > 1:60337/1‑149 (MQ=255) gATCGGTCGTTTCGGTCACGATGTTCAGACCCGCAGCCTGCTCGTCCAGCTCTTCAATGTTGATGCTCGGGGCGATAAAGCCGTGTTCCAGCATCAGCAGAGAGTAGATAGCTTCCTGTACGCCAGCAGCGCCCAGAGAGTGAc < 2:267648/144‑1 (MQ=255) gATCGGTCGTTTCGGTCACGATGTTCAGACCCGCAGCCTGCTCGTCCAGCTCTTCAATGTTGATGCTCAGGGCGATAAAGCCGTGTTCCAGCATCAGCAGAGAGTAGATAGCTTCCTGTACGCCAGCCGCGCCCAGAGAGTGAc > 2:236215/1‑144 (MQ=255) gTTTCGGTCACGATGTTCAGACCCGCAGCCTGCTCGTCCAGCTCTTCAATGTTGATGCTCGGGGCGATAAAGCCGTGTTCCAGCATCAGCAGCGAGGAGATAACTTCCGTTACGCCAGCAGCGCCCAGAGAGTGACCGGTCATGGCttt > 1:153161/1‑149 (MQ=255) cGGTCACGATGTTCAGGCCCGCAGCCTGCTCGTCGAGCTCTTCAATGTTGATGCTCGGGGCGATAAAGCCGTGTTCCAGCATCAGCAGACAGTAGATAGCTTCCTGTACGCCAGCAGCGCCCAGAGAGTGACCGGTCATGGCTTTGGtt < 1:264541/149‑1 (MQ=255) cGGTCACGATGTTCAGACCCGCAGCCTGCTCGTCCAGCTCTTCAATGTTGATGCTCGGGGCGATAAAGCCGTGTTCCAGCATCAGCAGAGAGTAGATAGCTTCCTGTACGCCAGCAGCGCCCAGAGAGTGACCGGTCATGGCtt < 2:97875/144‑1 (MQ=255) cGGTCACGATGTTCAGACCCGCAGCCTGCTCGTCCAGCTCTTCAATGTTGATGCTCGGGGCGATAAAGCCGTGTTCCAGCATCAGCAGAGAGTAGATAGCTTCCTGTACGCCAGCAGCGCCCAGAGAGTGACCGGTCATGGCTTTGGtt < 1:311665/149‑1 (MQ=255) tCACGATGTTCAGGCCCGCAGCCTGCTCGTCCAGCTCTTCAATGTTGATGCTCGGGGCGATAAAGCCGTGTTCCAGAATTAGGAGAGAGTAGATAACTTCCGTTACTCCAGCAGCGCCCAGAGAGTGACCGGTCATGGCTTTgg < 2:343510/144‑1 (MQ=255) gATGTTCAGACCCGCAGCCTGCTCGTCCAGCTCTTCAATGTTGATGCTCGGGGCGATAAAGCCGTGTTCCAGCATCAGCAGAGAGTAGATAGCTTCCTGTACGCCAGCAGCGCCCAGAGAGTGACCGGTCATGGCTTTGGTTGCAGaaa < 1:369382/149‑1 (MQ=255) gCAGCCTGCTCGTCCAGCTCTTCAATGTTGATGCTCGGGGCGATAAAGCCGTGTTCCAGCATCCGCAGAGAGTAGATAGCTTCCTGTACGCCAGCAGCGCCCAGAGAGTGCCCGGTCATGGCTTTGGTTGCAGAAATCGCCGGGCTCtt > 1:186829/1‑149 (MQ=255) cAGCCTGCTCGTCCAGCTCTTCAATGTTGATGCTCGGGGCGATAAAGCCGTGTTCCAGCATCAGCAGAGAGTAGATAGCTTCCTGTACGCCAGCAGCGCCCAGAGAGGGACCGGTCATGGCTTTGGTTGCAGAAATCGCCGGGCCCTTa > 1:193037/1‑149 (MQ=255) tGCTCGTCCAGCTCTTCAATGTTGATGCTCGGGGCGATAAAGCCGTGTTCCAGCATCAGCAGAGAGTAGATAGCTTCCTGTACGCCAGCAGCGGCCAGAGGGTGACGGGTCATGGCTTTGGTTGCAGAAATCGCCGGGCTCTTa > 2:354491/1‑144 (MQ=255) cTCGTCCAGCACGTCAATGTTGATGCTCGGGGCGATAAAGCCGTGTTCCAGCATCAGCAGAGAGTAGATAGCTTCCTGTACGCCAGCAGCGCCCAGAGAGTGACCGGTCATGGCTTTGGTTGCAGAAATCGCCGGGCTCTTATCGCCGa < 1:358053/149‑1 (MQ=255) tCGTCCAGCTCTTCAATGTTGATGCTCGGGGCGATAAAGCCGTGTTCCAGCATCAGCAGAGAGTAGATAGCTTCCTGTACGCCAGCAGCGCCCAGAGAGTGACCGGTCATGGCTTTGGTTGCAGAAATCGCCGGGCGCTTAtgg > 2:56817/1‑142 (MQ=255) gTCCAGCTCTTCAATGTTGATGCTCGGGGCGATAAAGCCGTGTTCCAGCATCAGCAGAGAGTAGATAGCTTCCTGTACGCCAGCAGCGCCCAGAGAGTGACCGGTCATGGCTTTGGTTGCAGAAATCGCCGGGCTCTTATCGcc > 2:258611/1‑144 (MQ=255) ccAGCTCTTCAATGTTGATGCTCGGGGCGATAAAGCCGTGTTCCAGCATCAGCAGAGAGTAGATAGCTTCCTGTACGCCAGCAGCGCCCAGAGAGTGACCGGTCATGGCTTTGGTTGCAGAAATCGCCGGGCTCTTATCGCCGa < 2:135232/144‑1 (MQ=255) cAGCTCTTCAATGTTGATGCTCGGGGCGATAAAGCCGTGTTCCAGCATCAGCAGAGAGTAGATAGCTTCCTGTACGCCAGCAGCGCCCAGAGAGTGACCGGTCAGGGCTTTGGTTGCAGAAATCGCCGGGCTCTTATCGCCGaa < 2:164646/144‑1 (MQ=255) aGCTCTTCAATGTTGATGCTCGGGGCGATAAAGCCGTGTTCCAGCATCAGCAGAGAGTAGATAGCTTCCTGTACGCCAGCAGCGCCCAGAGAGGGACCGGTCATGGCACTGGTGGCAGAAATCGCCGGGCTCTTATCGCCGCACACTTc > 1:307088/1‑149 (MQ=255) ctctGCAATGTTGATGCTCGGGGCGATAAAGCCGTGTTCCAGCATCAGCAGAGAGTAGATAGCTTCCTGTACGCCAGCAGCGCCCAGAGAGTGACCGGTCATGGCTTTGGTTGCAGAAATCGCCGGGCTCTTATCGCCGAACACTTCAc < 1:317041/149‑1 (MQ=255) ctTCAATGTTGATGCTCGGGGCGATAAAGCCGTGTTCCAGCATCAGCAGAGAGTAGATAGCTTCCTGTACGCCAGCAGCGCCCAGAGAGTGACCGGTCATGGCTTTGGTTGCAGAAATCGCCGGGCTCTTATCGCCGAACACtt < 2:44702/144‑1 (MQ=255) ttCAATGTTGATGCTCGGGGCGATAAAGCCGTGTTCCAGCATCAGCAGAGAGTAGATAGCTTCCTGTCCGCCAGCAGCGCCCAGAGAGTGACCGGTCATGGCTTTGGTTGCAGAAATCGCCGGGCTCTTATCGCCGAACACTTCACGGa > 1:111798/1‑149 (MQ=255) ttcacgcccGGGGCGATACAGCCGGGTTCCAGCACCAGCAGAGAGGAGATAGCTTCCTGTACGCCAGCAGCGCCCAGAGAGTGACCGGTCATGGCTTTGGTTGCAGAAAGCGCCGGGCTCTTATCGCCGAACACTTCACGGATAGCTGc < 1:2718/141‑1 (MQ=255) gcccGGGGCGATAAAGCCGTGTTCCAGCATCAGCAGAGAGTAGATAGCTTCCCGTACGCCAGCAGCGCCCAGAGAGTGACCGGTCATGGCTTTGGTTGCAGAAAACGCCGGGCTCTTATCGCCGAACACTTCACGGATAGCTGCCAGct < 1:39530/146‑1 (MQ=255) gCTCGGGGCGATAAAGCCGTGTTCCAGCATCAGCAGAGAGTAGATAGCTTCCCGTACGCCAGCAGCGCCCAGAGAGTGACCGGTCATGGCTTTGGTTGCAGAAATCGCCGGGCTCTTATCGCCGAACACTTCACGGATAGCTGCCAGct < 1:302324/149‑1 (MQ=255) tCGGGGCGATAAAGCCGTGTTCCAGCATCAGCAGAGAGTAGATAGCTTCCGGTACTCCAGCAGCGCCCAGAGAGTGACCGGGCATGGCTTTGGTTGCAGAAATCGCCGGGCTCTTATCGCCGAACACTTCACGGGTAGCGGCCa > 2:232498/1‑144 (MQ=255) gCATCAGCAGAGAGTAGATAGCTTCCTGTACGCCAGCAGCGCCCAGAGAGTGACCGGTCATGGCTTTGGTTGCAGAAATCGCCGGGCTCTTATCGCCGAACACTTCACGGATAGCTGCCAGCTCTTTCACGTCGCCAACCGGAGTCGaa < 1:165204/149‑1 (MQ=255) agTAGATAGCTTCCTGTACGCCAGCAGCGCCCAGAGACTGACCGGTCATGGCTTTGGTTGCAGAAATCGCCGGGCTCTTATCGCCGAACACTTCACGGATAGCTGCCAGCTCTTTCACGTCGCCAACCGGAGTCGAAGTACCGt < 2:83470/144‑1 (MQ=255) aTAGCTGCCTGTACGCCAGCAGCGCCCAGAGAGTGACCGGTCATGGCTTTGGTTGCAGAAATCGCCGGGCTCTTATCGCCGAACACTTCACGGATAGCTGCCAGCTCTTTCACGTCGCCAACCGGAGTCGAAGTACCGTGGGAGTTCAg < 1:101795/149‑1 (MQ=255) tCCTGTACGCCAGCAGCGCCCAGGGTGTGACCGGGCATGGCTTTGGTTGCAGAAATCGCCGGGCTCTTATCGCCGAACACTTCACGGATAGCTACCAGCTCTTTCACGTCGCCAACCGGAGTCGAAGTACCGTGGGAGTTCAgg < 2:11785/144‑1 (MQ=255) tGTACGCCAGCAGCGCCCAGAGAGTGACCGGTCATGGCTTTGGTTGCAGAAATCGCCGGGCTCTTATCGCCGAAACCTTCACGGATAGCTTGGAGCTCTTTCACGTCGCCAACCGGAGTCGAAGTACCGTGGGAGTTCAGGTaa > 2:305620/1‑144 (MQ=255) | GACATAACGGTGGTCAGTTCGCGATCGGTCGTTTCGGTCACGATGTTCAGACCCGCAGCCTGCTCGTCCAGCTCTTCAATGTTGATGCTCGGGGCGATAAAGCCGTGTTCCAGCATCAGCAGAGAGTAGATAGCTTCCTGTACGCCAGCAGCGCCCAGAGAGTGACCGGTCATGGCTTTGGTTGCAGAAATCGCCGGGCTCTTATCGCCGAACACTTCACGGATAGCTGCCAGCTCTTTCACGTCGCCAACCGGAGTCGAAGTACCGTGGGAGTTCAGGTAA > CP009273/2433925‑2434206 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |