PageRenderTime 25ms CodeModel.GetById 18ms app.highlight 4ms RepoModel.GetById 1ms app.codeStats 0ms

/matlab/get_training_sequences_intergenic_long.m

http://github.com/jcnossen/hmmgenefinder
Objective C | 14 lines | 13 code | 1 blank | 0 comment | 2 complexity | 63f980838867caf8b1e1b9061cb618b9 MD5 | raw file
Possible License(s): GPL-2.0
 1% Returns sequences for training HMM (used to export data for GHMM)
 2% Input:
 3%       seq - sequence with annotated genes
 4% Output:
 5%       res - cell array with intergnenic regions
 6function [res] = get_training_sequences_intergenic_long(seq)
 7    % Just in case, but better do it when loading the genome
 8    seq.Sequence = upper(seq.Sequence);
 9    
10    seq = HMM_Intergenic.get_intergenic(seq, true, true);
11    seq = HMM_Intergenic.select_intergenic(seq, 10, Inf);
12    seq = HMM.select_genes(seq, true, false, false);
13    res = HMM.get_training_sequences(seq);
14end