This is a test version of Biostars. For the public version, visit https://www.biostars.org.
MetaBat only produce one bin

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
• 7,422 views • link
updated by donkang34 • donkang34 • written by erikras1223 • erikras1223

4 answers

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.

• 0 views • link
written by jlfroula • jlfroula

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

• 0 views • link
written by h.l.wong h.l.wong

Nice to meet you. i Am Lilian Keita, i saw your profile today at www.biostars.org when i was browsing looking for honest partner can we be friends? i will like you to reply me through my email I will tell you more about my self when i get a reply from you OK.Best Regards. Miss Lilian Email. liliankeita17@gmail.com

• 0 views • link
written by emekannaji12 • emekannaji12

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

• 0 views • link
written by scarborough.matthew • scarborough.matthew

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

• 0 views • link
written by erikras1223 • erikras1223

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

• 0 views • link
written by donkang34 • donkang34

Log in to answer this question.

More posts like this