Hi,
I am handling different formats of Multiple sequence alignment files.
I can get the format of file by
but when I try to get sequence of different protein, I cannot get it.
I have some 10 accession ids, but unable to fetch it.
I am trying to run the script with above command in Linux, terminal
OR
I get address of objects.
I want to fetch aligned sequence of each accession/protein.
I want to print on screen not to pass to some out file.
Any help shall be appreciated.
Thanks.
I am handling different formats of Multiple sequence alignment files.
I can get the format of file by
Code:
my $format_guess = (new Bio::Tools::GuessSeqFormat( -file => $ARGV[0]))->guess;
I have some 10 accession ids, but unable to fetch it.
Code:
./perl_cript.pl file_multple_alg
Code:
my $in = Bio::AlignIO->new(-file => $ARGV[0] , -format => $format_guess);
Code:
while ( my $aln = $in->next_aln() ) { print $aln,"\n"; }
I want to fetch aligned sequence of each accession/protein.
I want to print on screen not to pass to some out file.
Any help shall be appreciated.
Thanks.
Comment