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 = []
for i in range(nbams):
reader = Process(target=read_pairs_from_bam,
args=(inq, i, hra.bams[i], hra, args.mapq,
args.slice),
daemon=False)
reader.start()
readers.append(reader)
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=[]
for i in range(nbams):
reader = Process(target=read_pairs_from_bam, args=(inq,i,hra.bams[i],hra,args.mapq,args.slice),daemon=False)
reader.start()
readers.append(reader)
q = JoinableQueue(maxsize=0)
histogram_queue = JoinableQueue(maxsize=0)
final_histogram_queue = JoinableQueue(maxsize=0)
