Use keep_original in pysam.tabix_index to not remove original data

This commit is contained in:
mvdbeek
2021-04-29 22:37:41 +02:00
parent ede2f0f6a8
commit 6a209388ee
+1 -1
View File
@@ -823,7 +823,7 @@ class VcfGz(BaseVcf, binary.Binary):
index_file = dataset.metadata.spec['tabix_index'].param.new_file(dataset=dataset)
try:
pysam.tabix_index(dataset.file_name, index=index_file.file_name, preset='vcf', force=True)
pysam.tabix_index(dataset.file_name, index=index_file.file_name, preset='vcf', keep_original=True, force=True)
except Exception as e:
raise Exception('Error setting VCF.gz metadata: %s' % (util.unicodify(e)))
dataset.metadata.tabix_index = index_file