Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Getting output to standard out (instead of writing to a VCF) #95

Open
pjshort opened this issue Jun 19, 2018 · 2 comments
Open

Getting output to standard out (instead of writing to a VCF) #95

pjshort opened this issue Jun 19, 2018 · 2 comments

Comments

@pjshort
Copy link

pjshort commented Jun 19, 2018

I am trying to annotate a VCF file with deepSEA results and would like to have the output VCF go to standard out, so I can pipe straight to another script.

I tried the following:

kipoi postproc score_variants DeepSEA/variantEffects --dataloader_args '{"fasta_file": "Homo_sapiens.GRCh37.74.dna.chr.fasta"}' -i min.vcf

using the following simple vcf:

##fileformat=VCFv4.1
##contig=<ID=1,length=249250621,assembly=b37>
##reference=file:Homo_sapiens_assembly19.fasta
#CHROM POS ID REF ALT
1 10002 - T A
1 10002 - T C
1 10002 - T G

This results in the following error:

Traceback (most recent call last):
File "/Users/ps14/anaconda3/envs/kipoi-DeepSEA__variantEffects/bin/kipoi", line 11, in
sys.exit(main())
File "/Users/ps14/anaconda3/envs/kipoi-DeepSEA__variantEffects/lib/python3.5/site-packages/kipoi/main.py", line 74, in main
command_fn(args.command, sys.argv[2:])
File "/Users/ps14/anaconda3/envs/kipoi-DeepSEA__variantEffects/lib/python3.5/site-packages/kipoi/cli/postproc.py", line 658, in cli_main
command_fn(args.command, raw_args[1:])
File "/Users/ps14/anaconda3/envs/kipoi-DeepSEA__variantEffects/lib/python3.5/site-packages/kipoi/cli/postproc.py", line 118, in cli_score_variants
dir_exists(os.path.dirname(args.output_vcf), logger)
File "/Users/ps14/anaconda3/envs/kipoi-DeepSEA__variantEffects/lib/python3.5/posixpath.py", line 148, in dirname
i = p.rfind(sep) + 1
AttributeError: 'NoneType' object has no attribute 'rfind'

If you could help me figure out how to print to standard out I would really appreciate it! Thank you!

@pjshort pjshort changed the title Getting output to standard out (instead of saving to a VCF) Getting output to standard out (instead of writing to a VCF) Jun 19, 2018
@Avsecz
Copy link
Contributor

Avsecz commented Jun 19, 2018

Good idea. We can add this functionality. Should be fairly easy.

@pjshort
Copy link
Author

pjshort commented Jun 19, 2018

Thank you, @Avsecz!

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

No branches or pull requests

3 participants