How to do batch processing with mpileup in a Linux shell for loop?
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
• 1,704 views
•
link
updated
by
Ram
•
written
by
elvissober •
•
written
by
elvissober •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
GATK SelectVariants behavior with gvcfs
written by Gianfilippo •
-
Set bam read group from existing tag
written by jbalberge
-
bcftools pileup call (multiple indels at same location)
written by microfuge
-
Set minimum variant frequency while calling variants with samtools
written by kspata
-
Variant Calling (Samtools Mpileup): Error
written by venuraherath •
-
Error in bcftools view -bvcg command
written by hafiz.talhamalik
-
How to use Ruby BioSamtools for SNP finding on .bam and .sam files?
written by elvissober •
-
Calling SNPs with 1X depth (samtools/bcftools)
written by Coryza •
-
Error in SNP calling, samtools
written by Adrian Pelin
-
Why Do We Need To Set Maximal Read Depth In Samtools
written by a.a.kechin •
Hello elvissober!
We believe that this post does not fit the main topic of this site.
Please look up bash loops. Do not open multiple questions addressing the same question.
For this reason we have closed your question. This allows us to keep the site focused on the topics that the community can help with.
If you disagree please tell us why in a reply below, we'll be happy to talk about it.
Cheers!