processor ID: 23388
==========================================================================

|------------------------------------------------------------------|
|                                                                  |
|              *** Running a seeded analysis ***                   |
|                                                                  |
|------------------------------------------------------------------|

command line: /home/grobertson/Fi/Code/20090926_v1.3/Fi_v1.3/bin/fi -fseq /projects/remc_bigdata/Karsan/HS1235_wt/compare/FA/HS1235-compare-lowP_10518-selfUnion-pm500bp.hg18.fa.3stage_dm_w4.fa -fout /projects/remc_bigdata/Karsan/HS1235_wt/compare/seed_RBPJ_MEME-m1/HS1235-compare-lowP_10518-selfUnion-pm500bp.hg18.fa.3stage_dm_w4.fa.pwmCSL-MEME-m1.mx.weight0.de_novo -minN 3000 -ev 5000 -pv 0.0002 -fpwm0 /home/grobertson/Fi/PWMs/CSL-MEME-m1.mx -em 80 -fEM 0.5 -posWt 0 -extTrim EXTRIM -verbose 1 -pgf 0 -fbm /home/grobertson/Fi/Code/20090926_v1.3/KmerFreq/hg18_kmer_1to9_freq.txt -bOrder 3 

maximal buffer length:					15000
maximal number of sequences set:			100000
maximal number of bases per seq read:			20000
maximal number of sites in a motif:			150000
input (ChIP) sequence file:				/projects/remc_bigdata/Karsan/HS1235_wt/compare/FA/HS1235-compare-lowP_10518-selfUnion-pm500bp.hg18.fa.3stage_dm_w4.fa
number of sequences in input file:			10518
average sequence length:				1082
total number of nucleotides:				11387432
max number of generations:				1
population size:					10
use a user-specified pwm as the seed			/home/grobertson/Fi/PWMs/CSL-MEME-m1.mx
fraction (number) input sequences subject to EM		1.00 (10518)
scale factor for converting (double)pwm to (int)pwm	 200
number of EM steps:					80
EM convergence criterion:				1.000000e-04
run EM on the starting pwm /home/grobertson/Fi/PWMs/CSL-MEME-m1.mx 10 times, each with a different maxp:
0.10*numSeq 0.20*numSeq 0.30*numSeq 0.40*numSeq 0.50*numSeq 0.60*numSeq 0.70*numSeq 0.80*numSeq 0.90*numSeq 1.00*numSeq 
no spaced dyads are generated and used. pop=10 gen=1 (no GA).

motif prior probability type (see documentation):	0
pwm score p-value cutoff for declaring binding site:	2.000000e-04

Approximate the null llr log{p(s|M)/p(s|B)} score distribution using
the llr scores of random/background sequences, where M is the EM-derived
motif model and B is the 3-th order Markov backgroun model.

The background sequences are simulated using the [a,c,g,t] frequencies
in the input data. The number sets of background sequences generated: 	10

pseudo count:						0.0005
minimal infomation for trimming/extending: 		0.40 0.50 0.60
minimal no. sites for each motif:			3000
base extension and trimming? 				no
sliding window for comparing pwm similarity:	6
PWM similarity cutoff:					0.300
log(E-value) cutoff:					5000.00
number of adjacent bases included in binding site output:	10

job started: Sat Jul 10 21:34:28 2010

=========================================================================

GADEM cycle[  1] generation[  1] number of unique motif(s): 2
   spacedDyad: yGTGGGAA             motifConsensus: TTTGGrAG              0.60 fitness: 1298.56
   spacedDyad: yGTGGGAA             motifConsensus: TTTTTAAA              1.00 fitness: 1881.39

GADEM cycle[  2] generation[  1] number of unique motif(s): 1
   spacedDyad: yGTGGGAA             motifConsensus: AATGGAAT              0.30 fitness:  669.71

GADEM cycle[  3] generation[  1] number of unique motif(s): 1
   spacedDyad: yGTGGGAA             motifConsensus: AAAkAAAA              0.90 fitness: -2989.26

GADEM cycle[  4] generation[  1] number of unique motif(s): 1
   spacedDyad: yGTGGGAA             motifConsensus: AATrmAAA              0.90 fitness: 3438.31

GADEM cycle[  5] generation[  1] number of unique motif(s): 2
   spacedDyad: yGTGGGAA             motifConsensus: AATrkAAA              0.50 fitness: 3146.19
   spacedDyad: yGTGGGAA             motifConsensus: AArAkAAA              0.80 fitness: 3960.62

GADEM cycle[  6] generation[  1] number of unique motif(s): 1
   spacedDyad: yGTGGGAA             motifConsensus: AAkGkrAA              0.40 fitness: 3559.32


finished: Sun Jul 11 05:53:40 2010

approximated processor time in seconds: 29952.000000