Question: Error in vcf, samtools
0
gravatar for wu.zhiqiang.1020
4.7 years ago by
United States
wu.zhiqiang.102020 wrote:

../software/samtools-1.1/samtools mpileup -uf ./S.int.fas ./S_new_F4_pe_mem.sort.bam | ../software/bcftools-1.1/bcftools view -vcg - > S_new_F4_pe_mem.sort.bam.raw.bcf 

[mpileup] 1 samples in 1 input files

<mpileup> Set max per-file depth to 8000

[E::cg] unknown typ

What's the reason? please help.

I used the same version of samtools and bcftools

 

 

snp • 1.7k views
ADD COMMENTlink modified 3.7 years ago by Biostar ♦♦ 20 • written 4.7 years ago by wu.zhiqiang.102020
0
gravatar for Pierre Lindenbaum
4.7 years ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum120k wrote:

you don't need to use bcftools view anymore: http://www.htslib.org/workflow/#mapping_to_variant

samtools mpileup -ugf <ref.fa> <sample1.bam> <sample2.bam> <sample3.bam> | bcftools call -vmO z -o <study.vcf.gz>
ADD COMMENTlink written 4.7 years ago by Pierre Lindenbaum120k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1396 users visited in the last hour