MetaBat only produce one bin
4
0
Entering edit mode
8.5 years ago
erikras1223 ▴ 10

I have just started using MetaBat and have tried follow the tutorial found hear https://bitbucket.org/berkeleylab/metabat/wiki/Home . I have supplied it with my assembled contigs file using megahit. I also have a bam file but am not aware of how to use it as input for metabat.

Here are the parameters used, note I did both --sensitive and --specific although I am only putting one command:

metabat -i contigs.fa -o bin1 --sensitive -v --saveDistance saved.gprob

The result I get is one bin file of contig sequences that contains almost exactly what the contig.fa file holds. I have used this same contigs.fa file in Anvio using its default clustering algorithm and can find at least Wolbachia (the most abundant bacteria in my sample) and few other bacteria at low abundance. Please let me know if I am missing something obvious I would love to use metabat since it seems pretty fast and efficent. Thanks in advance

metabat clustering TNF BIN • 6.1k views
ADD COMMENT
1
Entering edit mode
7.9 years ago
jlfroula ▴ 10

Yes...you need to run runMetaBat.sh unless you already have a "depth.txt" file. There are two ways to run metabat. The easiest is
runMetaBat.sh <assembly.fa> < *.bam >

The other way is generate the depth.txt file by doing . jgi_summarize_bam_contig_depths --outputDepth depth.txt Bams/*.bam
then
metabat -i assembly.fasta -a depth.txt -o bin

The "depth.txt" file shows the average read depths per contig for each sample. So if you had 5 samples, you'd have hopefully aligned each of the 5 read samples to the one assembly and get 5 bam files. These 5 bam files are input to either runMetaBat.sh or jgi_summarize_bam_contig_depths.

ADD COMMENT
0
Entering edit mode

Thanks jlfroula, do you need to sort the bam files first?

ADD REPLY
0
Entering edit mode
7.9 years ago

Hello! I am having the same problem as you. Did you ever find a fix?

ADD COMMENT
1
Entering edit mode

I believe i had to run metabat.sh not metabat.

ADD REPLY
0
Entering edit mode

Alright-- thanks for responding!

ADD REPLY
0
Entering edit mode
7.4 years ago
donkang34 • 0

We recommend you to try MetaBAT 2 which performs better with a few samples. And you don't need to tune parameters any more.

ADD COMMENT

Login before adding your answer.

Traffic: 1381 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