How Do I Align More Than Five Long Dna Scaffolds (Contigs) To The Reference Gene?
2
7
Entering edit mode
7.7 years ago
biostar ▴ 150

Hi everyone, I tried to align one of my reference genes with its contigs using Mauve, MAFFT, emboss needle etc, but could not get it aligned locally (I also needed to show the aligned sequences) . Is there any program I can use to do the local pairwise alignment of multiple sequences?? I just need to map the contigs to my reference gene and find the UTR regions and introns. Any inputs would be appreciated on this.. Thanks!!

alignment pairwise reference • 6.1k views
ADD COMMENT
2
Entering edit mode

you may want to try MUMmer - plotting your gene(s) on Y and contigs on X. the tool would plot diagonal bands (where the common elements are) which can shed some light on how things are arranged and assembled.

enter image description here

ADD REPLY
0
Entering edit mode

Thanks Pavel! I already tried this, but I also need to show the alignment of all sequences locally. I might need to do it manually using SeAl..

ADD REPLY
0
Entering edit mode

hi,

Which command line did you use to generate a plot like this? Thanks

ADD REPLY
3
Entering edit mode
7.7 years ago

I would suggest lastz:

http://www.bx.psu.edu/~rsharris/lastz/

I have used it to align whole chromosomes between divergent species. It has a flexible output and is efficient.

ADD COMMENT
0
Entering edit mode

Thank you Zev, but this program somehow generates error on installation..

ADD REPLY
1
Entering edit mode

Found this tool on NCBI, where I can input reference sequence and target sequences (Multiple) and gives nice result of pairwise alignment (local): http://blast.ncbi.nlm.nih.gov/Blast.cgi?PAGE_TYPE=BlastSearch&BLAST_SPEC=blast2seq&LINK_LOC=align2seq

ADD REPLY
3
Entering edit mode
6.1 years ago
Darked89 4.2k

Late but hopefully not too late. One can map anything from genome to genome (and obviously genes to contigs) up to FASTQ short reads using LAST: http://last.cbrc.jp/

LAST has tools to convert the output from MAF to SAM, so results are viewable in IGV. And using another path, it can output dotplots as well (see the web page).

ADD COMMENT

Login before adding your answer.

Traffic: 2333 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