parser.add_argument('-o', '--outfile', help="Output file name.") parser.add_argument('-H', '--histogram_file', help="Save score histograms to FILE") parser.add_argument('-j', '--nthreads', help="Number of threads.", type=int, default=16) args = parser.parse_args() if args.top: import psutil hra = HiriseAssembly() hra.load_assembly(args.infile) hra.merge_masked_regions(debug=args.debug) print(len(hra.layout_lines)) if len(hra.layout_lines) == 0: print("#make trivial layout") hra.make_trivial_layout_lines(debug=args.debug) ces.set_exp_insert_size_dist_fit_params(hra.model_params) model = ces.model nbams = len(hra.bams) inq = JoinableQueue(maxsize=0) readers = []
from __future__ import division from __future__ import print_function from builtins import range from past.utils import old_div from hirise_assembly import HiriseAssembly if __name__=="__main__": import sys import argparse parser = argparse.ArgumentParser(formatter_class=argparse.ArgumentDefaultsHelpFormatter) parser.add_argument('-d','--debug',default=False ,action="store_true",help="Turn on debugging ouput") parser.add_argument('-L','--layout',default=False ,help="A file containing a layout of contigs.") parser.add_argument('-i','--infile',default=False ,help="Filename for serialised assembly input file.") parser.add_argument('-o','--outfile',default=False,help="Filename for writing a list of segments on the raw contigs to mask for being promiscuous in linking.") args = parser.parse_args() if args.infile: asf = HiriseAssembly() asf.load_assembly(args.infile) if args.outfile: f=open(args.outfile,"wt") for contig in asf.contigs_iter(): f.write("{}\t{}\n".format(contig,asf.contig_length(contig)))