SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
BWA and masking the reference genome jullee Bioinformatics 2 09-18-2014 03:25 AM
bwa reference genome alignment nxtgenkid10 Genomic Resequencing 2 08-29-2013 04:12 AM
Alignment to a set of custom reference sequences along with standard genome reference eeyun Bioinformatics 4 05-08-2013 04:06 PM
Please help: imperfect reference genome/get consensus on genome/read alignment? KAP Bioinformatics 1 08-19-2011 07:14 AM
alignment and reference genome assembling slny Bioinformatics 3 03-16-2011 06:49 AM

Reply
 
Thread Tools
Old 01-20-2015, 03:41 PM   #1
hellsingwyk
Member
 
Location: HK

Join Date: Sep 2012
Posts: 18
Default How to rebuild query genome after BWA alignment to reference genome?

Dear all,

I recently sequence whole genome of Mycobacterium tuberculosis(unknown strain) with around 10X to 25X. Is it possible to rebuild query genome(genome of my target sample) after bwa alignment with MTB reference ?

And What kind of other tools or setting I need to used?

Since there are other unknown strains which would be sequenced for phylogenetic analysis, their genome sequence is necessary to be retrieved.

Thank you very much for all your kindly attention!
hellsingwyk is offline   Reply With Quote
Old 01-20-2015, 04:17 PM   #2
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 6,978
Default

Since you have good coverage you may be able to de novo assemble your strain's genome from the data you have and then compare it to known strains of TB.

SPAdes is a good choice for bacterial genome assembly assuming you have the right kind of data (illumina).
GenoMax is offline   Reply With Quote
Old 01-20-2015, 04:26 PM   #3
hellsingwyk
Member
 
Location: HK

Join Date: Sep 2012
Posts: 18
Default

Thank you very much for your reply!!

However, is there any possible way to obtain target genome after alignment?
hellsingwyk is offline   Reply With Quote
Old 01-20-2015, 04:40 PM   #4
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 6,978
Default

Quote:
Originally Posted by hellsingwyk View Post
Thank you very much for your reply!!

However, is there any possible way to obtain target genome after alignment?
You can create a consensus sequence based on the alignments. It would only be valid for those regions where your strain matches the reference.

http://samtools.sourceforge.net/mpileup.shtml

If you use new samtools (v.1.x): http://www.htslib.org/doc/samtools-1.1.html
GenoMax is offline   Reply With Quote
Old 01-20-2015, 05:09 PM   #5
hellsingwyk
Member
 
Location: HK

Join Date: Sep 2012
Posts: 18
Default

Dear GenoMax,

Thank you very much for your helpful reply!!!
I suppose my sequencing depth is suffice to cover whole reference genome.
And will that consensus sequence contain the snp and indel from my query genome?
hellsingwyk is offline   Reply With Quote
Reply

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 07:35 PM.


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