Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • samtools mpileup possible bug

    Hello,

    I noticed that last sample (last 3 columns) in mpileup format behave strangely in positions with low base quality (< 13). It looks like:

    0

    where it should be:

    0 * *

    Does anyone have similar issue? Is this a bug?

    It confuses my downstream scripts, because there are fewer columns in those positions than expected.

    I am using Version: 1.2 (using htslib 1.2.1)

  • #2
    It's a longstanding "feature" of samtools mpileup when coverage is 0 because of excluded bases rather than no mapped reads. I would imagine by now there are many scripts that check for this case and they might break if the mpileup were changed. I would simply modify your script to handle this case.

    Comment


    • #3
      Originally posted by dgscofield View Post
      It's a longstanding "feature" of samtools mpileup when coverage is 0 because of excluded bases rather than no mapped reads. I would imagine by now there are many scripts that check for this case and they might break if the mpileup were changed. I would simply modify your script to handle this case.
      I see, would you happen to know if this "feature" occurs only in last individual, or in your experience - it can occur anywhere? Many thanks for your reply!

      Comment


      • #4
        i believe it can occur anywhere, and i think you may find that the columns may be empty, but their positions are still delimited with tabs. so that instead of

        ... <tab> 0 <tab> * <tab> * <tab> ...

        it is

        .... <tab> 0 <tab> <tab> <tab> ...

        There is an additional wrinkle if you use the -s flag to include a read mapping quality column in the mpileup. In the case of "0 * *" entries, they are always three columns rather than the expected 4, while the "0 blank blank" entries are the expected four columns "0 blank blank blank".

        Comment


        • #5
          The additional -s wrinkle has been fixed since samtools 1.0.

          Thanks for the reminder to revisit the main problem reported here.

          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
          34 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
          82 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