action='store_true',
default=False,
help="Add pseudocounts in the generated MM"
)
parser.add_argument('-k','--order',dest='order',action='store',type=int,
default=3,
help='It sets the order of the Markov Model used.')
args = parser.parse_args()
######## 1 TT SETS ######
training_filename='training_'+args.infile
testing_filename='testing_'+args.infile
m.training_testing_sets_separation(args.infile,3,training_filename,testing_filename)
######## 2 +/- ###########
(background, signal)=m.background_signal_separation(training_filename)
######## 3 kmer dict #####
if (args.pseudocounts):
back_dict=m.build_hash_pseudocount(background,args.order)
sign_dit=m.build_hash_pseudocount(signal,args.order)
else:
back_dict=m.build_hash(background,args.order)
sign_dit=m.build_hash(signal,args.order)
######### 4 print ########
m.print_hash(sign_dit,back_dict,args.outfile)
sys.exit()
