Fetch terminal genes C. elegans
end_dist=30000
chrom_I=15072434
chrom_II=15279421
chrom_III=13783801
chrom_IV=17493829
chrom_V=20924180
chrom_X=17718942
tabix gene.gff.gz I:1-$end_dist > terminal_genes.txt
tabix gene.gff.gz II:1-$end_dist >> terminal_genes.txt
tabix gene.gff.gz III:1-$end_dist >> terminal_genes.txt
tabix gene.gff.gz IV:1-$end_dist >> terminal_genes.txt
tabix gene.gff.gz V:1-$end_dist >> terminal_genes.txt
tabix gene.gff.gz X:1-$end_dist >> terminal_genes.txt
tabix gene.gff.gz I:`expr $chrom_I - $end_dist`-$chrom_I >> terminal_genes.txt
tabix gene.gff.gz I:`expr $chrom_II - $end_dist`-$chrom_II >> terminal_genes.txt
tabix gene.gff.gz I:`expr $chrom_III - $end_dist`-$chrom_III >> terminal_genes.txt
tabix gene.gff.gz I:`expr $chrom_IV - $end_dist`-$chrom_IV >> terminal_genes.txt
tabix gene.gff.gz I:`expr $chrom_V - $end_dist`-$chrom_V >> terminal_genes.txt
tabix gene.gff.gz I:`expr $chrom_X - $end_dist`-$chrom_X >> terminal_genes.txt
cat terminal_genes.txt | grep 'Gene' | cut -f 9 | egrep -o 'Alias=([^,]+)'