Entering edit mode
4.6 years ago
cocchi.e89
▴
290
It may be silly but I am not getting how to pipe these 2 bcftools
commands where I need the file generated by the first command to be an atgument of the second one:
bcftools norm -m -any case.chr21.vcf.gz -Ou | bcftools isec -n~10 case.chr21.altered.vcf.gz
the plain command would be:
bcftools isec -n~10 case.chr21.vcf.gz case.chr21.altered.vcf.gz
I also tried a shell version:
bcftools isec -n~10 <(bcftools norm -m -any case.chr21.vcf.gz -Ou) case.chr21.altered.vcf.gz
but it doesn't work
Thank you so much for any help
May I know how to integrate
bcftools index
in the pipe? I encountered similar issue where I intend to usebcftools index
as the second command with inputs from the first command, something like|bcftools index
but the|bcftools index
command didn't run. Thank you!