bxlab / metaWRAP

MetaWRAP - a flexible pipeline for genome-resolved metagenomic data analysis
MIT License
390 stars 191 forks source link

How to split reassembled bins #549

Open Qilong50281 opened 2 months ago

Qilong50281 commented 2 months ago

Dear Developers

Thank you very much for the development of metawrap. I have completed the analysis according to the tutorial, and I have a question about Step 8.

Assume that I have 3 samples (for example, sample A, sample B, sample C). At the beginning, all R1 and R2 are concatenated, and the BIN_REASSEMBLY folder is generated in Step 8, and the reassembled_bins in it is used for taxonomic naming in Step 9.

Then my question is, if BIN_REASSEMBLY/reassembled_bins/bin.1.strict.fa is the combined assembly result containing all these three samples, how to split bin.1.strict.fa into results like bin.1.strict_sample A.fa bin.1.strict_sample B.fa bin.1.strict_sample C.fa according to the marker or index?

Because I found that the classified microorganisms (based on Step 9) represented by bin1 have differences in abundance in the heat map of bin_abundance_heatmap.png (via Step 7), and I want to compare it with the reference genome to see the differences at the nucleotide level.

Hope to get your advice. Thank you in advance for taking the time to check

Best regards, Chris

yqy6611 commented 1 month ago

Just assemble and bin them one sample after another.