[Bioperl-l] Newbie alert: Aligning and printing first position
Diego M. Vadell
dvadell at clustering.com.ar
Wed Mar 30 14:51:38 UTC 2011
Hi list,
I have an alignment ( done with Bio::Tools::Run::Alignment::Clustalw )
that I can see ( printing $aln->match_line() along with the sequences ):
********* .: (match_line)
.......TDDSFDTDVSSS.........[...]............. ($seq1)
SDKIIHLTDDSFDTDVLKADGA[...]EFLDANLA ($seq2)
It looks right ( that "[...]" is mine, it's larger ). Now I want to get
the position of the match ( the first "T" ), which should be "7". I thought I
needed to use something like:
my $seq = $aln->get_seq_by_pos(1);
my $loc = $seq->location_from_column(1);
say $loc;
And $loc would be 7, but it's undef. If I try to iterate from 1 to 10, I start
getting Bio::Location::Simple objects at the 8th place. So, should I loop
untill I location_from_column() gives me an object? Is it another way of
finding the position of the match?
As you can see, I'm just starting with bioinformatics. Any pointer will be
appreciated.
Thanks in advance,
-- Diego
More information about the Bioperl-l
mailing list