HenrikBengtsson / aroma.seq

🔬 R package: aroma.seq: High-Throughput Sequence Analysis using the Aroma Framework
https://github.com/HenrikBengtsson/aroma.seq
0 stars 1 forks source link

Analogue to *.gz files, add basic support for *.bgz files (bgzip) #44

Open HenrikBengtsson opened 7 years ago

HenrikBengtsson commented 7 years ago

The samtools software does not support gzip'ed files, because it can't index gzip-ed files. For instance, if we try to use a gzip'ed FASTA file, we get:

[fai_load] build FASTA index.
Cannot index files compressed with gzip, please use bgzip

Because it supports bgzip, the aroma.seq pathname parser should be updated to recognize *.bgz files too. Until then, we can only work with uncompressed FASTA files.

Related to https://github.com/HenrikBengtsson/Costello-PSCN-Seq/issues/5#issuecomment-310870978