| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NZ_CP009273 | 1,620,338 | A→G | S434P (TCA→CCA) | dcp ← | peptidyl‑dipeptidase Dcp |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 1,620,338 | 0 | A | G | 90.5% | 57.7 / ‑0.9 | 21 | S434P (TCA→CCA) | dcp | peptidyl‑dipeptidase Dcp |
| Reads supporting (aligned to +/- strand): ref base A (1/1); new base G (12/7); total (13/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 = 9.61e-01 | |||||||||||
ACCAAATTCATGGAATAAGGTTATGACATCATCCCAGAGTAACAACGCAGGCTCACCGGCAGCGGGTTTCTGATAATTGCAGACGTTATAAATTACCGGATGTGTTTTATTAAGCGTTGATTGCTCAACAAAATTGCCCATCCATGCACCGCCGCTTTTTGAATCACGGGCGAAGAAATCACCGTAAAATAACGCCAGCCCCACGCCATTATGATCAAAAATTTCCCACACACGAACGTCAGGATGGTAG > NZ_CP009273/1620219‑1620468 | aCCAAATTCATGGAATAAGGTTATGACATCATCCCAGAGTAACAACGCAGGCTCACCGGCAGCGGGTTTCTGATAATTGCAGACGTTATAAATTACCGGATGTGTTTTATTAAGCGTTGGTTGCTCAACAAAATTGcc < 1:244994/138‑1 (MQ=255) tGACATCATCCCAGAGTAACAACGCAGGCTCACCGGCAGCGGGTTTCTGATAATTGCAGACGTTATAAATTACCGGATGTGTTTTATTAAGCGTTGGTTGCt < 1:159050/102‑1 (MQ=255) tGACATCATCCCAGAGTAACAACGCAGGCTCACCGGCAGCGGGTTTCTGATAATTGCAGACGTTATAAATTACCGGATGTGTTTTATTAAGCGTTGGTTGCt > 2:159050/1‑102 (MQ=255) tGACATCATCCCAGAGTAACAACGCAGGCTCACCGGCAGCGGGTTTCTGATAATTGCAGACGTTATAAATTACCGGATGTGTTTTATTAAGCGTTGGTTGCTCAACAAAATTGcc < 2:123484/115‑1 (MQ=255) tGACATCATCCCAGAGTAACAACGCAGGCTCACCGGCAGCGGGTTTCTGATAATTGCAGACGTTATAAATTACCGGATGTGTTTTATTAAGCGTTGGTTGCTCAACAAAATTGcc > 1:123484/1‑115 (MQ=255) tGACATCATCCCAGAGTAACAACGCAGGCTCACCGGCAGCGGGTTTCTGATAATTGCAGACGTTATAAATTACCGGATGTGTTTTATTAAGCGTTGGTTGCTCAACAAAATTGCCCATCCATGCACCGCCGCTTTTTg > 2:399192/1‑138 (MQ=255) atcCCAGAGTAACAACGCAGGCTCACCGGCAGCGGGTTTCTGATAATTGCAGACGTTATAAATTACCGGATGTGTTTTATTAAGCGTTGGTTGCTCAACAAAATTGCCCATCCATGCACCGCCGCTTTTTGAATCACgg > 1:8053/1‑139 (MQ=255) cGCAGGCTCACCGGCAGCGGGTTTCTGATAATTGCAGACGTTATAAATTACCGGATGTGTTTTATTAAGCGTTGGTTGCTCAACAAAATTGCCCATCCATGCACCGCCGCTTTTTGAATCACGGGCGAAGAAATCACCg > 2:96363/1‑139 (MQ=255) gCTCACCGGCAGCGGGTTTCTGATAATTGCAGACGTTATAAATTACCGGATGTGTTTTATTAAGCGTTGGTTGCTCAACAAAATTGCCCATCCATGCACCGCCGCTTTTTGAATCACGGGCGAAGAAATCACCGTaaaa > 2:347670/1‑139 (MQ=255) ccGGCAGCGGGTTTCTGATAATTGCAGACGTTATAAATTACCGGATGTGTTTTATTAAGCGTTGGTTGCTCAACAAAATTGCCCATCCATGCACCGCCGCTTTTTGAATCACGGGCGAAGAAATCACCGTAAAATAACg > 2:639724/1‑139 (MQ=255) cGGGTTTCTGATAATTGCAGACGTTATAAATTACCGGATGTGTTTTATTAAGCgttggttg < 1:657477/61‑1 (MQ=255) cGGGTTTCTGATAATTGCAGACGTTATAAATTACCGGATGTGTTTTATTAAGCgttggttg > 2:657477/1‑61 (MQ=255) ggTTTCTGATAATTGCAGACGTTATAAATTACCGGATGTGTTTTATTAAGCGTTGATTGCTCAACAAAATTGCCCATCCATGCACCGCCGCTTTTTGAATCACGGGCGAAGAAATCACCg > 1:406411/1‑120 (MQ=255) ggTTTCTGATAATTGCAGACGTTATAAATTACCGGATGTGTTTTATTAAGCGTTGATTGCTCAACAAAATTGCCCATCCATGCACCGCCGCTTTTTGAATCACGGGCGAAGAAATCACCg < 2:406411/120‑1 (MQ=255) aCCGGATGTGTTTTATTAAGCGTTGGTTGCTCAACAAAATTGCCCATCCATGCACCGCCGCTTTTTGAATCACGGGCGAAGAAATCACCGTAAAATAACGCCAGCCCCACGCCATTATGATCAAAAATTTCCCacacac > 1:20256/1‑139 (MQ=255) cGGATGTGTTTTATTAAGCGTTGGTTGCTCAACAAAATTGCCCATCCATGCACCGCCGCTTTTTGAATCACGGGCga > 2:409866/1‑77 (MQ=255) cGGATGTGTTTTATTAAGCGTTGGTTGCTCAACAAAATTGCCCATCCATGCACCGCCGCTTTTTGAATCACGGGCga < 1:409866/77‑1 (MQ=255) attaAGCGTTGGTTGCTCAACAAAATTGCCCATCCATGCACCGCCGCTTTTTGAATCACGGGCGAAGAAATCACCGTAAAATAACGCCAGCCCCACGCCATTATGATCAAAAATTTCCCACACACGAACGTCAGGa > 1:627005/1‑136 (MQ=255) attaAGCGTTGGTTGCTCAACAAAATTGCCCATCCATGCACCGCCGCTTTTTGAATCACGGGCGAAGAAATCACCGTAAAATAACGCCAGCCCCACGCCATTATGATCAAAAATTTCCCACACACGAACGTCAGGa < 2:627005/136‑1 (MQ=255) aGCGTTGGTTGCTCAACAAAATTGCCCATCCATGCACCGCCGCTTTTTGAATCACGGGCGAAGAAATCACCGTAAAATAACGCCAGCCCCACGCCATTATGATCAAAAATTTCCCACACACGAACGTCAGGATGGTAg < 1:294597/138‑1 (MQ=255) aGCGTTGGTTGCTCAACAAAATTGCCCATCCATGCACCGCCGCTTTTTGAATCACGGGCGAAGAAATCACCGTAAAATAACGCCAGCCCCACGCCATTATGATCAAAAATTTCCCACACACGAACGTCAGGATGGTAg > 2:294597/1‑138 (MQ=255) cGTTGGTTGCTCAACAAAATTGCCCATCCATGCACCGCCGCTTTTTGAATCACGGGCGAAGAAATCACCGTAAAATAACGCCAGCCCCACGCCATTATGATCAAAAATTTCCCACACACGAACGt > 2:791450/1‑125 (MQ=255) cGTTGGTTGCTCAACAAAATTGCCCATCCAGGCACCGCCGCTTTTTGAATCACGGGCGAAGAAATCACCGTAAAATAACGCCAGCCCCACGCCATTATGATCAAAAATTTCCCACACACGAACGt < 1:791450/125‑1 (MQ=255) | ACCAAATTCATGGAATAAGGTTATGACATCATCCCAGAGTAACAACGCAGGCTCACCGGCAGCGGGTTTCTGATAATTGCAGACGTTATAAATTACCGGATGTGTTTTATTAAGCGTTGATTGCTCAACAAAATTGCCCATCCATGCACCGCCGCTTTTTGAATCACGGGCGAAGAAATCACCGTAAAATAACGCCAGCCCCACGCCATTATGATCAAAAATTTCCCACACACGAACGTCAGGATGGTAG > NZ_CP009273/1620219‑1620468 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |