Question: Getting protein fasta from RNAseq BAM file (human) paying attention on repeats
0
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.

Help
Access

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