trimmomatic竟然“无法解析”出gz压缩文本了,你遇到这个问题了吗?老老实实换回fastq-dump吧

在处理TCR beta序列时,使用fasterq-dump转换序列后为节省空间进行了.gz压缩。但在使用trimmomatic进行去接头操作时遇到报错,提示fastq文件解析错误。尝试了修改文件后缀和使用未压缩文件,问题依旧存在。最终发现只有使用未压缩的原始文件才能使trimmomatic正常运行,作者决定回归使用fastq-dump。如果你也遇到类似问题,可能需要寻找解决方案或考虑更换工具。
摘要由CSDN通过智能技术生成

下载了一个TCR beta的序列,用fasterq-dump转换序列,这个比fastq-dump 稍微快一些,不过没有了压缩功能,所以为了占用少的硬盘空间,在解析后需要自己压缩成gz或bz2。但是神奇一幕发生了,我手动压缩成.gz格式后,用trimmomatic进行去接头,但是死活就是运行不下去,并出现如下报错:

trimmomatic PE -threads 2 -phred33 SRR4449813.sra_1.fq.gz SRR4449813.sra_2.fq.gz paired_1_R1_paired.fq.gz unpaired_1_R1_unpaired.fq.gz paired_1_R2_paired.fq.gz unpaired_1_R2_unpaired.fq.gz ILLUMINACLIP:/home/decen/miniconda3/pkgs/trimmomatic-0.38-1/share/trimmomatic-0.38-1/adapters/TruSeq3-PE.fa:2:30:10 LEADING:3 TRAILING:3 SLIDINGWINDOW:4:15 MINLEN:36
TrimmomaticPE: Started with arguments:
 -threads 2 -phred33 SRR4449813.sra_1.fq.gz SRR4449813.sra_2.fq.gz paired_1_R1_paired.fq.gz unpaired_1_R1_unpaired.fq.gz paired_1_R2_paired.fq.gz unpaired_1_R2_unpaired.fq.gz ILLUMINACLIP:/home/decen/miniconda3/pkgs/trimmomatic-0.38-1/share/trimmomatic-0.38-1/adapters/TruSeq3-PE.fa:2:30:10 LEADING:3 TRAILING:3 SLIDINGWINDOW:4:15 MINLEN:36
Using 
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值