Question: Software For Gene Recontruction From Wgs Data
1
gravatar for Pawel Szczesny
5.6 years ago by
Pawel Szczesny3.2k
Poland
Pawel Szczesny3.2k wrote:

Many of you are familiar with EMIRGE, a software for reconstruction 16S rRNA genes from shotgun data. I've been trying to use EMIRGE for my own research, however the software is buggy and not that well supported. I've seen Megraft but authors claim it's not going to work well with data from Illumina.

Are there any other packages worth trying? I like a lot the idea of targeted reconstruction (saves time and increases sensitivity), but there no software packages to choose from to solve that task.

software ngs • 1.6k views
ADD COMMENTlink modified 5.6 years ago by Lee Katz2.9k • written 5.6 years ago by Pawel Szczesny3.2k
2
gravatar for Lee Katz
5.6 years ago by
Lee Katz2.9k
Atlanta, GA
Lee Katz2.9k wrote:

I think I've done something like what you want. Just create a multifasta file of genes that you would like to recover and then map your reads to it using something like Smalt, Bowtie, or BWA. Create a consensus sequence of the output sam files using samtools and vcfutils.pl.

I think that Nesoni could do this all-in-one using the shrimp mapper and samtools.

ADD COMMENTlink written 5.6 years ago by Lee Katz2.9k

Thanks. This might work for genes which are relatively uniform in terms of selective pressure. For 16S rRNA I don't see how could I map the reads and accurately recover consensus sequence given the structure of the gene (highly conserved regions mixed with hypervariable). But still, it's a good start.

ADD REPLYlink written 5.6 years ago by Pawel Szczesny3.2k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1902 users visited in the last hour