Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • How to convert complete bam to mpileup

    Hey,
    While doing mpileup, I could see that lot of reads are not counted in mpileup output although I could see them in tview,
    I am using following command:

    Code:
    samtools mpileup -d 10000 -f Ref.fasta accepted_hits.bam > mpileup
    (BTW: I have this problem in all coverage, so its not because of -d flag)

    I understood from this link that in samtools tview underlined reads are orphan therefore not counted in mpileup. Is there any way NOT to consider any parameter but to get the complete count at each site

    Thank you very much
    seq_g

  • #2
    Originally posted by Seq_g View Post
    Hey,
    I understood from this link that in samtools tview underlined reads are orphan therefore not counted in mpileup. Is there any way NOT to consider any parameter but to get the complete count at each site

    Thank you very much
    Hi- A work around this is to remove the flags that mpileup discards from the input bam. For example orphan reads should be those with flag 8 which means mate unmapped (see also http://picard.sourceforge.net/explain-flags.html). This could be done with a fairly simple script in your favourite language (maybe there are off-the-shelf tools to edit sam flags?).

    In my opinion it would be very useful if mpileup would let the user choose which reads include and exclude (the other day I found out that also reads marked as duplicate are automatically discarded).

    Dario

    Comment

    Latest Articles

    Collapse

    • seqadmin
      Quality Control Essentials for Next-Generation Sequencing Workflows
      by seqadmin




      Like all molecular biology applications, next-generation sequencing (NGS) workflows require diligent quality control (QC) measures to ensure accurate and reproducible results. Proper QC begins at nucleic acid extraction and continues all the way through to data analysis. This article outlines the key QC steps in an NGS workflow, along with the commonly used tools and techniques.

      Nucleic Acid Quality Control
      Preparing for NGS starts with isolating the...
      Yesterday, 01:58 PM
    • seqadmin
      An Introduction to the Technologies Transforming Precision Medicine
      by seqadmin


      In recent years, precision medicine has become a major focus for researchers and healthcare professionals. This approach offers personalized treatment and wellness plans by utilizing insights from each person's unique biology and lifestyle to deliver more effective care. Its advancement relies on innovative technologies that enable a deeper understanding of individual variability. In a joint documentary with our colleagues at Biocompare, we examined the foundational principles of precision...
      01-27-2025, 07:46 AM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by seqadmin, 02-07-2025, 09:30 AM
    0 responses
    19 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 02-05-2025, 10:34 AM
    0 responses
    37 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 02-03-2025, 09:07 AM
    0 responses
    37 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 01-31-2025, 08:31 AM
    0 responses
    41 views
    0 likes
    Last Post seqadmin  
    Working...
    X