| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 867,319 | A→G | Y251C (TAC→TGC) | iaaA → | Isoaspartyl peptidase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 867,319 | 0 | A | G | 90.5% | 60.5 / ‑1.1 | 21 | Y251C (TAC→TGC) | iaaA | Isoaspartyl peptidase |
| Reads supporting (aligned to +/- strand): ref base A (1/1); new base G (13/6); total (14/7) | |||||||||||
| 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.12e-01 | |||||||||||
TGGCGATAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTACGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAATGTCGCGCTACCGTTTAACACCG > NC_000913/867191‑867444 | tGGCGATAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTGCGGCGGATTa < 2:773528/139‑1 (MQ=255) cTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTGCGGCGGATTAAGTCTCGCGGaa > 1:194668/1‑139 (MQ=255) ccGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTGCGGCGGATTAAGTCTCGCgg < 2:435629/126‑1 (MQ=255) ccGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTGCGGCGGATTAAGTCTCGCgg > 1:435629/1‑126 (MQ=255) gATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCGGgt > 2:681666/1‑139 (MQ=255) aTAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGaaaa > 1:328061/1‑139 (MQ=255) tgtgGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTgcg > 2:123914/1‑139 (MQ=255) tgGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCt > 1:272671/1‑139 (MQ=255) gCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCggg > 1:331023/1‑96 (MQ=255) gCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCggg < 2:331023/96‑1 (MQ=255) gCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCtt > 1:125609/1‑139 (MQ=255) gCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCtt > 2:237302/1‑139 (MQ=255) tCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGa < 2:125609/139‑1 (MQ=255) cATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCggg > 2:107554/1‑80 (MQ=255) cATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCggg < 1:107554/80‑1 (MQ=255) cgccgcGTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAATGTCGCGCTACCGTTTAac > 2:655072/1‑139 (MQ=255) cgccgcGTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAATGTCGCGCTACCGTTTAac > 2:696202/1‑139 (MQ=255) gccgcGTTAATGGATTACGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCt > 2:162244/1‑102 (MQ=255) gccgcGTTAATGGATTACGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCt < 1:162244/102‑1 (MQ=255) ccgcGTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAATGTCGCGCTACCGTTTAacac > 1:438808/1‑139 (MQ=255) gcgTTAATGGATTGCGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAATGTCGCGCTACCGTTTAACACCg < 1:237302/139‑1 (MQ=255) | TGGCGATAGTCCCTTAGTGGGTGCCGGATGCTACGCCAATAACGCCAGTGTGGCGGTTTCTTGTACCGGCACGGGCGAAGTCTTCATCCGCGCGCTGGCGGCATATGACATCGCCGCGTTAATGGATTACGGCGGATTAAGTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAATGTCGCGCTACCGTTTAACACCG > NC_000913/867191‑867444 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |