Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
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:
-
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,Tags: None
Latest Articles
Collapse
-
by seqadmin
The complexity of cancer is clearly demonstrated in the diverse ecosystem of the tumor microenvironment (TME). The TME is made up of numerous cell types and its development begins with the changes that happen during oncogenesis. “Genomic mutations, copy number changes, epigenetic alterations, and alternative gene expression occur to varying degrees within the affected tumor cells,” explained Andrea O’Hara, Ph.D., Strategic Technical Specialist at Azenta. “As...-
Channel: Articles
07-08-2024, 03:19 PM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, Yesterday, 06:46 AM
|
0 responses
9 views
0 likes
|
Last Post
by seqadmin
Yesterday, 06:46 AM
|
||
Started by seqadmin, 07-24-2024, 11:09 AM
|
0 responses
26 views
0 likes
|
Last Post
by seqadmin
07-24-2024, 11:09 AM
|
||
Started by seqadmin, 07-19-2024, 07:20 AM
|
0 responses
159 views
0 likes
|
Last Post
by seqadmin
07-19-2024, 07:20 AM
|
||
Started by seqadmin, 07-16-2024, 05:49 AM
|
0 responses
127 views
0 likes
|
Last Post
by seqadmin
07-16-2024, 05:49 AM
|
Leave a comment: