Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 4,438,531 | T→G | D110A (GAC→GCC) | fbp ← | fructose‑1,6‑bisphosphatase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 4,438,531 | 0 | T | G | 96.8% | 118.1 / ‑4.0 | 31 | D110A (GAC→GCC) | fbp | fructose‑1,6‑bisphosphatase |
Reads supporting (aligned to +/- strand): ref base T (1/0); new base G (13/17); total (14/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.52e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01 |
GCGGCGGTAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGTCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAA > REL606/4438461‑4438603 | gcggcgGTAGATGGAGAAAATGGTACCGACAGAGACGTTAAGATCGCTGTTGGACGAGCCATCCAGGAGGGCCATc > 2:1310528/1‑76 (MQ=255) cggcggTAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCa < 2:172599/76‑1 (MQ=255) cggcggTAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCa < 2:607995/76‑1 (MQ=255) ggTAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGcac < 1:1311214/76‑1 (MQ=255) tAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGcacca > 2:83300/1‑76 (MQ=255) tAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGcacca > 2:326712/1‑76 (MQ=255) tAGAGGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGcacca < 1:1666195/76‑1 (MQ=255) tGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGTGCCATCAGCACCACGTACTTTGCgtg > 1:178919/1‑76 (MQ=255) tGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCgtg > 2:444990/1‑76 (MQ=255) ccGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTcac < 1:1069102/76‑1 (MQ=255) aCAGAAACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGc < 2:1823558/76‑1 (MQ=255) aCGTTAACATCGGTGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCaa < 2:190417/76‑1 (MQ=255) aCGTTAACATCGATGTTGGACGAGCCATCCAGGGGGTCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCaa > 1:1378142/1‑76 (MQ=255) gTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAg < 1:409832/76‑1 (MQ=255) ttAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAga < 1:665978/76‑1 (MQ=255) aaCATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAgacg > 2:616375/1‑76 (MQ=255) aaCATCGATGTCGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAgacg > 2:525640/1‑76 (MQ=255) aCATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAgacga > 2:1046763/1‑76 (MQ=255) cATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAgacgac < 1:1694044/76‑1 (MQ=255) aTCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACa < 2:1056966/76‑1 (MQ=255) aTCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACa < 2:1572039/76‑1 (MQ=255) cGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAAt > 2:1567841/1‑76 (MQ=255) ttGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCAt < 1:1111475/76‑1 (MQ=255) tGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCAtc > 2:1782923/1‑76 (MQ=255) ggACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCAtct > 2:116965/1‑76 (MQ=255) gCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCAtcttcttct > 1:888984/1‑76 (MQ=255) ccATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATcttcttctt < 1:427901/76‑1 (MQ=255) ccATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATcttcttctt < 1:219716/76‑1 (MQ=255) cATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCttcttcttc < 1:1530123/76‑1 (MQ=255) ccAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCaga > 2:420624/1‑76 (MQ=255) ggggggCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGc > 2:1110881/1‑76 (MQ=255) ggggggCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGc < 2:815824/76‑1 (MQ=255) ggggGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCa < 1:934713/76‑1 (MQ=255) ggggCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCaa > 1:187148/1‑76 (MQ=255) | GCGGCGGTAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGTCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAA > REL606/4438461‑4438603 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |