spec/gff3_assemble_spec.rb in bio-gff3-0.8.2 vs spec/gff3_assemble_spec.rb in bio-gff3-0.8.3

- old
+ new

@@ -75,21 +75,21 @@ it "should translate CDS 7838:7980 (in frame 0, + strand)" do recs = @cdslist['cds:MhA1_Contig1133.frz3.gene4'] component = @componentlist['cds:MhA1_Contig1133.frz3.gene4'] cds0 = recs[0] cds0.seqname.should == 'MhA1_Contig1133' - seq = @gff.assemble(@contigsequence,component.start,[cds0]) + seq = @gff.assemble(@contigsequence,component.start,[cds0],:trim=>false) seq.size.should == 143 seq.should == "ATGCGTCCTTTAACAGATGAAGAAACTGAAAAGTTTTTCAAAAAACTTTCAAATTATATTGGTGACAATATTAAACTTTTATTGGAAAGAGAAGATGGAGAATATGTTTTTCGTTTACATAAAGACAGAGTTTATTATTGCAG" aaseq = @gff.assembleAA(@contigsequence,component.start,[cds0]) aaseq.should == "MRPLTDEETEKFFKKLSNYIGDNIKLLLEREDGEYVFRLHKDRVYYC" end # MhA1_Contig1133 WormBase CDS 8065 8308 . + 1 ID=cds:MhA1_Contig1133.frz3.gene4;Parent=transcript:MhA1_Contig1133.frz3.gene4 it "should translate CDS 8065:8308 (in frame 1, + strand)" do recs = @cdslist['cds:MhA1_Contig1133.frz3.gene4'] component = @componentlist['cds:MhA1_Contig1133.frz3.gene4'] cds1 = recs[1] - seq = @gff.assemble(@contigsequence,component.start,[cds1]) + seq = @gff.assemble(@contigsequence,component.start,[cds1],:trim=>false) seq.size.should == 244 seq.should == "TGAAAAATTAATGCGACAAGCAGCATGTATTGGACGTAAACAATTGGGATCTTTTGGAACTTGTTTGGGTAAATTCACAAAAGGAGGGTCTTTCTTTCTTCATATAACATCATTGGATTATTTGGCACCTTATGCTTTAGCAAAAATTTGGTTAAAACCACAAGCTGAACAACAATTTTTATATGGAAATAATATTGTTAAATCTGGTGTTGGAAGAATGAGTGAAGGGATTGAAGAAAAACAA" seq = @gff.assemble(@contigsequence,component.start,[cds1],:phase => true) seq.size.should == 243 seq.should == "GAAAAATTAATGCGACAAGCAGCATGTATTGGACGTAAACAATTGGGATCTTTTGGAACTTGTTTGGGTAAATTCACAAAAGGAGGGTCTTTCTTTCTTCATATAACATCATTGGATTATTTGGCACCTTATGCTTTAGCAAAAATTTGGTTAAAACCACAAGCTGAACAACAATTTTTATATGGAAATAATATTGTTAAATCTGGTGTTGGAAGAATGAGTGAAGGGATTGAAGAAAAACAA"