Is It Possible To Align A New Sequence With An Existing Maf Multiple Alignment?
1
2
Entering edit mode
8.3 years ago
Ryan Thompson ★ 3.5k

I am trying to annotate the genome of a recently-sequenced primate species. I noticed that the UCSC browser folks have performed a 46-species whole-genome multiple alignment based on hg19, and I'd like to take advantage of that rather than starting a multi-species genome alignment from scratch. Is there a way for me to take my species' genome and align it into the existing 46-species alignment without restarting the whole thing? I am prepared to accept that the result may not be identical to what I would get if I started from scratch. Essentially my input is a set of MAF files (one for each human chromosome) and a fasta file of all the chromosomes and unlocalized scaffolds in my primate species.

alignment multiple-alignment • 3.3k views
ADD COMMENT
1
Entering edit mode

Hi! I am interested to know how you solved this issue. Thanks!

ADD REPLY
3
Entering edit mode
8.3 years ago
SES 8.5k

You can add sequences to an existing alignment using MAFFT, and they have written out a very good guide for adding sequences to an existing alignment. There is also a server for trying this out, but I'm not sure you can upload all 47 genomes to their server. Note that while this answers part of your question, I'm not sure this is the best tool for whole genome alignments. If MAFFT is not up to the task, you may want to also look into Mauve, MUMmer, or Encredo and Pecan.

ADD COMMENT
0
Entering edit mode

Hi, i read the the guides for MAFFT, but it requires the sequences or alignments with fasta format, is it any tool can using the MAF format?

ADD REPLY

Login before adding your answer.

Traffic: 2211 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6