args = parser.parse_args()
if args.infile:
asf = HiriseAssembly()
asf.load_assembly(args.infile)
# asf = HiriseAssembly()
# asf.load_playout(args.infile)
breaks=[]
scores={}
for l in open(args.breaks):
if l[0]=="#": continue
#Scaffold102239 741 1379 5097 -5.421961663655971
scaffold,a,b,slen,score = l.strip().split()
a=int(a) #start
b=int(b) #end
# c=int(c) #lowpoint
score=float(score)
c=int((a+b)/2)
breaks.append((scaffold,a,b,c))
scores[scaffold,a,b,c]=score
# asf.add_breakpoints_ranges(breaks,debug=args.debug,scores=scores)
asf.add_breakpoints_ranges(breaks,debug=args.debug)
asf.validate()
if args.outfile: asf.save_assembly( args.outfile )
