millironx/XAM.jl
1
0
Fork 0
mirror of https://github.com/MillironX/XAM.jl.git synced 2024-11-23 02:09:55 +00:00
Code Issues Releases Packages Activity

Update BAM.sequence to return BioSequences.LongDNA{4}

Browse source 
Same principles apply as updating SAM.sequence
This commit is contained in:
Thomas A. Christensen II 2022-07-07 09:19:43 -05:00
parent 1660a24a0f
commit a7a76db3d4
Signed by: millironx
GPG key ID: 139C07724802BC5D
1 changed files with 3 additions and 3 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

6
src/bam/record.jl
View file

@ -91,7 +91,7 @@ function Base.empty!(record::Record)
record.tlen = 0
#Note: data will be overwritten and indexed using data_size.
return record
end
@ -505,7 +505,7 @@ function hastemplength(record::Record)
end
"""
sequence(record::Record)::BioSequences.LongDNASeq
sequence(record::Record)::BioSequences.LongDNA{4}
Get the segment sequence of `record`.
"""
@ -522,7 +522,7 @@ function sequence(record::Record)
x = unsafe_load(src, i)
data[i] = (x & 0x0f0f0f0f0f0f0f0f) << 4 | (x & 0xf0f0f0f0f0f0f0f0) >> 4
end
return BioSequences.LongDNASeq(data, 1:seqlen, false)
return BioSequences.LongDNA{4}(data, UInt(seqlen))
end
function hassequence(record::Record)