jandot-ruby-ensembl-api 0.9.2

data/test/unit/release_45/variation/test_activerecord.rb

@@ -0,0 +1,32 @@
+#
+# = test/unit/test_seq.rb - Unit test for Ensembl::Core
+#
+# Copyright::   Copyright (C) 2008
+#               Jan Aerts <http://jandot.myopenid.com>
+# License::     Ruby's
+#
+# $Id:
+require 'pathname'
+libpath = Pathname.new(File.join(File.dirname(__FILE__), ['..'] * 4, 'lib')).cleanpath.to_s
+$:.unshift(libpath) unless $:.include?(libpath)
+
+require 'test/unit'
+require 'yaml'
+require 'ensembl'
+
+include Ensembl::Variation
+
+DBConnection.connect('homo_sapiens')
+
+class Simple < Test::Unit::TestCase
+  def test_allele
+    allele = Allele.find(1)
+    assert_equal('T', allele.allele)
+    assert_equal(0.04, allele.frequency)
+  end
+  
+  def test_allele_group
+    allele_group = AlleleGroup.find(1)
+    assert_equal('ABDR-1', allele_group.name)
+  end
+end

data/test/unit/release_50/core/run_tests.rb

@@ -0,0 +1,12 @@
+require 'test/unit'
+
+# Either run the tests that use Bos taurus
+require 'test_project'
+require 'test_transform'
+require 'test_slice'
+
+# Or the ones using Homo sapiens
+#require 'test_relationships'
+#require 'test_project_human'
+#require 'test_sequence'
+#require 'test_transcript'

data/test/unit/release_50/core/test_project.rb

