找了半天居然没有找到phred64to33的轮子

不知道从哪里拿到了一个上古测序数据,居然是用phred64来代表测序质量分数,可怕,导致我后面的流程报错,所以得先转换过来。结果Google了好久,都没有发现现成的轮子,可能是太简单了吧这个需求,本来想偷懒的。


现在只好自己写了

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

感觉代码写的不够优雅,希望有朋友可以帮忙改进一下。
ascii

 

Comments are closed.