Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,049,646 | A→T | 13.5% | K718M (AAG→ATG) ‡ | quiA → | quinate dehydrogenase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,049,646 | 0 | A | T | 13.5% | 99.0 / 8.6 | 37 | K718M (AAG→ATG) ‡ | quiA | quinate dehydrogenase |
Reads supporting (aligned to +/- strand): ref base A (16/16); new base T (2/3); total (18/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.74e-01 |
TCGCCGGCATTCCATGCCAGGCACCGCCGTTTGGCACCCTGACCGCCATCGACATGAAGACCCGCCAGGTGGCCTGGCAGGTGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAACGAAATCTGGAAGGCTCGCCTGCCCGTGGGCAGCCAG > NC_002947/4049505‑4049791 | tcgCCGGCATTCCATGCCAGGCACCGCCGTTTGGCACCCTGACCGCCATCGACATGAAGACCCGCCAGGTGGCCTGGCAGGTGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCTTGATCGGca < 2:441955/150‑2 (MQ=255) tGCCAGGCACCGCCGTTTGGCACCCTGACCGCCATCGACATGAAGACCCGCCAGGTGGCCTGGCAGGTGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTggc < 2:167911/150‑1 (MQ=255) ccgccgTTTGGCACCCTGACCGCCATCGACATGAAGACCCGCCAGGTGGCCTGGCAGGTGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTg > 1:359341/1‑150 (MQ=255) cgccgTTTGGCACCCTGACCGCCATCGACACGAAGACCCGCCAGGTGGCCTGGCAGGTGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGt > 1:482874/1‑150 (MQ=255) cGCCATCGACATGAAGCCCCGCCAGGTGGCCTGGCAGGTGCCGGTAGGCATCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGgt < 2:26552/150‑1 (MQ=255) aTCGACATGAAGACCCGCCAGGTGGCCTGGCAGGTGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGgtgtcc > 2:620504/1‑148 (MQ=255) aCATGAAGACCCGCCAGGTGGCCTGGCAGGTGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCTTGATCGGcaggt > 1:301682/1‑98 (MQ=255) aCATGAAGACCCGCCAGGTGGCCTGGCAGGTGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCTTGATCGGcaggt < 2:301682/102‑5 (MQ=255) acgAAGACCCGCCAGGTGGCCTGGCAGGTGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCa < 2:482874/148‑1 (MQ=255) gAAGACCCGCCAGGTGGCCTGGCAGGTGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCTTGATCGGca > 2:683400/1‑94 (MQ=255) gAAGACCCGCCAGGTGGCCTGGCAGGTGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCTTGATCGGca < 1:683400/95‑2 (MQ=255) aGACCCGCCAGGTGGCCTGGCAGGTGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTa < 2:247283/150‑1 (MQ=255) gACCCGCCAGGTGGCCTGGCAGGTGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTAc > 1:248126/1‑150 (MQ=255) gTGGCCTGGCAGGTGCCGGTAGGCACCGTTGAACAAAACGGCAAAACCGTTATCGGTATTGGCCTTGCGCTGGAGATCGGCCTGCCAACCCTTGGCGGAACCCTGTCGACCCAAGTCGGCCTGGTGGTCATCGCAGGTACCCAGGACTTc > 1:203533/1‑150 (MQ=255) gCCTGGCAGGTGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTa < 2:380747/136‑1 (MQ=255) gCCTGGCAGGTGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTa > 1:380747/1‑136 (MQ=255) aGGTGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTAc > 1:655968/1‑130 (MQ=255) aGGTGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTAc < 2:655968/130‑1 (MQ=255) ccGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCa > 2:50162/1‑120 (MQ=255) ccGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCa < 1:50162/120‑1 (MQ=255) tAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGAca < 2:248126/150‑1 (MQ=255) gACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAAc < 1:593045/150‑1 (MQ=255) acCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCaa > 2:472862/1‑140 (MQ=255) acCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCaa < 1:472862/140‑1 (MQ=255) ccGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAACGaaa > 1:735405/1‑150 (MQ=255) cGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAACGAAATCTGGa < 2:359341/150‑1 (MQ=255) gCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAACGAAATCTGGaa > 1:632501/1‑150 (MQ=255) gCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAACAAAATCTTGaa > 1:106087/1‑150 (MQ=255) ggCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAACGAAATCTGGAAGGCt > 1:437482/1‑150 (MQ=255) ccGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGAcag > 2:312540/1‑117 (MQ=255) ccGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGAcag < 1:312540/117‑1 (MQ=255) aCCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCAc < 2:558062/40‑1 (MQ=255) aCCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCAc > 1:558062/1‑40 (MQ=255) aCCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAACGAAATCTGGAAGGCTCGCCTGCCCGTgg < 2:203533/150‑1 (MQ=255) aTCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCAt < 2:566761/62‑1 (MQ=255) aTCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCAt > 1:566761/1‑62 (MQ=255) aTCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAACGAAATCTGGAAGGCTCGCCTGCCCGTGGGCAGCCAg < 2:437629/150‑1 (MQ=255) | TCGCCGGCATTCCATGCCAGGCACCGCCGTTTGGCACCCTGACCGCCATCGACATGAAGACCCGCCAGGTGGCCTGGCAGGTGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAACGAAATCTGGAAGGCTCGCCTGCCCGTGGGCAGCCAG > NC_002947/4049505‑4049791 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |