| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 4,558,839 | G→T | 24.9% | P360P (CCC→CCA) | PP_4041 ← | glycoside hydrolase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 4,558,839 | 0 | G | T | 24.9% | 44.8 / 5.3 | 20 | P360P (CCC→CCA) | PP_4041 | glycoside hydrolase |
| Reads supporting (aligned to +/- strand): ref base G (9/6); new base T (3/2); total (12/8) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.35e-01 | |||||||||||
AGTGCACGAACTGCTGCCGGCCGCCACGCACTTCCCAGATGCCTTCGTCCGGCTCGCGCCAGATACGTTCGACATAAGGCAATATCAAGCGCGAAATGGCGGCGCTGCGCGGGTGGCGGGGCAACCCGCCCTTGATCGCCTGGGTCATCGCATCGGCCAGTTCGCCATAGATGTCCAGTTGCCTCTGCGCCGAGGCGGCGTTACCCACCCGCACAGGCTGAGAATGCTCATAACCGGCCAGCCACGGCAAGGTG > NC_002947/4558721‑4558974 | aGTGCACGAACTGCTGCCGGCCGCCACGCACTTCCCAGATGCCTTCGTCCGGCTCGCGCCAGATACGTTCGACATAAGGCAATATCAAGCGCGAAATGGCGGCGCTGCGCGGGTGGCGGGGCAACCCGCCCTTGAt > 1:36766/1‑136 (MQ=255) gAACTGCTGCCGGCCGCCACGCACTTCCCAGATGCCTTCGTCCGGCTCGCGCCAGATACGTTCGACATAAGGCAATATCAAGCGCGAAATGGCGGCGCTGCGCGGGTGGCGGGGCAACCCGCCCTTGATCGCCTgg > 1:39776/1‑136 (MQ=255) cccAGATGCCTTCGTCCGGCTCGCGCCAGATACGTTCGACATAAGGCAATATCAAGCGCGAAATGGCGGCGCTGCGCGGGTGGCGGGGCAACCCGCCCTTGATCGCCTGGGTCATCGCATCGGCCAGTTCGCCATa < 2:39776/136‑1 (MQ=255) cAGATACGTTCGACATAAGGCAATATCAAGCGCGAAATGGCGGCGCTGCGCGGGTGGCGGGGCAACCCGCCCTTGATCGCCTGGGTCATCGCATCGGCCAGTTCGCCATAGATGTCCAGTTGCCTCTGCGCCGAgg > 2:50891/1‑136 (MQ=255) gACATAAGGCAATATCAAGCGCGAAATGGCGGCGCTGCGCGGGTGGCGTGGCAACCCGCCCTTGATCGCCTGGGTCATCGCATCGGCCAGTTCGCCATAGATGTCCAGTTGCCTCTGCGCCGAGGCGGCGTTaccc > 1:540487/1‑136 (MQ=255) gACATAAGGCAATATCAAGCGCGAAATGGCGGCGCTGCGCGGGTGGCGGGGCAACCCGCCCTTGATCGCCTGGGTCATCGCATCGGCCAGTTCGCCATAGATGTCCAGTTGCCTCTGCGCCGAGGCGGCGTTaccc > 1:519354/1‑136 (MQ=255) tAAGGCAATATCAAGCGCGAAATGGCGGCGCTGCGCGGGTGGCGTGGCAACCCGCCCTTGATCGCCTGGGTCATCGCATCGGCCAGTTCGCCATAGATGTCCAGtt > 1:301034/1‑106 (MQ=255) tAAGGCAATATCAAGCGCGAAATGGCGGCGCTGCGCGGGTGGCGTGGCAACCCGCCCTTGATCGCCTGGGTCATCGCATCGGCCAGTTCGCCATAGATGTCCAGtt < 2:5331/106‑1 (MQ=255) tAAGGCAATATCAAGCGCGAAATGGCGGCGCTGCGCGGGTGGCGTGGCAACCCGCCCTTGATCGCCTGGGTCATCGCATCGGCCAGTTCGCCATAGATGTCCAGtt < 2:301034/106‑1 (MQ=255) tAAGGCAATATCAAGCGCGAAATGGCGGCGCTGCGCGGGTGGCGTGGCAACCCGCCCTTGATCGCCTGGGTCATCGCATCGGCCAGTTCGCCATAGATGTCCAGtt > 1:5331/1‑106 (MQ=255) tAAGGCAATATCAAGCGCGAAATGGCGGCGCTGCGCGGGTGGCGGGGCAACCCGCCCTTGATCGCCTGGGTCATCGCATCGGCCAGTTCGCCATAGATGTCCAGtt > 1:438571/1‑106 (MQ=255) tAAGGCAATATCAAGCGCGAAATGGCGGCGCTGCGCGGGTGGCGGGGCAACCCGCCCTTGATCGCCTGGGTCATCGCATCGGCCAGTTCGCCATAGATGTCCAGtt < 2:438571/106‑1 (MQ=255) tAAGGCAATATCAAGCGCGAAATGGCGGCGCTGCGCGGGTGGCGGGGCAACCCGCCCTTGATCGCCTGGGTCATCGCATCGGCCAGTTCGCCATAGATGTCCAGTTGCCTCTGCGCCGAGGCGGCGTTacccaccc < 1:435861/136‑1 (MQ=255) tAAGGCAATATCAAGCGCGAAATGGCGGCGCTGCGCGGGTGGCGGGGCAACCCGCCCTTGATCGCCTGGGTCATCGCATCGGCCAGTTCGCCATAGATGTCCAGTTGCCTCTGCGCCGAGGCGGCGTTacccaccc < 2:377268/136‑1 (MQ=255) tCAAGCGCGAAATGGCGGCGCTGCGCGGGTGGCGGGGCAACCCGCCCTTGATCGCCTGGGTCATCGCATCGGCCAGTTCGCCATAGATGTCCAGTTGCCTCTGCGCCGAGGCGGCGTTACCCACCCGCACAGGCTg > 1:265178/1‑136 (MQ=255) cAAGCGCGAAATGGCGGCGCTGCGCGGGTGGCGGGGCAACCCGCCCTTGATCGCCTGGGTCATCGCATCGGCCAGTTCGCCATAGATGTCCAGTTGCCTCTGCGCCGAGGCGGCGTTACCCACCCGCACAGGCTga > 1:473064/1‑136 (MQ=255) cAAGCGCGAAATGGCGGCGCTGCGCGGGTGGCGGGGCAACCCGCCCTTGATCGCCTGGGTCATCGCATCGGCCAGTTCGCCATAGATGTCCAGTTGCCTCTGCGCCGAGGCGGCGTTACCCACCCGCACAGGCTga > 1:541744/1‑136 (MQ=255) gAAATGGCGGCGCTGCGCGGGTGGCGGGGCACCCCGCCCTTGATCGCCTGGGTCATCGCATCGGCCAGTTCGCCATAGATGTCCAGTTGCCTCTGCGCCGAGGCGGCGTTACCCACCCGCACAGGCTGAGAATGCt > 1:480207/1‑136 (MQ=255) tGGCGGGGCAACCCGCCCTTGATCGCCTGGGTCATCGCATCGGCCAGTTCGCCATAGATGTCCAGTTGCCTCTGCGCCGAGGCGGCGTTACCCACCCGCACAGGCTGAGAATGCTCATAACCGGCCAGCCACGGCa < 1:386521/136‑1 (MQ=255) gCGGGGCAACCCGCCCTTGATCGCCTGGGTCATCGCATCGGCCAGTTCGCCATAGATGTCCAGTTGCCTCTGCGCCGAGGCGGCGTTACCCACCCGCACAGGCTGAGAATGCTCATAACCGGCCAGCCACGGCAAg < 2:541744/136‑1 (MQ=255) gggCAACCCGCCCTTGAGCGCCTGGGTCATCGCATCGGCCAGTTCGCCATAGATGTCCAGTTGCCTCTGCGCCGAGGCGGCGTTACCCACCCGCACAGGCTGAGAATGCTCATAACCGGCCAGCCACGGCAAGgtg < 1:323843/136‑1 (MQ=255) | AGTGCACGAACTGCTGCCGGCCGCCACGCACTTCCCAGATGCCTTCGTCCGGCTCGCGCCAGATACGTTCGACATAAGGCAATATCAAGCGCGAAATGGCGGCGCTGCGCGGGTGGCGGGGCAACCCGCCCTTGATCGCCTGGGTCATCGCATCGGCCAGTTCGCCATAGATGTCCAGTTGCCTCTGCGCCGAGGCGGCGTTACCCACCCGCACAGGCTGAGAATGCTCATAACCGGCCAGCCACGGCAAGGTG > NC_002947/4558721‑4558974 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |