Hello, I am trying to filter low quality data off my mpileup files using Varscan but keep getting stuck with the error message: Invalid format for pileup at line 1. I don't think my vcf tools are corrupted, since they have made it through bedtools instersect.
The command lines I have used are the following:
samtools mpileup -uf genome.fa 1.bam 2.bam | bcftools view -bvcg -> *.bcf
I also used -vcg -> *.vcf and even ommitted -u. For VarScan I am typing:
java -jar VarScan.jar mpielup *.vcf --output-vcf
Any tip to as why is my vcf file not recognized by VarScan will be appreciated.