Samfile error when applying strainphlan to output .bz2 files of metaphlan

Hi,
I am trying to use strainphlan to build a phylogenetic tree, but when I apply sample2markers.py to the out put of metaphlan, I got errors:
Processing sample: /mnt/vstor/CSE_CSDS_VXC204/rxl761/consensus_markers/tmpxbb3zv62/11808.3205.SALB.W2V1.bam
Fri Feb 16 14:33:35 2024: Loading the bam file and extracting information…
Traceback (most recent call last):
File “/mnt/vstor/CSE_CSDS_VXC204/rxl761/anaconda3/envs/metaphlan/bin/sample2markers.py”, line 8, in
sys.exit(main())
File “/mnt/vstor/CSE_CSDS_VXC204/rxl761/anaconda3/envs/metaphlan/lib/python3.7/site-packages/metaphlan/utils/sample2markers.py”, line 399, in main
sampletomarkers.run_sample2markers()
File “/mnt/vstor/CSE_CSDS_VXC204/rxl761/anaconda3/envs/metaphlan/lib/python3.7/site-packages/metaphlan/utils/sample2markers.py”, line 281, in run_sample2markers
filtered=‘_filtered’ if len(self.clades) > 0 else ‘’)
File “/mnt/vstor/CSE_CSDS_VXC204/rxl761/anaconda3/envs/metaphlan/lib/python3.7/site-packages/metaphlan/utils/sample2markers.py”, line 119, in build_consensus_markers
results = self.get_consensuses_for_sample(i)
File “/mnt/vstor/CSE_CSDS_VXC204/rxl761/anaconda3/envs/metaphlan/lib/python3.7/site-packages/metaphlan/utils/sample2markers.py”, line 141, in get_consensuses_for_sample
sam_file = pysam.AlignmentFile(input_bam)
File “pysam/libcalignmentfile.pyx”, line 748, in pysam.libcalignmentfile.AlignmentFile.cinit
File “pysam/libcalignmentfile.pyx”, line 997, in pysam.libcalignmentfile.AlignmentFile._open
ValueError: file has no sequences defined (mode=‘r’) - is it SAM/BAM format? Consider opening with check_sq=False

I used conda to install metaphlan, strainphlan and samtools, and my version for metaphlan is 4.0.6, strainphlan is 4.0.6, and samtools is 1.9

Hi burdungy,

it seems the file you’re providing as input (11808.3205.SALB.W2V1.bam) is “empty”. You can verify it by running samtools view on the file. The bam file was generated by metaphlan? (usually it’s .sam.bz2 instead of .bam) Do you have the profile output from metaphlan? Does it contain any species? Either something went wrong with the metaphlan run or there’s not enough bacterial reads in the sequencing file.

Best
Michal