| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 2,749,880 | T→C | 20.0% | H218R (CAT→CGT) | PP_2406 ← | quinate/shikimate dehydrogenase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 2,749,880 | 0 | T | C | 20.0% | 36.1 / 4.5 | 20 | H218R (CAT→CGT) | PP_2406 | quinate/shikimate dehydrogenase |
| Reads supporting (aligned to +/- strand): ref base T (9/7); new base C (2/2); total (11/9) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.70e-01 | |||||||||||
GCCGCCCGCTGAACAGCTCGAATGCCTTCACCGCCTGGAACACCGCCATGTTGCTGCCATCCAGCGTGCGGCAGCCCAGTGCCCGGGCCGCGCGCAGCAGTTCGGTCTCCAGCGGGAAGTAGATGATCTCTGCCACCCACAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGG > NC_002947/2749731‑2750015 | gccgccCGCTGAACAGCTCGAATGCCTTCACCGCCTGGAACACCGCCATGTTGCTGCCATCCAGCGTGCGGCAGCCCAGTGCCCGGGCCGCGCGCAGCAGTTCGGTCTCCAGCGGGAAGTAGATGATCTCTGCCACCCACAGGCGGGGAt < 2:106475/150‑1 (MQ=255) cTCGAATGCCTTCACCGCCTGGAACACCGCCATGTTGCTGCCATCCAGCGTGCGGCAGCCCAGTGCCCGGGCCGCGCGCAGCAGTTCGGTCTCCAGCGGGAAGTAGATGATCTCTGCCACCCACAGGCGGGGATGCAACAACTCCACCgg > 1:436959/1‑150 (MQ=255) aCCGCCTGGAACACCGCCATGTTGCTGCCATCCAGCGTGCGGCAGCCCAGTGCCCGGGCCGCGCGCAGCAGTTCGGTCTCCAGCGGGAAGTAGATGATCTCTGCCACCCACAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGccc > 1:385569/1‑150 (MQ=255) aaCACCGCCATGTTGCTGCCATCCAGCGTGCGGCAGCCCAGTGCCCGGGCCGCGCGCAGCAGTTCGGTCTCCAGCGGGAAGTAGATGATCTCTGCCACCCACAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTg < 1:38285/150‑1 (MQ=255) tgcCATCCAGCGTGCGGCAGCCCAGTGCCCGGGCCGCGCGCAGCAGTTCGGTCTCCAGCGGGAAGTAGATGATCTCTGCCACCCACAGGCGGGGCTGCAACAACTCGGCCGGCAGTGGCGTGCGTGGTAGCTTGGGCATGCCCACCGGcc > 2:201365/1‑149 (MQ=255) aTCCAGCGTGCGGCAGCCCAGTGCCCGGGCCGCGCGCAGCAGTTCGGTCTCCAGCGGGAAGTAGATGATCTCTGCCACCCACAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAg < 2:9/126‑1 (MQ=255) aTCCAGCGTGCGGCAGCCCAGTGCCCGGGCCGCGCGCAGCAGTTCGGTCTCCAGCGGGAAGTAGATGATCTCTGCCACCCACAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGGAg > 1:9/1‑126 (MQ=255) cAGCGTGCGGCAGCCCAGTGCCCGGGCCGCGCGCAGCAGTTCGGTCTCCAGCGGGAAGTAGATGATCTCTGCCACCCACAGGCGGGGACGcaacaa < 1:270146/96‑1 (MQ=255) cAGCGTGCGGCAGCCCAGTGCCCGGGCCGCGCGCAGCAGTTCGGTCTCCAGCGGGAAGTAGATGATCTCTGCCACCCACAGGCGGGGACGcaacaa < 1:270094/96‑1 (MQ=255) cAGCGTGCGGCAGCCCAGTGCCCGGGCCGCGCGCAGCAGTTCGGTCTCCAGCGGGAAGTAGATGATCTCTGCCACCCACAGGCGGGGACGcaacaa > 2:270094/1‑96 (MQ=255) cAGCGTGCGGCAGCCCAGTGCCCGGGCCGCGCGCAGCAGTTCGGTCTCCAGCGGGAAGTAGATGATCTCTGCCACCCACAGGCGGGGACGcaacaa > 2:270146/1‑96 (MQ=255) gTGCGGCAGCCCAGTGCCCGGGCCGCGCGCAGCAGTTCGGTCTCCAGCGGGAAGTAGATGATCTCTGCCACCCACAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCgt < 2:19157/140‑1 (MQ=255) gTGCGGCAGCCCAGTGCCCGGGCCGCGCGCAGCAGTTCGGTCTCCAGCGGGAAGTAGATGATCTCTGCCACCCACAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCgt > 1:19157/1‑140 (MQ=255) ggCAGCCCAGTGCCCGGGCCGCGCGCAGCAGTTCGGTCTCCAGCGGGAAGTAGATGATCTCTGCCACCCACAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGcc > 2:212362/1‑150 (MQ=255) tGCCCGGGCCGCGCGCAGCAGTTCGGTCTCCAGCGGGAAGTAGATGATCTCTGCCACCCACAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCATCCCGTCCGCTTc < 1:17848/150‑1 (MQ=255) ccGCGCGCAGCAGTTCGGTCTCCAGCGGGAAGTAGATGATCTCTGCCACCCACAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAgg < 1:212362/150‑1 (MQ=255) gCGGGAAGTAGATGATCTCTGCCACCCACAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCa > 1:309945/1‑150 (MQ=255) gAAGTAGATGATCTCTGCCACCCACAGGCGGGCATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGgc < 2:385569/150‑1 (MQ=255) tgatCTCTGCCACCCACAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGGATCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCcg > 2:668950/1‑150 (MQ=255) ccacccacAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAAt > 1:137977/1‑150 (MQ=255) cccacAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATggg < 2:137977/150‑1 (MQ=255) | GCCGCCCGCTGAACAGCTCGAATGCCTTCACCGCCTGGAACACCGCCATGTTGCTGCCATCCAGCGTGCGGCAGCCCAGTGCCCGGGCCGCGCGCAGCAGTTCGGTCTCCAGCGGGAAGTAGATGATCTCTGCCACCCACAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGG > NC_002947/2749731‑2750015 |
| 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: ‑ |