Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • PacBio bam mpileup result in Segfault

    Hi, I'm new to this forum and I guess I would like to post this question here.

    I had an aligned bam file with reads from 3 PacBio SMRT runs using blasr.
    When I try to call snp between the bam and another bam file from illumina reads from the same library using "samtools mpileup -u". It raised an error with segmentation fault(core dump).

    The error seems to originate from the pacbio bam file since doing mpileup with the illumina bam file alone, it finishes without error.

    The pacbio bam file is sorted and when I try to index the bam, it also finishes without error. Later I tried mpileup without -u option, the text pileup finishes OK.

    Could anyone give me some hints or directions to debug this issue?

    Best,
    jack

  • #2
    How much memory does you machine have?
    What version of samtools are you using?
    What is your command line?
    Can you "samtools view" the region you think it might be at when it bombs?

    Comment


    • #3
      Originally posted by Richard Finney View Post
      How much memory does you machine have?
      What version of samtools are you using?
      What is your command line?
      Can you "samtools view" the region you think it might be at when it bombs?
      I'm using samtools Version: 0.1.17 (r973:277)

      The following command is run under the headnode of our server. I monitored the mem usage to be less than 1% until segfault. The head node has 70GB of memory.
      Code:
      samtools mpileup -u -f ../hg19/hg19-EBV-fos-ecoli/hg19.fa directmap/PacBio_primary.bam > bcfpileup.bcf
      [mpileup] 1 samples in 1 input files
      <mpileup> Set max per-file depth to 8000
      Segmentation fault (core dumped)
      I did the samtools view into that region and there appears to be normal just like the rest of the reads.

      Comment


      • #4
        Anyone have any idea on this issue?

        Or could someone tell me how to convert text mpileup into bcf file?

        Comment


        • #5
          I'm always frustrated to see somebody answer a how-to question by suggesting somebody use a different method, but that's what I'm going to do here.

          samtools mpileup is not designed to work with pacbio, and will probably give no result at best, and more likely something wrong.

          It is much better to use the quiver software that you can download from github.


          Originally posted by jackshen View Post
          Anyone have any idea on this issue?

          Or could someone tell me how to convert text mpileup into bcf file?

          Comment

          Latest Articles

          Collapse

          • 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
          • seqadmin
            The Impact of AI in Genomic Medicine
            by seqadmin



            Artificial intelligence (AI) has evolved from a futuristic vision to a mainstream technology, highlighted by the introduction of tools like OpenAI's ChatGPT and Google's Gemini. In recent years, AI has become increasingly integrated into the field of genomics. This integration has enabled new scientific discoveries while simultaneously raising important ethical questions1. Interviews with two researchers at the center of this intersection provide insightful perspectives into...
            02-26-2024, 02:07 PM

          ad_right_rmr

          Collapse

          News

          Collapse

          Topics Statistics Last Post
          Started by seqadmin, 03-14-2024, 06:13 AM
          0 responses
          33 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 03-08-2024, 08:03 AM
          0 responses
          72 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 03-07-2024, 08:13 AM
          0 responses
          81 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 03-06-2024, 09:51 AM
          0 responses
          68 views
          0 likes
          Last Post seqadmin  
          Working...
          X