+ date=2019-10-01
+ data_source=data_sources
+ organism='Homo sapiens'
++ dirname ../prep/study_metadata_layout.sh
+ sdir=../prep
+ study=DRP004845
+ org=human
+ dsource=sra
+ preds_file=../prep/Human.SRA_curated_predict.meta.w_rids.rtypes.tsv
+ perl -e '$study="DRP004845"; $study=~/(..)$/; $lo=$1; `mkdir -p $lo/$study`; print "$lo/$study\n";'
++ cat DRP004845.dir
+ dir=45/DRP004845
+ cat ../prep/samples.tsv.header /dev/fd/63
++ fgrep '	DRP004845	' ../prep/samples.tsv
+ cut_cols=1-3
+ [[ human == \m\o\u\s\e ]]
+ [[ human == \h\u\m\a\n ]]
+ gzip
+ cut -f 1-41 45/DRP004845/samples.tsv.DRP004845
+ cut -f 1-3,42-150 45/DRP004845/samples.tsv.DRP004845
+ 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
+ cat /dev/fd/63 /dev/fd/62
++ echo 'rail_id	external_id	study	project	organism	file_source	metadata_source	date_processed'
++ tail -n+2 45/DRP004845/samples.tsv.DRP004845
++ 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
++ head -1 ../prep/Human.SRA_curated_predict.meta.w_rids.rtypes.tsv
++ fgrep '	DRP004845	' ../prep/Human.SRA_curated_predict.meta.w_rids.rtypes.tsv
++ perl -e '$s="DRP004845"; $s=~/(..)$/; $lo=$1; print "$lo\n";'
+ lo=45
+ gzip
+ cat /dev/fd/63 /dev/fd/62
++ head -1 ../prep/seqtk.w_study.tsv.cut
++ fgrep '	DRP004845	' ../prep/seqtk.w_study.tsv.cut
+ rm 45/DRP004845/samples.tsv.DRP004845
+ rm DRP004845.dir