mspire 0.5.0 → 0.6.1

data/spec/ms/ident/protein_group_spec.rb

@@ -0,0 +1,68 @@
+require 'spec_helper'
+
+require 'ms/ident/protein_group'
+
+PeptideHit = Struct.new(:aaseq, :charge, :proteins) do
+  def inspect # easier to read output
+    "<PeptideHit aaseq=#{self.aaseq} charge=#{self.charge} proteins(ids)=#{self.proteins.map(&:id).join(',')}>"
+  end
+  def hash ; self.object_id end
+end
+ProteinHit = Struct.new(:id) do
+  def inspect # easier to read output
+    "<Prt #{self.id}>"
+  end
+  def hash ; self.object_id end
+end
+
+describe 'creating minimal protein groups from peptide hits' do
+  before do
+    @pep_hits = [ ['AABBCCDD', 2], 
+      ['BBCC', 2],
+      ['DDEEFFGG', 2],
+      ['DDEEFFGG', 3],
+      ['HIYA', 2],
+    ].map {|ar| PeptideHit.new(ar[0], ar[1], []) }
+    @prot_hits_hash = { 
+      'big_guy' => @pep_hits,
+      'little_guy' => [@pep_hits.last],
+      'medium_guy1' => @pep_hits[0,4],
+      'medium_guy2' => @pep_hits[0,4],
+      'subsumed_by_medium' => @pep_hits[2,2],
+    }
+    @prot_hits = @prot_hits_hash.keys.map {|id| ProteinHit.new(id) }
+  end
+
+  it 'is a greedy algorithm' do
+    @prot_hits.each {|prthit| @prot_hits_hash[prthit.id].each {|pep| pep.proteins << prthit } }
+    # big_guy has all the peptides, so it takes them all
+    protein_groups = MS::Ident::ProteinGroup.peptide_hits_to_protein_groups(@pep_hits)
+    protein_groups.first.size.should == 1# the group
+    protein_groups.first.first.id.should == 'big_guy'
+  end
+
+  it 'removes proteins accounted for only as little pieces of larger proteins' do
+    @prot_hits[1..-1].each {|prthit| @prot_hits_hash[prthit.id].each {|pep| pep.proteins << prthit } }
+    protein_groups = MS::Ident::ProteinGroup.peptide_hits_to_protein_groups(@pep_hits)
+    # no subsumed_by_medium
+    protein_groups.any? {|prot_group| prot_group.any? {|v| v.id == 'subsumed_by_medium' }}.should == false
+  end
+
+  it 'allows alternate sorting algorithms for greediness' do
+    @prot_hits.each {|prthit| @prot_hits_hash[prthit.id].each {|pep| pep.proteins << prthit } }
+    prot_groups = MS::Ident::ProteinGroup.peptide_hits_to_protein_groups(@pep_hits) do |prot_and_peptide_hits|
+      # deliberate using a counterintuitive sorting method to give little guys
+      # a chance
+      -prot_and_peptide_hits.last.size
+    end
+    # because the little proteins are given priority, they 'survive'. Bigger
+    # proteins may also survive if they have at least one unique peptide
+    # to add to the mix.  This demonstrates how proteins can be weighted in
+    # different ways based on their peptide hits.
+    seen = []
+    prot_groups.each {|pg| pg.each {|prot| seen << prot.id } }
+    # big guy is completely accounted for in the now prioritized little guy
+    # and medium guys, etc.
+    seen.sort.should == @prot_hits_hash.keys[1..-1].sort
+  end
+end

data/spec/ms/mass_spec.rb

@@ -0,0 +1,8 @@
+require 'spec_helper'
+
+describe 'MS::Mass' do
+  it 'can access elemental masses by string or symbol' do
+    MS::Mass::MONO['h'].should be_a_kind_of Float
+    MS::Mass::MONO[:h].should be_a_kind_of Float
+  end
+end

data/spec/ms/mzml/index_list_spec.rb

