Mercurial > repos > vipints > rdiff

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

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
version 0.3 uploaded
author vipints
date Thu, 14 Feb 2013 23:38:36 -0500
parents
children
line wrap: on
line source

function [genes]=compute_testing_region(CFG,genes)

%Iterate over gene

for i=1:size(genes,2)
    [SPLICINGEVENTS,SEQUENCE,EXONSEQUENCE]=splicingsequence(genes(i));
    genes(i).splicingevents=SPLICINGEVENTS;
    genes(i).sequence=SEQUENCE;
    genes(i).exonsequence=EXONSEQUENCE;
    
    [UNIQUE_NEW_EXONS,GRAPHNODES,ORDER_OF_GRAPHNODE,EIRS_IN_SEQ]=transform_single_end_reads(genes(i),CFG.sequenced_length-CFG.bases_to_clip*2);
    genes(i).unique_new_exons=UNIQUE_NEW_EXONS;
    genes(i).graphnodes=GRAPHNODES;
    genes(i).order_of_graphnodes=ORDER_OF_GRAPHNODE;
    genes(i).eirs_in_seq=EIRS_IN_SEQ;
 
end