Convert mothur classify.seqs output to BIOM format
1
0
Entering edit mode
5.0 years ago
pkozyrev ▴ 10

I got a lot of problems trying to convert output of mothur classify.seqs to BIOM format so that I can import data into phyloseq R package. I have FASTA file with paired-end reads already being merged with PEAR program.

I aligned my sequences to SILVA database with this command

mothur > align.seqs(fasta=sample.fasta, reference=silva.bacteria.fasta, processors=4, flip=t)

and then I classified it using the same database

mothur > classify.seqs(fasta=sample.align, reference=silva.bacteria.fasta, taxonomy=silva.bacteria.silva.tax, cutoff=80)

After this procedure I got different files: "sample.summary", "sample.silva.wang.taxonomy" and "sample.silva.wang.tax.summary". But I don't know how to import them into phyloseq R package.

I've read here https://github.com/joey711/phyloseq/issues/245 about shared file from which I can create BIOM file, but I don't have .list or .group files for make.shared command.

Can someone help?

Sincerely yours,

Petr

metagenomics mothur R phyloseq • 2.8k views
ADD COMMENT
2
Entering edit mode

1) Why would you use pear? when mothur has make.contigs, 2) When you say you have fasta file do you mean you have fastq file? 3) what is the error message you get when you use this command

make.biom(shared=final.tx.1.subsample.1.pick.shared, contaxonomy=final.tx.1.cons.taxonomy)

ADD REPLY
0
Entering edit mode

Thank you for your quick response!

1) Ah, sorry, I didn't mention that I omit some steps in mothur pipeline because I wanted to compare only RDPTools classifier and mothur classifier performance

2) As far as I understood, mothur align.seqs and mothur classify.seqs commands work with fasta files so I converted my fastq file with merged reads to fasta

3) The problem is that I don't have shared file because I don't have .list and .group files to make it

Petr

ADD REPLY
0
Entering edit mode

crosspost http://seqanswers.com/forums/showthread.php?t=75774

I guess you got your answer, need to make otus..

ADD REPLY
0
Entering edit mode

Yeah, it isn't a good practice. My bad.

Thank you for your help anyway!

ADD REPLY
1
Entering edit mode
5.0 years ago
h.mon 34k

To make a .list file, see cluster, cluster.split or phylotype.

To make a .group file, see make.group.

ADD COMMENT
0
Entering edit mode

Thank you!

I've already understood.

ADD REPLY

Login before adding your answer.

Traffic: 2180 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6