Question: Get positions in original sequence file after alignment/trimming
gravatar for palabari
6 months ago by
palabari10 wrote:

I would like to have suggestions into how I could solve this problem:

I have downloaded several cDNA from Ensembl, such as: ENSG00000139618. I have aligned them with other proteins, trimmed unnecesary bits and deleted some regions. Then, I run these alignments in some other program and I get a list of interesting aminoacid positions. The problem is, that I would like to map these positions (based on the trimmed alignment) to the equivalent original protein sequence, from where I got my initial sequences (ENST00000380152.7).

I am not too sure how to approach this. I was thinking of aligning my trimmed, translated ENSG00000139618 with the reference where I want my final positions (ENST00000380152.7). However, I do not know any software that would return my positions based on the reference.

Here I will try to show an example:

  #interesting positions are P 7 and R 17

   #interesting positions are located elsewhere and I would like those positions, which are now: P2, R14

It is not only deleting the gaps in the trimmed alignment, because there might be also insertions in the original sequence that were trimmed.

I would be grateful for any thoughts you might have on this.

ADD COMMENTlink modified 4 months ago by shenwei3565.2k • written 6 months ago by palabari10
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: 1917 users visited in the last hour