@@ -0,0 +1,215 @@
+#
+# = test/unit/test_project.rb - Unit test for Ensembl::Core
+#
+# Copyright::   Copyright (C) 2007
+#               Jan Aerts <http://jandot.myopenid.com>
+# License::     Ruby's
+#
+# $Id:
+require 'pathname'
+libpath = Pathname.new(File.join(File.dirname(__FILE__), ['..'] * 4, 'lib')).cleanpath.to_s
+$:.unshift(libpath) unless $:.include?(libpath)
+
+require 'test/unit'
+require 'ensembl'
+
+include Ensembl::Core
+
+DBConnection.connect('bos_taurus', 50)
+
+class CoordinateMappingsTestSimple < Test::Unit::TestCase
+  # First see if the relationships work
+  def test_assemblies
+    # Contig AAFC03055291 should only be a component of chromosome 20
+    contig_coord_system = CoordSystem.find_by_name('contig')
+    aafc03055291 = SeqRegion.find_by_name_and_coord_system_id('AAFC03055291', contig_coord_system.id)
+    assert_equal(1, aafc03055291.assembled_seq_regions.length)
+
+    # Chromosome 20 has 2970 components
+    chr_coord_system = CoordSystem.find_by_name('chromosome')
+    chr20 = SeqRegion.find_by_name_and_coord_system_id('20', chr_coord_system.id)
+    assert_equal(2970, chr20.component_seq_regions.length)
+
+    # Chromosome 20 has 2970 contigs
+    assert_equal(2970, chr20.component_seq_regions('contig').length)
+
+    # Positions of the link between Chr20 and AAFC03055291
+    # * Contig AAFC03055291 starts at position 13970982 on chromosome Chr20
+    assert_equal(13970982, aafc03055291.assembly_links_as_component('chromosome')[0].asm_start)
+    end
+end
+
+class Sequences < Test::Unit::TestCase
+  def setup
+    @seq_region = SeqRegion.find(92594)
+  end
+
+  def test_simple
+    assert_equal('AGCTATTTTATGACTT', @seq_region.seq.slice(4,16))
+  end
+
+  def test_subseq
+    assert_equal('AGCTATTTTATGACTT', @seq_region.subseq(5,20))
+  end
+end
+
+#class SliceProjectFromComponentToAssembly < Test::Unit::TestCase
+#  # |------------------------------------------> chromosome
+#  #     ^                 ^
+#  #     |                 |
+#  #     |-----------------> scaffold
+#  def test_project_from_whole_component_to_assembly
+#    source_slice = Slice.fetch_by_region('contig','AAFC03055291')
+#    target_slices = source_slice.project('chromosome')
+#
+#    # Start and stop of chr4_105 on Chr4
+#    assert_equal(13970982, target_slices[0].start)
+#    assert_equal(13982069, target_slices[0].stop)
+#  end
+#
+#  # |------------------------------------------> chromosome
+#  #         ^        ^
+#  #         |        |
+#  #     |-----------------> scaffold
+#  def test_project_from_component_to_assembly_with_positions
+#    source_slice = Slice.fetch_by_region('scaffold','Chr4.003.105', 42, 2007)
+#    target_slices = source_slice.project('chromosome')
+#
+#    # Position 42 on chr4_105 is position 96652152, position 2007 is 96654117
+#    assert_equal(96652152, target_slices[0].start)
+#    assert_equal(96654117, target_slices[0].stop)
+#  end
+#
+#  # |------------------------------------------> scaffold
+#  #         ^        ^
+#  #         |        |
+#  #       ----------------> contig
+#  #      /
+#  #   |--
+#  def test_project_from_component_to_assembly_with_positions_and_cmp_start_not_1
+#    source_slice = Slice.fetch_by_region('contig', 'AAFC03020247', 42, 2007)
+#    target_slices = source_slice.project('scaffold')
+#
+#    # Position 42 on AAFC03020247 is position 6570 on ChrUn.003.3522, position 2007 is 8565
+#    assert_equal(6570, target_slices[0].start)
+#    assert_equal(8535, target_slices[0].stop)
+#  end
+#
+#  # |------------------------------------------> scaffold
+#  #         ^        ^
+#  #         |        |
+#  #     <-----------------| contig
+#  def test_project_from_component_to_assembly_with_strand
+#    source_slice_fw = Slice.fetch_by_region('contig', 'AAFC03020247')
+#    target_slices_fw = source_slice_fw.project('scaffold')
+#
+#    assert_equal(1, target_slices_fw[0].strand)
+#
+#    source_slice_rev = Slice.fetch_by_region('contig', 'AAFC03061502')
+#    target_slices_rev = source_slice_rev.project('scaffold')
+#
+#    assert_equal(-1, target_slices_rev[0].strand)
+#  end
+#end
+
+#class SliceProjectFromComponentToAssemblyUsingTopLevel < Test::Unit::TestCase
+#  # |------------------------------------------> chromosome
+#  #     ^                 ^
+#  #     |                 |
+#  #     |-----------------> scaffold
+#  def test_project_from_whole_component_to_assembly
+#    source_slice = Slice.fetch_by_region('scaffold','Chr4.003.105')
+#    target_slices = source_slice.project('toplevel')
+#
+#    # Start and stop of chr4_105 on Chr4
+#    assert_equal(96652111, target_slices[0].start)
+#    assert_equal(97251689, target_slices[0].stop)
+#  end
+#
+#  # |------------------------------------------> chromosome
+#  #         ^        ^
+#  #         |        |
+#  #     |-----------------> scaffold
+#  def test_project_from_component_to_assembly_with_positions
+#    source_slice = Slice.fetch_by_region('scaffold','Chr4.003.105', 42, 2007)
+#    target_slices = source_slice.project('toplevel')
+#
+#    # Position 42 on chr4_105 is position 96652152, position 2007 is 96654117
+#    assert_equal(96652152, target_slices[0].start)
+#    assert_equal(96654117, target_slices[0].stop)
+#  end
+#end
+
+class SliceProjectFromAssemblyToComponentForwardStrands < Test::Unit::TestCase
+  def setup
+    @source_slice_single_contig = Slice.fetch_by_region('chromosome', '20', 175000, 180000)
+    @target_slices_single_contig = @source_slice_single_contig.project('contig')
+
+    @source_slice_two_contigs = Slice.fetch_by_region('chromosome','20', 175000, 190000)
+    @target_slices_two_contigs = @source_slice_two_contigs.project('contig')
+
+    @source_slice_contigs_with_strand = Slice.fetch_by_region('chromosome', '20', 160000, 190000)
+    @target_slices_contigs_with_strand = @source_slice_contigs_with_strand.project('contig')
+
+    @source_slice_contigs_with_strand_ends_in_gaps = Slice.fetch_by_region('chromosome', '20', 170950, 196000)
+    @target_slices_contigs_with_strand_ends_in_gaps = @source_slice_contigs_with_strand_ends_in_gaps.project('contig')
+  end
+
+  #     |-----------------> contig
+  #         ^          ^
+  #         |          |
+  # |------------------------------------------> chromosome
+  def test_project_from_assembly_to_single_component
+    # Position 175000 on chr20 is position 4030 on contig, position 180000 is 9030
+    assert_equal('AAFC03028970', @target_slices_single_contig[0].seq_region.name)
+    assert_equal(4030, @target_slices_single_contig[0].start)
+    assert_equal(9030, @target_slices_single_contig[0].stop)
+  end
+
+  #     |----->   |--------> contig
+  #         ^          ^
+  #         |          |
+  # |------------------------------------------> chromosome
+  def test_project_from_assembly_to_two_components
+    # This chromosomal region is covered by contigs AAFC03028970, a gap and AAFC03028962
+    #  * Position 175000 on chr 20 is position 4030 on contig AAFC03028970
+    #  * Position 190000 on chr 20 is position 35 on contig AAFC03028962
+    assert_equal(3, @target_slices_two_contigs.length)
+    assert_equal('contig:Btau_4.0:AAFC03028970:4030:17365:1', @target_slices_two_contigs[0].display_name)
+    assert_equal(Gap, @target_slices_two_contigs[1].class)
+    assert_equal('contig:Btau_4.0:AAFC03028962:1:35:1', @target_slices_two_contigs[2].display_name)
+  end
+
+  #     |----->   <-------|   |------->  |------->  contig
+  #         ^                                 ^
+  #         |                                 |
+  # |--------------------------------------------------> chromosome
+  def test_project_from_assembly_to_contigs_with_strand
+    # This chromosomal region is covered by 4 contigs and 3 gaps
+    # One of the contigs are on the reverse strand.
+    assert_equal(7, @target_slices_contigs_with_strand.length)
+    assert_equal('contig:Btau_4.0:AAFC03028964:90:9214:1', @target_slices_contigs_with_strand[0].display_name)
+    assert_equal(Gap, @target_slices_contigs_with_strand[1].class)
+    assert_equal('contig:Btau_4.0:AAFC03028959:1:1746:-1', @target_slices_contigs_with_strand[2].display_name)
+    assert_equal(Gap, @target_slices_contigs_with_strand[3].class)
+    assert_equal('contig:Btau_4.0:AAFC03028970:1:17365:1', @target_slices_contigs_with_strand[4].display_name)
+    assert_equal(Gap, @target_slices_contigs_with_strand[5].class)
+    assert_equal('contig:Btau_4.0:AAFC03028962:1:35:1', @target_slices_contigs_with_strand[6].display_name)
+  end
+
+  #      <--|  |----->        contig
+  #   ^                 ^
+  #   |                 |
+  # |--------------------------------------------------> chromosome
+  def test_project_from_assembly_to_contigs_with_strand_and_ending_in_gaps
+    # This chromosomal region is covered by 2 contigs and 2 gaps at the end: GaCoGaCoGa
+    assert_equal(5, @target_slices_contigs_with_strand_ends_in_gaps.length)
+    assert_equal(Gap, @target_slices_contigs_with_strand_ends_in_gaps[0].class)
+    assert_equal('contig:Btau_4.0:AAFC03028970:1:17365:1', @target_slices_contigs_with_strand_ends_in_gaps[1].display_name)
+    assert_equal(Gap, @target_slices_contigs_with_strand_ends_in_gaps[2].class)
+    assert_equal('contig:Btau_4.0:AAFC03028962:1:5704:1', @target_slices_contigs_with_strand_ends_in_gaps[3].display_name)
+    assert_equal(Gap, @target_slices_contigs_with_strand_ends_in_gaps[4].class)
+  end
+
+
+end

