Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • ERROR: [bcf_sync] incorrect number of fields..

    Has anyone encountered this error ('[bcf_sync] incorrect number of fields (0 != 5)') and found a solution already:


    vcfutils.pl splitchr -l 50000000 ucsc.hg19.fasta.fai | xargs -i echo samtools mpileup -I -C50 -m3 -F0.0002 -DSuf ucsc.hg19.fasta -r {} -b bam.list | bcftools view -bcvg - /> part-{}.bcf

    [bcf_sync] incorrect number of fields (0 != 5) at 0:0
    [afs] 0:0.000
    xargs: echo: terminated by signal 13

    Appreciate your advice.

    vsri

  • #2
    I get the same error message with a command like:

    nohup samtools mpileup -uf ref.fa s1.bam s2.bam

    but it runs fine without the 'nohup'. Have you tried the command sequence without the pipes?

    Comment


    • #3
      I ran into this problem as well, using mpileup/bcftools (samtools-0.1.18) to call variants from 19 bam files.
      I willl try to run it wiiltout nohup.

      Comment


      • #4
        Still have this problem when doing bcftools view for a bcf file (mpileup generated):[bcf_sync] incorrect number of fields (6 != 5) at 7:1396330564
        Any help?

        Comment


        • #5
          This is weird. I've used Samtools a lot before on bacterial references and never had this problem, yet on some human exomes it happens. Mpileup runs fine and then I repeatedly get

          [bcf_sync] incorrect number of fields (0 != 5) at 0:0

          Reinstalling/updating samtools and bcftools doesn't seem to help.

          Comment


          • #6
            usually this happens if your bcf file is truncated.

            Comment


            • #7
              Actually, I think it wasn't truncated, but something added by nohup.

              When I don't use nohup to run the mpileup command bcftools seems to work fine.

              Comment


              • #8
                I still got the error message without nohup: "[bcf_sync] incorrect number of fields (6 != 5) at 7:1396330564"

                Comment


                • #9
                  I just had the same error and by adding '-u' option to 'samtools mplieup' part of the pipeline , this solved the issue in my case.
                  Last edited by salturki; 04-23-2012, 10:18 AM.

                  Comment


                  • #10
                    Make sure you're using the correct reference genome

                    In case anyone else encounters this same problem ("[bcf_sync] incorrect number of fields (6 != 5)..."), here's another possible solution:

                    I was having this problem with multiple files, but then I realized I was using the wrong reference fasta as the argument to the mpileup -f option. It wasn't the same reference fasta that was used to create the bam files. When I switched to the correct one, the problem went away.

                    Comment


                    • #11
                      I also had the problem when using nohup. Without nohup. it goes well.

                      Comment


                      • #12
                        The problem occurs when nohup is used with a Unix pipe.

                        I've just sent this reply to the samtools-help mailing list.

                        nohup merges stdout and stderr of the command that it runs. bcftools then
                        doesn't know what to do with the input data since it also contains messages
                        such as "[mpileup] 1 samples in 1 input files".

                        To avoid the problem, you can run the command like this:

                        nohup samtools mpileup ... 2> stderr.txt | bcftools ... - > out.raw.bcf &

                        But I strongly suggest you try the 'screen' program. It is less fragile and
                        much more powerful than nohup. This is the briefest introduction I could find: http://www.askbjoernhansen.com/2006/...ix_screen.html

                        Comment


                        • #13
                          I'm getting this error with this command
                          samtools mpileup -D -S -u -f hg19.fasta child.bam father.bam mother.bam | bcftools view -vcgT trioauto -s sample.txt - > samtoolsCLR.vcf
                          The resulting vcf file just has the INFO and FORMAT ## lines, then:
                          #CHROM POS ID REF ALT QUAL FILTER INFO FORMAT
                          [bcf_sync] incorrect number of fields (0 != 5) at 0:0
                          [afs] 0:0.000
                          I am using the correct fa reference file.
                          My samples file looks like:
                          child.bam
                          father.bam
                          mother.bam

                          Any advice appreciated.

                          Comment


                          • #14
                            Solved my own problem, simple typo in samples file I missed.

                            Comment


                            • #15
                              I had the same problem. I think the problen is that bcftools works with bcf-file but you don't specify for samtools to output bcf. So you should add to samtools mpileup command argument -g.

                              Hope this was helpful.

                              Comment

                              Latest Articles

                              Collapse

                              • seqadmin
                                Strategies for Sequencing Challenging Samples
                                by seqadmin


                                Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
                                03-22-2024, 06:39 AM
                              • seqadmin
                                Techniques and Challenges in Conservation Genomics
                                by seqadmin



                                The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

                                Avian Conservation
                                Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
                                03-08-2024, 10:41 AM

                              ad_right_rmr

                              Collapse

                              News

                              Collapse

                              Topics Statistics Last Post
                              Started by seqadmin, Yesterday, 06:37 PM
                              0 responses
                              10 views
                              0 likes
                              Last Post seqadmin  
                              Started by seqadmin, Yesterday, 06:07 PM
                              0 responses
                              9 views
                              0 likes
                              Last Post seqadmin  
                              Started by seqadmin, 03-22-2024, 10:03 AM
                              0 responses
                              49 views
                              0 likes
                              Last Post seqadmin  
                              Started by seqadmin, 03-21-2024, 07:32 AM
                              0 responses
                              67 views
                              0 likes
                              Last Post seqadmin  
                              Working...
                              X