| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 135,481 | A→G | G34G (GGT→GGC) | speD ← | S‑adenosylmethionine decarboxylase proenzyme |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 135,481 | 0 | A | G | 96.0% | 77.4 / ‑4.0 | 25 | G34G (GGT→GGC) | speD | S‑adenosylmethionine decarboxylase proenzyme |
| Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (14/10); total (15/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 = 8.49e-01 | |||||||||||
TCCTGGCGGGCGATGTTAAGAATATTAGCCCCGATAATGGAACAGGTTTCTGACAGGATTTCGGTCAGACGGTTGGCATTATAGAGTTCATCGATATAAGCAATATAACCGTCGCGCTCTTCGGCAGTTTTGGCGTAGCAGATATCGTAAATACAAAAACTCAGACTTTTGGTCAGATTATTAAAGCCATGCAGTTTCAGTTTTTTCAATTTCTTATCTTCTCCTTAGGACGGCTG > NC_000913/135374‑135609 | tCCTGGCGGGCGATGTTAAGAATATTAGCCCCGATAATGGAACAGGTTTCTGACAGGATTTCGGTCAGACGGTTGGCATTATAGAGTTCATCGATATAAGCAATATAGCCGTCGCGCTCTTCGGCAGTTTTgg > 1:8630/1‑133 (MQ=255) ggcgATGTCAAGAATATTAGCCCCGATAATGGAACAGGTTTCTGACAGGATTTCGGTCAGACGGTTGGCATTATAGAGTTCATCGATATAAGCAATATAGCCGTCGCGCTCTTCGGCAGTTTTGGCGTAGCAg < 1:219936/133‑1 (MQ=255) atTAGCCCCGATAATGGAACAGGTTTCTGACAGGATTTCGGTCAGACGGTTGGCATTATAGAGTTCATCGATATAAGCAATATAGCCGTCGCGCTCTTCGGCAGTTTTGGCGTAGCAGATATCGTAAATACaa > 2:6885/1‑133 (MQ=255) tAGCCCCGATAATGGAACAGGTTTCTGACAGGATTTCGGTCAGACGGTTGGCATTATAGAGTTCATCGATATAAGCAATATAGCCGTCGCGCTCTTCGGCAGTTTTGGCGTAGCAGATATCGTAAATACaaaa > 2:377184/1‑133 (MQ=255) tAGCCCCGATAATGGAACAGGTTTCTGACAGGATTTCGGTCAGACGGTTGGCATTATAGAGTTCATCGATATAAGCAATATAGCCGTCGCGCTCTTCGGCAGTTTTGGCGTAGCAGATATCGTAAATACaaaa > 2:403513/1‑133 (MQ=255) tAGCCCCGATAATGGAACAGGTTTCTGACAGGATTTCGGTCAGACGGTTGGCATTATAGAGTTCATCGATATAAGCAATATAGCCGTCGCGCTCTTCGGCAGGTTTGGCGGAGCGGATATCGTAGATACaaaa > 1:129744/1‑133 (MQ=255) cagATAATGCAACAGATTTCTGATAGTACTTCCGTCAGACGTTTGGCATTATAGAGTTCATCGATATAAGCAATATAGCCGGCGCGCTCTTCGGCAGTTTTGGCGTAGCAGATATCGTAAATACAAAAACTCa < 2:246842/131‑1 (MQ=255) gATAATGGAACAGGTTTCTGACAGGATTTCGGTCAGACGGTTGGCATTATAGAGTTCATCGATATAAGCAATATAGCCGTCGCGCTCTTCGGCAGTTTTGGCGTAGCAGATATCGTAAATACAAAAACTCAGa > 1:9073/1‑133 (MQ=255) aTAATGGAACAGGTTTCTGACAGGATTTCGGTCAGACGGTTGGCATTATAGAGTTCATCGATATAAGCAATATAGCCGTCGCGCTCTTCGGCAGTTTTGGCGTAGCAGATATCGTAAATACAAAAACTCAGAc > 2:370995/1‑133 (MQ=255) aaCAGGTTTCTGACAGGATTTCGGTCAGACGGTTGGCATTATAGAGTTCATCGATATAAGCAATATAGCCGTCGCGCTCTTCGGCAGTTTTGGCGTAGCAGATATCGTAAATACAAAAACTCAGACTTTTGGt < 1:403513/133‑1 (MQ=255) aaCAGGTTTCTGACAGGATTTCGGTCAGACGGTTGGCATTATAGAGTTCATCGATATAAGCAATATAGCCGTCGCGCTCTTCGGCAGTTTTGGCGTAGCAGATATCGTAAATACAAAAACTCAGACTTTTGGt < 1:370995/133‑1 (MQ=255) aaCAGGTTTCTGACAGGATTTCGGTCAGACGGTTGGCATTATAGAGTTCATCGATATAAGCAATATAGCCGTCGCGCTCTTCGGCAGTTTTGGCGTAGCAGATATCGTAAATACAAAAACTCAGACTTTTGGt < 1:6885/133‑1 (MQ=255) aGGTTTCTGACAGGATTTCGGTCAGACGGTTGGCATTATAGAGTTCATCGATATAAGCAATATAGCCGTCGCGCTCTTCGGCAGTTTTGGCGTAGCAGATATCGTAAATACAAAAACTCAGACTTTTGGTCAg < 2:9073/133‑1 (MQ=255) aGGTTTCTGACAGGATTTCGGTCAGACGGTTGGCATTATAGAGTTCATCGATATAAGCAATATAGCCGTCGCGCTCTTCGGCAGTTTTGGCGTAGCAGATATCGTAAATACAAAAACTCAGACTTTTGGTCAg < 1:7939/133‑1 (MQ=255) gATTTCGGTCAGACGGTTGGCATTATAGAGTTCATCGATATAAGCAATATAGCCGTCGCGCTCTTCGGCAGTTTTGGCGTAGCAGATATCGTAAATACAAAAACTCAGACTTTTGGTCAGATTATTAAAGCCa > 2:323740/1‑133 (MQ=255) aTTTCGGTCAGACGGTTGGCATTATAGAGTTCATCGATATAAGCAATATAGCCGTCGCGCTCTTCGGCAGTTTTGGCGTAGCAGATATCGTTAATACAAAAACTAAGACTTTTGGTCCGATTATTAAAGCCAt > 2:219895/1‑133 (MQ=255) aTTTCGGTCAGACGGTTGACATTATAGAGTTCATCGATATAAGCAATATAGCCGTCGCGCTCTTCGGCAGTTTTGGCGTAGCAGATATCGTAAATACAAAAACTCAGACTTTTGGTCAGATTATTAAAGCCAt < 1:219895/133‑1 (MQ=255) ttGGCATTATAGAGTTCATCGATATAAGCAATATAACCGTCGCGCTCTTCGGCAGTTTTGGCGTAGCAGATATCGTAAATACAAAAACTCAGACTTTTGGTCAGATTATTAAAGCCATGCAGTTTCAGttttt > 1:370024/1‑133 (MQ=255) tCATCGATATCAGCATTATAGCCGCCGCCCTCGTCGGCAGTCTTGGCGTAGCAGATATCGTAAATACAAAAACTCAGACTTTTGGTCAGATTATTAAAGCCATGCAGTTTCAGTTTTTTCAATTTCTTAtctt < 2:129744/133‑1 (MQ=255) aTCGATATAAGCAATATAGCCGTCGCGCTCTTCGGCAGTTTTGGCGTAGCAGATATCGTAAATACAAAAACTCAGACTTTTGGTCAGATTATTAAAGCCATGCAGTTTCAGTTTTTTCAATTTCTTAtcttct > 1:356851/1‑133 (MQ=255) cGATATAAGCAATATAGCCGTCGCGCTCTTCGGCAGTTTTGGCGTAGCAGATATCGTAAATACAAAAACTCAGACTTTTGGTCAGATTATTAAAGCCATGCAGTTTCAGTTTTTTCAATTTCTTATCTTCTcc > 2:26441/1‑133 (MQ=255) cGATATAAGCAATATAGCCGTCGCGCTCTTCGGCAGTTTTGGCGTAGCAGATATCGTAAATACAAAAACTCAGACTTTTGGTCAGATTATTAAAGCCATGCAGTTTCAGTTTTTTCAATTTCTTATCTTCTcc > 1:409355/1‑133 (MQ=255) tataAGCAATATAGCCGTCGCGCTCTTCGGCAGTTTTGGCGTAGCAGATATCGTAAATACAAAAACTCAGACTTTTGGTCAGATTATTAAAGCCATGCAGTTTCAGTTTTTTCAATTTCTTATCTTCTCCTTa > 1:96479/1‑133 (MQ=255) tataGCCGTCGCGCTCTTCGGCAGTTTTGGCGTAGCAGATATCGTAAATACAAAAACTCAGACTTTTGGTCAGATTATTAAAGCCATGCAGTTTCAGTTTTTTCAATTTCTTATCTTCTCCTTAGGACGGCtg > 1:44077/1‑133 (MQ=255) tataGCCGTCGCGCTCTTCGGCAGTTTTGGCGTAGCAGATATCGTAAATACAAAAACTCAGACTTTTGGTCAGATTATTAAAGCCATGCAGTTTCAGTTTTTTCAATTTCTTATCTTCTCCTTAGGACGGCtg < 2:96479/133‑1 (MQ=255) | TCCTGGCGGGCGATGTTAAGAATATTAGCCCCGATAATGGAACAGGTTTCTGACAGGATTTCGGTCAGACGGTTGGCATTATAGAGTTCATCGATATAAGCAATATAACCGTCGCGCTCTTCGGCAGTTTTGGCGTAGCAGATATCGTAAATACAAAAACTCAGACTTTTGGTCAGATTATTAAAGCCATGCAGTTTCAGTTTTTTCAATTTCTTATCTTCTCCTTAGGACGGCTG > NC_000913/135374‑135609 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |