| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | W3110S.gb | 3,622,543 | C→T | 72.7% | A362A (GCG→GCA) | metE ← | 5‑methyltetrahydropteroyltriglutamate‑ homocysteine S‑methyltransferase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | W3110S.gb | 3,622,543 | 0 | C | T | 72.7% | 48.8 / 21.2 | 33 | A362A (GCG→GCA) | metE | 5‑methyltetrahydropteroyltriglutamate‑ homocysteine S‑methyltransferase |
| Reads supporting (aligned to +/- strand): ref base C (5/4); new base T (8/16); total (13/20) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.25e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01 | |||||||||||
AGAGTGACGACGTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTT > W3110S.gb/3622477‑3622602 | agagTGACGACGTGCATGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAgtgc > 1:2615316/1‑69 (MQ=255) acgacgTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGTGCATcgcg < 1:1070548/69‑1 (MQ=255) acgacgTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGTGCATcgcg < 1:2749140/69‑1 (MQ=255) acgacgTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGTGCATcg > 1:2497309/1‑67 (MQ=255) acgacgTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGTGCATcg > 1:590713/1‑67 (MQ=255) acgacgTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATcg > 1:3406367/1‑67 (MQ=255) acgTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGTGCATCGCgcag < 1:3382245/69‑1 (MQ=255) ccTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGTGCATCGCgcag < 1:2272902/64‑1 (MQ=255) tGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGTGCATCGCGCAGCAGTGc < 1:3439249/68‑1 (MQ=255) tCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGTGCATCGCGCAGCAGTGCCAGt < 1:1965220/68‑1 (MQ=255) tCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGTGCATCGCGCAGCAGTGCCAGt > 1:3517111/1‑68 (MQ=255) tCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGTGCATCGCGCAGCAGTGCCAGt < 1:1186986/68‑1 (MQ=255) tCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGtt < 1:62888/69‑1 (MQ=255) tCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGtt < 1:1524555/69‑1 (MQ=255) tCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTCCCAGtt < 1:2169879/69‑1 (MQ=255) gcTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGcc > 1:389825/1‑58 (MQ=255) tctGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTg < 1:2407459/68‑1 (MQ=255) ctGCCAGAGCTGCCGTGTCACCACTGTTCAGTGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGt < 1:1539353/68‑1 (MQ=255) ctGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTa > 1:2335207/1‑69 (MQ=255) ctGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTa > 1:1442017/1‑69 (MQ=255) ctGCCAGAGCTGCCGTGTCACAACTGTTCAGCGCATCGCGCAGCAGTg > 1:520623/1‑48 (MQ=255) agagCTGCCGTGTCACCACTGTTCAGTGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCg < 1:2571662/69‑1 (MQ=255) cTGCCGTGTCACCACTGTTCAGTGCATCGCGCAGCAGTGCCAGTTCa > 1:2452896/1‑47 (MQ=255) gCCGTGTCACCACTGTTCAGTGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCa < 1:3082525/69‑1 (MQ=255) gCCGTGTCACCACTGTTCAGTGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCa < 1:998178/69‑1 (MQ=255) tgtCACCACTGTTCAGTGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAAcc < 1:3148421/69‑1 (MQ=255) tgtCACCACTGTTCAGTGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAAcc < 1:1581944/69‑1 (MQ=255) accacTGTTCAGTGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCaa > 1:3387458/1‑62 (MQ=255) accacTGTTCAGTGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCaa > 1:864983/1‑62 (MQ=255) tgtGTTCAGTGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCtttt < 1:1561397/67‑1 (MQ=255) cAGTGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCtttt < 1:1039896/63‑1 (MQ=255) cAGTGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCtttt > 1:3650984/1‑63 (MQ=255) cAGTGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCtttt < 1:2404714/63‑1 (MQ=255) cAGTGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCtttt > 1:1924890/1‑63 (MQ=255) | AGAGTGACGACGTGCCTGAATCGGGGCGCTCCACTCTGCCAGAGCTGCCGTGTCACCACTGTTCAGCGCATCGCGCAGCAGTGCCAGTTCATGGCATTTTTGTAGGGCGAAGGCAAACCAGCTTTT > W3110S.gb/3622477‑3622602 |
| 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: ‑ |