data/test/unit/release_50/core/test_project_human.rb

@@ -0,0 +1,58 @@
+#
+# = test/unit/test_project.rb - Unit test for Ensembl::Core
+#
+# Copyright::   Copyright (C) 2007
+#               Jan Aerts <http://jandot.myopenid.com>
+# License::     Ruby's
+#
+# $Id:
+require 'pathname'
+libpath = Pathname.new(File.join(File.dirname(__FILE__), ['..'] * 4, 'lib')).cleanpath.to_s
+$:.unshift(libpath) unless $:.include?(libpath)
+
+require 'test/unit'
+require 'ensembl'
+
+include Ensembl::Core
+
+DBConnection.connect('homo_sapiens', 50)
+
+class AssemblyExceptions < Test::Unit::TestCase
+  def test_chr_x
+    source_slice = Slice.fetch_by_region('chromosome','X', 2709497, 2709520)
+    assert_equal('tagttatagattaaaagaagttaa', source_slice.seq)
+  end
+
+  def test_slice_overlapping_PAR_and_allosome
+    source_slice = Slice.fetch_by_region('chromosome','Y',2709500,2709540)
+    target_slices = source_slice.project('contig')
+    assert_equal('contig::AC006209.25.1.141759:23323:23343:-1', target_slices[0].display_name)
+    assert_equal('contig::AC006040.3.1.186504:57272:57291:1', target_slices[1].display_name)
+  end
+
+  def test_seq_slice_overlapping_PAR
+    seq = ''
+    File.open('../../data/seq_y.fa').reject{|l| l=~/^>/}.each do |line|
+      line.chomp!
+      seq += line
+    end
+    seq.downcase!
+
+    source_slice = Slice.fetch_by_region('chromosome', 'Y', 2709497, 2709542)
+    assert_equal(seq.downcase, source_slice.seq)
+  end
+
+  # The MHC haplotypes for human are not implemented yet, so we raise an error
+  # in the code.
+  def test_seq_slice_overlapping_HAP
+    seq = ''
+    File.open('../../data/seq_c6qbl.fa').reject{|l| l=~/^>/}.each do |line|
+      line.chomp!
+      seq += line
+    end
+    seq.downcase!
+
+    source_slice = Slice.fetch_by_region('chromosome', 'c6_QBL', 33451191, 33451690)
+    assert_raise(NotImplementedError) {source_slice.seq}
+  end
+end

