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!!
2 answers
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.
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
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).
Log in to answer this question.
-
Beginner problems: Aligning "genomes" that come in contigs
written by ABS •
-
Aligning two genomes
written by Nicolas Carter •
-
Intron finder for 3' UTR sequences (non-model organism)?
written by LAB_2019 •
-
How to align new DNA sequences with existing multiple DNA sequence alignment
written by liaoyunshi •
-
reads distribution from the reference
written by Alex •
-
local alignment with a sets of sequences in a fasta file to a reference sequence.
written by kanshenglong •
-
How Can I Find Utr Regions Of This Gene?
written by biostar •
-
Any Free Software Similar To Sequencher?
written by biostar •
-
Codon Alignment Via Python?
written by a1ultima
-
Mostly Used (Fast) Global Pairwise Alignment Tool
written by Pdubois •
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.
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..
hi,
Which command line did you use to generate a plot like this? Thanks