Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • lre1234
    Senior Member
    • Aug 2011
    • 110

    #1

    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.
  • Brian Bushnell
    Super Moderator
    • Jan 2014
    • 2709

    #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

    • Richard Finney
      Senior Member
      • Feb 2009
      • 701

      #3
      How many samples?

      Comment

      • lre1234
        Senior Member
        • Aug 2011
        • 110

        #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

        • fanli
          Senior Member
          • Jul 2014
          • 197

          #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

          • SEQadmin2
            Multiomic Liquid Biopsy: What DNA/RNA Panels Reveal That DNA Alone Can’t
            by SEQadmin2



            Tumors contain complex mixtures of cell populations that evolve in response to treatment and host-specific factors. Tissue biopsies provide only static snapshots of a tumor’s state. But what if we could follow that evolution over time using a minimally invasive technique?

            Liquid biopsies can detect tumor-derived material circulating in the bloodstream, providing a way to study treatment response, detect emerging resistance, and discover new biomarkers. Circulating
            ...
            10-14-2026, 09:55 AM
          • SEQadmin2
            New Genomics Technologies Take Aim at Long-Standing Limits
            by SEQadmin2


            Researchers using sequencing and genomics tools often have to make trade-offs. They can choose between speed or scale, short reads or long-range information, or targeted panels or a view of the whole transcriptome. New technologies that have been released this year are built to address those tough choices.

            We asked six companies the same four questions to learn about their latest products. The new technologies bring a lot to the table, including rethinking sequencing
            ...
            09-28-2026, 10:25 AM

          ad_right_rmr

          Collapse

          News

          Collapse

          Topics Statistics Last Post
          Started by SEQadmin2, 10-22-2026, 11:13 AM
          0 responses
          28 views
          0 reactions
          Last Post SEQadmin2  
          Started by SEQadmin2, 10-22-2026, 10:48 AM
          0 responses
          21 views
          0 reactions
          Last Post SEQadmin2  
          Started by SEQadmin2, 10-15-2026, 10:52 AM
          0 responses
          38 views
          0 reactions
          Last Post SEQadmin2  
          Started by SEQadmin2, 10-14-2026, 05:39 AM
          0 responses
          44 views
          0 reactions
          Last Post SEQadmin2  
          Working...