Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • nonunique read mapping and snp calling

    Hi all,
    We are having a little debate here that we would like some other insights into.

    The problem is that we are doing targeted sequencing of a few hundred genes, and plan on looking for snps within those genes. the issue is that a subset of the genes, have some pseudogenes, and/or repetitive regions within them, that would prevent the reads from mapping uniquely to those genes. So we have two choices:

    first would be to allow reads to map to multiple locations in the genome, which would include those repetitive regions within our target genes, and then do SNP calling specifically within the target regions.
    or
    second, be conservative with the alignments, and allow for reads to only map uniquely in the genome, and then do the SNP calling.

    From our discussions, the advantage of the first approach would be that there may be some useful snps within those regions that would be able to see, which may be important for our disease or phenotype. We seem to be on the fence as to which approach to do and which would be better.

    Just wanted to know if anyone may have a thought about this.

  • #2
    I suggest throwing out all reads that do not map uniquely before doing variant-calling. If that excludes too much data, then you should probably use longer reads and a longer insert size for pairs.

    Comment


    • #3
      How many samples?

      Comment


      • #4
        I suggest throwing out all reads that do not map uniquely before doing variant-calling. If that excludes too much data, then you should probably use longer reads and a longer insert size for pairs.
        we are using 2x150 PE reads and a 300bp insert size (I think, but need to double check). I guess, my question would be why throwout the reads? What would (in your opinion) be the issue with doing the snp calling in those regions?

        How many samples?
        ~1000 samples

        Comment


        • #5
          Let's say that some of the non-uniquely mapped reads actually come from the pseudogenes. Then you'd probably get false positive SNP calls in your target regions

          Agree w/ BB, the standard is to throw out non-unique mappers.

          Comment

          Latest Articles

          Collapse

          • seqadmin
            Essential Discoveries and Tools in Epitranscriptomics
            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 on Modified Bases...
            Today, 07:01 AM
          • seqadmin
            Current Approaches to Protein Sequencing
            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...
            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
          37 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 04-10-2024, 10:19 PM
          0 responses
          41 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 04-10-2024, 09:21 AM
          0 responses
          35 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 04-04-2024, 09:00 AM
          0 responses
          54 views
          0 likes
          Last Post seqadmin  
          Working...
          X