Hi all,
I uploaded a script that I used for my own research just in case anyone else feel it is useful.
This script takes two input file: one is gene_exp.diff generated by cuffdiff; the other is multifasta file (genome file) that your reads were mapped to. And the output is differential expression gene sequences (in fasta format). It takes advantage of the gene positions information of gene_exp.diff to extract the sequence in genome file.
usage: perl extractGene.pl [-d gene_exp.diff] [-g genome.fasta] [--help]
I uploaded a script that I used for my own research just in case anyone else feel it is useful.
This script takes two input file: one is gene_exp.diff generated by cuffdiff; the other is multifasta file (genome file) that your reads were mapped to. And the output is differential expression gene sequences (in fasta format). It takes advantage of the gene positions information of gene_exp.diff to extract the sequence in genome file.
usage: perl extractGene.pl [-d gene_exp.diff] [-g genome.fasta] [--help]
Comment