Changes between Version 18 and Version 19 of SOPs/variant_calling

Timestamp:

07/14/17 12:14:04 (8 years ago)

Author:

gbell

Comment:

--

SOPs/variant_calling

@@ -61 +61 @@
 In our experience, "-B" (disable BAQ) or "-E" (recalculate BAQ) works better than the default method, which can remove some obvious variants.
 {{{
-samtools mpileup -d100000 -t DP,AD,DV,ADF,ADR,INFO/AD,SP -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
+samtools mpileup -d100000 -t DP,AD,ADF,ADR,INFO/AD,SP -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
 #in bcftools call: -m for multi-allelic, -c for bi-alleleic
 }}}
 Call variants (multiple sample vs. reference) using a set of BAM files
 {{{
-samtools mpileup -d100000 -t DP,AD,DV,ADF,ADR,INFO/AD,SP -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
+samtools mpileup -d100000 -t DP,AD,ADF,ADR,INFO/AD,SP -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
 }}}