SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



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

Reply
 
Thread Tools
Old 05-05-2013, 08:37 PM   #1
rpbaptista
Junior Member
 
Location: Athens, GA

Join Date: Mar 2013
Posts: 2
Question 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

rpbaptista is offline   Reply With Quote
Old 04-28-2014, 10:56 PM   #2
tw7649116
Junior Member
 
Location: Yesan,Korea

Join Date: Mar 2014
Posts: 8
Default

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!
tw7649116 is offline   Reply With Quote
Old 04-29-2014, 07:51 AM   #3
rpbaptista
Junior Member
 
Location: Athens, GA

Join Date: Mar 2013
Posts: 2
Default

Hi tw7649116,

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

hope it helps you

rpbaptista is offline   Reply With Quote
Old 04-29-2014, 01:54 PM   #4
Matt Kearse
Member
 
Location: New Zealand

Join Date: Mar 2014
Posts: 20
Default

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.
Matt Kearse is offline   Reply With Quote
Old 04-29-2014, 04:56 PM   #5
tw7649116
Junior Member
 
Location: Yesan,Korea

Join Date: Mar 2014
Posts: 8
Default

Thank you very much!I'll try!!
tw7649116 is offline   Reply With Quote
Old 04-29-2014, 11:44 PM   #6
tw7649116
Junior Member
 
Location: Yesan,Korea

Join Date: Mar 2014
Posts: 8
Default

Quote:
Originally Posted by rpbaptista View Post
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!
tw7649116 is offline   Reply With Quote
Reply

Tags
alignment, consensus sequence, fasta, mummer, reference guided assembly

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 06:33 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO