Question: Getting protein fasta from RNAseq BAM file (human) paying attention on repeats
gravatar for pauli.pauli
3 months ago by
pauli.pauli0 wrote:

Hi everybody,

I have several RNAseq BAM files (mapped with BWA and GATK IndelRealigner) of cancer patients of which I should compare if there any copynumber in repeats appeared (whole transcriptome). I need to find a pipeline to find repeats in with our repeat detection software which needs FASTA and protein sequences as input.

Now I would like to have each patient's personal transcriptome as Fasta format which I could only feed to our software. What do I do now? For example, can I use Blastx, which compares directly to a reference proteome? Does this pay attention to the reading-frame? How do I get a Fasta output with each protein? In this pipeline, it is important to pay special attention to repeats that deviate from the reference.

Many thanks for your help!

rna-seq blastx bam repeats fasta • 198 views
ADD COMMENTlink modified 3 months ago by RamRS21k • written 3 months ago by pauli.pauli0
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: 1117 users visited in the last hour