DaehwanKimLab / tophat

Spliced read mapper for RNA-Seq
http://ccb.jhu.edu/software/tophat
Boost Software License 1.0
90 stars 46 forks source link

Do not double open genome fasta in gtf_to_fasta #46

Closed dondelelcaro closed 6 years ago

dondelelcaro commented 7 years ago

No longer open the genome fasta in the GTFToFasta() and again in make_transcriptome.

Remove genome_fhandle initializer and member as it is not used anywhere else.

[This enables the following to work: gtf_to_fasta <(gzip -dc Homo_sapiens.GRCh38.84.gtf.gz ) <( gzip -dc Homo_sapiens.GRCh38.dna.toplevel.fa.gz) gtf_fasta.fa which fails with the current version of gtf_to_fasta.]