@@ -0,0 +1,122 @@
+require 'spec_helper'
+
+require 'ms/mzml'
+require 'ms/mzml/index_list'
+
+describe 'non-indexed uncompressed peaks, mzML file' do
+  subject { TESTFILES + "/ms/mzml/openms.noidx_nocomp.12.mzML" }
+
+  describe 'MS::Mzml getting the IndexList' do
+    before do
+      @io = File.open(subject)
+    end
+    after do
+      @io.close
+    end
+
+    it 'works for an un-indexed file' do
+      mzml = MS::Mzml.new(@io)
+      index_list = mzml.index_list
+      spectrum_idx = index_list[:spectrum]
+      spectrum_idx.name.should == :spectrum
+      spectrum_idx.should == [5755, 6616, 9244, 10077, 12706, 13539, 16168, 17001, 19630, 20463, 23092, 23926]
+      spectrum_idx.ids.should == [
+        "scan=10929", "scan=10930", "scan=10931", 
+        "scan=10932", "scan=10933", "scan=10934", 
+        "scan=10935", "scan=10936", "scan=10937", 
+        "scan=10938", "scan=10939", "scan=10940"
+      ]
+      # right now, the behavior is to expect a chromatogram index, even if
+      # there are no chromatograms... is that the right behavior???
+      chromatogram_idx = index_list[:chromatogram]
+      chromatogram_idx.name.should == :chromatogram
+      chromatogram_idx.should == []
+      chromatogram_idx.ids.should == []
+    end
+  end
+end
+
+
+describe 'indexed, compressed peaks, mzML file' do
+  subject { TESTFILES + "/ms/mzml/j24z.idx_comp.3.mzML"  }
+
+
+  describe 'MS::Mzml getting the IndexList' do
+    before do
+      @io = File.open(subject)
+    end
+    after do
+      @io.close
+    end
+
+    correct_index_list = lambda do |index_list|
+      spectrum_idx = index_list[:spectrum]
+      spectrum_idx.name.should == :spectrum
+      spectrum_idx.should == [4398, 219667, 227895]
+      spectrum_idx.ids.should == [
+        "controllerType=0 controllerNumber=1 scan=1", 
+        "controllerType=0 controllerNumber=1 scan=2", 
+        "controllerType=0 controllerNumber=1 scan=3"
+      ]
+      chromatogram_idx = index_list[:chromatogram]
+      chromatogram_idx.name.should == :chromatogram
+      chromatogram_idx.should == [239184]
+      chromatogram_idx.ids.should == ["TIC"]
+    end
+
+    it 'works for an indexed file' do
+      mzml = MS::Mzml.new(@io)
+      index_list = mzml.index_list
+      correct_index_list.call(index_list)
+    end
+
+    it 'can create the index manually, if requested' do
+      mzml = MS::Mzml.new(@io)
+      index_list =  mzml.create_index_list
+      correct_index_list.call(index_list)
+    end
+
+    describe 'an MS::Mzxml::IndexList' do
+      before do
+        @inner_io = File.open(subject)
+        mzml = MS::Mzml.new(@inner_io)
+        @index_list = mzml.index_list
+      end
+      after do
+        @inner_io.close
+      end
+      it 'can access indices like an array' do
+        @index_list[0].should_not be_nil
+        @index_list[1].should_not be_nil
+        @index_list.map(&:name).should == [:spectrum, :chromatogram]
+      end
+      it 'can access indices like a hash' do
+        @index_list[:spectrum].should_not be_nil
+        @index_list[:chromatogram].should_not be_nil
+      end
+    end
+
+    describe 'an MS::Mzxml::Index' do
+      before do
+        @inner_io = File.open(subject)
+        mzml = MS::Mzml.new(@inner_io)
+        @spec_index = mzml.index_list[:spectrum]
+      end
+      after do
+        @inner_io.close
+      end
+      it 'can create a scan to index hash' do
+        @spec_index.create_scan_to_index.should == {1=>0, 2=>1, 3=>2}
+      end
+      it 'returns nil if there are no scans defined in the id (scan=XX)' do
+        @spec_index.ids.map! { "no scan here" }
+        @spec_index.create_scan_to_index.should be_nil
+      end
+      it 'returns false if there are non-unique scan numbers' do
+        @spec_index.ids.map! { "scan=1" }
+        @spec_index.create_scan_to_index.should == false
+      end
+    end
+  end
+end
+

data/spec/ms/mzml/plms1_spec.rb

