Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
Samtools-mpileup aslihan Bioinformatics 3 01-18-2017 10:53 PM
samtools mpileup: -B or not to -B? yenhuahuang1 Bioinformatics 4 03-19-2013 03:21 AM
Many Ns in samtools mpileup mikheyev Bioinformatics 5 03-22-2012 09:56 AM
samtools mpileup hv4 Bioinformatics 2 10-29-2011 01:21 PM
samtools mpileup bair Bioinformatics 2 03-03-2011 03:38 AM

Thread Tools
Old 12-06-2010, 05:26 PM   #1
Junior Member
Location: Vancouver

Join Date: Dec 2010
Posts: 1
Default samtools mpileup

Hi everyone,

When I try mpileup in samtools-0.1.12.
samtools-0.1.12/samtools mpileup -C50 -gf ref.fasta file.bam  > file.bcf
There is segfault:
[mpileup] 1 samples in 1 input files
<mpileup> Set max per-sample depth to 8000
Segmentation fault

Also, I tried
samtools-0.1.12/samtools mpileup -ugf file.fasta file.bam | samtools-0.1.12/bcftools/bcftools view -bvcg - > var.raw.bcf

samtools-0.1.12/bcftools/bcftools view var.raw.bcf | samtools-0.1.12/bcftools/ varFilter -D100 > var.flt.vcf

awk '{print $1}' var.flt.vcf | sort | uniq -c
The output after awk was:
12111 10
1640 11
1 ##fileformat=VCFv4.0

Can anyone give me some hints on what the problem is? Thank you!
Ling1 is offline   Reply With Quote
Old 05-24-2012, 10:09 AM   #2
Junior Member
Location: Cuernavaca, México

Join Date: Jul 2009
Posts: 6

Hi I have the same problem, did you found the solution?? any hint would help
yina is offline   Reply With Quote

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 08:58 PM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2019, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO