Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 2,982,740 | C→G | 6.9% | R161P (CGC→CCC) | PP_2608 ← | shikimate dehydrogenase family protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 2,982,740 | 0 | C | G | 6.9% | 178.5 / 4.3 | 58 | R161P (CGC→CCC) | PP_2608 | shikimate dehydrogenase family protein |
Reads supporting (aligned to +/- strand): ref base C (29/25); new base G (2/2); total (31/27) | |||||||||||
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.49e-01 |
GCGGACAACGGCAGCTCTGCGCGGGTGCCCATGCCGACTGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCTCGAAGCCATGCTTGTGCGCAGCGCCGA > NC_002947/2982596‑2982884 | gcggACAACGGCAGCTCTGCGCGGGTGCCCATGCCGACTGGGGACGCGTTGGCCACTAGGTTGAAGTCTTCCAGCCCGCTGATCTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCAtgcgtgcg < 1:456536/150‑1 (MQ=255) aCAACGGCAGCTCTGCGCGGGTGCCCATGCCGACTGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACGCATGGGTGCGGtat < 1:415585/150‑3 (MQ=255) ggCAGCTCTGGGCGGGTGCCCATGCAGAATGGGGACGCGTTGGTCACCAGGTCGAAGTCGTCCAGCCGGCTGAACTGTGTGTAGACTGTAAGGCAGGGGAAGGCGTTGCCGAGCAGTTCGCATAGCGCACCCATGCGTGCGTGGATgggg > 2:444476/1‑150 (MQ=255) gCAGCTCTGCGCGGGTGCCCATGCCGACTGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGt > 1:1078559/1‑150 (MQ=255) tctGCGCGGGTGCCCATGCCGACTGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCag > 1:1344262/1‑150 (MQ=255) tctGCGCGGGTGCCCATGCCGACTGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGATGGGGTCGCag > 2:239321/1‑150 (MQ=255) tGCGCGGGTGCCCATGCCGACTGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCagag > 2:164489/1‑150 (MQ=255) cgcgGGTGCCCATGCCGACGGGGGACGCGTTGGCCACCAGGTCGACGTCTTCCAGCCCACTGAACTGTGTGCCGCCCGTTAGCCCGGGGAAGCCGTTGCCGGGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAggg < 1:320401/150‑1 (MQ=255) cATGCCGACTGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGc < 1:137426/150‑1 (MQ=255) tGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGGGATCCCTGc < 2:910875/150‑1 (MQ=255) gggACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCt > 1:819483/1‑150 (MQ=255) gACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCa < 2:819483/150‑1 (MQ=255) cgcgTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGc < 1:741085/150‑1 (MQ=255) gcgTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATg < 1:979061/129‑1 (MQ=255) gcgTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATg > 2:979061/1‑129 (MQ=255) ttGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGt < 1:1055277/150‑1 (MQ=255) ccaccaGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGATGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCAt < 1:239321/150‑1 (MQ=255) ccaGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAgg > 1:1007202/1‑150 (MQ=255) aGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCa > 1:1421964/1‑150 (MQ=255) tCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATg < 2:1344262/150‑1 (MQ=255) ttCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCTAGCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGcc > 1:785406/1‑150 (MQ=255) aCTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAgg < 2:178443/91‑1 (MQ=255) aCTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAgg > 1:178443/1‑91 (MQ=255) tgtgtgCTGACCGTAAGCCCGGGGAAGCCGTTGCCGCGCAGTTGGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGGTATGCTGGCGATCCCTGCTTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCg > 1:180017/1‑150 (MQ=255) tgtgtgCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACGCATGGGTGCGGtat < 1:286596/70‑3 (MQ=38) tgtgtgCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACGCATGGGTGCGGtat > 2:286596/1‑68 (MQ=38) tgtgCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCt < 1:1185925/100‑1 (MQ=255) tgtgCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCt > 2:1185925/1‑100 (MQ=255) gtgCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATc < 2:1078559/150‑1 (MQ=255) tgCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCa < 1:157610/108‑1 (MQ=255) tgCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCa > 2:157610/1‑108 (MQ=255) tGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCCTCGACAATGGCACTGCCAACACCGCCGCAACCGATCACCa > 2:337350/1‑150 (MQ=255) aCCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATg < 2:395293/117‑1 (MQ=255) aCCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATg > 1:395293/1‑117 (MQ=255) aCCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGt > 1:741337/1‑150 (MQ=255) cGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGc > 1:132470/1‑150 (MQ=255) aGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCAATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGt > 1:341544/1‑150 (MQ=255) gggAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAg < 1:1009435/99‑1 (MQ=255) gggAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAg > 2:1009435/1‑99 (MQ=255) ggAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGctcg > 2:1127751/1‑147 (MQ=255) aaGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTACGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCATACATCAGCCATTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATAACCAGTGCCCGTTTGCccgcc > 1:261139/1‑150 (MQ=255) aGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCccgccg < 2:1272079/150‑1 (MQ=255) gCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCgg > 1:886220/1‑150 (MQ=255) cGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCt < 1:444476/150‑1 (MQ=255) gCCGAGCAGTTCGCATACCGCACCAATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCTCGaa < 2:341544/150‑1 (MQ=255) cccccAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCTCGAAGc > 2:347168/5‑150 (MQ=255) cGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGTGTCGCAGAGGGTTATGCTgg < 1:965844/60‑1 (MQ=255) cGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTgg > 2:965844/1‑60 (MQ=255) gCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCTCGAAGCCATGCTtgtg < 1:347168/150‑1 (MQ=255) gCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCTCGAAGCCATGCTtgtg < 1:98512/150‑1 (MQ=255) aTACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGc < 1:421337/58‑1 (MQ=255) aTACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGc > 2:421337/1‑58 (MQ=255) ccGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCTCGAAGCCATGCTtgt > 1:1163698/1‑144 (MQ=255) ccGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCTCGAAGCCATGCTtgt < 2:1163698/144‑1 (MQ=255) gCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCaa < 2:1332236/73‑1 (MQ=255) gCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCaa > 1:1332236/1‑73 (MQ=255) cACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCTCGAAGCCATGCTTGTGCGCAgcgc > 2:1291089/1‑150 (MQ=255) ccATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCTCGAAGCCATGCTTGTGCGCAGCGCCGa < 2:419053/150‑1 (MQ=255) aTGGGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACt > 2:1282935/1‑74 (MQ=255) | GCGGACAACGGCAGCTCTGCGCGGGTGCCCATGCCGACTGGGGACGCGTTGGCCACCAGGTCGAAGTCTTCCAGCCCACTGAACTGTGTGCTGACCGTAAGCCCGGGGAAGCCGTTGCCGAGCAGTTCGCATACCGCACCCATGCGTGCGGTGCTGGGGTCGCAGAGGGTTATGCTGGCGATCCCTGCCTCAGCCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCTCGAAGCCATGCTTGTGCGCAGCGCCGA > NC_002947/2982596‑2982884 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |