gtc2vcf icon indicating copy to clipboard operation
gtc2vcf copied to clipboard

Failed to open file "Ou" : No such file or directory

Open crbeam opened this issue 1 year ago • 1 comments

Hello, freeseek,

I cannot seem to convert my .gtc files to a vcf file using the following code:

bpm_manifest_file="InfiniumOmni2-5-8v1-5_A1.bpm" csv_manifest_file="InfiniumOmni2-5-8v1-5_A1.csv" egt_cluster_file="InfiniumOmni2-5-8v1-5_A1_ClusterFile.egt" ref="$HOME/GRCh37/human_g1k_v37.fasta" out_prefix="batch1_vcf" bcftools +gtc2vcf --no-version -Ou --bpm $bpm_manifest_file --csv $csv_manifest_file --egt $egt_cluster_file --gtcs $path_to_gtc_folder --fasta-ref $ref --extra $out_prefix.tsv | bcftools sort -Ou -T ./bcftools-sort.XXXXXX | bcftools norm --no-version -Ob -c x -f $ref | tee $out_prefix.bcf | bcftools index --force --output $out_prefix.bcf.csi

The error that I receive is [E: :hts_open_format] Failed to open file "Ou" : No such file or directory Reading BPM file InfiniumOmni2-5-8v1-5_A1.bpm Could not read Ou Failed to read from standard input: unknown file type index: "-" is in a format that cannot be usefully indexed

I've tried adapting the command by reading through the other issues that have come up, but have had no luck creating a bcf file that has > 0 bytes. May I ask for assistance in resolving this issue? I should mention that the manifest and cluster files provided by illumina are in the same directory in which I am running this command.

Thank you, Chris

crbeam avatar Oct 19 '22 21:10 crbeam

Please disregard. A dash (-) before "Ou" in the bcftools sort command was missing. Thanks for a nice set of tools. Chris

crbeam avatar Oct 19 '22 23:10 crbeam