| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,478,087 | A→G | 46.2% | N159N (AAT→AAC) | glpQ ← | periplasmic glycerophosphodiester phosphodiesterase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,478,087 | 0 | A | G | 46.2% | ‑2.3 / 18.5 | 26 | N159N (AAT→AAC) | glpQ | periplasmic glycerophosphodiester phosphodiesterase |
| Reads supporting (aligned to +/- strand): ref base A (3/11); new base G (6/6); total (9/17) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.18e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GATGATGGAACCACGGCGCTTTGATTTCTGGATAAATACCGATATTTTTCCCGGTAGAGTGATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACC > minE/1478026‑1478149 | gatgatGGAACCACGGCGCTTTGATTTCTGGATAAATACCGATATTTTTCCCGGTAGAGTGATTTAAc < 1:1322798/68‑1 (MQ=255)gatgatGGAACCACGGCGCTTTGATTTCTGGATAAATACCGATATTTTTCCCGGTAGAGTGATTTAAc < 1:2369341/68‑1 (MQ=255)gatgatGGAACCACGGCGCTTTGATTTCTGGATAAATACCGATATTTTTCCCGGTAGAGTGATTTAAc < 1:762283/68‑1 (MQ=255)gatgatGGAACCACGGCGCTTTAATTTCTGGATAGATACCGATATTTTTCCCGGTAGAGTGGTTTAAcc < 1:3253435/69‑1 (MQ=255)gatgatGGAACCACGGCGCTTTAATTTCTGGATAGATACCGATATTTTTCCCGGTAGAGTGGTTTAAcc < 1:1785215/69‑1 (MQ=255)gatgatGGAACCACGGCGCTTTAATTTCTGGATAGATACCGATATTTTTCCCGGTAGAGTGGTTTAAcc < 1:1853494/69‑1 (MQ=255)gatgatGGAACCACGGCGCTTTAATTTCTGGATAGATACCGATATTTTTCCCGGTAGAGTGGTTTAAcc < 1:1987551/69‑1 (MQ=255) tgGAACCACGGCGCTTTAATTTCTGGATAGATACCGATAATTTTCCCGGTAGAGTGGTTTAAcc < 1:2303889/64‑1 (MQ=255) aCGGCGCTTTGATTTCTGGATAAATACCGATATTTTTCCCGGTAGAGTGATTTAACCCCTGAACAAAt > 1:467058/1‑68 (MQ=255) tGGATAGATACCGATATTTTTCCCGGTAGAGTGGTTTAACCCCTGAACAAATTCAATCTCTTCTTCaa < 1:1642105/68‑1 (MQ=255) tAAATACCGATATTTTTCCCGGTAGAGTGATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAg < 1:270726/66‑1 (MQ=255) tAAATACCGATATTTTTCCCGGTAGAGTGATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAg < 1:1033842/66‑1 (MQ=255) aaTACCGATATTTTTCCCGGTAGAGTGATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGtgtg < 1:946279/69‑1 (MQ=255) aaTACCGATATTTTTCCCGGTAGAGTGATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGtgtg < 1:802439/69‑1 (MQ=255) aaTACCGATATTTTTCCCGGTAGAGTGATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGtgtg < 1:240556/69‑1 (MQ=255) aaTACCGATATTTTTCCCGGTAGAGTGATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGtgtg < 1:128732/69‑1 (MQ=255) tttCCCGGTAGAGTGATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCAc > 1:289940/1‑60 (MQ=255) tttCCCGGTAGAGTCATTTAACCCCTGAACAAATTCAATCTCTTCTTCCAAGGTGTGCACCCGGAAGt > 1:2927750/1‑68 (MQ=255) tAGAGTGGTTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTca > 1:2403084/1‑61 (MQ=255) tAGAGTGGTTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCAGATTTa > 1:841378/1‑68 (MQ=255) tAGAGTGGTTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCAGATTTa > 1:2670187/1‑68 (MQ=255) tAGAGTGGTTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCAGATTTa > 1:560868/1‑68 (MQ=255) tAGAGTGGTTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCAGATTTa > 1:3040089/1‑68 (MQ=255) tAGAGTGGTTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCAGATTTa > 1:2927689/1‑68 (MQ=255) agagTGATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTAcc < 1:2460895/69‑1 (MQ=255) agagTGATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTAcc < 1:1201387/69‑1 (MQ=255) | GATGATGGAACCACGGCGCTTTGATTTCTGGATAAATACCGATATTTTTCCCGGTAGAGTGATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACC > minE/1478026‑1478149 |
| 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: ‑ |