Entering edit mode
                    5.5 years ago
        cocchi.e89
        
    
        ▴
    
    300
    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 indexin the pipe? I encountered similar issue where I intend to usebcftools indexas the second command with inputs from the first command, something like|bcftools indexbut the|bcftools indexcommand didn't run. Thank you!