data/test/unit/release_50/core/test_relationships.rb

@@ -0,0 +1,66 @@
+#
+# = test/unit/test_transfers.rb - Unit test for Ensembl::Core
+#
+# Copyright::   Copyright (C) 2007
+#               Jan Aerts <http://jandot.myopenid.com>
+# License::     Ruby's
+#
+# $Id:
+require 'pathname'
+libpath = Pathname.new(File.join(File.dirname(__FILE__), ['..'] * 4, 'lib')).cleanpath.to_s
+$:.unshift(libpath) unless $:.include?(libpath)
+
+require 'test/unit'
+
+require 'ensembl'
+
+include Ensembl::Core
+
+DBConnection.connect('homo_sapiens', 50)
+
+# Let's see if we can 'find' things
+class SimpleRecordsTest < Test::Unit::TestCase
+  def setup
+    @sry_gene = Gene.find(34927)
+    @sry_transcript = Transcript.find(60290)
+  end
+
+  def test_coord_system
+    coord_system = CoordSystem.find(17)
+    assert_equal('chromosome', coord_system.name)
+  end
+
+  def test_coord_system_toplevel
+    coord_system = CoordSystem.find_toplevel
+    assert_equal('chromosome', coord_system.name)
+  end
+
+  def test_coord_system_seqlevel
+    coord_system = CoordSystem.find_seqlevel
+    assert_equal('contig', coord_system.name)
+  end
+
+  def test_display_label
+    assert_equal('SRY', @sry_gene.display_label)
+    assert_equal('SRY', @sry_gene.display_name)
+    assert_equal('SRY', @sry_gene.label)
+    assert_equal('SRY', @sry_gene.name)
+
+    assert_equal('SRY-001', @sry_transcript.display_label)
+    assert_equal('SRY-001', @sry_transcript.display_name)
+    assert_equal('SRY-001', @sry_transcript.label)
+    assert_equal('SRY-001', @sry_transcript.name)
+  end
+end
+
+class RelationshipsTest < Test::Unit::TestCase
+  def test_go_terms
+    gene = Gene.find(34928)
+    assert_equal(["GO:0005576", "GO:0042742"], gene.go_terms.sort)
+  end
+  
+  def test_hgnc
+    gene = Gene.find_by_stable_id('ENSG00000169740')
+    assert_equal('ZNF32', gene.hgnc)
+  end
+end

