| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 4,160,540 | A→G | 11.8% | E412E (GAA→GAG) | PP_3662 → | hypothetical protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 4,160,540 | 0 | A | G | 11.8% | 83.7 / 3.6 | 34 | E412E (GAA→GAG) | PP_3662 | hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base A (15/15); new base G (2/2); total (17/17) | |||||||||||
| 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.22e-01 | |||||||||||
CGTTCGACCCGACCCACCAGGCCATGGCCGACCTGGACTTGCGCCGCGAGCTGCCGCCCCATGAACTGGCGGCCAACCTGCGCCGGGCGTTTTCCGGGATTGTCGCAGGCAACGTCAAGGACAAAGGCATCCGCCTGATCGAGGAACACGGGCCGTATCAGATCCACGGTGACAGCGCCGTACTGGACCCGCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTG > NC_002947/4160395‑4160685 | cGTTCGACCCGACCCACCAGGCCATGGCCGACCTGGACTTGCGCCGCGAGCTGCCGCCCCATGAACTGGCGGCCAACCTGCGCCGGGCGTTTTCCGGGATTGTCGCAGGCAACGTCAAGGACAAAGGCATCCGCCTGATCGAGGAACACg > 2:246832/1‑150 (MQ=255) ccgacccACCAGGCCATGGCCGACCTGGACTTGCGCCGCGAGCTGCCGCCCCATGAACTGGCGGCCAACCTGCGCCGGGCGTTTTCCGGGATTGTCGCAGGCAACGTCAAGGACAAAGGCATCCGCCTGATCGAGGAACACGGGCCGTAt < 2:553063/150‑1 (MQ=255) gacccACCAGGCCATGGCCGACCTGGACTTGCGCCGCGAGCTGCCGCCCCATGAACTGGCGGCCAACCTGCGCCGGGCGTTTTCCGGGATTGTCGCAGGCAACGTCAAGGACAAAGGCATCCGCCTGATCGAGGAACACGGGCCGTATCa > 2:350865/1‑150 (MQ=255) ccaccaGGCCATGGCCGACCTGGACTTGCGCCGCGAGCTGCCGCCCCATGAACTGGCGGCCAACCTGCGCCGGGCGTTTTCCGGGATTGTCGCAGGCAACGTCAAGGACAAAGGCATCCGCCTGATCGAGGAACACGGGCCGTATCAGAt > 1:386496/1‑150 (MQ=255) caGGCCATGGCCGACCTGGACTTGCGCCGCGAGCTGCCGCCCCATGAACTGGCGGCCAACCTGCGCCGGGCGTTTTCCGGGATTGTCGCAGGCAACGTGAAGGACAAAGGCATCCGCCTGATCGAGGAACACGGGCCGTATCAGATCCAc < 1:350865/150‑1 (MQ=255) ggCCATGGCCGACCTGGACTTGCGCCGCGAGCTGCCGCCCCATGAACTGGCGGCCAACCTGCGCCGGGCGTTTTCCGGGATTGTCGCAGGCAACGTCAAGGACAAAGGCATCCGCCTGATCGAGGAACACGGGCCGTATCAGATCCACgg > 2:245853/1‑150 (MQ=255) gCCATGGCCGACCTGGACTTGCGCCGCGAGCTGCCGCCCCATGAACTGGCGGCCAACCTGCGCCGGGCGTTTTCCGGGATTGTCGCAGGCAACGTCAAGGACAAAGGCATCCGCCTGATCGAGGAACACGGGCCGTATCAGATCCACGGt < 2:385773/150‑1 (MQ=255) aTGGCCGACCTGGACTTGCGCCGCGAGCTGCCGCCCCATGAACTGGCGGCCAACCTGCGCCGGGCGTTTTCCGGGATTGTCGCAGGCAACGTCAAGGACAAAGGCATCCGCCTGATCGAGGAACACGGGCCGTATc > 2:485852/1‑136 (MQ=255) aTGGCCGACCTGGACTTGCGCCGCGAGCTGCCGCCCCATGAACTGGCGGCCAACCTGCGCCGGGCGTTTTCCGGGATTGTCGCAGGCAACGTCAAGGACAAAGGCATCCGCCTGATCGAGGAACACGGGCCGTATc < 1:485852/136‑1 (MQ=255) ccGACCTGGACTTGCGCCGCGAGCTGCCGCCCCATGAACTGGCGGCCAACCTGCGCCGGGCGTTTTCCGGGATTGTCGCAGGCAACGTCAAGGACAAAGGCATCCGCCTGATCGAGGAACACGGGCCGTATCAGATCCACGGTGACAgcg < 2:11130/150‑1 (MQ=255) cTGGACTTGCGCCGCGAGCTGCCGCCCCATGAACTGGCGGCCAACCTGCGCCGGGCGTTTTCCGGGATTGTCGCAGGCAACGTCAAGGACAAAGGCATCCGCCTGATCGAGGAACACGGGCCGTATCAGATCCACGGTGACAGCGCCGTa < 1:337546/150‑1 (MQ=255) ttGCGCCGCGAGCTGCCGCCCCATGAACTGGCGGCCAACCTGCGCCGGGCGTTTTCCGGGATTGTCGCAGGCAACGTCAAGGACAAAGGCATCCGCCTGATCGAGGAACACGGGCCGTATCAGATCCACGGTGACAGCGCCGTACTGGAc < 2:386496/150‑1 (MQ=255) gAGCTGCCGCCCCATGAACTGGCGGCCAACCTGCGCCGGGCGTTTTCCGGGATTGTCGCAGGCAACGTCAAGGACAAAGGCATCCGCCTGATCGAGGAACACGGGCCGTATCAGATCCACGGTGACAGCGCCGTACTGGACCCGCTGGGc < 1:246832/150‑1 (MQ=255) ccccATGAACTGGCGGCCAACCTGCGCCGGGCGTTTTCCGGGATTGTCGCAGGCAACGTCAAGGACAAAGGCATCCGCCTGATCGAGGAACACGGGCCGTATCAGATCCACGGTGACAGCGCCGTACTGGACCCGCTGGGCCGGctgctg > 2:127520/1‑150 (MQ=255) cATGAACTGGCGGCCAACCTGCGCCGGGCGTTTTCCGGGATTGTCGCAGGCAACGTCAAGGACAAAGGCATCCGCCTGATCGAGGAACACGGGCCGTATCAGATCCACGGTGACAGCGCCGTACTGGACCCGCTGAGCCGGCTGCTGCAg < 1:245853/150‑1 (MQ=255) gcggcCAACCTGCGCCGGGCGTTTTCCGGGATGGTCGCAGGCAACGTCAAGGACAAAGGCATCCGCCTGATCGAGGAACACGGGCCGTATCAGATCCACGGTGACAGCGCCGTACTGGACCCGCTGGGCCGGCTGCTGCAGGCATtcgtc > 2:23759/1‑150 (MQ=255) aaCCTGCGCCGGGCGTTTTCCGGGATTGTCGCAGGCAACGTCAAGGACAAAGGCATCCGCCTGATCGAGGAACACGGGCCGTATCAGATCCACGGTGACAGCGCCGTACTGGACCCGCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAg < 1:55085/150‑1 (MQ=255) ggCGTTTTCCGGGATTGTCGCAGGCAACGTCAAGGACAAAGGCATCCGCCTGATCGAGGAACACGGGCCGTATCAGATCCACGGTGACAGCGCCGTACTGGACCCGCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGaa < 1:334750/150‑1 (MQ=255) cGGGATTGTCGCAGGCAACGTCAAGGACAAAGGCATCCGCATGATCGAGGAACACGGGCCGTATCAGATCCACGGTGACAGAGCCGTACTGGACCCGCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGAGGAAGCTGTCcgg > 1:583499/1‑150 (MQ=255) gggATTGTCGCAGGCAACGTCAAGGACAAAGGCATCCACCTGATCGAGGAACACGGGCCGTATCAGATCCACGGTGACAGCGCCGTACTGGACCCGCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCcggc > 1:570964/1‑150 (MQ=255) tCGCAGGCAACGTCAAGGACAAAGGCATCCGCCTGATCGAGGAGCACGGGCCGTATCAGATCCACGGTGACAGCGCCGTACTGGACCCGCTGGGCCGgct > 1:260548/1‑100 (MQ=255) tCGCAGGCAACGTCAAGGACAAAGGCATCCGCCTGATCGAGGAGCACGGGCCGTATCAGATCCACGGTGACAGCGCCGTACTGGACCCGCTGGGCCGgct > 1:260540/1‑100 (MQ=255) tCGCAGGCAACGTCAAGGACAAAGGCATCCGCCTGATCGAGGAGCACGGGCCGTATCAGATCCACGGTGACAGCGCCGTACTGGACCCGCTGGGCCGgct < 2:260540/100‑1 (MQ=255) tCGCAGGCAACGTCAAGGACAAAGGCATCCGCCTGATCGAGGAGCACGGGCCGTATCAGATCCACGGTGACAGCGCCGTACTGGACCCGCTGGGCCGgct < 2:260548/100‑1 (MQ=255) tCGCAGGCAACGTCAAGGACAAAGGCATCCGCCTGATCGAGGAACACGGGCCGTATCAGATCCACGGTGACAGCGCCGTACTGGACCCGCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGgcggcgcgg > 1:151643/1‑150 (MQ=255) aaCGTCAAGGACAAAGGCATCCGCCTGATCGAGGAACACGGGCCGTATCAGATCCACGGTGACAGCGCCGTACTGGACCCGCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTa > 1:148066/1‑150 (MQ=255) aaCGTCAAGGACAAAGGCATCCGCCTGATCGAGGAACACGGGCCGTATCAGATCCACGGTGACAGCGCCGTACTGGACCCGCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTa > 1:331443/1‑150 (MQ=255) gTCAAGGACAAAGGCATCCGCCTGATCGAGGAACACGGGCCGTATCAGATCCACGGTGACAGCGCCGTACTGGACCCGCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCg > 1:2302/1‑150 (MQ=255) ggCATCCACCTGATCGAGGAACACGGGCCGTATCAGATCCACGGTGACAGCGCCGTACTGGACCCGCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTg < 2:570964/150‑1 (MQ=255) cTGATCGAGGAACACGGGCCGTATCAGATCCACGGTGACAGCGCCGTACTGGACCCGCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGa < 1:127520/150‑1 (MQ=255) gATCGAGGAACACGGGCCGTATCAGATCCACGGTGACAGCGCCGTACTGGACCCGCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCGATCATGGTGCCACCTGata < 2:151643/150‑3 (MQ=255) aTCGAGGAACACGGGCCGTATCAGATCCACGGTGACAGCGCCGTACTGGACCCGCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGAccc > 2:85289/1‑150 (MQ=255) gAGGAACACGGGCCGTATCAGATCCACGGTGACAGCGCCGTACTGGACCCGCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCgct < 2:55119/150‑1 (MQ=255) aGGAACACGGGCCGTATCAGATCCACGGTGACAGCGCCGTACTGGACCCGCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCgctg > 1:322774/1‑150 (MQ=255) | CGTTCGACCCGACCCACCAGGCCATGGCCGACCTGGACTTGCGCCGCGAGCTGCCGCCCCATGAACTGGCGGCCAACCTGCGCCGGGCGTTTTCCGGGATTGTCGCAGGCAACGTCAAGGACAAAGGCATCCGCCTGATCGAGGAACACGGGCCGTATCAGATCCACGGTGACAGCGCCGTACTGGACCCGCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTG > NC_002947/4160395‑4160685 |
| 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: ‑ |