ms-sequest 0.1.1 → 0.1.2
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.
- data/History +4 -0
- data/Rakefile +4 -4
- data/VERSION +1 -1
- data/lib/ms/sequest/sqt.rb +1 -2
- data/lib/ms/sequest/srf.rb +10 -7
- data/lib/ms/sequest/srf/pepxml.rb +8 -8
- data/spec/ms/sequest/srf/sqt_spec.rb +1 -0
- metadata +8 -19
data/History
CHANGED
data/Rakefile
CHANGED
@@ -11,13 +11,13 @@ Jeweler::Tasks.new do |gem|
|
|
11
11
|
gem.email = "jtprince@gmail.com"
|
12
12
|
gem.authors = ["John T. Prince"]
|
13
13
|
gem.rubyforge_project = 'mspire'
|
14
|
-
gem.add_runtime_dependency "ms-ident", ">= 0.
|
15
|
-
gem.add_runtime_dependency "ms-core", ">= 0.0.
|
14
|
+
gem.add_runtime_dependency "ms-ident", ">= 0.1.1"
|
15
|
+
gem.add_runtime_dependency "ms-core", ">= 0.0.19"
|
16
16
|
#gem.add_runtime_dependency "ms-msrun", ">= 0.3.4"
|
17
|
-
gem.add_runtime_dependency "trollop", "~> 1.16"
|
17
|
+
gem.add_runtime_dependency "trollop", "~> 1.16.2"
|
18
18
|
gem.add_development_dependency "jeweler", "~> 1.5.2"
|
19
19
|
gem.add_development_dependency "ms-testdata", ">= 0.1.1"
|
20
|
-
gem.add_development_dependency "spec-more", ">= 0"
|
20
|
+
gem.add_development_dependency "spec-more", ">= 0.0.4"
|
21
21
|
end
|
22
22
|
Jeweler::RubygemsDotOrgTasks.new
|
23
23
|
|
data/VERSION
CHANGED
@@ -1 +1 @@
|
|
1
|
-
0.1.
|
1
|
+
0.1.2
|
data/lib/ms/sequest/sqt.rb
CHANGED
data/lib/ms/sequest/srf.rb
CHANGED
@@ -14,8 +14,7 @@ require 'ms/sequest/params'
|
|
14
14
|
module Ms ; end
|
15
15
|
module Ms::Sequest ; end
|
16
16
|
|
17
|
-
class Ms::Sequest::Srf
|
18
|
-
include Ms::Ident::Search
|
17
|
+
class Ms::Sequest::Srf < Ms::Ident::Search
|
19
18
|
|
20
19
|
class NoSequestParamsError < ArgumentError
|
21
20
|
end
|
@@ -102,7 +101,7 @@ class Ms::Sequest::Srf
|
|
102
101
|
# # will attempt to read peptide hit information (equivalent to .out
|
103
102
|
# # files), otherwise, just reads the dta information.
|
104
103
|
def initialize(filename=nil, opts={})
|
105
|
-
@
|
104
|
+
@peptide_hits = []
|
106
105
|
|
107
106
|
@dta_files = []
|
108
107
|
@out_files = []
|
@@ -270,7 +269,7 @@ class Ms::Sequest::Srf
|
|
270
269
|
outfile.charge = ind[2]
|
271
270
|
|
272
271
|
pep_hits = @out_files[i].hits
|
273
|
-
@
|
272
|
+
@peptide_hits.push( *pep_hits )
|
274
273
|
pep_hits.each do |pep_hit|
|
275
274
|
pep_hit[15] = @base_name
|
276
275
|
pep_hit[16] = ind[0]
|
@@ -286,7 +285,7 @@ class Ms::Sequest::Srf
|
|
286
285
|
filter_by_precursor_mass_tolerance! if params
|
287
286
|
|
288
287
|
if opts[:link_protein_hits]
|
289
|
-
(@
|
288
|
+
(@peptide_hits, @proteins) = merge!([self.peptide_hits]) do |_protein, _peptides|
|
290
289
|
Ms::Sequest::Srf::Out::Protein.new(_protein.reference, _peptides)
|
291
290
|
end
|
292
291
|
end
|
@@ -615,7 +614,7 @@ end
|
|
615
614
|
# deltacn should also be changed to reflect it.
|
616
615
|
# mh - the theoretical mass + h
|
617
616
|
# proteins are created as SRF prot objects with a reference and linked to their
|
618
|
-
#
|
617
|
+
# peptide_hits (from global hash by reference)
|
619
618
|
# ppm = 10^6 * ∆m_accuracy / mass_measured [ where ∆m_accuracy = mass_real – mass_measured ]
|
620
619
|
# This is calculated for the M+H mass!
|
621
620
|
# num_other_loci is the number of other loci that the peptide matches beyond
|
@@ -728,7 +727,7 @@ end
|
|
728
727
|
Ms::Sequest::Srf::Out::Protein = Struct.new( *%w(reference peptides).map(&:to_sym) )
|
729
728
|
|
730
729
|
class Ms::Sequest::Srf::Out::Protein
|
731
|
-
include Ms::Ident::Protein
|
730
|
+
#include Ms::Ident::Protein
|
732
731
|
|
733
732
|
## we shouldn't have to do this because this is inlcuded in SpecID::Protein, but
|
734
733
|
## under some circumstances it won't work without explicitly calling it.
|
@@ -741,6 +740,10 @@ class Ms::Sequest::Srf::Out::Protein
|
|
741
740
|
end
|
742
741
|
$VERBOSE = tmp
|
743
742
|
|
743
|
+
def first_entry
|
744
|
+
reference.split(/[\s\r]/)[0]
|
745
|
+
end
|
746
|
+
|
744
747
|
# "<Ms::Sequest::Srf::Out::Protein reference=\"#{@reference}\">"
|
745
748
|
|
746
749
|
undef_method :inspect
|
@@ -155,24 +155,24 @@ class Ms::Sequest::Srf
|
|
155
155
|
msms_pipeline_analysis.merge!(:summary_xml => summary_xml_filename, :pepxml_version => opt[:pepxml_version]) do |msms_run_summary|
|
156
156
|
# prep the sample enzyme and search_summary
|
157
157
|
msms_run_summary.merge!(
|
158
|
-
:base_name => File.join(opt[:mz_dir], srf.base_name_noext),
|
158
|
+
{ :base_name => File.join(opt[:mz_dir], srf.base_name_noext),
|
159
159
|
:ms_manufacturer => opt[:ms_manufacturer],
|
160
160
|
:ms_model => opt[:ms_model],
|
161
161
|
:ms_ionization => opt[:ms_ionization],
|
162
162
|
:ms_mass_analyzer => opt[:ms_mass_analyzer],
|
163
163
|
:ms_detector => opt[:ms_detector],
|
164
164
|
:raw_data => opt[:raw_data].first,
|
165
|
-
:raw_data_type => opt[:raw_data].first,
|
165
|
+
:raw_data_type => opt[:raw_data].first, }
|
166
166
|
) do |sample_enzyme, search_summary, spectrum_queries|
|
167
167
|
sample_enzyme.merge!(params.sample_enzyme_hash)
|
168
168
|
sample_enzyme.name = opt[:enzyme] if opt[:enzyme]
|
169
169
|
search_summary.merge!(
|
170
|
-
:base_name=> srf.resident_dir + '/' + srf.base_name_noext,
|
170
|
+
{ :base_name=> srf.resident_dir + '/' + srf.base_name_noext,
|
171
171
|
:search_engine => 'SEQUEST',
|
172
172
|
:precursor_mass_type => params.precursor_mass_type,
|
173
173
|
:fragment_mass_type => params.fragment_mass_type,
|
174
174
|
:out_data_type => opt[:out_data_type],
|
175
|
-
:out_data => opt[:out_data],
|
175
|
+
:out_data => opt[:out_data], }
|
176
176
|
) do |search_database, enzymatic_search_constraint, modifications_ar, parameters_hash|
|
177
177
|
search_database.merge!(:local_path => db_filename, :seq_type => opt[:db_seq_type], :database_name => opt[:db_name], :orig_database_url => opt[:db_orig_url], :database_release_date => opt[:db_release_date], :database_release_identifier => opt[:db_release_id])
|
178
178
|
|
@@ -184,9 +184,9 @@ class Ms::Sequest::Srf
|
|
184
184
|
end
|
185
185
|
|
186
186
|
enzymatic_search_constraint.merge!(
|
187
|
-
:enzyme => opt[:enzyme] ? opt[:enzyme] : params.enzyme,
|
187
|
+
{ :enzyme => opt[:enzyme] ? opt[:enzyme] : params.enzyme,
|
188
188
|
:max_num_internal_cleavages => params.max_num_internal_cleavages,
|
189
|
-
:min_number_termini => params.min_number_termini,
|
189
|
+
:min_number_termini => params.min_number_termini, }
|
190
190
|
)
|
191
191
|
modifications_ar.replace(modifications_obj.modifications)
|
192
192
|
parameters_hash.merge!(params.opts)
|
@@ -240,10 +240,10 @@ class Ms::Sequest::Srf
|
|
240
240
|
end
|
241
241
|
end
|
242
242
|
Ms::Ident::Pepxml::SpectrumQuery.new(
|
243
|
-
:spectrum => [srf.base_name_noext, *i_ar].join('.'), :start_scan => i_ar[0], :end_scan => i_ar[1],
|
243
|
+
{ :spectrum => [srf.base_name_noext, *i_ar].join('.'), :start_scan => i_ar[0], :end_scan => i_ar[1],
|
244
244
|
:precursor_neutral_mass => dta_file.mh - h_plus, :assumed_charge => i_ar[2],
|
245
245
|
:retention_time_sec => ret_time,
|
246
|
-
:search_results => [sr],
|
246
|
+
:search_results => [sr], }
|
247
247
|
)
|
248
248
|
end
|
249
249
|
spectrum_queries.replace(spec_queries)
|
metadata
CHANGED
@@ -2,7 +2,7 @@
|
|
2
2
|
name: ms-sequest
|
3
3
|
version: !ruby/object:Gem::Version
|
4
4
|
prerelease:
|
5
|
-
version: 0.1.
|
5
|
+
version: 0.1.2
|
6
6
|
platform: ruby
|
7
7
|
authors:
|
8
8
|
- John T. Prince
|
@@ -10,7 +10,7 @@ autorequire:
|
|
10
10
|
bindir: bin
|
11
11
|
cert_chain: []
|
12
12
|
|
13
|
-
date: 2011-
|
13
|
+
date: 2011-07-05 00:00:00 -06:00
|
14
14
|
default_executable:
|
15
15
|
dependencies:
|
16
16
|
- !ruby/object:Gem::Dependency
|
@@ -21,7 +21,7 @@ dependencies:
|
|
21
21
|
requirements:
|
22
22
|
- - ">="
|
23
23
|
- !ruby/object:Gem::Version
|
24
|
-
version: 0.
|
24
|
+
version: 0.1.1
|
25
25
|
type: :runtime
|
26
26
|
version_requirements: *id001
|
27
27
|
- !ruby/object:Gem::Dependency
|
@@ -32,7 +32,7 @@ dependencies:
|
|
32
32
|
requirements:
|
33
33
|
- - ">="
|
34
34
|
- !ruby/object:Gem::Version
|
35
|
-
version: 0.0.
|
35
|
+
version: 0.0.19
|
36
36
|
type: :runtime
|
37
37
|
version_requirements: *id002
|
38
38
|
- !ruby/object:Gem::Dependency
|
@@ -43,7 +43,7 @@ dependencies:
|
|
43
43
|
requirements:
|
44
44
|
- - ~>
|
45
45
|
- !ruby/object:Gem::Version
|
46
|
-
version:
|
46
|
+
version: 1.16.2
|
47
47
|
type: :runtime
|
48
48
|
version_requirements: *id003
|
49
49
|
- !ruby/object:Gem::Dependency
|
@@ -76,7 +76,7 @@ dependencies:
|
|
76
76
|
requirements:
|
77
77
|
- - ">="
|
78
78
|
- !ruby/object:Gem::Version
|
79
|
-
version:
|
79
|
+
version: 0.0.4
|
80
80
|
type: :development
|
81
81
|
version_requirements: *id006
|
82
82
|
description: reads .SRF, .SQT and supports conversions
|
@@ -161,16 +161,5 @@ rubygems_version: 1.6.2
|
|
161
161
|
signing_key:
|
162
162
|
specification_version: 3
|
163
163
|
summary: An mspire library supporting SEQUEST, Bioworks, SQT, etc
|
164
|
-
test_files:
|
165
|
-
|
166
|
-
- spec/ms/sequest/params_spec.rb
|
167
|
-
- spec/ms/sequest/pepxml/modifications_spec.rb
|
168
|
-
- spec/ms/sequest/pepxml_spec.rb
|
169
|
-
- spec/ms/sequest/sqt_spec.rb
|
170
|
-
- spec/ms/sequest/sqt_spec_helper.rb
|
171
|
-
- spec/ms/sequest/srf/pepxml_spec.rb
|
172
|
-
- spec/ms/sequest/srf/search_spec.rb
|
173
|
-
- spec/ms/sequest/srf/sqt_spec.rb
|
174
|
-
- spec/ms/sequest/srf_spec.rb
|
175
|
-
- spec/ms/sequest/srf_spec_helper.rb
|
176
|
-
- spec/spec_helper.rb
|
164
|
+
test_files: []
|
165
|
+
|