def OpenVCFfileAsArray(fpath): VariantsArr = allel.vcf_to_recarray(fpath) return VariantsArr
import binascii import sys import timeit import os from Crypto.PublicKey import RSA from hashlib import sha512 from hashlib import sha256 import hmac from fastpbkdf2 import pbkdf2_hmac mt = MerkleTools() RSAbitsNum = 2048 hash_type = "sha256" ra = allel.vcf_to_recarray('../files/NB72462M.vcf') Iter = int(sys.argv[1]) ResiltFileName = "MHTexperimentResults1IterNov11fast128Iter" + str(Iter) ExperimentResults = open(ResiltFileName, "w") RootSigFile = "root_sig" SaltsFileName = "salts" TreeName = "mtlevels" salts = [] lvs = [] for x in ra: lvs.append(str(x)) hash_function = getattr(hashlib, hash_type) def GenerateMHT(leavess):