| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,477,520:1 | +G | 100% | coding (1044/1077 nt) | glpQ ← | periplasmic glycerophosphodiester phosphodiesterase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,477,520 | 1 | . | G | 88.9% | 88.0 / 4.8 | 27 | coding (1044/1077 nt) | glpQ | periplasmic glycerophosphodiester phosphodiesterase |
| Reads supporting (aligned to +/- strand): ref base . (3/0); new base G (24/0); total (27/0) | |||||||||||
| 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.91e-01 | |||||||||||
| Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
| Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch. | |||||||||||
GCAGTAACATTACTCTTTATTAAGAAATTTTACTGCCTTATCA‑GGGAAATCAGTAAACAGCCCATTTACACCCGCT > minE/1477478‑1477553 | gCAGTAACATTACTCTTTATTAAGAAATTTTACTGCCTTATCA‑GGGAAATCAGTAAACAGCCCAtt > 1:706406/1‑66 (MQ=255)gCAGTAACATTACTCTTTATTAAGAAATTTTACTGCCTTATCA‑GGGAAATCAGTAAACAGCCCAtt > 1:1424947/1‑66 (MQ=255) tACTCTTTATTAAGAAATTTTACTGCCTTATCAGGGGAAATCAGTa > 1:333567/1‑46 (MQ=255) tACTCTTTATTAAGAAATTTTACTGCCTTATCAGGGGAAATCAGTAAACa > 1:907305/1‑50 (MQ=255) tACTCTTTATTAAGAAATTTTACTGCCTTATCAGGGGAAATCAGTAAACAg > 1:1760185/1‑51 (MQ=255) tACTCTTTATTAAGAAATTTTACTGCCTTATCAGGGGAAATCAGTAAACAGCCCATTTACACCCg > 1:1075874/1‑65 (MQ=255) tACTCTTTATTAAGAAATTTTACTGCCTTATCAGGGGAAATCAGTAAACAGCCCATTTACACCCg > 1:595367/1‑65 (MQ=255) tACTCTTTATTAAGAAATTTTACTGCCTTATCAGGGGAAATCAGTAAACAGCCCATTTACACCCGc > 1:1002525/1‑66 (MQ=255) tACTCTTTATTAAGAAATTTTACTGCCTTATCAGGGGAAATCAGTAAACAGCCCATTTACACCCGc > 1:1002224/1‑66 (MQ=255) tACTCTTTATTAAGAAATTTTACTGCCTTATCAGGGGAAATCAGTAAACAGCCCATTTACACCCGc > 1:1064399/1‑66 (MQ=255) tACTCTTTATTAAGAAATTTTACTGCCTTATCAGGGGAAATCAGTAAACAGCCCATTTACACCCGc > 1:12352/1‑66 (MQ=255) tACTCTTTATTAAGAAATTTTACTGCCTTATCAGGGGAAATCAGTAAACAGCCCATTTACACCCGCt > 1:2005809/1‑67 (MQ=255) tACTCTTTATTAAGAAATTTTACTGCCTTATCAGGGGAAATCAGTAAACAGCCCATTTACACCCGCt > 1:962625/1‑67 (MQ=255) tACTCTTTATTAAGAAATTTTACTGCCTTATCAGGGGAAATCAGTAAACAGCCCATTTACACCCGCt > 1:933605/1‑67 (MQ=255) tACTCTTTATTAAGAAATTTTACTGCCTTATCAGGGGAAATCAGTAAACAGCCCATTTACACCCGCt > 1:364590/1‑67 (MQ=255) tACTCTTTATTAAGAAATTTTACTGCCTTATCAGGGGAAATCAGTAAACAGCCCATTTACACCCGCt > 1:2097087/1‑67 (MQ=255) tACTCTTTATTAAGAAATTTTACTGCCTTATCAGGGGAAATCAGTAAACAGCCCATTTACACCCGCt > 1:2010258/1‑67 (MQ=255) tACTCTTTATTAAGAAATTTTACTGCCTTATCAGGGGAAATCAGTAAACAGCCCATTTACACCCGCt > 1:1336749/1‑67 (MQ=255) tACTCTTTATTAAGAAATTTTACTGCCTTATCAGGGGAAATCAGTAAACAGCCCATTTACACCCGCt > 1:1881706/1‑67 (MQ=255) tACTCTTTATTAAGAAATTTTACTGCCTTATCAGGGGAAATCAGTAAACAGCCCATTTACACCCGCt > 1:1821716/1‑67 (MQ=255) tACTCTTTATTAAGAAATTTTACTGCCTTATCAGGGGAAATCAGTAAACAGCCCATTTACACCCGCt > 1:1784098/1‑67 (MQ=255) tACTCTTTATTAAGAAATTTTACTGCCTTATCAGGGGAAATCAGTAAACAGCCCATTTACACCCGCt > 1:128007/1‑67 (MQ=255) tACTCTTTATTAAGAAATTTTACTGCCTTATCAGGGGAAATCAGTAAACAGCCCATTTACACCCGCt > 1:15032/1‑67 (MQ=255) tACTCTTTATTAAGAAATTTTACTGCCTTATCAGGGGAAATCAGTAAACAGCCCATTTACACCCGCt > 1:1446751/1‑67 (MQ=255) tACTCTTTATTAAGAAATTTTACTGCCTTATCAGGGGAAATCAGTAAACAGCCCATTTACACCCGCt > 1:1296228/1‑67 (MQ=255) tACTCTTTATTAAGAAATTTTACTGCCTTATCAGGGGAAATCAGTAAAAAGCCCATTTACACCCGCt > 1:1132159/1‑67 (MQ=255) tACTCTTTATTAAGAAATTTTACTGCCTTATCA‑GGGgaaat > 1:1748761/1‑39 (MQ=255) | GCAGTAACATTACTCTTTATTAAGAAATTTTACTGCCTTATCA‑GGGAAATCAGTAAACAGCCCATTTACACCCGCT > minE/1477478‑1477553 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |