species="Oryza_sativa_Japonica_Group_x_Oryza_sativa_Indica_Group"
species_dir="oryza"
awk 'BEGIN{print "Type\tShape\tChr\tStart\tEnd\tcolor"}{print "Top\ttriangle\t"$1"\t"$2"\t"$3"\tff7f00"}' ./${species}.top100.txt > ./data/${species}.top100.mark.txt
awk 'ARGIND==1{a[$1]}ARGIND==2{if($3 in a){print $0}}' ~/project/Genome/${species_dir}/${species_dir}_karyotype.txt ./data/${species}.top100.mark.txt > ./data/${species}.top100.mark2.txt
Rscript plot_mark.r ./data/${species}.top100.mark2.txt ~/project/Genome/${species_dir}/${species_dir}_karyotype.txt ${species}.top100 ./ 240 220
