comparison read2mut.py @ 81:612c110305db draft

planemo upload for repository https://github.com/Single-Molecule-Genetics/VariantAnalyzerGalaxy/tree/master/tools/variant_analyzer commit ee4a8e6cf290e6c8a4d55f9cd2839d60ab3b11c8-dirty
author mheinzl
date Wed, 27 Jul 2022 07:08:10 +0000
parents 8336a4f2b647
children c2e8932b4d8d
comparison
equal deleted inserted replaced
80:8336a4f2b647 81:612c110305db
425 counts_mut = 0 425 counts_mut = 0
426 chimeric_tag_list = [] 426 chimeric_tag_list = []
427 chimeric_tag = {} 427 chimeric_tag = {}
428 if (key1 in pure_tags_dict_short.keys()) or (key1 in pure_tags_dict_ref.keys()): # ref or alt 428 if (key1 in pure_tags_dict_short.keys()) or (key1 in pure_tags_dict_ref.keys()): # ref or alt
429 429
430 if key1 not in np.array(['#'.join(str(i) for i in z) 430 # if key1 not in np.array(['#'.join(str(i) for i in z)
431 for z in zip(mut_array[:, 0], mut_array[:, 1], mut_array[:, 2], mut_array[:, 3])]): 431 # for z in zip(mut_array[:, 0], mut_array[:, 1], mut_array[:, 2], mut_array[:, 3])]):
432 continue 432 # continue
433 433
434 change_tier_after_print = [] 434 change_tier_after_print = []
435 i = np.where(np.array(['#'.join(str(i) for i in z) 435 i = np.where(np.array(['#'.join(str(i) for i in z)
436 for z in zip(mut_array[:, 0], mut_array[:, 1], mut_array[:, 2], mut_array[:, 3])]) == key1)[0][0] 436 for z in zip(mut_array[:, 0], mut_array[:, 1], mut_array[:, 2], mut_array[:, 3])]) == key1)[0][0]
437 ref = mut_array[i, 2] 437 ref = mut_array[i, 2]