INFO  @ 22 Mar 2026 18:16:10: [35 MB] # read alignment files... 
INFO  @ 22 Mar 2026 18:16:10: [35 MB] # read treatment tags... 
INFO  @ 22 Mar 2026 18:16:10: [55 MB] Detected format is: BAM 
INFO  @ 22 Mar 2026 18:16:10: [55 MB] * Input file is gzipped. 
INFO  @ 22 Mar 2026 18:16:13: [66 MB]  1000000 reads parsed 
INFO  @ 22 Mar 2026 18:16:16: [73 MB]  2000000 reads parsed 
INFO  @ 22 Mar 2026 18:16:19: [74 MB]  3000000 reads parsed 
INFO  @ 22 Mar 2026 18:16:22: [81 MB]  4000000 reads parsed 
INFO  @ 22 Mar 2026 18:16:25: [81 MB]  5000000 reads parsed 
INFO  @ 22 Mar 2026 18:16:28: [94 MB]  6000000 reads parsed 
INFO  @ 22 Mar 2026 18:16:32: [94 MB]  7000000 reads parsed 
INFO  @ 22 Mar 2026 18:16:35: [95 MB]  8000000 reads parsed 
INFO  @ 22 Mar 2026 18:16:38: [98 MB]  9000000 reads parsed 
INFO  @ 22 Mar 2026 18:16:41: [107 MB] 9959723 reads have been read. 
INFO  @ 22 Mar 2026 18:16:41: [107 MB] tag size is determined as 100 bps 
INFO  @ 22 Mar 2026 18:16:41: [107 MB] # tag size = 100 
INFO  @ 22 Mar 2026 18:16:41: [107 MB] # total tags in alignment file: 9959723 
INFO  @ 22 Mar 2026 18:16:41: [107 MB] # Build Peak Model... 
INFO  @ 22 Mar 2026 18:16:41: [107 MB] #2 looking for paired plus/minus strand peaks... 
INFO  @ 22 Mar 2026 18:16:43: [108 MB] #2 Total number of paired peaks: 125 
INFO  @ 22 Mar 2026 18:16:44: [108 MB] #2 Model building with cross-correlation: Done 
INFO  @ 22 Mar 2026 18:16:44: [108 MB] # finished! 
INFO  @ 22 Mar 2026 18:16:44: [108 MB] # predicted fragment length is 57 bps 
INFO  @ 22 Mar 2026 18:16:44: [108 MB] # alternative fragment length(s) may be 57,76,109,156,201,258,290,318,357,385,418,456,503,506,514,555 bps 
INFO  @ 22 Mar 2026 18:16:44: [108 MB] # Generate R script for model : /home/jforment/biovice/internal_projects/260320_mblazquez_chipseq/02-mapping/DRIF3.predictd/predictd_model.R