readfasta.jl

using CodecZlib
using BioSequences

function readfasta(filepath)
    open(GzipDecompressionStream, filepath) do stream
        reader = FASTA.Reader(stream)
        for record in reader
            id  = FASTA.identifier(record)
            seq = FASTA.sequence(record)
            @show id, seq
        end
    end
end

readfasta("data.fa.gz")