SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
BCF Tools JohnK Bioinformatics 4 11-28-2013 07:27 AM
Dindel last stage error nothing in the var.VCF file pankajsvats Bioinformatics 3 02-02-2012 06:00 AM
DESeq estimates of mu,var bruce01 Bioinformatics 0 01-10-2012 06:53 AM
Clarification on BCF file format phatjoe Bioinformatics 2 01-09-2012 05:22 PM
SAMTOOLS: bcf-fix.pl problem michalkovac Bioinformatics 0 12-06-2010 03:54 AM

Reply
 
Thread Tools
Old 05-19-2011, 07:23 PM   #1
dongshenglulv
Member
 
Location: Shanghai

Join Date: May 2011
Posts: 15
Default *.bam -> var.bcf No result

Hi all,

When I tried to call variants by 25 BAM files like this,
$ samtools mpileup -P ILLUMINA -ugf hg19.fasta *.bam | bcftools view -bcvg - > var.raw.bcf ; bcftools view var.raw.bcf | vcfutils.pl varFilter -D 2000 > var.flt.vcf

Then the standard output was
[mpileup] 25 samples in 25 input files
<mpileup> Set max per-file depth to 320
[afs] 0:0.000 1:0.000 2:0.000 3:0.000 4:0.000 5:0.000 6:0.000 7:0.000 8:0.000 9:0.000 10:0.000 11:0.000 12:0.000 13:0.000 14:0.000 15:0.000 16:0.000 17:0.000 18:0.000 19:0.000 20:0.000 21:0.000 22:0.000 23:0.000 24:0.000 25:0.000 26:0.000 27:0.000 28:0.000 29:0.000 30:0.000 31:0.000 32:0.000 33:0.000 34:0.000 35:0.000 36:0.000 37:0.000 38:0.000 39:0.000 40:0.000 41:0.000 42:0.000 43:0.000 44:0.000 45:0.000 46:0.000 47:0.000 48:0.000 49:0.000 50:0.000

and the program just stopped, the VCF file just only had the header information. I don't why, could you help me?

thanks
dongshenglulv is offline   Reply With Quote
Old 05-22-2011, 07:24 AM   #2
DZhang
Senior Member
 
Location: East Coast, US

Join Date: Jun 2010
Posts: 177
Default

Hi,

According to my experience, this usually means that your parameters for SNP calling is not optimal so there is NO SNP called. You may loosen the criteria to try to get some callings, and then tighten them up.

Douglas
www.contigexpress.com
DZhang is offline   Reply With Quote
Old 05-22-2011, 07:40 AM   #3
dongshenglulv
Member
 
Location: Shanghai

Join Date: May 2011
Posts: 15
Default

The only parameter for criteria here is -D. I just tried this command,
$ samtools mpileup -P ILLUMINA -ugf hg19.fasta *.bam | bcftools view -cvg - > var.raw.vcf

Then the standard output was the same, and the file var.raw.vcf just had the header information without the body.
Is there any suggestion you can give to me? Thanks very much
dongshenglulv is offline   Reply With Quote
Old 05-22-2011, 07:45 AM   #4
DZhang
Senior Member
 
Location: East Coast, US

Join Date: Jun 2010
Posts: 177
Default

Hi,

I am not entirely sure but I recall there is a step from bcf to vcf. Please check samtools manual, especially the part related to bcftools.

Douglas
www.contigexpress.com
DZhang is offline   Reply With Quote
Old 08-04-2011, 11:24 AM   #5
Hkins552
Member
 
Location: Baltimore MD

Join Date: Jun 2011
Posts: 18
Default

Have you resolved this problem?
Hkins552 is offline   Reply With Quote
Old 08-04-2011, 01:02 PM   #6
swbarnes2
Senior Member
 
Location: San Diego

Join Date: May 2008
Posts: 912
Default

Try running just mpileup to generate a pileup. Does that look normal?

I'm not sure that all those 0's are normal in the output to the terminal.
swbarnes2 is offline   Reply With Quote
Old 09-01-2011, 02:20 PM   #7
Lynette
Junior Member
 
Location: Canada

Join Date: Apr 2008
Posts: 2
Default

I had the same problem as you (had no output from mpileup except headers) from solid paired end data and when I added the -A option I got results!
Lynette is offline   Reply With Quote
Old 09-01-2011, 03:15 PM   #8
swbarnes2
Senior Member
 
Location: San Diego

Join Date: May 2008
Posts: 912
Default

Quote:
Originally Posted by Lynette View Post
I had the same problem as you (had no output from mpileup except headers) from solid paired end data and when I added the -A option I got results!
Hmm, real SNPs should be visible without relying solely on anomalous read pairs.

The other thing to try is run mpileup with -B. That will disengage the BAQ calculations. It's possible that you will get more false positives, but I've seen the BAQ calculations hide true positives.
swbarnes2 is offline   Reply With Quote
Old 09-02-2011, 07:03 AM   #9
Lynette
Junior Member
 
Location: Canada

Join Date: Apr 2008
Posts: 2
Default

Bit of background: I'm using solid 5500 data, aligning with bfast, picard remove duplicates, srma (for local-realignment), picard fixmate, and then gatk q-score recalibration. I've been used the -B option and still got no output until I used the -A option

I'm not too if one of those steps is messing up the mate information so I tried picard fixmate after gatk q-score recalibration and still no results from mpileup until I used the -A option. It's very weird when I look at the picard insert size metrics of this sample I have like half (wrong) RF orientation and half (correct) FR orientation.
Lynette is offline   Reply With Quote
Reply

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 10:53 PM.


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