data/test/unit/release_50/core/test_sequence.rb

@@ -0,0 +1,175 @@
+#
+# = test/unit/test_seq.rb - Unit test for Ensembl::Core
+#
+# Copyright::   Copyright (C) 2007
+#               Jan Aerts <http://jandot.myopenid.com>
+# License::     Ruby's
+#
+# $Id:
+require 'pathname'
+libpath = Pathname.new(File.join(File.dirname(__FILE__), ['..'] * 4, 'lib')).cleanpath.to_s
+$:.unshift(libpath) unless $:.include?(libpath)
+
+require 'test/unit'
+require 'yaml'
+require 'ensembl'
+
+include Ensembl::Core
+
+DBConnection.connect('homo_sapiens', 50)
+
+#class SequenceForSlice < Test::Unit::TestCase
+#  def test_forward_strand_seqlevel
+#    slice = Slice.new(SeqRegion.find(170931),5,15)
+#    seq = 'gcagtggtgtg'
+#    assert_equal(seq, slice.seq)
+#  end
+#
+#  def test_reverse_strand_seqlevel
+#    slice = Slice.new(SeqRegion.find(170931),5,15, -1)
+#    seq = 'cacaccactgc'
+#    assert_equal(seq, slice.seq)
+#  end
+#
+#  def test_forward_strand_not_seqlevel_single_target
+#    slice = Slice.new(SeqRegion.find(226044),69437100,69437110)
+#    seq = 'gtctatttaca'
+#    assert_equal(seq, slice.seq)
+#  end
+#
+#  def test_reverse_strand_not_seqlevel_single_target
+#    slice = Slice.new(SeqRegion.find(226044),69437100,69437110,-1)
+#    seq = 'tgtaaatagac'
+#    assert_equal(seq, slice.seq)
+#  end
+#
+#  def test_forward_strand_not_seqlevel_composite_target
+#    seq = ''
+#    File.open('../../data/seq_forward_composite.fa').reject{|l| l=~/^>/}.each do |line|
+#      line.chomp!
+#      seq += line
+#    end
+#    seq.downcase!
+#    slice = Slice.new(SeqRegion.find(226044),69387650,69487649)
+#    assert_equal(seq, slice.seq)
+#  end
+#
+#  def test_reverse_strand_not_seqlevel_composite_target
+#    slice = Slice.new(SeqRegion.find(226044),69437061,69437160,-1)
+#    assert_equal('nnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnatgtaaatagacaactaacaagatagctgagattgtttccttatccagaca', slice.seq)
+#  end
+#
+#end
+#
+#class SequenceForUnsplicedFeature < Test::Unit::TestCase
+#  def test_forward_strand_seqlevel
+#    marker_feature = MarkerFeature.find(1323757)
+#    marker_seq = 'ggcttacttggaaaggtctcttccaacccaatattattcaaatactttcaattttcttctaatgtttttagtgtgtgtgtgtgtgtgtgtgtgtgtgtgtgtgtgtgtgtggttttgttttatttcttttaattctctgatacatttagaatttcttttattattttattttattttattatttatttatttatttttgagacagagttttgctc'
+#    assert_equal(marker_seq, marker_feature.seq)
+#  end
+#
+#  def test_reverse_strand_seqlevel
+#    gene = Gene.find_by_name('ANKRD56')
+#
+#    ankrd56_seq = 'atggcccgagagctgagccaggaggcactactggactttctgtgccaggctgggggccgcgtgaccaacgctgccttgctgagccacttcaagagctttctccgagaccccgacgcgtcccccagccagcaccagcaccgccgcgagctcttcaagggcttcgtcaactcggtcgccgcagtgcgccaggaccccgacggcaccaagtacgtggtgctcaagaggagatacagggaccttttgggggaggaggggctgcagcgaccccgcgagccgcccgcggccgcccccagtgcagggggagctgcgccctgctccccgcgaggcgcgcgccggggggagccgccccagcagcagcccaggcggcggcggcgcgagaaggagccggaggaggagccagcaggtgcagcagccagagccgccgacgcagcttgcaatggactcccgggcagcgactcccgtagggcgcccgggaagggcggcggatcgaagggcagtcccggacagaggccgccggtgcccgcagctgcagcggcaggggcccaggcgagagcgagctgcgcggcggcgaagacgcagggccgctgctgctgggaatgcctccagaacaacctggctgtactgccgggagagctcggcgcactcccgcactcggccaccgcggaggagaagccggcacgggctctgcctgcccaggatgaccgcggggcttccagggagcgggaagaaggcgcgctagctgagccggcgcctgtgcctgcagtggctcactcgcctcccgccaccgtcgaggctgcgacaagcagggcttccccgcctgctctcctgcccggccccgctccccgcggagaccggccggagctgctgacccccagctccctgcattattcgaccctgcagcagcagcagcagcgcactcgagagtgggtggccaggcacccgcaggtgcccgaggcccgtgatcagggccctatccgcgcctggtcggtgctgccagacaacttcctccagctgcccttggaacccggctccacggagcctaattcagagccgccagacccctgtctttcctcgcactctctctttcctgttgttccggatgagtcctgggaatcctgggcggggaacccttcattgactgtctttcgcagcattcgttgtcagctgtccctccaagatctggatgactttgtggaccaggagagtgatggcagtgaggagagcagcagtgggcccaaagactccccgggggcttctgaagaggggctgcaggttgtcttgggaaccccagatagggggaagctcaggaatccagctgggggcctttctgtatctcggaaggagggcagccccagccggagccctcagggtctcagaaacagaggggatggtcacatctctcagcaggtccctgcaggggctaatggccttgcaggccaccccctgaagcctttgccttggccagttcctaagttaaggaggtccctcaggaggagctctctggcagggagagccaaattgtcctcctctgatgaggagtacctcgatgagggcttgctgaaaagaagtcggcgcccacctcgatccaggaagccctccaaggcaggaacggcacccagcccaagggttgatgcaggtttatcactaaaacttgcagaggttaaggctgttgtggccgagcggggttggcgacacagcctgtgggtccccagtggggaggggtctgcagccttggccccccacagaacttctgagcacaaatcatccctggttccactagatgccagggagcatgagtggattgtgaagcttgccagtggctcctggattcaggtgtggactttgttctgggaggaccctcaactggccttgcacaaagactttttgactgggtacactgcgttgcactggatagccaaacatggtgacctcagggcccttcaggacttggtgtctggagcaaagaaggcagggattgtccttgatgtaaacgtgaggtccagttgtggatataccccgctgcaccttgcagccattcacggccaccagggggtcatcaaattgctagtgcaaaggttggcttctcgggtaaatgtcagggacagcagtgggaagaagccatggcagtatctaaccagtaatacctctggggaaatatggcagctgttgggagctcctcggggcaagcccattttccctgtctatcccttagttggaagttcttcccctaccagaaaggccaagagcaaggaaatatctagaagtgtcacccgaaaaacttccttcgctgcactactcaaaagtcagcacaacaagtggaaactggccaaccagtatgagaaattccacagtccaagggaaagagaagagtatagtgactga'
+#    assert_equal(ankrd56_seq, gene.seq)
+#  end
+#
+#  def test_reverse_strand_not_seqlevel
+#    gene = Gene.find_by_name('DRD3')
+#    drd3_gene_seq = ''
+#    File.open('../../data/seq_drd3_gene.fa').reject{|l| l=~/^>/}.each do |line|
+#      line.chomp!
+#      drd3_gene_seq += line
+#    end
+#    drd3_gene_seq.downcase!
+#    assert_equal(drd3_gene_seq, gene.seq)
+#
+#  end
+#
+#  def test_exon
+#    exon = Exon.find(719588)
+#    assert_equal('atggcatctctgagccagctgagtggccacctgaactacacctgtggggcagagaactccacaggtgccagccaggcccgcccacatgcctactatgccctctcctactgcgcgctcatcctggccatcgtcttcggcaatggcctggtgtgcatggctgtgctgaaggagcgggccctgcagactaccaccaactacttagtagtgagcctggctgtggcagacttgctggtggccaccttggtgatgccctgggtggtatacctggag', exon.seq)
+#  end
+#
+#end
+
+class SequenceForSlicedFeature < Test::Unit::TestCase
+  def test_transcript_foward
+    transcript = Transcript.find(73491) # UB2R1 = CDC34
+    ub2r1_transcript_seq = ''
+    File.open('../../data/seq_ub2r1_transcript.fa').reject{|l| l=~/^>/}.each do |line|
+      line.chomp!
+      ub2r1_transcript_seq += line
+    end
+    ub2r1_transcript_seq.downcase!
+    assert_equal(ub2r1_transcript_seq, transcript.seq)
+
+  end
+
+  def test_transcript_reverse
+    transcript = Transcript.find(107548)
+    rnd3_transcript_seq = ''
+    File.open('../../data/seq_rnd3_transcript.fa').reject{|l| l=~/^>/}.each do |line|
+      line.chomp!
+      rnd3_transcript_seq += line
+    end
+    rnd3_transcript_seq.downcase!
+    assert_equal(rnd3_transcript_seq, transcript.seq)
+  end
+
+end
+
+class SequenceForCDS < Test::Unit::TestCase
+  def setup
+    # Transcript tr_fw is ENST00000215574
+    @tr_fw = Transcript.find(73491)
+    # Transcript tr_rev is ENST00000315489
+    @tr_rev = Transcript.find(73411)
+  end
+
+  def test_cds_fw
+    ub2r1_coding_seq = ''
+    File.open('../../data/seq_ub2r1_coding.fa').reject{|l| l=~/^>/}.each do |line|
+      line.chomp!
+      ub2r1_coding_seq += line
+    end
+    ub2r1_coding_seq.downcase!
+    assert_equal(ub2r1_coding_seq, @tr_fw.cds_seq)
+  end
+
+  def test_cds_rev
+    cso19_coding_seq = ''
+    File.open('../../data/seq_cso19_coding.fa').reject{|l| l=~/^>/}.each do |line|
+      line.chomp!
+      cso19_coding_seq += line
+    end
+    cso19_coding_seq.downcase!
+    assert_equal(cso19_coding_seq, @tr_rev.cds_seq)
+  end
+
+  def test_five_prime_utr_fw
+    assert_equal('GGCAAGCGCCGGTGGGGCGGCGGCGCCAGAGCTGCTGGAGCGCTCGGGGTCCCCGGGCGGCGGCGGCGGCGCAGAGGAGGAGGCAGGCGGCGGCCCCGGTGGCTCCCCCCCGGACGGTGCGCGGCCCGGCCCGTCTCGCGAACTCGCGGTGGTCGCGCGGCCCCGCGCTGCTCCGACCCCGGGCCCCTCCGCCGCCGCC'.downcase, @tr_fw.five_prime_utr_seq)
+  end
+
+  def test_five_primer_utr_rev
+    assert_equal('ACTCGATCCCGGCCCCACTTCCAGGCCAGTGTCCGGCCGACCAGCCTGCCTTGGGCCAGGGCCCCACGACTCCCTGCTGCGGGACAAGAGGCCGTCTGTGCGGCTGTGGTCGTGGGAGGGTGTGGTGAGGCCGTGAAGGTGGGGACGGTGCCTGGGCCTGTGGCCGCCAGAGCTGCTGCGGCTCAGAAGGTAGCACCAGGCCCCGTGGGTGCTGTGGGGGCCATCGCCTGCCCACC'.downcase, @tr_rev.five_prime_utr_seq)
+  end
+
+  def test_three_prime_utr_fw
+    assert_equal('CACCACCAGAATAAACTTGCCGAGTTTACCTCACTAGGGCCGGACCCGTGGCTCCTTAGACGACAGACTACCTCACGGAGGTTTTGTGCTGGTCCCCGTCTCCTCTGGTTGTTTCGTTTTGGCTTTTTCTCCCTCCCCATGTCTGTTCTGGGTTTTCACGTGCTTCAGAGAAGAGGGGCTGCCCCACCGCCACTCACGTCACTCGGGGCTCGGTGGACGGGCCCAGGGTGGGAGCGGCCGGCCCACCTGTCCCCTCGGGAGGGGAGCTGAGCCCGACTTCTACCGGGGTCCCCCAGCTTCCGGACTGGCCGCACCCCGGAGGAGCCACGGGGGCGCTGCTGGGAACGTGGGCGGGGGGCCGTTTCCTGACACTACCAGCCTGGGAGGCCCAGGTGTAGCGGTCCGAGGGGCCCGGTCCTGCCTGTCAGCTCCAGGTCCTGGAGCCACGTCCAGCACAGAGTGGACGGATTCACCGTGGCCGACTCTTTTCCCTGCTTTGGTTTGTTTGAAATCTAAATAAAACTACTTTATG'.downcase, @tr_fw.three_prime_utr_seq)
+  end
+
+  def test_three_prime_utr_rev
+    assert_equal('GCCAGGCCTGTTTCCCACACAGTGAACGGGGGGCTGGCGGGCCTCTCCTGGGGCCTGGCTTTGGCTGGAGGGGCTAGAAGGCGAGAGGGGCTGGGAAAGGAGTCTTCCCTCCCCTTTCCTGAACTGTCCAGGTCCTTCTAGGACACCTGATCCTTCCAGTCCCTGGGGGCTGTGACCCATGGCCCTGTGCAGGGTGCAGGGTGAGTCTCCTTCCAGTGCCCCAGTTTCTTCCAGGTCACGCCAGGCCAGGCCAGGCCAGGTGGGGAAAGGGACACCTTCCGGCCCTCCCCAGTAGCTGGCTGGAGATGGAGCTTCCTGTGTCCCAGAACTCGGCGTCCAGCTCCACTAGGGCCTGGATCCCCATCACAGCTTGGGTTAGCCCCGGTCCCAGCCCAAACTCAGGCTGGAGGCAGCCCCGAGGCCTGTGCCTTTCCCACTCCACCTTCTACAGTTGCTTAGCCAATAAACCTTTCCTGGGCTGGAG'.downcase, @tr_rev.three_prime_utr_seq)
+  end
+
+  def test_protein_fw
+    assert_equal('MARPLVPSSQKALLLELKGLQEEPVEGFRVTLVDEGDLYNWEVAIFGPPNTYYEGGYFKARLKFPIDYPYSPPAFRFLTKMWHPNIYETGDVCISILHPPVDDPQSGELPSERWNPTQNVRTILLSVISLLNEPNTFSPANVDASVMYRKWKESKGKDREYTDIIRKQVLGTKVDAERDGVKVPTTLAEYCVKTKAPAPDEGSDLFYDDYYEDGEVEEEADSCFGDDEDDSGTEES*', @tr_fw.protein_seq)
+  end
+
+  def test_protein_rev
+    assert_equal('MGTLSCDSTPRLATAPLGRRVTEGQIPETGLRKSCGTATLENGSGPGLYVLPSTVGFINHDCTRVASPAYSLVRRPSEAPPQDTSPGPIYFLDPKVTRFGRSCTPAYSMQGRAKSRGPEVTPGPGAYSPEKVPPVRHRTPPAFTLGCRLPLKPLDTSAPAPNAYTMPPLWGSQIFTKPSSPSYTVVGRTPPARPPQDPAEIPGPGQYDSPDANTYRQRLPAFTMLGRPRAPRPLEETPGPGAHCPEQVTVNKARAPAFSMGIRHSKRASTMAATTPSRPAGHRLPGRCC*', @tr_rev.protein_seq)
+  end
+end
+