This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Get positions from FASTA alignment in bed format

Hello, I have created multiple species aligments and I have found regions of interests in these alignments. They are certain aminoacids in the alignments that I have extracted based on the human position. 

>human
GMRGGVATSCHVAW
>mouse
GMRGGVATSMHVAW

Interest regions: Human, position 9, "C"

I would like to obtain the bed coordinates (in hg19) for these positions but I am not too sure how to "map" my protein alignments to the hg19 genome and then add my positions (based in my original alignment) to get their localization in bed format. Could anyone suggest me an approach to perform this? Thank you very much.

fasta bed alignment bam
• 1,172 views • link
written by kalahari • kalahari

Do you have more information about these protein sequences, like transcript or gene?

• 1 views • link
written by finswimmer finswimmer

0 answers

No answers yet.

Log in to answer this question.

More posts like this