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
            Recent Innovations in Spatial Biology
            by seqadmin


            Spatial biology is an exciting field that encompasses a wide range of techniques and technologies aimed at mapping the organization and interactions of various biomolecules in their native environments. As this area of research progresses, new tools and methodologies are being introduced, accompanied by efforts to establish benchmarking standards and drive technological innovation.

            3D Genomics
            While spatial biology often involves studying proteins and RNAs in their...
            Yesterday, 07:30 PM
          • seqadmin
            Advancing Precision Medicine for Rare Diseases in Children
            by seqadmin




            Many organizations study rare diseases, but few have a mission as impactful as Rady Children’s Institute for Genomic Medicine (RCIGM). “We are all about changing outcomes for children,” explained Dr. Stephen Kingsmore, President and CEO of the group. The institute’s initial goal was to provide rapid diagnoses for critically ill children and shorten their diagnostic odyssey, a term used to describe the long and arduous process it takes patients to obtain an accurate...
            12-16-2024, 07:57 AM

          ad_right_rmr

          Collapse

          News

          Collapse

          Topics Statistics Last Post
          Started by seqadmin, 12-30-2024, 01:35 PM
          0 responses
          23 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 12-17-2024, 10:28 AM
          0 responses
          41 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 12-13-2024, 08:24 AM
          0 responses
          55 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 12-12-2024, 07:41 AM
          0 responses
          41 views
          0 likes
          Last Post seqadmin  
          Working...
          X