Seqanswers Leaderboard Ad

Collapse

Announcement

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

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

  • #2
    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.

    Comment


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

      Comment


      • #4
        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

        Comment


        • #5
          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.

          Comment


          • #6
            Alternatively:

            samtools view -b <in.bam> <chr>:<pos>-<pos> | samtools pileup - | grep <pos> | awk '{print $4}'
            Mendelian Disorder: A blogshare of random useful information for general public consumption. [Blog]
            Breakway: A Program to Identify Structural Variations in Genomic Data [Website] [Forum Post]
            Projects: U87MG whole genome sequence [Website] [Paper]

            Comment

            Latest Articles

            Collapse

            • seqadmin
              Recent Advances in Sequencing Technologies
              by seqadmin







              Innovations in next-generation sequencing technologies and techniques are driving more precise and comprehensive exploration of complex biological systems. Current advancements include improved accessibility for long-read sequencing and significant progress in single-cell and 3D genomics. This article explores some of the most impactful developments in the field over the past year.

              Long-Read Sequencing
              Long-read sequencing has...
              Yesterday, 01:49 PM
            • seqadmin
              Genetic Variation in Immunogenetics and Antibody Diversity
              by seqadmin



              The field of immunogenetics explores how genetic variations influence immune responses and susceptibility to disease. In a recent SEQanswers webinar, Oscar Rodriguez, Ph.D., Postdoctoral Researcher at the University of Louisville, and Ruben Martínez Barricarte, Ph.D., Assistant Professor of Medicine at Vanderbilt University, shared recent advancements in immunogenetics. This article discusses their research on genetic variation in antibody loci, antibody production processes,...
              11-06-2024, 07:24 PM

            ad_right_rmr

            Collapse

            News

            Collapse

            Topics Statistics Last Post
            Started by seqadmin, Yesterday, 09:29 AM
            0 responses
            14 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, Yesterday, 09:06 AM
            0 responses
            11 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, Yesterday, 08:03 AM
            0 responses
            11 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 11-22-2024, 07:36 AM
            0 responses
            65 views
            0 likes
            Last Post seqadmin  
            Working...
            X