Getting the consensus file from a reference based alignment

rpbaptista · 05-05-2013, 08:37 PM

Hi,

I have a Mapped genome fasta file with chromosomes scaffolds and a denovo contigs sequences all together in a Multifasta file.
I aligned them to a multifasta reference genome (each sequence is a chromosome) using Nucmer (from Mummer), and some of my gaps were filled as I expected for the divergences , but I couldnt take the consensus fasta file from this Nucmer alignment, does someone have an idea about how could I make it? Or other genome aligner that works like Nucmer that gave me the same result?
(My genome is highly similar to this Reference)
thanks

tw7649116 · 04-28-2014, 10:56 PM

Dear rpbaptista:
I have the same problem now. How can you solve the problem?
I used SOAPdenovo get the scafseq,how can I use the reference to get the consensus sequence?
Thank you！

rpbaptista · 04-29-2014, 07:51 AM

Hi tw7649116,

Try to use ABACAS. http://abacas.sourceforge.net/

hope it helps you

Matt Kearse · 04-29-2014, 01:54 PM

The free (unlicensed) version of Geneious can do what you want. It also has an option to fill the consensus with either the reference sequence bases or Ns when either the coverage or quality of the aligned sequences are below a specified threshold.

tw7649116 · 04-29-2014, 04:56 PM

Thank you very much！I'll try!!

tw7649116 · 04-29-2014, 11:44 PM

Quote:

Originally Posted by rpbaptista

Hi tw7649116,

Try to use ABACAS. http://abacas.sourceforge.net/

hope it helps you

I have many sequences need to be aligned to get the consensus sequence.
Is there any command line software can do this?
Thank you!

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
How to build reference genome for targeted region(from a BED file) for alignment?	serenaliao	Epigenetics	6	08-07-2015 02:44 PM
vcf to consensus call reference instead of N	spike1985	Bioinformatics	5	09-17-2012 09:29 AM
how to create junctions file based off reference	zorph	Bioinformatics	0	11-08-2011 09:42 AM
Please help: imperfect reference genome/get consensus on genome/read alignment?	KAP	Bioinformatics	1	08-19-2011 08:14 AM
Help needed to generate Consensus from alignment	sarbashis	Illumina/Solexa	0	08-10-2011 10:58 PM

05-05-2013, 08:37 PM	#1
rpbaptista Junior Member Location: Athens, GA Join Date: Mar 2013 Posts: 2	Getting the consensus file from a reference based alignment Hi, I have a Mapped genome fasta file with chromosomes scaffolds and a denovo contigs sequences all together in a Multifasta file. I aligned them to a multifasta reference genome (each sequence is a chromosome) using Nucmer (from Mummer), and some of my gaps were filled as I expected for the divergences , but I couldnt take the consensus fasta file from this Nucmer alignment, does someone have an idea about how could I make it? Or other genome aligner that works like Nucmer that gave me the same result? (My genome is highly similar to this Reference) thanks

04-28-2014, 10:56 PM	#2
tw7649116 Junior Member Location: Yesan,Korea Join Date: Mar 2014 Posts: 8	Dear rpbaptista: I have the same problem now. How can you solve the problem? I used SOAPdenovo get the scafseq,how can I use the reference to get the consensus sequence? Thank you！

04-29-2014, 07:51 AM	#3
rpbaptista Junior Member Location: Athens, GA Join Date: Mar 2013 Posts: 2	Hi tw7649116, Try to use ABACAS. http://abacas.sourceforge.net/ hope it helps you

04-29-2014, 01:54 PM	#4
Matt Kearse Member Location: New Zealand Join Date: Mar 2014 Posts: 20	The free (unlicensed) version of Geneious can do what you want. It also has an option to fill the consensus with either the reference sequence bases or Ns when either the coverage or quality of the aligned sequences are below a specified threshold.

04-29-2014, 04:56 PM	#5
tw7649116 Junior Member Location: Yesan,Korea Join Date: Mar 2014 Posts: 8	Thank you very much！I'll try!!