Changes between Version 10 and Version 11 of SOPs/variant_calling

Timestamp:

02/13/15 14:10:52 (10 years ago)

Author:

gbell

Comment:

--

SOPs/variant_calling

@@ -56 +56 @@
 In our experience, "-B" is needed to disable BAQ computation, which would otherwise remove many obvious variants.
 {{{
-samtools mpileup -d100000 -Buf /nfs/genomes/sgd_2010/bwa/sacCer3.fa A_reads.bt2.sorted_unique.bam | bcftools view -bvcg - >| A_reads.bt2.sorted_unique.raw.bcf
+samtools mpileup -d100000 -Buf /nfs/genomes/sgd_2010/bwa/sacCer3.fa A_reads.bt2.sorted_unique.bam | bcftools call -O b -v -c - >| A_reads.bt2.sorted_unique.raw.bcf
 }}}
 Call variants (multiple sample vs. reference) using a set of BAM files
 {{{
-samtools mpileup -d100000 -Buf /nfs/genomes/sgd_2010/bwa/sacCer3.fa *_reads.bt2.sorted_unique.bam | bcftools view -bvcg - >| ALL_reads.bt2.sorted_unique.raw.bcf
+samtools mpileup -d100000 -Buf /nfs/genomes/sgd_2010/bwa/sacCer3.fa *_reads.bt2.sorted_unique.bam | bcftools call -O b -v -c - >| ALL_reads.bt2.sorted_unique.raw.bcf
 }}}