SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
ERROR: [bcf_sync] incorrect number of fields.. vinodhsri Bioinformatics 21 08-11-2017 01:42 PM
GATK multisample vcf to BayeScan input ndeshpan Bioinformatics 2 03-24-2017 07:34 AM
Question: bcftools: Wrong number of PL fields? nals=4 npl=6 schumric Illumina/Solexa 3 01-09-2017 10:37 AM
VCF file / relationship between QUAL and DP fields cambyses Bioinformatics 1 11-26-2014 10:22 AM
Filtering multisample vcf files on DP DavyK Bioinformatics 3 01-17-2013 09:39 AM

Reply
 
Thread Tools
Old 04-04-2018, 08:46 AM   #1
Greg
Member
 
Location: British Columbia

Join Date: Oct 2009
Posts: 31
Default bbmap callvariants.sh multisample vcf format, wrong number of fields

Trying to use BBMap callvariants on multiple samples. The resulting VCF does not have the sample names and yields an error with any vcf tools:

Code:
vcf-validator tmp.vcf 
Wrong number of fieldsin tmp.vcf; expected 15, got 16. The offending line was: [... the first line of variants]
The VCF:
Code:
#CHROM  POS     ID      REF     ALT     QUAL    FILTER  INFO    FORMAT  FORMAT  FORMAT  FORMAT  FORMAT  FORMAT  FORMAT
BBMap command:
Code:
callvariants.sh multisample=t list=SamList out=tmp.vcf ref=genome.fa 
#ive also tried using the output command, no change
callvariants.sh multisample=t list=SamList out=tmp.vcf ref=genome.fa vcf0=S1,S2,S3,S4,S5,S5
#As well as each sample using the "in" command
callvariants.sh multisample=t in=S1.sam,S2.sam,S3.sam,S4.sam,S5.sam,S6.sam out=tmp.vcf ref=genome.fa overwrite=t vcf0=AO,GG,OG,SF,TH,Tum
I think BBMap could be a very powerful set of tools and I am expecting I am missing something pretty basic!
Greg is offline   Reply With Quote
Old 04-04-2018, 09:10 AM   #2
Brian Bushnell
Super Moderator
 
Location: Walnut Creek, CA

Join Date: Jan 2014
Posts: 2,707
Default

Hi Greg,

Sorry about that, there was a bug that slipped in around v37.87 with regards to multisample VCF names. It's fixed in v37.96 which I will release this week.

-Brian
Brian Bushnell is offline   Reply With Quote
Old 04-04-2018, 09:27 AM   #3
Greg
Member
 
Location: British Columbia

Join Date: Oct 2009
Posts: 31
Default

Okay cool, thanks for the quick response!
Greg is offline   Reply With Quote
Old 04-24-2018, 04:45 AM   #4
denisDDS
Junior Member
 
Location: France

Join Date: Feb 2013
Posts: 6
Default

Same error here: downloaded BBMap 2 weeks ago, finishing my snp calling pipeline today, same problem and discovering the patch went out 10 days ago...

Thx for the fix
denisDDS is offline   Reply With Quote
Old 04-24-2018, 08:42 AM   #5
Brian Bushnell
Super Moderator
 
Location: Walnut Creek, CA

Join Date: Jan 2014
Posts: 2,707
Default

Oh, yep - it's fixed now, sorry for not mentioning that!
Brian Bushnell is offline   Reply With Quote
Old 04-25-2018, 12:00 AM   #6
denisDDS
Junior Member
 
Location: France

Join Date: Feb 2013
Posts: 6
Default

I have another problem: when using vcf0=, followed by files locations separated by comma, callvariants.sh create subdirectory instead of just creating files at specified location.
denisDDS is offline   Reply With Quote
Reply

Tags
bbmap, vcf

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 06:22 PM.


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