Announcement

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

  • Michael.James.Clark
    replied
    Alternatively:

    samtools view -b <in.bam> <chr>:<pos>-<pos> | samtools pileup - | grep <pos> | awk '{print $4}'

    Leave a comment:


  • cristae8
    replied
    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:


  • tanmoy63
    replied
    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:


  • PeteH
    replied
    The coverageBed function in bedtools (http://code.google.com/p/bedtools/) should be able to do what you want.

    Leave a comment:


  • kwatts59
    replied
    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:


  • cristae8
    started a topic *** Per base coverage ***

    *** 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

  • seqadmin
    Multiomics Techniques Advancing Disease Research
    by seqadmin


    New and advanced multiomics tools and technologies have opened new avenues of research and markedly enhanced various disciplines such as disease research and precision medicine1. The practice of merging diverse data from various ‘omes increasingly provides a more holistic understanding of biological systems. As Maddison Masaeli, Co-Founder and CEO at Deepcell, aptly noted, “You can't explain biology in its complex form with one modality.”

    A major leap in the field has
    ...
    02-08-2024, 06:33 AM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, 02-21-2024, 08:52 AM
0 responses
30 views
0 likes
Last Post seqadmin  
Started by seqadmin, 02-20-2024, 08:57 AM
0 responses
20 views
0 likes
Last Post seqadmin  
Started by seqadmin, 02-14-2024, 09:19 AM
0 responses
55 views
0 likes
Last Post seqadmin  
Started by seqadmin, 02-12-2024, 03:37 PM
0 responses
458 views
0 likes
Last Post seqadmin  
Working...
X