533
edits
Changes
From Genome Analysis Wiki
Jump to navigationJump to search→Initial set of variant calls
samtools mpileup -Iuf ref/human_g1k_v37_chr20.fa bams/SAMPLE*bam | bcftools view -bvcg - > result/chr20.mpileup.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/result/chr20.mpileup.bcf | sed s/AF1/AF/g | grep -v AF=1 > result/chr20.mpileup.vcf
