美文网首页RNA-seq转录组软件学习
快捷判断测序数据为Phred33还是64

快捷判断测序数据为Phred33还是64

作者: 热衷组培的二货潜 | 来源:发表于2018-08-09 14:44 被阅读0次

脚本来源https://www.plob.org/article/11672.html

PS: 矿工写的一系列文章都有非常详细

$ cat fq_qual_type.sh 
less $1 | head -n 1000 | awk '{if(NR%4==0) printf("%s",$0);}' \
| od -A n -t u1 -v \
| awk 'BEGIN{min=100;max=0;} \
{for(i=1;i<=NF;i++) {if($i>max) max=$i; if($i<min) min=$i;}}END \
{if(max<=126 && min<59) print "Phred33"; \
else if(max>73 && min>=64) print "Phred64"; \
else if(min>=59 && min<64 && max>73) print "Solexa64"; \
else print "Unknown score encoding"; \
print "( " min ", " max, ")";}'

测试脚本

$ sh fq_qual_type.sh R1.fq 
Phred33
( 44, 75 )

相关文章

网友评论

    本文标题:快捷判断测序数据为Phred33还是64

    本文链接:https://www.haomeiwen.com/subject/aezqbftx.html