Extract fasta sequences from CDS #fasta
exon_file="N2_GFF_WS245/CDS.gff.gz"
reference="reference/c_elegans.PRJNA13758.WS245.genomic.fa"
paste <(gunzip -kfc ${exon_file}) <(bedtools getfasta -fi ${reference} -bed ${exon_file} -tab -fo -) | grep 'WormBase'