htslib 0.1.0 → 0.2.2
Sign up to get free protection for your applications and to get access to all the features.
- checksums.yaml +4 -4
- data/README.md +30 -11
- data/lib/hts/bam/aux.rb +1 -0
- data/lib/hts/bam/cigar.rb +5 -7
- data/lib/hts/bam/flag.rb +5 -0
- data/lib/hts/bam/header.rb +1 -0
- data/lib/hts/bam/record.rb +7 -4
- data/lib/hts/bam.rb +128 -35
- data/lib/hts/bcf/format.rb +4 -4
- data/lib/hts/bcf/header.rb +1 -0
- data/lib/hts/bcf/header_record.rb +11 -0
- data/lib/hts/bcf/info.rb +3 -2
- data/lib/hts/bcf/record.rb +31 -30
- data/lib/hts/bcf.rb +183 -30
- data/lib/hts/faidx.rb +62 -21
- data/lib/hts/hts.rb +54 -2
- data/lib/hts/libhts/bgzf.rb +5 -5
- data/lib/hts/libhts/constants.rb +42 -20
- data/lib/hts/libhts/cram.rb +312 -279
- data/lib/hts/libhts/hfile.rb +29 -11
- data/lib/hts/libhts/hts.rb +158 -25
- data/lib/hts/libhts/sam.rb +695 -94
- data/lib/hts/libhts/sam_funcs.rb +131 -588
- data/lib/hts/libhts/tbx.rb +1 -1
- data/lib/hts/libhts/tbx_funcs.rb +25 -0
- data/lib/hts/libhts/vcf.rb +460 -234
- data/lib/hts/libhts/vcf_funcs.rb +277 -424
- data/lib/hts/libhts.rb +1 -0
- data/lib/hts/tbx.rb +39 -8
- data/lib/hts/version.rb +1 -1
- metadata +5 -3
data/lib/hts/libhts.rb
CHANGED
data/lib/hts/tbx.rb
CHANGED
@@ -8,7 +8,7 @@ module HTS
|
|
8
8
|
class Tbx < Hts
|
9
9
|
include Enumerable
|
10
10
|
|
11
|
-
attr_reader :file_name
|
11
|
+
attr_reader :file_name, :index_name, :mode, :nthreads
|
12
12
|
|
13
13
|
def self.open(*args, **kw)
|
14
14
|
file = new(*args, **kw) # do not yield
|
@@ -22,24 +22,55 @@ module HTS
|
|
22
22
|
file
|
23
23
|
end
|
24
24
|
|
25
|
-
def initialize(file_name, threads: nil)
|
25
|
+
def initialize(file_name, mode = "r", index: nil, threads: nil, build_index: false)
|
26
26
|
if block_given?
|
27
27
|
message = "HTS::Tbx.new() dose not take block; Please use HTS::Tbx.open() instead"
|
28
28
|
raise message
|
29
29
|
end
|
30
30
|
|
31
|
-
@file_name = file_name
|
32
|
-
|
33
31
|
# NOTE: Do not check for the existence of local files, since file_names may be remote URIs.
|
34
32
|
|
35
|
-
@
|
33
|
+
@file_name = file_name
|
34
|
+
@index_name = index
|
35
|
+
@mode = mode
|
36
|
+
@nthreads = threads
|
36
37
|
@hts_file = LibHTS.hts_open(@file_name, @mode)
|
37
38
|
|
38
39
|
raise Errno::ENOENT, "Failed to open #{@file_name}" if @hts_file.null?
|
39
40
|
|
40
|
-
if threads
|
41
|
-
|
42
|
-
|
41
|
+
set_threads(threads) if threads
|
42
|
+
|
43
|
+
# return if @mode[0] == "w"
|
44
|
+
raise "Not implemented" if @mode[0] == "w"
|
45
|
+
|
46
|
+
# build_index(index) if build_index
|
47
|
+
@idx = load_index(index)
|
48
|
+
|
49
|
+
super # do nothing
|
50
|
+
end
|
51
|
+
|
52
|
+
def build_index
|
53
|
+
raise "Not implemented yet"
|
54
|
+
end
|
55
|
+
|
56
|
+
def load_index(index_name = nil)
|
57
|
+
if index_name
|
58
|
+
LibHTS.tbx_index_load2(@file_name, index_name)
|
59
|
+
else
|
60
|
+
LibHTS.tbx_index_load3(@file_name, nil, 2)
|
61
|
+
end
|
62
|
+
end
|
63
|
+
|
64
|
+
def tid(name)
|
65
|
+
LibHTS.tbx_name2id(@idx, name)
|
66
|
+
end
|
67
|
+
|
68
|
+
def seqnames
|
69
|
+
nseq = FFI::MemoryPointer.new(:int)
|
70
|
+
LibHTS.tbx_seqnames(@idx, nseq).then do |pts|
|
71
|
+
pts.read_array_of_pointer(nseq.read_int).map do |pt|
|
72
|
+
pt.read_string
|
73
|
+
end
|
43
74
|
end
|
44
75
|
end
|
45
76
|
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.
|
4
|
+
version: 0.2.2
|
5
5
|
platform: ruby
|
6
6
|
authors:
|
7
7
|
- kojix2
|
8
8
|
autorequire:
|
9
9
|
bindir: bin
|
10
10
|
cert_chain: []
|
11
|
-
date: 2022-
|
11
|
+
date: 2022-08-19 00:00:00.000000000 Z
|
12
12
|
dependencies:
|
13
13
|
- !ruby/object:Gem::Dependency
|
14
14
|
name: ffi
|
@@ -140,6 +140,7 @@ files:
|
|
140
140
|
- lib/hts/bcf.rb
|
141
141
|
- lib/hts/bcf/format.rb
|
142
142
|
- lib/hts/bcf/header.rb
|
143
|
+
- lib/hts/bcf/header_record.rb
|
143
144
|
- lib/hts/bcf/info.rb
|
144
145
|
- lib/hts/bcf/record.rb
|
145
146
|
- lib/hts/faidx.rb
|
@@ -158,6 +159,7 @@ files:
|
|
158
159
|
- lib/hts/libhts/sam.rb
|
159
160
|
- lib/hts/libhts/sam_funcs.rb
|
160
161
|
- lib/hts/libhts/tbx.rb
|
162
|
+
- lib/hts/libhts/tbx_funcs.rb
|
161
163
|
- lib/hts/libhts/thread_pool.rb
|
162
164
|
- lib/hts/libhts/vcf.rb
|
163
165
|
- lib/hts/libhts/vcf_funcs.rb
|
@@ -176,7 +178,7 @@ required_ruby_version: !ruby/object:Gem::Requirement
|
|
176
178
|
requirements:
|
177
179
|
- - ">="
|
178
180
|
- !ruby/object:Gem::Version
|
179
|
-
version: '
|
181
|
+
version: '3.1'
|
180
182
|
required_rubygems_version: !ruby/object:Gem::Requirement
|
181
183
|
requirements:
|
182
184
|
- - ">="
|