Bitwise parallel implementations of gc_content and hasambiguity

Parallel_Seqops.jl

import BioSequences: index, offset, bitindex

function parallel_gc_content( seq::BioSequence )
    gc   = 0
    s = bitindex(seq, 1)
    e = bitindex(seq, length(seq))
    for i in index(s):index(e)
        @inbounds data = seq.data[i]
        # mask nucleotides upstream of first(seq.part)
        if i == index(s)
            data &= 0xFFFFFFFFFFFFFFFF << offset(s)
        end
        # mask nucleotides downstream of last(seq.part)
        if i == index(e)
            data &= 0xFFFFFFFFFFFFFFFF >> (64 - (offset(e)+4))
        end
        # count unambiguous GC nucleotides
        gc += count_ones(data & 0x6666666666666666)
    end
    return length(seq) > 0 ? gc / length(seq) : 0.0
end

function parallel_hasambiguity( seq::BioSequence )
    nucs = 0
    s = bitindex(seq, 1)
    e = bitindex(seq, length(seq))
    for i in index(s):index(e)
        @inbounds data = seq.data[i]
        # mask nucleotides upstream of first(seq.part)
        if i == index(s)
            data &= 0xFFFFFFFFFFFFFFFF << offset(s)
        end
        # mask nucleotides downstream of last(seq.part)
        if i == index(e)
            data &= 0xFFFFFFFFFFFFFFFF >> (64 - (offset(e)+4))
        end
        # count
        block_ones = count_ones(data)
        if block_ones > 16
            return true
        end
        nucs += block_ones
    end
    return nucs > length(seq) ? true : false
end