Entering edit mode
8.3 years ago
elvissober
▴
20
I want to make VCFs on a set of .sams (.bams):
for I in ls (*.sam) # calling variants on .sam files
do
./samtools mpileup -R -d | bcftools view -g - > test.vcf $i # what to write here correctly?
done