| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 2,982,931 | C→A | 21.0% | E97D (GAG→GAT) ‡ | 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,931 | 0 | C | A | 21.0% | 40.5 / 2.6 | 19 | E97D (GAG→GAT) ‡ | PP_2608 | shikimate dehydrogenase family protein |
| Reads supporting (aligned to +/- strand): ref base C (7/8); new base A (2/2); total (9/10) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
CCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCTCGAAGCCATGCTTGTGCGCAGCGCCGAGGAAGCCGGCGCCGTCCACATTGTCGCCCAGCAGGCGTCCGTCGCGCTCGCGACGTATCACATTGATCGAGCCCAGGGCAGCCGCGCGCTCGCTCAGGCCGTCCACGCGGTTGGCCAGCGCCTGCTTGTAAGGCACGGTGACCACGCAACCGCGCAGGTTCTGCCAGCCACGCAGGGTATCCGCAAAG > NC_002947/2982789‑2983072 | ccAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCTCGAAGCCATGCTTGTGCGCAGCGCCGAGGAAGCCGGCGCCGTCCACATTGTCGCCCAGCAGGCGTCCGTCGCGCTCGCGAc > 1:251677/1‑150 (MQ=255) aTAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCTCGAAGCCATGCTTGTGCGCAGCGCCGAGGAAGCCGGCGCCGTCCACATTGTCGCCCAGCAGGCGTCCGTCGCGAGCGCGACgga > 1:349500/1‑144 (MQ=255) aTAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCTCGAAGCCATGCTTGTGCGCAGCGCCGAGGAAGCCGGCGCCGTCCACATTGTCGCCCAGCAGGCGTCCGTCGCGAGCGCGACgga < 2:349500/146‑3 (MQ=255) ggCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCTCGAAGCCATGCTTGTGCGCAGCGCCGAGGAAGCCGGCGCCGTCCACATTGTCGCCCAGCAGGCGTCCGTCGCGCTCGCGACGTATCACATTGATCGa < 2:378595/150‑1 (MQ=255) gCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCATTTGCCCGCCGGCTCGAAGCCATGCTTGTGCGCAGCGCCGAGGAAGCCGGCGCCGTCCACATTGTCGCCCAGCAGGCGTACGTCGCGAGCGCGACggac < 2:444119/137‑4 (MQ=255) gCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCATTTGCCCGCCGGCTCGAAGCCATGCTTGTGCGCAGCGCCGAGGAAGCCGGCGCCGTCCACATTGTCGCCCAGCAGGCGTACGTCGCGAGCGCGACggac > 1:444119/1‑134 (MQ=255) tGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCTCGAAGCCATGCTTGTGCGCAGCGCCGAGGAAGCCGGCGCCGTCCACATCGTCGCCCAGCAGGCGTCCGTCGCGCTCGCGACGTATCACATTGATCGAGCCCa > 1:382300/1‑150 (MQ=255) acacCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCTCGAAGCCATGCTTGTGCGCAGCGCCGAGGAAGCCGGCGCCGTCCACATTGTCGCCCAGCAGGCGTCCGTCGCGCTCGCGACGTATCACATTGATCGAGCCCAGGGCa > 2:380942/1‑150 (MQ=255) aCCAGTGCCCGTTTGCCCGCCGGCTCGAAGCCATGCTTGTGCGCAGCGCCGAGGAAGCCGGCGCCGTCCACATTGTCGCCCAGCAGGCGTCCGTCGCGCTCGCGACGTATCACATTGATCGAGCCCAGGGCAGCCGCGCGCTCGCTCAgg < 1:112212/150‑1 (MQ=255) gTTTGCCCGCCGGCTCGAAGCCATGCTTGTGCGCAGCGCCGAGGAAGCCGGCGCCGTCCACATTGTCGCCCAGCAGGCGTCCGTCGCGCt < 1:554061/90‑1 (MQ=255) gTTTGCCCGCCGGCTCGAAGCCATGCTTGTGCGCAGCGCCGAGGAAGCCGGCGCCGTCCACATTGTCGCCCAGCAGGCGTCCGTCGCGCt > 2:554061/1‑90 (MQ=255) tttGCCCGCCGGCTCGAAGCCATGCTTGTGCGCAGCGCCGAGGAAGCCGGCGCCGTCCACATTGTCGCCCAGCAGGCGTCCGTCGCGCTCGCGACGTATCACATTGATCGAGCCCAGGGCAGCCGCGCGCTCGCTCAGGCCGTCCAcgcg > 1:491619/1‑150 (MQ=255) cGAAGCCATGCTTGTGCGCAGCGCCGAGGAAGCCGGCGCCGTCCCCATTGTCGCCCAGCAGGCGTCCGTCGCGCTCGCGACGTATCACATTGATCGAGCCCAGGGCAGCCGCGCGCTCGCTCAGGCCGTCCACGCGGTTGGCCAGCGCCt < 1:204554/150‑1 (MQ=255) gCTTGTGCGCAGCGCCGAGGAAGCCGGCGCCGTCCACATTGTCGCCCAGCAGGCGTCCGTCGCGCTCGCGACGTATCACATTGATCGAGCCCAGGGCAGCCGCGCGCTCGCTCAGGCCGTCCACGCGGTTGGCCAGCGCCTGCTTGTAAg < 2:251677/150‑1 (MQ=255) cgcAGCGCCGAGGAAGCCGGCGCCGTCCACATTGTCGCCCAGCAGGCGTCCGTCGCGCTCGCGACGTATCACATTGATCGAGCCCAGGGCAGCCGCGCGCTCGCTCAGGCCGTCCACGCGGTTGGCCAGCGCCTGCTTGTAAGGCACGGt < 2:491619/150‑1 (MQ=255) gcgcCGTCCACATTGTCGCCCAGCAGGCGTCCGTCGCGCTCGCGACGTATCACATTGATCGAGCCCAGGGCAGCCGCGCGCTCGCTCAgg < 2:184955/90‑1 (MQ=255) gcgcCGTCCACATTGTCGCCCAGCAGGCGTCCGTCGCGCTCGCGACGTATCACATTGATCGAGCCCAGGGCAGCCGCGCGCTCGCTCAgg > 1:184955/1‑90 (MQ=255) gtccgtcGCGCTCGCGACGTATCACATTGATCGAGCCCAGGGCAGCCGCGCGCTCGCTCAGGCCGTCCACGCGGTTGGCCAGCGCCTGCTTGTAAGGCACGGTGACCACGCAACCGCGCAGGTTCTGCCAGCCACGCAGGGTATCCGCaa > 1:379901/1‑150 (MQ=255) ccgtcGCGCTCGCGACGTATCACATTGATCGAGCCCAGGGCAGCCGCGCGCTCGCTCAGGCCGTCCACGCGGTTGGCCAGCGCCTGCTTGTAAGGCACGGTGACCACGCAACCGCGCAGGTTCTGCCAGCCACGCAGGGTATCCGCAAAg < 2:382300/150‑1 (MQ=255) | CCAGTGCATAGGCAATGGCACTGCCAACACCGCCGCAACCGATCACCAGTGCCCGTTTGCCCGCCGGCTCGAAGCCATGCTTGTGCGCAGCGCCGAGGAAGCCGGCGCCGTCCACATTGTCGCCCAGCAGGCGTCCGTCGCGCTCGCGACGTATCACATTGATCGAGCCCAGGGCAGCCGCGCGCTCGCTCAGGCCGTCCACGCGGTTGGCCAGCGCCTGCTTGTAAGGCACGGTGACCACGCAACCGCGCAGGTTCTGCCAGCCACGCAGGGTATCCGCAAAG > NC_002947/2982789‑2983072 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |