get_expr.sh 381 Bytes
Newer Older
Marine Bergot's avatar
Marine Bergot committed
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17
#!/bin/bash


WHERE=$1

for file in $WHERE/*.txt
        do
        TAD=$(basename $file _list_genes.txt)
	mkdir -p $WHERE/$TAD
	for gene in $(cat $file)
		do
		head -n 1 /work/gad/shared/analyse/Integration/RNA_SEQ/normalized_counts.csv >> $WHERE/$TAD/$gene.txt
		grep $gene /work/gad/shared/analyse/Integration/RNA_SEQ/normalized_counts.csv >> $WHERE/$TAD/$gene.txt
	done
done