Question: Aligning protein sequences to local RNA-seq data
gravatar for friedrichlab
2.7 years ago by
friedrichlab0 wrote:

I have a large .BAM file containing cleaned reads of RNA Seq data from a partnered sequencing project. My goal is to take existing protein sequences and align them to the RNA Seq data in order to identify corresponding sequence regions. I have looked into using multiple resources to do this (BLAST command line, genome workbench, picard), but have not found much success as my bioinformatics experience only goes so far as R, muscle, and the online BLAST tool. I was wondering if anyone had any recommendations or experience with something like this that they could offer.

Thank you in advance.

rna-seq alignment • 1.1k views
ADD COMMENTlink written 2.7 years ago by friedrichlab0

You should do the alignments the other way around. Use the reads to map against protein DB. DIAMOND is an option and will require significant hardware to be available.

Why are you going this route though? Is the data already aligned to a reference? If yes then you should be able to see what the genes are.

ADD REPLYlink written 2.7 years ago by genomax80k

Do you want to annotate the genome, or to annotate the transcriptome? Or do you want to do something else entirely? You are describing the task you want to perform, but not the goal you want to achieve. As this task seems a bit odd, maybe if you tell what is your ultimate goal, someone can be more helpful.

ADD REPLYlink modified 2.7 years ago • written 2.7 years ago by h.mon29k
Please log in to add an answer.


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