@@ -0,0 +1,62 @@
+require 'spec_helper'
+
+require 'ms/mzml'
+require 'ms/mzml/plms1'
+
+describe 'converting mzml to plms1' do
+  describe 'converting no spectra mzml' do
+    before do 
+      @mzml_file = File.open( TESTFILES + "/ms/mzml/openms.noidx_nocomp.12.mzML" )
+      @mzml = MS::Mzml.new(@mzml_file)
+    end
+
+    after do
+      @mzml_file.close
+    end
+
+    it 'can be converted into a plms1 object' do
+      scan_nums = (10929..10940).to_a
+      times =[6604.58, 6605.5, 6605.91, 6606.48, 6606.98, 6607.53, 6607.93, 6608.49, 6608.92, 6609.49, 6609.94, 6610.53]
+      plms1 = @mzml.to_plms1
+      plms1.spectra.respond_to?(:each).should be_true
+      plms1.times.should == times
+      plms1.scan_numbers.should == scan_nums
+      plms1.spectra.each do |spec|
+        p spec.size
+        p spec.class
+        p spec.mzs
+        p spec.intensities
+      end
+      plms1.write("tmp.tmp.bin")
+    end
+
+  end
+
+  describe 'converting normal mzml' do
+    before do 
+      @mzml_file = File.open( TESTFILES + "/ms/mzml/j24z.idx_comp.3.mzML" )
+      @mzml = MS::Mzml.new(@mzml_file)
+    end
+
+    after do
+      @mzml_file.close
+    end
+
+    it 'can be converted into a plms1 object' do
+      #scan_nums = ###
+      #times =[6604.58, 6605.5, 6605.91, 6606.48, 6606.98, 6607.53, 6607.93, 6608.49, 6608.92, 6609.49, 6609.94, 6610.53]
+      plms1 = @mzml.to_plms1
+      plms1.spectra.respond_to?(:each).should be_true
+      p plms1.times
+      p plms1.scan_numbers
+      plms1.spectra.each do |spec|
+        p spec.size
+        p spec.class
+        p spec.mzs.size
+        p spec.intensities.size
+      end
+      #plms1.write("tmp.tmp.bin")
+    end
+
+  end
+end

data/spec/ms/mzml_spec.rb

@@ -0,0 +1,50 @@
+require 'spec_helper'
+
+require 'ms/mzml'
+
+describe 'indexed, compressed peaks, mzML file' do
+
+  describe MS::Mzml do
+
+    describe 'reading a spectrum' do
+
+      before do
+        @file = TESTFILES + "/ms/mzml/j24z.idx_comp.3.mzML"
+        @io = File.open(@file)
+        @mzml = MS::Mzml.new(@io)
+      end
+      after do
+        @io.close
+      end
+
+      it '#spectrum (or #[]) returns a spectrum when queried by index (Integer)' do
+        spectrum = @mzml.spectrum(1) # getting the second spectrum
+        spectrum1 = @mzml[1]
+        spectrum.should == spectrum1
+        spectrum.should be_a(MS::Spectrum)
+        spectrum.should respond_to(:mzs)
+        spectrum.should respond_to(:intensities)
+        spectrum.mzs.size.should == 315
+        spectrum.intensities.size.should == 315
+        spectrum.mzs[2].should be_within(1e-10).of(164.32693481445312)
+      end
+
+      it '#spectrum (or #[]) returns a spectrum when queried by id (String)' do
+        spectrum = @mzml.spectrum("controllerType=0 controllerNumber=1 scan=2")
+        spectrum1 = @mzml["controllerType=0 controllerNumber=1 scan=2"]
+        spectrum.should == spectrum1
+        spectrum.should be_a(MS::Spectrum)
+      end
+
+      it 'goes through spectrum with #each or #each_spectrum' do
+        mz_sizes = [20168, 315, 634]
+        @mzml.each do |spec|
+          spec.mzs.size.should == mz_sizes.shift
+        end
+      end
+
+    end
+  end
+end
+
+

data/spec/ms/plms1_spec.rb

@@ -0,0 +1,38 @@
+require 'spec_helper'
+
+require 'ms/plms1'
+require 'ms/spectrum'
+
+describe 'plms1 - Prince Lab MS 1 specification' do
+  before do
+    @keyfile = TESTFILES + "/plms1/output.key" 
+    times = [0.55, 0.9]
+    scan_numbers = [1,2]
+    spectra = [
+      MS::Spectrum.new([[300.0, 301.5, 303.1], [10, 20, 35.5]]),
+      MS::Spectrum.new([[300.5, 302, 303.6], [11, 21, 36.5]])
+    ]
+    @plms1_obj = MS::Plms1.new(scan_numbers, times, spectra)
+    @outfile = @keyfile.sub(/\.key$/, ".tmp")
+  end
+
+  it 'has a detailed specification' do
+    specification = MS::Plms1::SPECIFICATION
+    specification.should be_an_instance_of String
+    (specification.size > 50).should == true
+  end
+
+  it 'writes a plms1 file' do
+    @plms1_obj.write(@outfile)
+    File.exist?(@outfile).should == true
+    IO.read(@outfile, :mode => 'rb').should == IO.read(@keyfile, :mode => 'rb')
+    File.unlink(@outfile) if File.exist?(@outfile)
+  end
+
+  it 'reads a plms1 file' do
+    obj = MS::Plms1.new.read(@keyfile)
+    [:scan_numbers, :times, :spectra].each do |val|
+      obj.send(val).should == @plms1_obj.send(val)
+    end
+  end
+end

