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.
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 ?
•
written
by
Arpitha Hebbar
2 answers
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.
i used it worked well....
Log in to answer this question.
-
Percentage coverage of reference genome by de novo genome
written by Lemonhope •
-
filtering alignned seqeunces
written by priscilladraj •
-
Align pacbio subreads against assembled pacbio reads
written by wes •
-
Obtaining only the aligned regions in BLAST
written by arya •
-
How to compare contigs with an existing protein from a species ?
written by carerjose •
-
How To Use Coordinates In Order To Extract Sequences In Fasta File?
written by srpatel144 •
-
Is There A Tutorial/Guide Explaining How To Create A Draft Genome From Assembled Contigs?
written by Arpitha Hebbar
-
How Can I Do Scaffolding With The Single End Data ?
written by Arpitha Hebbar
-
How To Do Scaffolding Using Mira For Ion Torrent Data?
written by Arpitha Hebbar
-
About Four Hundreds Contigs Mapping
written by wHo •
