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 | 98.2% | 225.5 / ‑4.1 | 54 | D110A (GAC→GCC) | fbp | fructose‑1,6‑bisphosphatase |
Reads supporting (aligned to +/- strand): ref base T (1/0); new base G (20/33); total (21/33) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.89e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.46e-01 |
CGGCGGTAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGTCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAAT > REL606/4438462‑4438604 | cggcggTAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCa > 2:82356/1‑76 (MQ=255) cggcggTAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCa < 2:1296040/76‑1 (MQ=255) cggcggTAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCa > 1:1561967/1‑76 (MQ=255) ggcggTAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAg > 2:377994/1‑76 (MQ=255) ggcggTAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAg < 2:1537609/76‑1 (MQ=255) gTAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGcacc > 1:1295547/1‑76 (MQ=255) gTAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGcacc > 1:63743/1‑76 (MQ=255) gTAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGcacc > 1:40317/1‑76 (MQ=255) tAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGcacca < 1:925672/76‑1 (MQ=255) tAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGcacca > 2:1216487/1‑76 (MQ=255) aGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGcaccac > 2:764807/1‑76 (MQ=255) aGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGcaccac > 1:846072/1‑76 (MQ=255) tGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTa < 2:1563851/76‑1 (MQ=255) gagaAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTAtt < 1:305734/76‑1 (MQ=255) agaAAATGGTGCCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTAttt < 1:1301161/76‑1 (MQ=255) agaAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTAttt < 2:777314/76‑1 (MQ=255) gaAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTAtttt < 1:774436/76‑1 (MQ=255) gaAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTAtttt < 1:971364/76‑1 (MQ=255) aaaaaTGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTAtttt < 2:353030/75‑1 (MQ=255) tGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCgtg > 2:1362187/1‑76 (MQ=255) tGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCgtg < 2:1589783/76‑1 (MQ=255) ggTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCgtgt > 1:1238976/1‑76 (MQ=255) aCCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTca < 2:692805/76‑1 (MQ=255) ccGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTcac < 2:883708/76‑1 (MQ=255) ccGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTcac < 2:1192608/76‑1 (MQ=255) ccGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTcac < 1:1318209/76‑1 (MQ=255) ccGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTcac < 2:1513499/76‑1 (MQ=255) cGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTcaca < 1:1149448/76‑1 (MQ=255) cGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTcaca < 1:1326869/76‑1 (MQ=255) gACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAg < 1:1504959/76‑1 (MQ=255) agagACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACTGCCt < 2:63743/76‑1 (MQ=255) agagACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCt < 2:1610927/76‑1 (MQ=255) aCGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCaa > 1:892874/1‑76 (MQ=255) aCGTTAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCaa < 1:12870/76‑1 (MQ=255) ttAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCTCCACGTATTTTGCGTGTTCACAGCCTTCAAAga < 2:798678/76‑1 (MQ=255) ttAACATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCTCCACGTATTTTGCGTGTTCACAGCCTTCAAAga < 2:798673/76‑1 (MQ=255) aaCATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAgacg < 1:1517098/76‑1 (MQ=255) aaCATCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAgacg < 1:736511/76‑1 (MQ=255) aTCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACa < 1:694149/76‑1 (MQ=255) tCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACaa > 2:451314/1‑76 (MQ=255) tCGATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACaa > 2:610406/1‑76 (MQ=255) gATGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAACGACGACAAtc > 1:1447161/1‑76 (MQ=255) tGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAAtctc < 1:1672602/76‑1 (MQ=255) tGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAAtctc > 2:943565/1‑76 (MQ=255) tGTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAAtctc > 2:122789/1‑76 (MQ=255) cgTTGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAAtctc > 1:578505/2‑76 (MQ=255) ttGGACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCAt > 2:977122/1‑76 (MQ=255) gACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCAtctt < 1:1246371/76‑1 (MQ=255) gACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCAtctt < 1:537776/76‑1 (MQ=255) gACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCAtctt < 2:954713/76‑1 (MQ=255) gACGAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCAtctt > 2:974255/1‑76 (MQ=255) gAGCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCAtcttctt < 1:969356/76‑1 (MQ=255) gCCATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCAtcttcttct < 2:974729/76‑1 (MQ=255) cATCCAGGGGGGCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCttcttcttc > 2:631512/1‑76 (MQ=255) ggggCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCaa < 2:1287661/76‑1 (MQ=255) ggTCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAAt > 1:206442/1‑76 (MQ=255) | CGGCGGTAGATGGAGAAAATGGTACCGACAGAGACGTTAACATCGATGTTGGACGAGCCATCCAGGGGGTCCATCAGCACCACGTATTTTGCGTGTTCACAGCCTTCAAAGACGACAATCTCATCTTCTTCTTCAGAGGCAAT > REL606/4438462‑4438604 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |