comparison rDiff/src/tools/compute_testing_region.m @ 0:0f80a5141704

version 0.3 uploaded
author vipints
date Thu, 14 Feb 2013 23:38:36 -0500
parents
children
comparison
equal deleted inserted replaced
-1:000000000000 0:0f80a5141704
1 function [genes]=compute_testing_region(CFG,genes)
2
3 %Iterate over gene
4
5 for i=1:size(genes,2)
6 [SPLICINGEVENTS,SEQUENCE,EXONSEQUENCE]=splicingsequence(genes(i));
7 genes(i).splicingevents=SPLICINGEVENTS;
8 genes(i).sequence=SEQUENCE;
9 genes(i).exonsequence=EXONSEQUENCE;
10
11 [UNIQUE_NEW_EXONS,GRAPHNODES,ORDER_OF_GRAPHNODE,EIRS_IN_SEQ]=transform_single_end_reads(genes(i),CFG.sequenced_length-CFG.bases_to_clip*2);
12 genes(i).unique_new_exons=UNIQUE_NEW_EXONS;
13 genes(i).graphnodes=GRAPHNODES;
14 genes(i).order_of_graphnodes=ORDER_OF_GRAPHNODE;
15 genes(i).eirs_in_seq=EIRS_IN_SEQ;
16
17 end