Alternatively:
samtools view -b <in.bam> <chr>:<pos>-<pos> | samtools pileup - | grep <pos> | awk '{print $4}'
Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
Thanks, everyone. bedTools seems to work the best and it's fully documented so very easy to use! Much thanks to Aaron (the author) & Pete.
Leave a comment:
-
you can use samtools pileup -c from on sorted indexed bam file of your alignment..this pileup consensus will give you read depth in each position and base composition in each also
Leave a comment:
-
The coverageBed function in bedtools (http://code.google.com/p/bedtools/) should be able to do what you want.
Leave a comment:
-
I loaded all my short read data into an SQL database so that this information can be pulled out with a quick query.
For example, to find how many reads cover position 72342 of chr 1, just type in the following query. This assumes 50bp reads.
select count(*) from short_reads
where (chromosome = 'chr1')
and (start_position >= 72293)
and (start_position <= 72342)
To perform this query for many SNPs, a PERL script can be written.
Leave a comment:
-
*** Per base coverage ***
Hi,
This seems like a silly question, but I can't quite find the answer for it anywhere. I have a list of single nucleotide position, ie,
chr 1 72342
chr 2 32912
chr X 184343
etc,
to which I would like to find their coverage, ie, how many reads cover that position and what are those bases.
This is almost like calling SNP, but not quite because I just like to know the coverage even if no SNP is present. For this reason, I can't just use a SNP caller. Does anyone have a solution for this?
Thanks so much.
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, 07-25-2024, 06:46 AM
|
0 responses
9 views
0 likes
|
Last Post
by seqadmin
07-25-2024, 06:46 AM
|
||
Started by seqadmin, 07-24-2024, 11:09 AM
|
0 responses
28 views
0 likes
|
Last Post
by seqadmin
07-24-2024, 11:09 AM
|
||
Started by seqadmin, 07-19-2024, 07:20 AM
|
0 responses
161 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: