manu check RAID GEM count and Fragment count of local FN>1

cat MANU_CHK.sh


echo "input Chr"
read C
echo "input Start"
read S
echo "input End"
read E

cat $C.bed|awk -v s=$S -v e=$E '{if($2>=s && $3<=e)print}'|cut -f5,6|sort|uniq -c|awk -v OFS="\t" '{if($1>1)print $2,$3,$1}'|awk 'BEGIN{S=0;G=0}{S=S+$3;G=G+1}END{print "total GEM count = "G"  total fragment count = "S}'
sh MANU_CHK.sh
input Chr
chr3R
input Start
3711265
input End
3953741
total GEM count = 1028  total fragment count = 2361
上一篇:Redis搭建集群


下一篇:Lab: Exploiting Ruby deserialization using a documented gadget chain:使用小工具链利用 Ruby 反序列化(Ruby 2.x Uni