Question: How To Align The Reference Sequence With The Contigs ?
gravatar for Arpitha Hebbar
7.2 years ago by
Arpitha Hebbar200 wrote:

i have a fasta file with 113 contigs and a reference sequence in fasta file. i have to align my contigs to the reference sequence and i have to see in which order my contigs are aligned to the reference sequence. pls can you help me which software i can use for this ?

ngs alignment contigs • 11k views
ADD COMMENTlink modified 7.2 years ago by sentausa640 • written 7.2 years ago by Arpitha Hebbar200
gravatar for Michael Dondrup
7.2 years ago by
Bergen, Norway
Michael Dondrup47k wrote:

You have to make up your mind: you can either use global-local alignment or local alignments. If you are quite sure, the contigs are not misassembled, you can use global-local alignment. You have only few contigs, therefore you can possibly use glsearch or if that doesn't work MUMer. On the other hand you won't necessarily find misassembled contigs that way, e.g. a contig that alignes only partially in two or more sections, and without testing for such contigs, you cannot know if your assembly is good. Therefore, use also local alignment tools, e.g. FASTA or BLAST.

ADD COMMENTlink modified 7.2 years ago • written 7.2 years ago by Michael Dondrup47k

For genome-wide contigs, misassemblies and translocations almost always happen, which will break glocal alignment; if contigs can really be aligned in one piece, proper local alignment can identify the correct alignment at a high chance. There is no need to try glocal alignment. As to programs, MUMmer is usually the first choice for bacterial (btw, MUMmer reports local hits). Cross_match also works well for me in small scale. With FASTA, you cannot align >40kb query, IIRC. BLAST is slow and usually produces very fragmented alignments under the default setting which are hard to process. I would not use these two. I have heard good things about Mauve, but have not tried myself. There are also bwa-sw and bwa-mem for mammalian genomes (my programs), but they are less used.

ADD REPLYlink modified 7 months ago by RamRS27k • written 7.2 years ago by lh332k
gravatar for sentausa
7.2 years ago by
sentausa640 wrote:

I've used the "move contigs" feature in Mauve to do similar thing for my bacterial contigs and it worked really well.

ADD COMMENTlink written 7.2 years ago by sentausa640

i used it worked well....

ADD REPLYlink written 7.1 years ago by Arpitha Hebbar200
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1157 users visited in the last hour