Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • Mansequencer
    Member
    • Jul 2010
    • 16

    #1

    Samtools pileup filesize limit

    Is there a file size limitation for samtools pileup command? I converted a 82GB sam file (output from bowtie) into a 25 GB compressed bam and then sorted it using samtools. Now when I try running samtools pileup on the sorted file (25GB), samtools quits saying memory leak. If samtools cannot handle large file, is there a way to split the sorted bam, analyze it with samtools pileup separately and then merge the output.
    Thanks.
  • NicoBxl
    not just another member
    • Aug 2010
    • 264

    #2
    what is your computer config ? how much RAM ?

    Comment

    • Mansequencer
      Member
      • Jul 2010
      • 16

      #3
      Hi Nicobxl,
      I am running samtools on a 16 core 64 bit linux server with 256 GB RAM. Going through the samtools manual, I did not notice any performance related command/options. Please let me know if you have any suggestions.
      Thanks,

      Comment

      • lh3
        Senior Member
        • Feb 2008
        • 686

        #4
        error message? I am sure 200GB BAMs have no problems. Pileup only uses tiny amount of memory.

        Comment

        • Mansequencer
          Member
          • Jul 2010
          • 16

          #5
          Hi Ih3,
          I am using this command to generate a txt file.
          samtools pileup -vcf filename.fas filename_sorted.bam > rawsnp.txt

          I ran samtools on a 1.7 GB file under 32 bit linux system and it worked fine. When I tried running the same file on the 64 bit server, samtools quits saying memory leak 1. continue anyway and then quits.

          Can you suggest some solution to this.
          Thanks,

          Comment

          • Jetse
            Member
            • Nov 2012
            • 38

            #6
            I also had a memory leak, I could solve this by sorting my bam files.

            (this was highest hit when googeling memory leak of mpileup and this is a common mistake)

            Comment

            Latest Articles

            Collapse

            • kinpict.com
              Reply to 16S/ITS Profiling vs. Shotgun Metagenomics: What Changes in Your Workflow?
              by kinpict.com
              The part about this topic feels surprisingly grounded, and the details on the page is the bit that made me slow down. I’d keep <a href="https://kinpict.com/">ai family photo generator</a> bookmarked for later.
              Today, 06:37 PM
            • SEQadmin2
              16S/ITS Profiling vs. Shotgun Metagenomics: What Changes in Your Workflow?
              by SEQadmin2





              NGS has revolutionized clinical microbiology by offering new ways to identify and surveil pathogens without prior knowledge of the target. One popular approach is 16S rRNA and ITS amplicon sequencing, which targets specific microbial markers. But studies are increasingly incorporating shotgun metagenomics, a technique that sequences all DNA from a sample without targeting specific microbial markers. This approach allows more detailed taxonomic classification,...
              10-30-2026, 06:46 AM

            ad_right_rmr

            Collapse

            News

            Collapse

            Topics Statistics Last Post
            Started by SEQadmin2, 11-12-2026, 07:59 AM
            0 responses
            27 views
            0 reactions
            Last Post SEQadmin2  
            Started by SEQadmin2, 11-05-2026, 09:25 AM
            0 responses
            35 views
            0 reactions
            Last Post SEQadmin2  
            Started by SEQadmin2, 11-03-2026, 11:54 AM
            2 responses
            56 views
            0 reactions
            Last Post kinpict.com  
            Started by SEQadmin2, 10-28-2026, 09:57 AM
            0 responses
            175 views
            0 reactions
            Last Post SEQadmin2  
            Working...