Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • samtools pileup manual description

    Hello all,

    I am trying to use a shell script that calls samtools pileup function. Because the pileup function is deprecated, the shell script produces nonsense output. I would like to translate the line that calls pileup into a line calling mpileup that produces the same output. Unfortunately, I can not find any detailed information on the pileup command line arguments. The offending line is listed below:

    samtools pileup -c -f $FA -r 0.0000007 -l $BED $BAM > $PILEUP

    in which $FA is the reference fasta file, $BED is a file formatted for the l option of pileup, $BAM is the bam file and $pileup is the pileup file to write to.

    Does anyone know what the -c and -r arguments refer to under the deprecated pileup function and if this call can be converted into a call to mpileup and produce the same output?

    Thanks,
    Michael Gormley

  • #2
    Hi

    I was troubleshooting this on another issue. Installing version 0.1.16 samtools gave me:

    Usage: samtools pileup [options] <in.bam>|<in.sam>

    Option: -s simple (yet incomplete) pileup format
    -S the input is in SAM
    -B disable BAQ computation
    -A use the original MAQ model for SNP calling (DEPRECATED)
    -2 output the 2nd best call and quality
    -i only show lines/consensus with indels
    -Q INT min base quality (possibly capped by BAQ) [13]
    -C INT coefficient for adjusting mapQ of poor mappings [0]
    -m INT filtering reads with bits in INT [0x704]
    -M INT cap mapping quality at INT [60]
    -d INT limit maximum depth for indels [1024]
    -t FILE list of reference sequences (force -S)
    -l FILE list of sites at which pileup is output
    -f FILE reference sequence in the FASTA format

    -c compute the consensus sequence
    -v print variants only (for -c)
    -g output in the GLFv3 format (DEPRECATED)
    -T FLOAT theta in maq consensus calling model (for -c) [0.83]
    -N INT number of haplotypes in the sample (for -c) [2]
    -r FLOAT prior of a difference between two haplotypes (for -c) [0.001]
    -G FLOAT prior of an indel between two haplotypes (for -c) [0.00015]
    -I INT phred prob. of an indel in sequencing/prep. (for -c) [40]


    hope that helps.

    Comment


    • #3
      I am on the same boat. I am trying to revise the shell script that uses depreciated pileup function you mentioned. Did you by any chance figure out how to replace it with mpileup command? Any help would be appreciated. Thanks.

      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, 03-27-2024, 06:37 PM
      0 responses
      12 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 03-27-2024, 06:07 PM
      0 responses
      11 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 03-22-2024, 10:03 AM
      0 responses
      53 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 03-21-2024, 07:32 AM
      0 responses
      69 views
      0 likes
      Last Post seqadmin  
      Working...
      X