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

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

现在只好自己写了

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

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

ascii码对照表

居然说字数不够300字,所以没办法申请原创,难道技术的高低是由废话的多少来评价的吗?

居然说字数不够300字,所以没办法申请原创,难道技术的高低是由废话的多少来评价的吗?

居然说字数不够300字,所以没办法申请原创,难道技术的高低是由废话的多少来评价的吗?

居然说字数不够300字,所以没办法申请原创,难道技术的高低是由废话的多少来评价的吗?

居然说字数不够300字,所以没办法申请原创,难道技术的高低是由废话的多少来评价的吗?

居然说字数不够300字,所以没办法申请原创,难道技术的高低是由废话的多少来评价的吗?

(0)

相关推荐