This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Aligning de-novo assembly to reference strain genome

Dear Biostars Community -

I generated a de-novo assembly of my viral strain. I wanted to analyze differences of my denovo result to a reference strain, so got a very useful report from Assemblytics in conjunction with using the nucmer program from MUMmer. However the output of the alignment is in the delta format. There is delta-filterprogram to manipulate the alignment. But what I really want to do is ailgn my de-novo to the reference and have the alignment in FASTA for MAF format. The size of genome is fairly modest ( < 200kb) I was thinking of using mafft. Are there any other alignment programs that might be more appropriate?

denovo-assembly multiple-sequence-alignment
• 4,002 views • link
updated by Ram Ram • written by aindap aindap

3 answers

You may use mVISTA http://genome.lbl.gov/vista/mvista/submit.shtml or ask Inna Dubchak directly

• 0 views • link
written by Petr Ponomarenko Petr Ponomarenko

I would use stampy or bowtie2 alignment tools

• 0 views • link
written by vmicrobio vmicrobio

I actually had the same issue with nucmer. I would say Mauve is the best tool for multi-genome/assembly alignment, since it's fast and has it's own viewer. If you wanted an alignment in FASTA or MAF, mafft does a good job when there are chromosome length scaffolds, but if there are contigs spanning portions of chromosomes, I found it to struggle quite a bit. kAlign did a much better job when you don't penalize end gaps for those "floating" contigs. Mugsy uses nucmer to produce a MAF file.

• 0 views • link
written by asur • asur

Log in to answer this question.

More posts like this