[ensembl-dev] Multiple alignment for more than 2 species

Dazong Zhou zhoudazong at gmail.com
Wed Sep 10 09:42:50 BST 2014


Hi,

I get problems on multiple sequence alignment for more than 2 species using
the sample code at:
http://www.ensembl.org/info/docs/api/compara/compara_tutorial.html (the
code below "Here is an example script with all of this:").

The sample runs well on default, but fails on adding additional species.
That is to say, the below code works when changing "human:mouse" to
"human:dog" but returns no data when changing to "human:mouse:dog" .


my $species = "human";
my $coord_system = "chromosome";
my $seq_region = "14";
my $seq_region_start = 75000000;
my $seq_region_end = 75010000;
my $alignment_type = "LASTZ_NET";
my $set_of_species = "human:mouse";
my $output_file = undef;
my $output_format = "clustalw";
my $help;

GetOptions(
    "help" => \$help,
    "species=s" => \$species,
    "coord_system=s" => \$coord_system,
    "seq_region=s" => \$seq_region,
    "seq_region_start=i" => \$seq_region_start,
    "seq_region_end=i" => \$seq_region_end,
    "alignment_type=s" => \$alignment_type,
    "set_of_species=s" => \$set_of_species,
    "output_format=s" => \$output_format,
    "output_file=s" => \$output_file);

Thanks,
Hunter
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.ensembl.org/pipermail/dev_ensembl.org/attachments/20140910/a0fed330/attachment.html>


More information about the Dev mailing list