+ date=2019-10-01
+ data_source=data_sources
+ organism='Homo sapiens'
++ dirname ../prep/study_metadata_layout.sh
+ sdir=../prep
+ study=ERP016351
+ org=human
+ dsource=sra
+ preds_file=../prep/Human.SRA_curated_predict.meta.w_rids.rtypes.tsv
+ perl -e '$study="ERP016351"; $study=~/(..)$/; $lo=$1; `mkdir -p $lo/$study`; print "$lo/$study\n";'
++ cat ERP016351.dir
+ dir=51/ERP016351
+ cat ../prep/samples.tsv.header /dev/fd/63
++ fgrep '	ERP016351	' ../prep/samples.tsv
+ cut_cols=1-3
+ [[ human == \m\o\u\s\e ]]
+ [[ human == \h\u\m\a\n ]]
+ gzip
+ cut -f 1-41 51/ERP016351/samples.tsv.ERP016351
+ cut -f 1-3,42-150 51/ERP016351/samples.tsv.ERP016351
+ gzip
+ [[ human == \g\t\e\x ]]
+ [[ human == \t\c\g\a ]]
+ proj_header='rail_id	external_id	study	project	organism	file_source	metadata_source	date_processed'
+ dsource_path=data_sources/sra
+ extra_proj_cols='Homo sapiens	data_sources/sra	data_sources/sra	2019-10-01'
+ gzip
++ echo 'rail_id	external_id	study	project	organism	file_source	metadata_source	date_processed'
+ cat /dev/fd/63 /dev/fd/62
++ tail -n+2 51/ERP016351/samples.tsv.ERP016351
++ cut -f 1-3
++ perl -ne 'chomp; $f=$_; ($rid,$run,$study)=split(/\t/,$f); print "$f\t$study\tHomo sapiens	data_sources/sra	data_sources/sra	2019-10-01\n";'
+ [[ -n ../prep/Human.SRA_curated_predict.meta.w_rids.rtypes.tsv ]]
+ gzip
+ cat /dev/fd/63 /dev/fd/62
++ fgrep '	ERP016351	' ../prep/Human.SRA_curated_predict.meta.w_rids.rtypes.tsv
++ head -1 ../prep/Human.SRA_curated_predict.meta.w_rids.rtypes.tsv
++ perl -e '$s="ERP016351"; $s=~/(..)$/; $lo=$1; print "$lo\n";'
+ lo=51
+ gzip
+ cat /dev/fd/63 /dev/fd/62
++ head -1 ../prep/seqtk.w_study.tsv.cut
++ fgrep '	ERP016351	' ../prep/seqtk.w_study.tsv.cut
+ rm 51/ERP016351/samples.tsv.ERP016351
+ rm ERP016351.dir