Post navigation ← Previous Next → 找了半天居然没有找到phred64to33的轮子 Posted on 2018年1月3日 by ulwvfje 不知道从哪里拿到了一个上古测序数据,居然是用phred64来代表测序质量分数,可怕,导致我后面的流程报错,所以得先转换过来。结果Google了好久,都没有发现现成的轮子,可能是太简单了吧这个需求,本来想偷懒的。 现在只好自己写了 zcat clean.1.fq.gz |head -12|perl -alne '{if($.%4==0){ print join"",map{chr(ord($_)-31)}(split//,$_) }else{print}}' 感觉代码写的不够优雅,希望有朋友可以帮忙改进一下。