This is a test version of Biostars. For the public version, visit https://www.biostars.org.
translating a bam file to protein sequences

we would like to translate the genomic sequences in a bam file to its protein sequence (all six frames).

We would like to do it while keeping the structure of the bam file ( or at least a sam file), if possible.

I was wondering if anyone know about a tool to do such a task.

thanks

bam protein translate samtools
• 2,399 views • link
written by Assa Yeroslaviz Assa Yeroslaviz

we would like to translate the genomic sequences in a bam file

Sequence present in each line/record?

• 0 views • link
written by GenoMax GenoMax

Yes, We might not need it for the whole bam file (see above), but I would like to know if there is an option.

• 0 views • link
written by Assa Yeroslaviz Assa Yeroslaviz

You can convert the reads into fasta format by using BBMap suite.

reformat.sh in=your.bam out=file.fa
reformat.sh in=your.bam out1=R1.fa out2=R2.fa

Then use sixpack or transeq from EMBOSS to get the translations in all frames.

• 0 views • link
written by GenoMax GenoMax

Thanks I'll try that. I was hoping for a tool which can keep the bam/sam structure, but I'll keep looking.

• 0 views • link
written by Assa Yeroslaviz Assa Yeroslaviz

That is something really uncommon because you have millions of short reads (nucleotides) in the BAM file.

If you a looking for a particular protein, it is better to do some assembly and look for your protein in the contigs.

• 0 views • link
written by JC JC

Yes, I know this is uncommon und I couldn't find any solution for that yet. But I still look for one, if possible. It must not be a huge bam file, it can be one where we have extracted reads based on a bed file for specific positions. It is mainly meant to be used in identifying SAP (single AA polymorphism), so we might split a possible bam file for only the positions we need.

• 0 views • link
written by Assa Yeroslaviz Assa Yeroslaviz

0 answers

No answers yet.

Log in to answer this question.

More posts like this