Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • lh3
    replied
    Use "-A" to override the default behaviour.

    Leave a comment:


  • kopi-o
    replied
    I've run into this issue as well. I think it would be helpful if there was a option to allow non-properly paired reads in mpileup.

    Leave a comment:


  • Hit
    replied
    Does anyone have any opinion on these questions? I just posted a similar question in an independent thread before seeing this one:

    Discussion of next-gen sequencing related bioinformatics: resources, algorithms, open source efforts, etc


    I also observed the same number of SNPs from mpileup when not using "-q 1" to filter any reads and when using "-q 1" to filter alignments with mapping quality score 0.

    Any suggestion would be appreciated!

    Leave a comment:


  • samtools mpileup only uses reads with IsProperPair set to 1?

    Hi,

    I have some questions about the samtools mpileup usage, what kind of filters does it apply before generating the consensus (or variant call)?

    1. it seems all reads with the IsProperPair flag unset (0) are all discarded. Unfortunately, none of the reads in my BAM files has this flag set, so running samtools mpileup gives a blank output.

    2. what about orphan reads (i.e., one read is mapped but not the mate, so it could not possibly be a "proper pair")? Yet, I'd like to incorporate info from these reads to make variant calling

    3. I'd suppose mpileup will discard reads with 0 mapping quality, but just want to confirm

    Thanks,

Latest Articles

Collapse

  • seqadmin
    Current Approaches to Protein Sequencing
    by seqadmin


    Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...
    04-04-2024, 04:25 PM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, 04-11-2024, 12:08 PM
0 responses
34 views
0 likes
Last Post seqadmin  
Started by seqadmin, 04-10-2024, 10:19 PM
0 responses
37 views
0 likes
Last Post seqadmin  
Started by seqadmin, 04-10-2024, 09:21 AM
0 responses
33 views
0 likes
Last Post seqadmin  
Started by seqadmin, 04-04-2024, 09:00 AM
0 responses
54 views
0 likes
Last Post seqadmin  
Working...
X