| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,352,737 | T→C | 100% | E259G (GAG→GGG) | thiM ← | hydoxyethylthiazole kinase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,352,737 | 0 | T | C | 97.5% | 143.4 / ‑3.7 | 40 | E259G (GAG→GGG) | thiM | hydoxyethylthiazole kinase |
| Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (0/39); minor base A (0/1); total (0/40) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01 | |||||||||||
| Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
| Rejected as polymorphism: Frequency below/above cutoff threshold. | |||||||||||
| Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
| Rejected as polymorphism: Polymorphic base substitution creates a homopolymer stretch. | |||||||||||
AGCGTTAATTCGTTTCATGCCTGCACCTCCTGCGTCAATTGCCAGAGCGCATCAAGGAAATGTGGAACAAA > minE/1352710‑1352780 | agCGTTAATTCGTTTCATGCCTGCACCCCCTGCGTCAATTGCCAGAGCGCTTCAAGGAAATGTGGAACaaa < 1:171741/71‑1 (MQ=255)agCGTTAATTCGTTTCATGCCTGCACCCCCTGCGTCAATTGCCAGAGCGCATCAAGGAAATGTGGAACaaa < 1:561675/71‑1 (MQ=255)agCGTTAATTCGTTTCATGCCTGCACCCCCTGCGTCAATTGCCAGAGCGCATCAAGGAAATGTGGAACaaa < 1:343983/71‑1 (MQ=255)agCGTTAATTCGTTTCATGCCTGCACCCCCTGCGTCAATTGCCAGAGCGCATCAAGGAAATGTGGAACaaa < 1:398097/71‑1 (MQ=255)agCGTTAATTCGTTTCATGCCTGCACCCCCTGCGTCAATTGCCAGAGCGCATCAAGGAAATGTGGAACaaa < 1:431827/71‑1 (MQ=255)agCGTTAATTCGTTTCATGCCTGCACCCCCTGCGTCAATTGCCAGAGCGCATCAAGGAAATGTGGAACaaa < 1:476111/71‑1 (MQ=255)agCGTTAATTCGTTTCATGCCTGCACCCCCTGCGTCAATTGCCAGAGCGCATCAAGGAAATGTGGAACaaa < 1:497217/71‑1 (MQ=255)agCGTTAATTCGTTTCATGCCTGCACCCCCTGCGTCAATTGCCAGAGCGCATCAAGGAAATGTGGAACaaa < 1:506836/71‑1 (MQ=255)agCGTTAATTCGTTTCATGCCTGCACCCCCTGCGTCAATTGCCAGAGCGCATCAAGGAAATGTGGAACaaa < 1:52485/71‑1 (MQ=255)agCGTTAATTCGTTTCATGCCTGCACCCCCTGCGTCAATTGCCAGAGCGCATCAAGGAAATGTGGAACaaa < 1:550964/71‑1 (MQ=255)agCGTTAATTCGTTTCATGCCTGCACCCCCTGCGTCAATTGCCAGAGCGCATCAAGGAAATGTGGAACaaa < 1:390088/71‑1 (MQ=255)agCGTTAATTCGTTTCATGCCTGCACCCCCTGCGTCAATTGCCAGAGCGCATCAAGGAAATGTGGAACaaa < 1:627428/71‑1 (MQ=255)agCGTTAATTCGTTTCATGCCTGCACCCCCTGCGTCAATTGCCAGAGCGCATCAAGGAAATGTGGAACaaa < 1:639909/71‑1 (MQ=255)agCGTTAATTCGTTTCATGCCTGCACCCCCTGCGTCAATTGCCAGAGCGCATCAAGGAAATGTGGAACaaa < 1:641718/71‑1 (MQ=255)agCGTTAATTCGTTTCATGCCTGCACCCCCTGCGTCAATTGCCAGAGCGCATCAAGGAAATGTGGAACaaa < 1:644015/71‑1 (MQ=255)agCGTTAATTCGTTTCATGCCTGCACCCCCTGCGTCAATTGCCAGAGCGCATCAAGGAAATGTGGAACaaa < 1:65594/71‑1 (MQ=255)agCGTTAATTCGTTTCATGCCTGCACCCCCTGCGTCAATTGCCAGAGCGCATCAAGGAAATGTGGAACaaa < 1:68615/71‑1 (MQ=255)agCGTTAATTCGTTTCATGCCTGCACCCCCTGCGTCAATTGCCAGAGCGCATCAAGGAAATGTGGAACaaa < 1:9278/71‑1 (MQ=255)agCGTTAATTCGTTTCATGCCTGCACCCCCTGCGTCAATTGCCAGAGCGCATCAAGGAAATGTGGAACaaa < 1:215457/71‑1 (MQ=255)agCGTTAATTCGTTTCATGCCTGCACCCCCTGCGTCAATTGCCAGAGCGCATCAAGGAAATGTGGAACaaa < 1:120978/71‑1 (MQ=255)agCGTTAATTCGTTTCATGCCTGCACCCCCTGCGTCAATTGCCAGAGCGCATCAAGGAAATGTGGAACaaa < 1:123097/71‑1 (MQ=255)agCGTTAATTCGTTTCATGCCTGCACCCCCTGCGTCAATTGCCAGAGCGCATCAAGGAAATGTGGAACaaa < 1:129715/71‑1 (MQ=255)agCGTTAATTCGTTTCATGCCTGCACCCCCTGCGTCAATTGCCAGAGCGCATCAAGGAAATGTGGAACaaa < 1:155079/71‑1 (MQ=255)agCGTTAATTCGTTTCATGCCTGCACCCCCTGCGTCAATTGCCAGAGCGCATCAAGGAAATGTGGAACaaa < 1:164738/71‑1 (MQ=255)agCGTTAATTCGTTTCATGCCTGCACCCCCTGCGTCAATTGCCAGAGCGCATCAAGGAAATGTGGAACaaa < 1:329500/71‑1 (MQ=255)agCGTTAATTCGTTTCATGCCTGCACCCCCTGCGTCAATTGCCAGAGCGCATCAAGGAAATGTGGAACaaa < 1:240812/71‑1 (MQ=255)agCGTTAATTCGTTTCATGCCTGCACCCCCTGCGTCAATTGCCAGAGCGCATCAAGGAAATGTGGAACaaa < 1:253520/71‑1 (MQ=255)agCGTTAATTCGTTTCATGCCTGCACCCCCTGCGTCAATTGCCAGAGCGCATCAAGGAAATGTGGAACaaa < 1:266383/71‑1 (MQ=255)agCGTTAATTCGTTTCATGCCTGCACCCCCTGCGTCAATTGCCAGAGCGCATCAAGGAAATGTGGAACaaa < 1:282500/71‑1 (MQ=255)agCGTTAATTCGTTTCATGCCTGCACCCCCTGCGTCAATTGCCAGAGCGCATCAAGGAAATGTGGAACaaa < 1:288399/71‑1 (MQ=255)agCGTTAATTCGTTTCATGCCTGCACCCCCTGCGTCAATTGCCAGAGCGCATCAAGGAAATGTGGAACaaa < 1:323092/71‑1 (MQ=255)agCGTTAATTCGTTTCATGCCTGCACCCCCTGCGTCAATTGCCAGAGCGCATCAAGGAAATGTGGAACaaa < 1:324129/71‑1 (MQ=255)agCGTTAATTCGTTTCATGCCTGCACCCCCTGCGTCAATTGCCAGAGCGCATAAAGGAAATGTGGAACaaa < 1:152647/71‑1 (MQ=255)agCGTTAATTCGTTTCATGCCTGCACCACCTGCGTCAATTGCCAGAGCGCATCAAGGAAATGTGGAACaaa < 1:395105/71‑1 (MQ=255) gCGTTAATTCGTTTCATGCCTGCACCCCCTGCGTCAATTGCCAGAGCGCATCAAGGAAATGTGGAACaaa < 1:583276/70‑1 (MQ=255) gCGTTAATTCGTTTCATGCCTGCACCCCCTGCGTCAATTGCCAGAGCGCATCAAGGAAATGTGGAACaaa < 1:197342/70‑1 (MQ=255) gCGTTAATTCGTTTCATGCCTGCACCCCCTGCGTCAATTGCCAGAGCGCATCAAGGAAATGTGGAACaaa < 1:120862/70‑1 (MQ=255) gCGTTAATTCGTTTCATGCCTGCACCCCCTGCGTCAATTGCCAGAGCGCATCAAGGAAATGTGGAACaaa < 1:380715/70‑1 (MQ=255) aTTCGTTTCATGCCTGCACCCCCTGCGTCAATTGCCAGAGCGCATCAAGGAAATGTGGAACaaa < 1:564879/64‑1 (MQ=255) ttCGTTTCATGCCTGCACCCCCTGCGTCAATTGCCAGAGCGCATCAAGGAAATGTGGAACaaa < 1:57720/63‑1 (MQ=255) | AGCGTTAATTCGTTTCATGCCTGCACCTCCTGCGTCAATTGCCAGAGCGCATCAAGGAAATGTGGAACAAA > minE/1352710‑1352780 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |