From Genome Analysis Wiki
Jump to navigationJump to search
7 bytes removed
, 10:12, 10 February 2012
Line 120: |
Line 120: |
| | | |
| mkdir result | | mkdir result |
− | samtools mpileup -r 20:1500000-2500000 -Iuf ref/human_g1k_v37_chr20.fa bams/*bam | bcftools view -bvcg - > result/chr20.bcf
| + | samtools mpileup -Iuf ref/human_g1k_v37_chr20.fa bams/SAMPLE*bam | bcftools view -bvcg - > result/chr20.mpileup.bcf |
| | | |
| bcftools view result/chr20.bcf | sed s/AF1/AF/g | grep -v AF=1 > result/chr20.mpileup.vcf | | bcftools view result/chr20.bcf | sed s/AF1/AF/g | grep -v AF=1 > result/chr20.mpileup.vcf |