def __init__(self, options): # Openning tabix file representing the reference genome self.fastafile = pysam.Fastafile(options.args['reference']) # Retrieving the sequence of a genomic region