Dear Qiime2 community,
I am trying to run itsxpress plugin installed in Qiime2 2021.4 to analyse my ITS data produced by MiSeq. First I import the data, then I trim the primers with cutadapt and finally I run itsxpress:
#Importing the data
qiime tools import --type 'SampleData[PairedEndSequencesWithQuality]' --input-path $INP --input-format CasavaOneEightSingleLanePerSampleDirFmt --output-path
demux-soil_ITS.qza
#Cutting the primers
qiime cutadapt trim-paired \
--i-demultiplexed-sequences demux-soil_ITS.qza \
--p-front-f GGGTGARTCATCRARTYTTTG \
--p-front-r CCTSCSCTTANTDATATGC \
--p-cores 12 \
--o-trimmed-sequences trimmed-seqs_ITS.qza \
--verbose
#itsxpress trimming:
qiime itsxpress trim-pair-output-unmerged --i-per-sample-sequences trimmed-seqs_ITS.qza --p-region ITS2 --p-taxa F --p-threads 12 --verbose --o-trimmed trim
med_demux_ITS.qza
However, I am facing with the error pasted below:
ERROR:root:Could not perform clustering with Vsearch. Error from Vsearch was:
vsearch v2.13.3_linux_x86_64, 251.7GB RAM, 12 cores
https://github.com/torognes/vsearch
Fatal error: Files compressed with gzip are not supported
Traceback (most recent call last):
File "/hpc-home/miniconda3/envs/qiime2-2021.4/lib/python3.8/site-packages/itsxpress/main.py", line 540, in cluster
p2.check_returncode()
File "/hpc-home/miniconda3/envs/qiime2-2021.4/lib/python3.8/subprocess.py", line 448, in check_returncode
raise CalledProcessError(self.returncode, self.args, self.stdout,
subprocess.CalledProcessError: Command '['vsearch', '--cluster_size', '/tmp/itsxpress_xgc13947/seq.fq.gz', '--centroids', '/tmp/itsxpress_xgc13947/rep.fa', '
--uc', '/tmp/itsxpress_xgc13947/uc.txt', '--strand', 'both', '--id', '0.995', '--threads', '12']' returned non-zero exit status 1.
Traceback (most recent call last):
File "/hpc-home/miniconda3/envs/qiime2-2021.4/lib/python3.8/site-packages/q2cli/commands.py", line 329, in __call__
results = action(**arguments)
File "<decorator-gen-285>", line 2, in trim_pair_output_unmerged
File "/hpc-home/miniconda3/envs/qiime2-2021.4/lib/python3.8/site-packages/qiime2/sdk/action.py", line 244, in bound_callable
outputs = self._callable_executor_(scope, callable_args,
File "/hpc-home/miniconda3/envs/qiime2-2021.4/lib/python3.8/site-packages/qiime2/sdk/action.py", line 390, in _callable_executor_
output_views = self._callable(**view_args)
File "/hpc-home/miniconda3/envs/qiime2-2021.4/lib/python3.8/site-packages/q2_itsxpress/_itsxpress.py", line 151, in trim_pair_output_unmerged
results = main(per_sample_sequences=per_sample_sequences,
File "/hpc-home/miniconda3/envs/qiime2-2021.4/lib/python3.8/site-packages/q2_itsxpress/_itsxpress.py", line 208, in main
sobj.cluster(threads=threads, cluster_id=cluster_id)
File "/hpc-home/miniconda3/envs/qiime2-2021.4/lib/python3.8/site-packages/itsxpress/main.py", line 543, in cluster
raise e
File "/hpc-home/miniconda3/envs/qiime2-2021.4/lib/python3.8/site-packages/itsxpress/main.py", line 540, in cluster
p2.check_returncode()
File "/hpc-home/miniconda3/envs/qiime2-2021.4/lib/python3.8/subprocess.py", line 448, in check_returncode
raise CalledProcessError(self.returncode, self.args, self.stdout,
subprocess.CalledProcessError: Command '['vsearch', '--cluster_size', '/tmp/itsxpress_xgc13947/seq.fq.gz', '--centroids', '/tmp/itsxpress_xgc13947/rep.fa', '
--uc', '/tmp/itsxpress_xgc13947/uc.txt', '--strand', 'both', '--id', '0.995', '--threads', '12']' returned non-zero exit status 1.
Plugin error from itsxpress:
Command '['vsearch', '--cluster_size', '/tmp/itsxpress_xgc13947/seq.fq.gz', '--centroids', '/tmp/itsxpress_xgc13947/rep.fa', '--uc', '/tmp/itsxpress_xgc139
47/uc.txt', '--strand', 'both', '--id', '0.995', '--threads', '12']' returned non-zero exit status 1.
See above for debug info.
It looks like the problem is with vsearch. How can I fix it?
Thank you very much for your help!