awk '{print $1}' Gene_Corr8_Symb.out | sort -u | awk '{print $1, NR}' > j1 sort Gene_Corr8_Symb.out | 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 > big_matrix.in ./big_matrix.pl big_matrix.in big_matrix.png zcat ~areverte/BLAST/another_annotation.txt.gz | \ awk '{print $2, $3}' | sort -u | join -a1 j1 - > Genes_616.lst rm j1 big_matrix.in