生信分析經常使用腳本（一）

1.計算reads PE長度腳本html 01.Cal_PE_Depth.shpost Read_count=`gzip -dc Reads1.fq.gz |wc -l |awk '{print $1/4}'` echo "read pair count : $Read_count" average_depth=`expr $Read_count \* 200 / 6000000 ` echo a

>>阅读原文<<