Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Samtools call variant

    Hi,

    I used samtools to call variant. I use the latest version 1.0. First, I have to generate bcf file. I run the following command.

    samtools mpileup -t DP -ugf my.fasta my.bam

    the result is always like this:

    [mpileup] 1 samples in 1 input files
    <mpileup> Set max per-file depth to 8000
    Aborted

    If I just run

    samtools mpileup -f my.fasta my.bam

    it works fine.

    I have no idea what happens.

    If I go back to use samtools-0.1.19 to run

    samtools mpileup -uD -f my.fasta my.bam

    everything is okay.

    Could anyone helps me? Thank you.

    Silas
    Last edited by silasysh; 09-17-2014, 01:31 PM.

  • #2
    You're not doing anything obviously wrong. If you can share the BAM file (privately at least) then it'd be a good idea to post this on the samtools-help email list. There have been a few issues with the newest release that still need to get resolved and perhaps you've just found another one.

    BTW, what happens if you just do the following with version 1.0:
    Code:
     samtools mpileup -ugf my.fasta my.bam
    I ask because that could help track down exactly where the problem is happening.

    Comment


    • #3
      Thank you.

      If I just do the following:

      samtools mpileup -ugf my.fasta my.bam

      the result is like what I posted before. Nothing more.

      I will try to post bam file to samtools-help email list. Thank you again.

      Silas

      Comment

      Latest Articles

      Collapse

      • seqadmin
        Recent Advances in Sequencing Technologies
        by seqadmin







        Innovations in next-generation sequencing technologies and techniques are driving more precise and comprehensive exploration of complex biological systems. Current advancements include improved accessibility for long-read sequencing and significant progress in single-cell and 3D genomics. This article explores some of the most impactful developments in the field over the past year.

        Long-Read Sequencing
        Long-read sequencing has...
        12-02-2024, 01:49 PM
      • seqadmin
        Genetic Variation in Immunogenetics and Antibody Diversity
        by seqadmin



        The field of immunogenetics explores how genetic variations influence immune responses and susceptibility to disease. In a recent SEQanswers webinar, Oscar Rodriguez, Ph.D., Postdoctoral Researcher at the University of Louisville, and Ruben Martínez Barricarte, Ph.D., Assistant Professor of Medicine at Vanderbilt University, shared recent advancements in immunogenetics. This article discusses their research on genetic variation in antibody loci, antibody production processes,...
        11-06-2024, 07:24 PM

      ad_right_rmr

      Collapse

      News

      Collapse

      Topics Statistics Last Post
      Started by seqadmin, 12-02-2024, 09:29 AM
      0 responses
      144 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 12-02-2024, 09:06 AM
      0 responses
      51 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 12-02-2024, 08:03 AM
      0 responses
      42 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 11-22-2024, 07:36 AM
      0 responses
      72 views
      0 likes
      Last Post seqadmin  
      Working...
      X