66
77import click
88
9- from bio2zarr import icf , vcf , verification
9+ from bio2zarr import vcf2zarr
1010
1111# TODO add support here for split vcfs. Perhaps simplest to take a
1212# directory provided as input as indicating this, and then having
@@ -40,11 +40,11 @@ def cli(vcfs, worker_processes, force):
4040 else :
4141 files = [f ]
4242 source_file = f
43- exploded = tmp_path / (f .name + ".exploded " )
43+ exploded = tmp_path / (f .name + ".icf " )
4444 if force and exploded .exists ():
4545 shutil .rmtree (exploded )
4646 if not exploded .exists ():
47- icf .explode (
47+ vcf2zarr .explode (
4848 exploded ,
4949 files ,
5050 worker_processes = worker_processes ,
@@ -53,21 +53,21 @@ def cli(vcfs, worker_processes, force):
5353 spec = tmp_path / (f .name + ".schema" )
5454 if force or not spec .exists ():
5555 with open (spec , "w" ) as specfile :
56- vcf .mkschema (exploded , specfile )
56+ vcf2zarr .mkschema (exploded , specfile )
5757
58- zarr = tmp_path / (f .name + ".zarr " )
58+ zarr = tmp_path / (f .name + ".vcz " )
5959 if force and zarr .exists ():
6060 shutil .rmtree (zarr )
6161 if not zarr .exists ():
62- vcf .encode (
62+ vcf2zarr .encode (
6363 exploded ,
6464 zarr ,
6565 spec ,
6666 worker_processes = worker_processes ,
6767 show_progress = True ,
6868 )
6969
70- verification . validate (source_file , zarr , show_progress = True )
70+ vcf2zarr . verify (source_file , zarr , show_progress = True )
7171
7272
7373if __name__ == "__main__" :
0 commit comments