Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • samtools mpileup segmentation fault

    hi:

    I clearly do not understand behavior of samtools mpileup. it is frustrating at times for not knowing the exact reasons for segmentation fault.

    1. is seg-fault in mpileup related to a bam file? If so how do I know why it segfaults?

    here is my segfault:
    samtools mpileup -d 100000 -f hg19.fasta -g bc1.bam | bcftools view -e -g -c -v -N -i -1 -I -

    (here I clearly gave -I -- to avoid indels)

    chr3 10183782 . T G 5.29 . // info truncated///
    [afs] 0:1003.293 1:3.825 2:8.882
    Segmentation fault


    2. I take a different bam file and enforce -I -- I get segfault with a different message. In a previous analysis this bam file went through bcftools pipe.

    [bcf_sync] incorrect number of fields (6 != 5) at 2:10191496
    [afs] 0:1438.765 1:13.635 2:23.601
    Segmentation fault


    Is this due to corrput bcftools or samtools.




    3. Cannot exactly understand the reason for segfault in both 1 and 2.


    Is this because I have some problem in bam file or bcftools got corrputed somehow?

    This bam is generated through (PGM) tmap aligner.



    thanks for your help.

    -Adrian

  • #2
    Notice that it´s necesary to sort the bam file before calling the pileup:

    samtools sort alignment.bam > alignment.sort

    hope it helps
    Last edited by yina; 05-24-2012, 10:25 AM.

    Comment


    • #3
      Rather than piping into bcftools, why don't you try just making a pileup file. Then you'll know whether or not bcftools is respnsible.

      One quick thing worth checking: remake the reference index with samtools faidx, and confirm that it gets made with no error. mpileup will check to see if the file exists, and it will make it if it doesn't see a file with the right name, but if the index file was made wrongly, mpileup won't know what to do.

      Comment


      • #4
        Hi,
        I think first you can create mpileup file without error and then try for bcftools.
        Follows these steps

        1.samtools faidx reference.fa (this will produce reference.fai file)
        2.samtools import reference.fai sample.sam sample.sam.bam
        3.samtools sort sample.sam.bam sample.sorted.bam
        4.samtools index sample.sorted.bam sample.sorted.bam.bai
        5.samtools mpileup reference.fa sample.sorted.bam.mpileup

        Please be careful in giving file extensions

        For tyring bcfltools you must need both sorted.bam & bam.bai file

        I hope this helps

        Comment


        • #5
          Originally posted by bioman1 View Post
          Hi,
          I think first you can create mpileup file without error and then try for bcftools.
          Follows these steps

          1.samtools faidx reference.fa (this will produce reference.fai file)
          2.samtools import reference.fai sample.sam sample.sam.bam
          3.samtools sort sample.sam.bam sample.sorted.bam
          4.samtools index sample.sorted.bam sample.sorted.bam.bai
          5.samtools mpileup reference.fa sample.sorted.bam.mpileup

          Please be careful in giving file extensions

          For tyring bcfltools you must need both sorted.bam & bam.bai file

          I hope this helps
          samtool import is deprecated, and you don't actually need an .bam index to run mpileup. And running sort like that will make a file called sample.sorted.bam.bam. And you forgot to specify a input .bam in your pileup command.

          Comment


          • #6
            Segmentation error FIX

            When you try the command, samtools faidx ref.fasta
            You may get this error message: [fai_build_core] different line length in sequence 'scaffold50'. Segmentation fault.

            Use sed '/^$/d' myFile > tt to reformat the linebreaks.

            You can reformat text files in the BASH using several SED commands such as cut/sort/replace/find.

            Comment


            • #7
              I'm extremely annoyed by this problem too.
              The problem seems to have been caused by bcftools when it's trying to convert the uncompressed bcf output from samtools mpileup. This might be caused by some unknown glitches in the bam file. It's extremely hard to pinpoint the problematic read though because usually it takes days to get to that point and when it hangs there's no record of where it stopped. Disabling indel calling by supplying -I seems to solve some but not all "incorrect number of field" error.
              The bottom line is: if there's a glitch in the bam file or the uncompressed bcf output, bcftools should just give a warning, make a record of it instead of failing altogether wasting several days of computation. Someone should improve the source code.

              Comment


              • #8
                Originally posted by melop View Post
                I'm extremely annoyed by this problem too.
                The problem seems to have been caused by bcftools when it's trying to convert the uncompressed bcf output from samtools mpileup. This might be caused by some unknown glitches in the bam file. It's extremely hard to pinpoint the problematic read though because usually it takes days to get to that point and when it hangs there's no record of where it stopped. Disabling indel calling by supplying -I seems to solve some but not all "incorrect number of field" error.
                The bottom line is: if there's a glitch in the bam file or the uncompressed bcf output, bcftools should just give a warning, make a record of it instead of failing altogether wasting several days of computation. Someone should improve the source code.
                I am also fed up of similar error
                [afs] 0:99774.498 1:110.438 2:115.064
                [bcf_sync] incorrect number of fields (4 != 5) at 15:16294867

                Have you found the reason for it . If you can share that that would be a great help

                Regards

                Comment


                • #9
                  Unfortunately I haven't found a solution for it. Perhaps try remaking the bam files would help?

                  Originally posted by figo1019 View Post
                  I am also fed up of similar error
                  [afs] 0:99774.498 1:110.438 2:115.064
                  [bcf_sync] incorrect number of fields (4 != 5) at 15:16294867

                  Have you found the reason for it . If you can share that that would be a great help

                  Regards

                  Comment

                  Latest Articles

                  Collapse

                  • seqadmin
                    Non-Coding RNA Research and Technologies
                    by seqadmin




                    Non-coding RNAs (ncRNAs) do not code for proteins but play important roles in numerous cellular processes including gene silencing, developmental pathways, and more. There are numerous types including microRNA (miRNA), long ncRNA (lncRNA), circular RNA (circRNA), and more. In this article, we discuss innovative ncRNA research and explore recent technological advancements that improve the study of ncRNAs.

                    Nobel Prize for MicroRNA Discovery
                    This week,...
                    10-07-2024, 08:07 AM
                  • seqadmin
                    Recent Developments in Metagenomics
                    by seqadmin





                    Metagenomics has improved the way researchers study microorganisms across diverse environments. Historically, studying microorganisms relied on culturing them in the lab, a method that limits the investigation of many species since most are unculturable1. Metagenomics overcomes these issues by allowing the study of microorganisms regardless of their ability to be cultured or the environments they inhabit. Over time, the field has evolved, especially with the advent...
                    09-23-2024, 06:35 AM

                  ad_right_rmr

                  Collapse

                  News

                  Collapse

                  Topics Statistics Last Post
                  Started by seqadmin, Today, 06:35 AM
                  0 responses
                  7 views
                  0 likes
                  Last Post seqadmin  
                  Started by seqadmin, Yesterday, 02:44 PM
                  0 responses
                  7 views
                  0 likes
                  Last Post seqadmin  
                  Started by seqadmin, 10-11-2024, 06:55 AM
                  0 responses
                  15 views
                  0 likes
                  Last Post seqadmin  
                  Started by seqadmin, 10-02-2024, 04:51 AM
                  0 responses
                  111 views
                  0 likes
                  Last Post seqadmin  
                  Working...
                  X