[Bioperl-l] how to not count gaps in the multiple sequence alignment

[Remo Sanges]

To get the location in the initial sequence starting from a column in a 
multiple alignment you can:

1) create a Bio::LocatableSeq compliant object by using the method 
each_seq_with_id on the SimpleAlign object

2) then using the method location_from_column on the created 
LocatableSeq object

HTH

ERemo


-- 
Remo Sanges
Bioinformatics - Animal Physiology and Evolution
Stazione Zoologica Anton Dohrn
Villa Comunale, 80121 Napoli - Italy
+39 081 5833428


On 11/2/11 5:29 PM, Jun Yin wrote:
> Hi,
>
> You need to calculate the coordinates of the protein coding gene in the alignment by yourself. After that, you can use the slice function to get the alignment block for the selected gene, e.g.
>
> $aln2 = $aln->slice(20, 30);
>
> Cheers,
> Jun
>
>
> ----- Original Message -----
> From: wenbin mei<wenbinmei at gmail.com>
> Date: Wednesday, November 2, 2011 5:51 am
> Subject: [Bioperl-l] how to not count gaps in the multiple sequence alignment
> To: bioperl-l at lists.open-bio.org
>
>> Hi,
>>
>> I need some help in coding. I have a multiple sequence alignment
>> which has
>> gaps. And also I have a reference genome sequence in the
>> alignment which I
>> know all the coordinates for the protein coding genes. I want to
>> extractall these protein coding genes alignment from the big
>> alignment. I am using
>> Bio SimpleAlign but the question is that due to the gaps in the
>> alignment,the coordinates has shifted in the alignment. I wonder
>> is there a way I can
>> not count the gaps and still be able to extract the protein
>> alignment. One
>> way I can do is remove the gaps in the reference first and then
>> extract the
>> sequence. But I don't like this way ... Thank you for help.
>>
>> -best,
>> wenbin
>> _______________________________________________
>> Bioperl-l mailing list
>> Bioperl-l at lists.open-bio.org
>> http://lists.open-bio.org/mailman/listinfo/bioperl-l
> _______________________________________________
> Bioperl-l mailing list
> Bioperl-l at lists.open-bio.org
> http://lists.open-bio.org/mailman/listinfo/bioperl-l