| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | W3110S.gb | 556,858 | A→T | 100% | L36Q (CTG→CAG) | folD ← | bifunctional 5,10‑methylene‑tetrahydrofolate dehydrogenase and 5,10‑methylene‑tetrahydrofolate cyclohydrolase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | W3110S.gb | 556,858 | 0 | A | T | 94.6% | 98.6 / ‑1.1 | 37 | L36Q (CTG→CAG) | folD | bifunctional 5,10‑methylene‑tetrahydrofolate dehydrogenase and 5,10‑methylene‑tetrahydrofolate cyclohydrolase |
| Reads supporting (aligned to +/- strand): ref base A (0/2); new base T (17/18); total (17/20) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.89e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
| Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
ACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCAGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCA > W3110S.gb/556791‑556921 | acaAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTc < 1:730032/71‑1 (MQ=255)acaAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTc > 1:608762/1‑71 (MQ=255)acaAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTc < 1:3496931/71‑1 (MQ=255)acaAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTc > 1:1788302/1‑71 (MQ=255)acaAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTc > 1:1972217/1‑71 (MQ=255)acaAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACAAGCACAACGGCCTGTc > 1:1434359/1‑71 (MQ=255) caAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTcc > 1:2107873/1‑71 (MQ=255) cTTGCGTTTGCTTGCGACATAAATTTGCGATG‑GGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGc > 1:3683708/1‑71 (MQ=255) cGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCg < 1:1324747/71‑1 (MQ=255) gTTTGCTTGCGACATAAATTTGCGATGCAGGGTTGCTACCCACCAGCACAACGGCCAGTCCTGGCGCCCGc < 1:33570/71‑1 (MQ=255) gTTTGCTTGCGACATAAATTTGCGATGCAGGGTTGCTACCCACCAGCACAACGGCCAGTCCTGGCGCCCGc < 1:2146000/71‑1 (MQ=255) tttGCTTGCGACATAAATTTTCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCg < 1:1088704/69‑1 (MQ=255) tttGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCg < 1:3120440/69‑1 (MQ=255) tgcttgcGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGt > 1:1089698/1‑71 (MQ=255) tgcttgcGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGt > 1:2021343/1‑71 (MQ=255) tgcttgcGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGt > 1:1320103/1‑71 (MQ=255) aTAAATTTTCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCaa < 1:1825399/71‑1 (MQ=255) aTAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCaa < 1:863265/71‑1 (MQ=255) aTAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCaa < 1:2293380/71‑1 (MQ=255) aTAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCaa > 1:2828673/1‑71 (MQ=255) tAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAAt < 1:1709035/71‑1 (MQ=255) aaTTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATa > 1:1367611/1‑70 (MQ=255) ttGCGATGCAGGGTTACTACCCAACAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATAcgc > 1:1676819/1‑70 (MQ=255) tGCGATGCAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACgcgc > 1:325398/1‑71 (MQ=255) cAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACt > 1:2625792/1‑71 (MQ=255) cAGGGTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACt > 1:2630236/1‑71 (MQ=255) gTTACTACCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACttttt < 1:1610573/71‑1 (MQ=255) ctacCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGc < 1:91118/71‑1 (MQ=255) ctacCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGc < 1:1636142/71‑1 (MQ=255) ctacCCACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGc < 1:2802443/71‑1 (MQ=255) cccACCAGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAAc > 1:1645585/1‑71 (MQ=255) caccaGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACtt < 1:119104/71‑1 (MQ=255) aGCACAACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCaga < 1:1637011/71‑1 (MQ=255) cacaACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCaa > 1:650815/1‑38 (MQ=255) caACGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAgcgc < 1:3490203/71‑1 (MQ=255) aaCGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCa < 1:2046457/71‑1 (MQ=255) aCGGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCa < 1:313438/70‑1 (MQ=255) | ACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAACGGCCAGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCA > W3110S.gb/556791‑556921 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |