millironx/XAM.jl
1
0
Fork 0
mirror of https://github.com/MillironX/XAM.jl.git synced 2024-11-14 22:33:14 +00:00
Code Issues Releases Packages Activity

Improve BAM.quality performance

Browse source 
Closes https://github.com/BioJulia/BioAlignments.jl/pull/33/files.
This commit is contained in:
Ciarán O'Mara 2020-07-01 17:48:39 +10:00 committed by Ciarán O'Mara
parent dce44b3f55
commit cf168acaa3
1 changed files with 4 additions and 4 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

6
src/bam/record.jl
View file

@ -543,15 +543,15 @@ function hasseqlength(record::Record)
end
"""
quality(record::Record)::Vector{UInt8}
quality(record::Record)
Get the base quality of `record`.
"""
function quality(record::Record)::Vector{UInt8}
function quality(record::Record)
checkfilled(record)
seqlen = seqlength(record)
offset = seqname_length(record) + n_cigar_op(record, false) * 4 + cld(seqlen, 2)
return [reinterpret(Int8, record.data[i+offset]) for i in 1:seqlen]
return record.data[(1+offset):(seqlen+offset)]
end
function hasquality(record::Record)