Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Calculating coverage of SNPs

    Hello all,

    I'm having a problem calculating the depth of coverage of a set of known variants (SNPs). What I want to know is 1) what percentage of known variant sites are captured in my sequencing, 2) to what depth are those variant sites sequenced to (ie. average depth, >5 reads, >10 reads, ect).

    The problem is I have hundreds of large bam files from separate experiments, and a very large set of SNPs, so I can't spare the space to output a samtools depth file and just iterate through. I was wondering if there is a tool that alone or in combination with basic unix commands will give me what I want. I have tried using coverageBed (bedtools v. 2.17.0), but I can't seem to get it to give me what I want.

    I have the SNPs in a vcf file (that bedtools errors out on), and a bed file that looks like this:

    M 647 648 A/C +
    M 934 935 G/A +
    M 941 942 T/C +
    M 1098 1099 T/C +
    M 1222 1223 G/A +
    M 1247 1248 T/C +
    M 1520 1521 G/A +
    M 1692 1693 C/T +
    M 1831 1832 A/G +
    M 1855 1856 G/A +

    But bedtools doesn't find any coverage or overlap of these SNPs (which is not true).

    Thanks in advance.
    April
    Last edited by aprilw; 08-28-2015, 04:31 AM.

Latest Articles

Collapse

  • seqadmin
    Techniques and Challenges in Conservation Genomics
    by seqadmin



    The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

    Avian Conservation
    Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
    03-08-2024, 10:41 AM
  • seqadmin
    The Impact of AI in Genomic Medicine
    by seqadmin



    Artificial intelligence (AI) has evolved from a futuristic vision to a mainstream technology, highlighted by the introduction of tools like OpenAI's ChatGPT and Google's Gemini. In recent years, AI has become increasingly integrated into the field of genomics. This integration has enabled new scientific discoveries while simultaneously raising important ethical questions1. Interviews with two researchers at the center of this intersection provide insightful perspectives into...
    02-26-2024, 02:07 PM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, 03-14-2024, 06:13 AM
0 responses
34 views
0 likes
Last Post seqadmin  
Started by seqadmin, 03-08-2024, 08:03 AM
0 responses
72 views
0 likes
Last Post seqadmin  
Started by seqadmin, 03-07-2024, 08:13 AM
0 responses
82 views
0 likes
Last Post seqadmin  
Started by seqadmin, 03-06-2024, 09:51 AM
0 responses
68 views
0 likes
Last Post seqadmin  
Working...
X