htslib 0.2.5 → 0.2.6
This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
- checksums.yaml +4 -4
- data/README.md +41 -43
- data/TUTORIAL.md +270 -0
- data/lib/hts/bam/auxi.rb +28 -2
- data/lib/hts/bam/cigar.rb +46 -6
- data/lib/hts/bam/flag.rb +43 -4
- data/lib/hts/bam/record.rb +2 -1
- data/lib/hts/bam.rb +21 -13
- data/lib/hts/bcf/format.rb +28 -24
- data/lib/hts/bcf/info.rb +28 -24
- data/lib/hts/bcf.rb +34 -23
- data/lib/hts/faidx/sequence.rb +64 -0
- data/lib/hts/faidx.rb +64 -15
- data/lib/hts/hts.rb +10 -7
- data/lib/hts/libhts/constants.rb +14 -14
- data/lib/hts/libhts/cram.rb +1 -1
- data/lib/hts/libhts/sam.rb +4 -4
- data/lib/hts/libhts/tbx.rb +2 -0
- data/lib/hts/libhts/tbx_funcs.rb +1 -1
- data/lib/hts/libhts/vcf.rb +1 -1
- data/lib/hts/tbx.rb +49 -0
- data/lib/hts/version.rb +1 -1
- metadata +5 -3
data/lib/hts/tbx.rb
CHANGED
@@ -59,6 +59,10 @@ module HTS
|
|
59
59
|
end
|
60
60
|
end
|
61
61
|
|
62
|
+
def index_loaded?
|
63
|
+
!@idx.null?
|
64
|
+
end
|
65
|
+
|
62
66
|
def tid(name)
|
63
67
|
LibHTS.tbx_name2id(@idx, name)
|
64
68
|
end
|
@@ -69,5 +73,50 @@ module HTS
|
|
69
73
|
pts.read_array_of_pointer(nseq.read_int).map(&:read_string)
|
70
74
|
end
|
71
75
|
end
|
76
|
+
|
77
|
+
def query(region, start = nil, end_ = nil, &block)
|
78
|
+
check_closed
|
79
|
+
raise "Index file is required to call the query method." unless index_loaded?
|
80
|
+
|
81
|
+
if start && end_
|
82
|
+
queryi(tid(region), start, end_, &block)
|
83
|
+
else
|
84
|
+
querys(region, &block)
|
85
|
+
end
|
86
|
+
end
|
87
|
+
|
88
|
+
private
|
89
|
+
|
90
|
+
def queryi(id, start, end_)
|
91
|
+
return to_enum(__method__, id, start, end_) unless block_given?
|
92
|
+
|
93
|
+
qiter = LibHTS.tbx_itr_queryi(@idx, id, start, end_)
|
94
|
+
|
95
|
+
raise "Failed to query region: #{id}:#{start}-#{end_}" if qiter.null?
|
96
|
+
|
97
|
+
r = LibHTS::KString.new
|
98
|
+
begin
|
99
|
+
yield r[:s] while LibHTS.tbx_itr_next(@hts_file, @idx, qiter, r) > 0
|
100
|
+
ensure
|
101
|
+
LibHTS.hts_itr_destroy(qiter)
|
102
|
+
end
|
103
|
+
self
|
104
|
+
end
|
105
|
+
|
106
|
+
def querys(region)
|
107
|
+
return to_enum(__method__, region) unless block_given?
|
108
|
+
|
109
|
+
qiter = LibHTS.tbx_itr_querys(@idx, region)
|
110
|
+
|
111
|
+
raise "Failed to query region: #{region}" if qiter.null?
|
112
|
+
|
113
|
+
r = LibHTS::KString.new
|
114
|
+
begin
|
115
|
+
yield r[:s].split("\t") while LibHTS.tbx_itr_next(@hts_file, @idx, qiter, r) > 0
|
116
|
+
ensure
|
117
|
+
LibHTS.hts_itr_destroy(qiter)
|
118
|
+
end
|
119
|
+
self
|
120
|
+
end
|
72
121
|
end
|
73
122
|
end
|
data/lib/hts/version.rb
CHANGED
metadata
CHANGED
@@ -1,14 +1,14 @@
|
|
1
1
|
--- !ruby/object:Gem::Specification
|
2
2
|
name: htslib
|
3
3
|
version: !ruby/object:Gem::Version
|
4
|
-
version: 0.2.
|
4
|
+
version: 0.2.6
|
5
5
|
platform: ruby
|
6
6
|
authors:
|
7
7
|
- kojix2
|
8
8
|
autorequire:
|
9
9
|
bindir: bin
|
10
10
|
cert_chain: []
|
11
|
-
date:
|
11
|
+
date: 2023-01-09 00:00:00.000000000 Z
|
12
12
|
dependencies:
|
13
13
|
- !ruby/object:Gem::Dependency
|
14
14
|
name: ffi
|
@@ -131,6 +131,7 @@ extra_rdoc_files: []
|
|
131
131
|
files:
|
132
132
|
- LICENSE.txt
|
133
133
|
- README.md
|
134
|
+
- TUTORIAL.md
|
134
135
|
- lib/hts/bam.rb
|
135
136
|
- lib/hts/bam/auxi.rb
|
136
137
|
- lib/hts/bam/cigar.rb
|
@@ -145,6 +146,7 @@ files:
|
|
145
146
|
- lib/hts/bcf/info.rb
|
146
147
|
- lib/hts/bcf/record.rb
|
147
148
|
- lib/hts/faidx.rb
|
149
|
+
- lib/hts/faidx/sequence.rb
|
148
150
|
- lib/hts/ffi_ext/README.md
|
149
151
|
- lib/hts/ffi_ext/pointer.rb
|
150
152
|
- lib/hts/ffi_ext/struct.rb
|
@@ -187,7 +189,7 @@ required_rubygems_version: !ruby/object:Gem::Requirement
|
|
187
189
|
- !ruby/object:Gem::Version
|
188
190
|
version: '0'
|
189
191
|
requirements: []
|
190
|
-
rubygems_version: 3.
|
192
|
+
rubygems_version: 3.4.1
|
191
193
|
signing_key:
|
192
194
|
specification_version: 4
|
193
195
|
summary: HTSlib bindings for Ruby
|