Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • No output from samtool's mpileup

    I am trying to run samtools for SNP detection, but neither mpileup or the deprecated pileup produces any output.

    My data are in SAM format, which I convert to BAM using samtools thusly:

    ./samtools view -bSt FAIFILE.fai INFILE.sam > OUTFILE.bam

    No errors, seem to get a functional BAM file.

    Now, following the manual,

    ./samtools mpileup -ugf REFSEQ.fasta INFILE.bam | bcftools view -bvcg - > var.raw.bcf

    This produces the following output to sdterr, and no output to screen
    [mpileup] 1 samples in 1 input files
    <mpileup> Set max per-sample depth to 8000

    I have aslo tried a similar approach with pileup:

    ./samtools pileup -vcf REFSEQ.fasta INFILE.bam

    This computes for LONG time, but still produces no output.

    I have no idea what is going wrong... Please help!

    Sasha

  • #2
    I had the very same problem, and found this:



    Not sure it will help you, but it did help me because I had exactly such a bam file (paired-end reads where the read pairs had not been matched up, due to some issues with upstream software.)

    Comment


    • #3
      I also tried to run the command like that, and I got the same thing. It didn't write anywhere.

      So I don't pipe. I run mpileup with -ugf to make a file, (and that step can take a very long time) and then I feed that file to bcftools to make the bcf once it's done.

      The second step is so much faster than the first that you don't save much time by piping.

      Comment


      • #4
        Thanks a lot for both suggestions. However, when I tried downgrading, all I got was a memory error, and no joy.
        Not piping with -ugf gives the following output, but no output file.

        [mpileup] 1 samples in 1 input files
        ?BCa??BCF#chromosom1chromosome2chromosome3 s7_2.bam%##samtoolsVersion=0.1.13 (r926:134)
        <mpileup> Set max per-sample depth to 8000

        I checked to see whether there was some problem with my SAM to BAM conversion by uploading my test file to Galaxy and running mpileup there, and everything seemed to work. So, could there be some weirdness in the way samtools are compiled on my system? I run Mac OS, but received no errors during compilation.

        Comment


        • #5
          Too bad.

          You could try to trawl the Samtools mailing list archives (e g http://sourceforge.net/mailarchive/f...=samtools-help) for further clues.

          Comment


          • #6
            I guess, you need to sort your bam file. i had similar problem and I got away with it by simply sorting it.

            Comment


            • #7
              Relevant for me was lh3's comment in another thread:: Use the -A option with mpileup in order to include non-properly-paired reads.

              Comment

              Latest Articles

              Collapse

              • seqadmin
                Exploring the Dynamics of the Tumor Microenvironment
                by seqadmin




                The complexity of cancer is clearly demonstrated in the diverse ecosystem of the tumor microenvironment (TME). The TME is made up of numerous cell types and its development begins with the changes that happen during oncogenesis. “Genomic mutations, copy number changes, epigenetic alterations, and alternative gene expression occur to varying degrees within the affected tumor cells,” explained Andrea O’Hara, Ph.D., Strategic Technical Specialist at Azenta. “As...
                07-08-2024, 03:19 PM
              • seqadmin
                Exploring Human Diversity Through Large-Scale Omics
                by seqadmin


                In 2003, researchers from the Human Genome Project (HGP) announced the most comprehensive genome to date1. Although the genome wasn’t fully completed until nearly 20 years later2, numerous large-scale projects, such as the International HapMap Project and 1000 Genomes Project, continued the HGP's work, capturing extensive variation and genomic diversity within humans. Recently, newer initiatives have significantly increased in scale and expanded beyond genomics, offering a more detailed...
                06-25-2024, 06:43 AM

              ad_right_rmr

              Collapse

              News

              Collapse

              Topics Statistics Last Post
              Started by seqadmin, 07-16-2024, 05:49 AM
              0 responses
              22 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 07-15-2024, 06:53 AM
              0 responses
              30 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 07-10-2024, 07:30 AM
              0 responses
              40 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 07-03-2024, 09:45 AM
              0 responses
              205 views
              0 likes
              Last Post seqadmin  
              Working...
              X