data/spec/ms/quant/qspec_spec.rb

@@ -0,0 +1,25 @@
+require 'spec_helper'
+
+require 'ms/quant/qspec'
+require 'csv'
+
+describe 'running qspec' do 
+  before do
+    @file = TESTFILES + '/ms/quant/max_quant_output.txt'
+    rows = IO.readlines(@file).map {|line| line.chomp.split("\t") }
+    p rows.map(&:size)
+    abort 'here'
+    #columns = IO.readlines(@file).map {|line| line.split("\t") }.transpose
+    headers = columns.map(&:first)
+    p headers
+  end
+
+  describe 'on spectral count data' do
+    it 'works' do
+      1.should == 1
+    end
+  end
+
+  describe 'on MaxQuant LFQ data' do
+  end
+end

data/spec/msplat_spec.rb

@@ -0,0 +1,24 @@
+require 'spec_helper'
+
+require 'msplat'
+
+class String
+  def const
+    self.split('::').inject(Kernel) {|scope, const_name| scope.const_get(const_name)}
+  end
+end
+
+describe "msplat" do
+  it "requires major classes/modules" do
+    %w(
+      MS::Mass
+      MS::Mass::AA
+    ).each do |str|
+      expect { str.const }.to be_true
+    end
+    # for sanity's sake
+    %w(MS::DoesntExit).each do |str|
+      expect { str.const }.to raise_error(NameError)
+    end
+  end
+end

data/spec/obo_spec.rb

@@ -0,0 +1,25 @@
+require 'spec_helper'
+
+require 'obo/ms'
+require 'obo/ims'
+require 'obo/unit'
+
+describe 'accessing a specific Obo::Ontology' do
+  it 'can access MS obo' do
+    Obo::MS.id_to_name['MS:1000004'].should == 'sample mass'
+    Obo::MS.name_to_id['sample mass'].should == 'MS:1000004'
+    Obo::MS.id_to_element['MS:1000004'].should be_a(Obo::Stanza)
+  end
+
+  it 'can access IMS obo' do
+    Obo::IMS.id_to_name['IMS:1000004'].should == 'image'
+    Obo::IMS.name_to_id['image'].should == 'IMS:1000004'
+    Obo::IMS.id_to_element['IMS:1000004'].should be_a(Obo::Stanza)
+  end
+
+  it 'can access Unit obo' do
+    Obo::Unit.id_to_name['UO:0000005'].should == 'temperature unit'
+    Obo::Unit.name_to_id['temperature unit'].should == 'UO:0000005'
+    Obo::Unit.id_to_element['UO:0000005'].should be_a(Obo::Stanza)
+  end
+end

data/spec/spec_helper.rb

@@ -0,0 +1,25 @@
+$LOAD_PATH.unshift(File.join(File.dirname(__FILE__), '..', 'lib'))
+$LOAD_PATH.unshift(File.dirname(__FILE__))
+
+require 'rubygems'
+require 'rspec'
+
+TESTFILES = File.dirname(__FILE__) + '/testfiles'
+
+# Requires supporting files with custom matchers and macros, etc,
+# in ./support/ and its subdirectories.
+Dir["#{File.dirname(__FILE__)}/support/**/*.rb"].each {|f| require f}
+
+RSpec.configure do |config|
+  config.color_enabled = true
+  config.tty = true
+  #config.formatter = :documentation  # :progress, :html, :textmate
+  config.formatter = :progress # :progress, :html, :textmate
+end
+
+
+def benchmark(*args, &block)
+  if ENV['BENCHMARK']
+    Benchmark.bmbm(*args, &block)
+  end
+end

data/spec/testfiles/ms/ident/peptide/db/uni_11_sp_tr.fasta

