phred64转33

zcat clean.1.fq.gz |head -12|perl -alne '{if(.%4==0){ print join"",map{chr(ord()-31)}(split//,$) }else{print}}'
引自：

http://www.bio-info-trainee.com/2946.html
©著作权归作者所有,转载或内容合作请联系作者
平台声明：文章内容（如有图片或视频亦包括在内）由作者上传并发布，文章内容仅代表作者本人观点，简书系信息发布平台，仅提供信息存储服务。

推荐阅读更多精彩内容