It appears that samtools is not using reads which are not marked as proper pairs. This was my problem, and it led me to this thread. I fixed it with bamtools resolve.
Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
Need more information!
Dear All,
Could you please let me know what is rangehg18.chr?
Originally posted by Richard Finney View PostYou can force samtools mpileup with pipe to bcftools to produce lots'o'info like this ...
samtools mpileup -R -d 1000 -l rangehg18.chr -uf hg18.fa in.bam | bcftools view -g - > out.vcf
Note the -g parameter, it forces calls and outputs for all locations.
Check why you're not getting calls in the details of the VCF file.
You can write your own caller based on the VCF information if you need to.Last edited by Soumya18; 02-15-2018, 09:35 PM.
Comment
-
The "-l parameter " (i.e. "dash el") specifies a file of (genomic postions) OR (genomic ranges in bed format) to produce output for. Regions OR positions not specified are dropped (not produced in output, i.e. "skipped").
typing "samtools mpileup" on the command line provides information on all the parameters appropriate to the "mpileup" feature of samtools.
bash-3.2$ samtools mpileup
Usage: samtools mpileup [options] in1.bam [in2.bam [...]]
Input options:
-6, --illumina1.3+ quality is in the Illumina-1.3+ encoding
... [ deleted stuff ]
-G, --exclude-RG FILE exclude read groups listed in FILE
-l, --positions FILE skip unlisted positions (chr pos) or regions (BED)
... [ deleted stuff ]
Comment
Latest Articles
Collapse
-
by seqadmin
The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist...-
Channel: Articles
04-22-2024, 07:01 AM -
-
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...-
Channel: Articles
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
59 views
0 likes
|
Last Post
by seqadmin
04-11-2024, 12:08 PM
|
||
Started by seqadmin, 04-10-2024, 10:19 PM
|
0 responses
57 views
0 likes
|
Last Post
by seqadmin
04-10-2024, 10:19 PM
|
||
Started by seqadmin, 04-10-2024, 09:21 AM
|
0 responses
51 views
0 likes
|
Last Post
by seqadmin
04-10-2024, 09:21 AM
|
||
Started by seqadmin, 04-04-2024, 09:00 AM
|
0 responses
55 views
0 likes
|
Last Post
by seqadmin
04-04-2024, 09:00 AM
|
Comment