# reads a column file with gene names of interest

awk '{print $1, NR}' $1 | sort > j1

#zcat Gene_pCorr8.out | sort | join j1 - | \
zcat Gene_Corr8.out | sort | join j1 - | \
awk '{print $3, $0}' | sort | join j1 - | \
awk '{print $3, $4, $1, $2, $6}' | sort -n +1 -n +3 | \
awk '$2<$4 {print $0}' | awk -f pick_largest.awk > $1.temp.in

./any_matrix.pl $1.temp.in $1_matrix.png
#./any_matrix_grey.pl $1.temp.in $1_matrix.png

#rm j1 $1.temp.in