@@ -0,0 +1,69 @@
+>sp|P31946|1433B_HUMAN 14-3-3 protein beta/alpha OS=Homo sapiens GN=YWHAB PE=1 SV=3
+MTMDKSELVQKAKLAEQAERYDDMAAAMKAVTEQGHELSNEERNLLSVAYKNVVGARRSS
+WRVISSIEQKTERNEKKQQMGKEYREKIEAELQDICNDVLELLDKYLIPNATQPESKVFY
+LKMKGDYFRYLSEVASGDNKQTTVSNSQQAYQEAFEISKKEMQPTHPIRLGLALNFSVFY
+YEILNSPEKACSLAKTAFDEAIAELDTLNEESYKDSTLIMQLLRDNLTLWTSENQGDEGD
+AGEGEN
+>sp|P31946-2|1433B_HUMAN Isoform Short of 14-3-3 protein beta/alpha OS=Homo sapiens GN=YWHAB
+MDKSELVQKAKLAEQAERYDDMAAAMKAVTEQGHELSNEERNLLSVAYKNVVGARRSSWR
+VISSIEQKTERNEKKQQMGKEYREKIEAELQDICNDVLELLDKYLIPNATQPESKVFYLK
+MKGDYFRYLSEVASGDNKQTTVSNSQQAYQEAFEISKKEMQPTHPIRLGLALNFSVFYYE
+ILNSPEKACSLAKTAFDEAIAELDTLNEESYKDSTLIMQLLRDNLTLWTSENQGDEGDAG
+EGEN
+>sp|P62258|1433E_HUMAN 14-3-3 protein epsilon OS=Homo sapiens GN=YWHAE PE=1 SV=1
+MDDREDLVYQAKLAEQAERYDEMVESMKKVAGMDVELTVEERNLLSVAYKNVIGARRASW
+RIISSIEQKEENKGGEDKLKMIREYRQMVETELKLICCDILDVLDKHLIPAANTGESKVF
+YYKMKGDYHRYLAEFATGNDRKEAAENSLVAYKAASDIAMTELPPTHPIRLGLALNFSVF
+YYEILNSPDRACRLAKAAFDDAIAELDTLSEESYKDSTLIMQLLRDNLTLWTSDMQGDGE
+EQNKEALQDVEDENQ
+>sp|Q04917|1433F_HUMAN 14-3-3 protein eta OS=Homo sapiens GN=YWHAH PE=1 SV=4
+MGDREQLLQRARLAEQAERYDDMASAMKAVTELNEPLSNEDRNLLSVAYKNVVGARRSSW
+RVISSIEQKTMADGNEKKLEKVKAYREKIEKELETVCNDVLSLLDKFLIKNCNDFQYESK
+VFYLKMKGDYYRYLAEVASGEKKNSVVEASEAAYKEAFEISKEQMQPTHPIRLGLALNFS
+VFYYEIQNAPEQACLLAKQAFDDAIAELDTLNEDSYKDSTLIMQLLRDNLTLWTSDQQDE
+EAGEGN
+>sp|P61981|1433G_HUMAN 14-3-3 protein gamma OS=Homo sapiens GN=YWHAG PE=1 SV=2
+MVDREQLVQKARLAEQAERYDDMAAAMKNVTELNEPLSNEERNLLSVAYKNVVGARRSSW
+RVISSIEQKTSADGNEKKIEMVRAYREKIEKELEAVCQDVLSLLDNYLIKNCSETQYESK
+VFYLKMKGDYYRYLAEVATGEKRATVVESSEKAYSEAHEISKEHMQPTHPIRLGLALNYS
+VFYYEIQNAPEQACHLAKTAFDDAIAELDTLNEDSYKDSTLIMQLLRDNLTLWTSDQQDD
+DGGEGNN
+>sp|P31947|1433S_HUMAN 14-3-3 protein sigma OS=Homo sapiens GN=SFN PE=1 SV=1
+MERASLIQKAKLAEQAERYEDMAAFMKGAVEKGEELSCEERNLLSVAYKNVVGGQRAAWR
+VLSSIEQKSNEEGSEEKGPEVREYREKVETELQGVCDTVLGLLDSHLIKEAGDAESRVFY
+LKMKGDYYRYLAEVATGDDKKRIIDSARSAYQEAMDISKKEMPPTNPIRLGLALNFSVFH
+YEIANSPEEAISLAKTTFDEAMADLHTLSEDSYKDSTLIMQLLRDNLTLWTADNAGEEGG
+EAPQEPQS
+>tr|D2KLI3|D2KLI3_HUMAN Cytochrome b OS=Homo sapiens GN=CYTB PE=3 SV=1
+MTPTRKTNPLMKLINHSFIDLPTPSNISAWWNFGSLLGACLILQITTGLFLAMHYSPDAS
+TAFSSIAHITRDVNYGWIIRYLHANGASMFFICLFLHIGRGLYYGSFLYSETWNIGIILL
+LATMATAFMGYVLPWGQMSFWGATVITNLLSAIPYIGTDLVQWIWGGYSVDSPTLTRFFT
+FHFILPFIIAALAALHLLFLHETGSNNPLGITSHSDKITFHPYYTIKDALGLLLFLLSLM
+TLTLFSPDLLGDPDNYTLANPLNTPPHIKPEWYFLFAYTILRSVPNKLGGVLALLLSILI
+LAMIPILHMSKQQSMMFRPLSQSLYWLLAADLLILTWIGGQPVSYPFTIIGQVASVLYFT
+TILILMPTISLIENKMLKWA
+>tr|D2KTA8|D2KTA8_HUMAN Putative uncharacterized protein FCAMR OS=Homo sapiens GN=FCAMR PE=4 SV=1
+MDGEATVKPGEQVPLWTHGWPPDDPSPSFAAGSSFALPQKRPHPRWLWEGSLPSRTHLRA
+MGTLRPSSPLCWREESSFAAPNSLKGSRLVSGEPGGAVTIQCHYAPSSVNRHQRKYWCRL
+GPPRWICQTIVSTNQYTHHRYRDRVALTDFPQRGLFVVRLSQLSPDDIGCYLCGIGSENN
+MLFLSMNLTISAGPASTLPTATPAAGELTMRSYGTASPVANRWTPGTTQTLGQGTAWDTV
+ASTPGTSKTTASAEGRRTPGATRPAAPGTGSWAEGSVKAPAPIPESPPSKSRSMSNTTEG
+VWEGTRSSVTNRARASKDRREMTTTKADRPREDIEGVRIALDAAKKVLGTIGPPALVSET
+LAWEILPQATPVSKQQSQGSIGETTPAAGMWTLGTPAADVWILGTPAADVWTSMEAASGE
+GSAAGDLDAATGDRGPQATLSQTPAVGPWGPPGKESSVKRTFPEDESSSRTLAPVSTMLA
+LFMLMALVLLQRKLWRRRTSQEAERVTLIQMTHFLEVNPQADQLPHVERKMLQDDSLPAG
+ASLTAPERNPGP
+>tr|D2KTA9|D2KTA9_HUMAN Putative uncharacterized protein XKR5 OS=Homo sapiens GN=XKR5 PE=4 SV=1
+MHARLLGLSALLQAAEQSARLYTVAYYFTTGRLLWGWLALAVLLPGFLVQALSYLWFRAD
+GHPGHCSLMMLHLLQLGVWKRHWDAALTSLQKELEAPHRGWLQLQEADLSALRLLEALLQ
+TGPHLLLQTYVFLASDFTDIVPGVSTLFSWSSLSWALVSYTRFMGFMKPGHLAMPWAALF
+CQQLWRMGMLGTRVLSLVLFYKAYHFWVFVVAGAHWLVMTFWLVAQQSDIIDSTCHWRLF
+NLLVGAVYILCYLSFWDSPSRNRMVTFYMVMLLENIILLLLATDFLQGASVDQPADHSWG
+PVWISDWQCLTGNLLQPAASKIHRHLAGLPKEVLWHCRR
+>tr|D3DSH8|D3DSH8_HUMAN HCG2036819, isoform CRA_a OS=Homo sapiens GN=hCG_2036819 PE=4 SV=1
+MLGWIQPSRQPQLRAAPPTRTPSAKRCILCNFLPGCWLVGDVAGSRQPSAPQTLRQRQHT
+RPPPQERGSGRRSPLREARRANPHFKSFPVLEARGLPCGARRTGPRRPVREMTLPSDPER
+ATLPNPRLGAPAVPRRGPRSHGGRR
+>tr|D3DX18|D3DX18_HUMAN Putative uncharacterized protein LOC128977 OS=Homo sapiens GN=LOC128977 PE=4 SV=1
+MADGSGWQPPRPCEAYRAEWKLCRSARHFLHHYYVHGERPACEQWQRDLASCRDWEERRN
+AEAQQSLCESERARVRAARKHILVWAPRQSPPPDWHLPLPQEKDE