Question: Extract SNPs flanking sequence from organisme without reference genome
gravatar for Konan
2.2 years ago by
Konan20 wrote:

Hi, I'm new in managing DNA sequences and I'm looking for help. I have fastq and vcf files from the sequencing of my samples. The plant species on which I work (Lagenaria siceraria) has no reference genome. What I want to do is to extract SNPs flanking sequences and make a blast in Plant RefSeq for determining the putative functions of my SNPs. All the posts Ive read are related to species that have model organism. Then, could someone help me please?


ADD COMMENTlink modified 2.2 years ago by Jennifer Hillman Jackson25k • written 2.2 years ago by Konan20
gravatar for Jennifer Hillman Jackson
2.2 years ago by
United States
Jennifer Hillman Jackson25k wrote:


Some form of an assembled reference fasta file will be needed to perform this function. This can be a transcriptome or genome assembly (complete or partial, for either) - with the sequences mapped to that assembled fasta file. Use a Custom reference genome/build when working with your own assembly.


Recent publication that outlines computational approaches when performing de-novo variant analysis. This may be a fit for your analysis, or portions of it may be helpful. A Galaxy tutorial/example is included:

Thanks, Jen, Galaxy team

ADD COMMENTlink written 2.2 years ago by Jennifer Hillman Jackson25k

Thanks a lot Jen, I will follow your Avice and let keep you informed!

ADD REPLYlink written 2.2 years ago by Konan20
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 180 users visited in the last hour