| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | W3110S.gb | 3,047,650 | T→C | 30.0% | E16E (GAA→GAG) | gcvP ← | glycine decarboxylase, PLP‑dependent, subunit (protein P) of glycine cleavage complex |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | W3110S.gb | 3,047,650 | 0 | T | C | 30.0% | 49.4 / 30.6 | 40 | E16E (GAA→GAG) | gcvP | glycine decarboxylase, PLP‑dependent, subunit (protein P) of glycine cleavage complex |
| Reads supporting (aligned to +/- strand): ref base T (11/17); new base C (7/5); total (18/22) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.15e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.62e-01 | |||||||||||
GTTTAACGATTGTGCACCAACGGCATTCAGCATTTCTTGCTGTTGCGCGGCGTCCGGTCCGATATGGCGTTCAATAAAAGCGCCGCTGTTTTCAAGCTGGCTTAACGTCTGTGTCATGAGCGATGGTTCCTGAAACG > W3110S.gb/3047581‑3047717 | gTTTAACGATTGTGCACCAACGGCATTCAGCATTTCTTGCTGTTGCGCGGCGTCCGGTCCGATATGGCGtt < 1:1246886/71‑1 (MQ=255)gTTTAACGATTGTGCACCAACGGCATTCAGCATTTCTTGCTGTTGCGCGGCGTCCGGTCCGATATGGCGtt < 1:1889933/71‑1 (MQ=255) tAACGATTGTGCACCAACGGCATTCAGCATTTCTTGCTGTTGCGCGGCGTCCGGTCCGATATGGCGTTCaa > 1:341085/1‑71 (MQ=255) tAACGATTGTGCACCAACGGCATTCAGCATTTCTTGCTGTTGCGCGGCGTCCGGTCCGATATGGCGTTCa < 1:1212398/70‑1 (MQ=255) tAACGATTGTGCACCAACGGCATTCAGCATTTCTTGCTGTTGCGCGGCGTCCGGTCCGATATGGCGTTCa < 1:3616922/70‑1 (MQ=255) tAACGATTGTGCACCAACGGCATTCAGCATTTCTTGCTGTTGCGCGGCGTCCGGTCCGATATGGCGCTCaa > 1:2708448/1‑71 (MQ=255) aTTGTGCACCAACGGCATTCAGCATTTCTTGCTGTTGCGCGGCGTCCGGTCCGATATGGCGTTCAATaaa < 1:1701548/70‑1 (MQ=255) ccAACGGCATTCAGCATTTCTTGCTGTTGCGCGGCGTCCGGTCCGATATGGCGCTCa < 1:1031008/57‑1 (MQ=255) aaCGGCATTCAGCATTTCTTGCTGTTGCGCGGCGTCCGGTCCGATATGGCGTTCAATaa > 1:1416329/1‑59 (MQ=255) ggCATTCAGCATTTCTTGCTGTTGCGCGGCGTCCGGTCCGATATGGCGTTCAATAAAAGCGCCGCTGtttt < 1:3632127/71‑1 (MQ=255) ttCAGCATTTCTTGCTGTTGCGCGGCGTCCGGTCCGATATGGCGTTCaa < 1:3248904/49‑1 (MQ=255) ttCAGCATTTCTTGCTGTTGCGCGGCGTCCGGTCCGATATGGCGTTCaa < 1:3243347/49‑1 (MQ=255) ttCAGCATTTCTTGCTGTTGCGCGGCGTCCGGTCCGATATGGCGTTCaa < 1:2629930/49‑1 (MQ=255) ttCAGCATTTCTTGCTGTTGCGCGGCGTCCGGTCCGATATGGCGTTCa < 1:509530/48‑1 (MQ=255) ttCAGCATTTCTTGCTGTTGCGCGGCGTCCGGTCCGATATGGCGTTCAATaaaa < 1:920418/54‑1 (MQ=255) tCAGCATTTCTTGCTGTTGCGCGGCGTCCGGTCCGATATGGCGTTCAATAAAAgc > 1:733583/1‑55 (MQ=255) tCAGCATTTCTTGCTGTTGCGCGGCGTCCGGTCCGATATGGCGTTCAATAAAAg < 1:2865907/54‑1 (MQ=255) tCAGCATTTCTTGCTGTTGCGCGGCGTCCGGTCCGATATGGCGTTCAATAAAAg < 1:3030606/54‑1 (MQ=255) tCAGCATTTCTTGCTGTTGCGCGGCGTCCGGTCCGATATGGCGCTCAAt > 1:146877/1‑49 (MQ=255) tttCTTGCTGTTGCGCGGCGTCCGGTCCGATATGGCGTTCAATaaa < 1:3126799/46‑1 (MQ=255) tttCTTGCTGTTGCGCGGCGTCCGGTCCGATATGGCGCTCAATAAAAgc < 1:588396/49‑1 (MQ=255) tttCTTGCTGTTGCGCGGCGTCCGGTCCGATATGGCGCTCAATAAAAgc < 1:3536303/49‑1 (MQ=255) tCTTGCTGTTGCGCGGCGTCCGGTCCGATATGGCGCTCAATAAAAGCGCCGCTGTTTTCAAGCTGGCTTaa > 1:210804/1‑71 (MQ=255) tCTTGCTGTTGCGCGGCGTCCGGTCCGATATGGCGCTCAATAAAAGCGCCGCTGTTTTCAAGCTGGCTTaa > 1:285243/1‑71 (MQ=255) tCTTGCTGTTGCGAGGCGTCCGGTCCGATATGGCGCTCAATAAAAGCGCCGCTGttt > 1:2349989/1‑57 (MQ=255) tGCTGTTGCGCGGCGTCCGGTCCGATATGGCGCTCAATAAAAGCGCCGCTGTTTTCAAGCTGGCTTAACGt > 1:2437537/1‑71 (MQ=255) cTGTTGCGCGGCGTCCGGTCCGATATGGCGTTCAATAAAAGCGCCGCTGTTTTCAAGCTGGCTTAACGTc < 1:2428163/70‑1 (MQ=255) tGTTGCGCGGCGTCCGGTCCGATATGGCGTTCAAt < 1:1848803/35‑1 (MQ=255) gTTGCGCGGCGTCCGGTCCGATATGGCGTTCAATa < 1:4312470/35‑1 (MQ=255) gTTGCGCGGCGTCCGGTCCGATATGGCGCTCAATAAAAGCGCCGCTGTTTTCAAGc < 1:2770684/56‑1 (MQ=255) ttGCGCGGCGTCCGGTCCGATATGGCGTTCAATaa < 1:2635758/35‑1 (MQ=255) gcgcGGCGTCCGGTCCGATATGGCGTTCAATAAAATCGCCGCTGTTTTCAAGCTGGCTTAACGTCTgtgt > 1:2677202/1‑70 (MQ=255) cggtccgATATGGCGTTCAATAAAAGCGCCGCTGTTTTCAAg > 1:541013/1‑42 (MQ=255) ccgATATGGCGTTCAATAAAAGCGCCGCTGTTTTCAAGCTGGCTTAACGTCTGTGTCa > 1:329688/1‑58 (MQ=255) ccgATATGGCGTTCAATAAAAGCGCCGCTGTTTTCAAGCTGGCTTAACGTCTGTGTCATGAGCGATGGTTc > 1:1744876/1‑71 (MQ=255) cgATATGGCGTTCAATAAAAGCGCCGCTGTTTTCAAGCTGGCTTAACGTCTGTGTCATGAGCGATGGTTc < 1:2962569/70‑1 (MQ=255) gATATGGCGCTCAATAAAAGCGCCGCTGTTTTCAAGCTGGCTTAACGTCTGTGTCATGAGCTATGGTTCCt < 1:1498867/71‑1 (MQ=255) atGGCGTTCAATAAAAGCGCCGCTGTTTTCAAGCTGGCt > 1:1059253/1‑39 (MQ=255) tGGCGTTCAATAAAAGCGCCGCTGTTTTCAAGCTGGCTTAACGTCTGTGTCATg > 1:116509/1‑54 (MQ=255) ggCGTTCAATAAAAGCGCCGCTGTTTTCAAGCTGGCTTAACGTCTGTGTCATGAGCGATGGTTCCTGAAAc > 1:3866117/1‑71 (MQ=255) ggCGCTCAATAAAAGCGCCGCTGTTTTCAAGCTGGCTTAACGTCTGTGTCATGAGCGATGGt > 1:418692/1‑62 (MQ=255) gcgcTCAATAAAAGCGCCGCTGTTTTCAAGCTGGCTTAACGTCTGTGTCATGAGCGATg > 1:4225303/1‑59 (MQ=255) gCGTTCAATAAAAGCGCCGCTGTTTTCAAGCTGGCTTAACGTCTGTGTCATGAGCAATGGTTCCTGAAACg > 1:2135973/1‑71 (MQ=255) | GTTTAACGATTGTGCACCAACGGCATTCAGCATTTCTTGCTGTTGCGCGGCGTCCGGTCCGATATGGCGTTCAATAAAAGCGCCGCTGTTTTCAAGCTGGCTTAACGTCTGTGTCATGAGCGATGGTTCCTGAAACG > W3110S.gb/3047581‑3047717 |
| 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: ‑ |