Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • How do I identify reads that overhang the end of a reference?

    Hi,
    I have 76nt single-end reads and I want to map them to a reference and identify the reads that overhang the end of the reference.
    What would be the best mapper to do the mapping with?
    If I use BWA as the mapper, I can find reads that have up to about 9nt overhangs (by comparing the reads' unclipped start site to the length of the reference the read maps to), but what I'm looking for is some way to find reads that overhang by about 40 nts, i.e. the first 37nts of the read match the end of the reference, but (obviously) the last 40nts don't.
    Do any next-gen mappers allow such mapping?

  • #2
    Well the SAM/BAM file format itself doesn't allow this directly (the start/end co-ordinates must be within the length of the reference), so any reads which did map off the ends would probably have to be represented with a large insert or soft clip at the end.

    I'm not sure of any specialized tool to do this. What I would consider is trying to map the first and last halves of your reads to the start/end of the reference. i.e. Some scripting to prepare the modified input files, then reuse a mapper like BWA, and interpret the results accordingly. Fiddly - so I hope someone else has done the hard work already!

    Comment


    • #3
      Dear all,

      I want to re-open this discussion as I have a similar question: how can I detect and extract overhanging reads, and this at both ends of my reference?
      I have overlapping MiSeq reads which I merged, resulting in reads with an average length of 300 - 350 bp. I'm only interested in what we call hybrid reads: reads partly mapping with the borders of my reference.
      Any advice on a mapping tool? Even use Blast?
      Thanks for the input.

      Comment


      • #4
        strob,

        If you run BBMap with default parameters, it will soft-clip ONLY reads that map off the ends of reference scaffolds. So, you would just have to filter the sam file for reads containing the "S" symbol in their cigar strings.

        You can adjust the fraction of the read that must be on the reference with the "minratio" flag (default 0.56). I suggest something like this:

        minratio=0.2 maxindel=0

        ...which will look for alignments with as little as 20% of the read over the reference and not look for alignments containing indels.

        Comment

        Latest Articles

        Collapse

        • seqadmin
          Exploring the Dynamics of the Tumor Microenvironment
          by seqadmin




          The complexity of cancer is clearly demonstrated in the diverse ecosystem of the tumor microenvironment (TME). The TME is made up of numerous cell types and its development begins with the changes that happen during oncogenesis. “Genomic mutations, copy number changes, epigenetic alterations, and alternative gene expression occur to varying degrees within the affected tumor cells,” explained Andrea O’Hara, Ph.D., Strategic Technical Specialist at Azenta. “As...
          07-08-2024, 03:19 PM
        • seqadmin
          Exploring Human Diversity Through Large-Scale Omics
          by seqadmin


          In 2003, researchers from the Human Genome Project (HGP) announced the most comprehensive genome to date1. Although the genome wasn’t fully completed until nearly 20 years later2, numerous large-scale projects, such as the International HapMap Project and 1000 Genomes Project, continued the HGP's work, capturing extensive variation and genomic diversity within humans. Recently, newer initiatives have significantly increased in scale and expanded beyond genomics, offering a more detailed...
          06-25-2024, 06:43 AM

        ad_right_rmr

        Collapse

        News

        Collapse

        Topics Statistics Last Post
        Started by seqadmin, 07-19-2024, 07:20 AM
        0 responses
        141 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 07-16-2024, 05:49 AM
        0 responses
        116 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 07-15-2024, 06:53 AM
        0 responses
        109 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 07-10-2024, 07:30 AM
        0 responses
        43 views
        0 likes
        Last Post seqadmin  
        Working...
        X