study=$1
#sra, gtex, or tcga
dsource=$2
#RNA-seq, tissue, celltype Predictions/curations file
#e.g. Human.SRA_curated_predict.meta.w_rids.tsv
#or
#Mouse.SRA_curated_predict.meta.w_rids.w_cell_types_nas.tsv
preds_file=$3

perl -e '$study="'$study'"; $study=~/(..)$/; $lo=$1; print "$lo/$study\n";' > ${study}.dir
dir=`cat ${study}.dir`

cat <(head -1 $preds_file) <(fgrep "	$study	" $preds_file) | gzip > $dir/${dsource}.recount_pred.${study}.MD.gz