Question: How to identify promoter sequences from plant NGS raw data
0
gravatar for anikng
5.0 years ago by
anikng0
Korea, Republic Of
anikng0 wrote:

Recently, I got genome sequencing data of a plant subspecies for which I HAVE WELL ANNOTATED REFERENCE SEQUENCE information. My objective is to comparative analysis of promoter regions of a few genes. How can I get the promoter region from NGS genome sequencing reads? Can anyone provide a correct method for that? 

 

Thanks..

ADD COMMENTlink modified 5.0 years ago by Istvan Albert ♦♦ 81k • written 5.0 years ago by anikng0
0
gravatar for Devon Ryan
5.0 years ago by
Devon Ryan92k
Freiburg, Germany
Devon Ryan92k wrote:

If the well-annotated species and yours are highly related, then you'll likely have luck just aligning your reads against it and using it's reference annotation to call variants in promoter regions. Otherwise, you'd be better served by assembling your species and trying to align the contigs (you'd then find differences in the results). Since these are subspecies we're talking about, I assume the former possibility will work well and prove easier than the latter.

ADD COMMENTlink written 5.0 years ago by Devon Ryan92k

Thanks @ Devon Ryan.. I will look into it..Also i worry whether i can retrieve sequence rather than direct analysis to reference sequence.

ADD REPLYlink written 5.0 years ago by anikng0

Sorry, I can't parse your last sentence.

ADD REPLYlink written 5.0 years ago by Devon Ryan92k

Sorry, what I mean is that in my case, I need to retrieve ~5kb upstream region of that gene (because I need that sequence for further MSA with other subspecies sequences I have). I am new to 'variant call in promoter region'. So, my impression is that using that method, I can identify the variation, but CANNOT fetch the sequence. @Devon Ryan

ADD REPLYlink written 5.0 years ago by anikng0

Ah, right, yes, calling the variants won't change the sequence. GATK has a tool that will convert a reference fasta file according to variants in a VCF file though.

ADD REPLYlink written 5.0 years ago by Devon Ryan92k

Thanks a lot @Devon Ryan...I  will look into it..

ADD REPLYlink written 5.0 years ago by anikng0
Please log in to add an answer.

Help
Access

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