Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Brian Bushnell
    replied
    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.

    Leave a comment:


  • strob
    replied
    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.

    Leave a comment:


  • maubp
    replied
    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!

    Leave a comment:


  • 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?

Latest Articles

Collapse

  • seqadmin
    The Impact of AI in Genomic Medicine
    by seqadmin



    Artificial intelligence (AI) has evolved from a futuristic vision to a mainstream technology, highlighted by the introduction of tools like OpenAI's ChatGPT and Google's Gemini. In recent years, AI has become increasingly integrated into the field of genomics. This integration has enabled new scientific discoveries while simultaneously raising important ethical questions1. Interviews with two researchers at the center of this intersection provide insightful perspectives into...
    02-26-2024, 02:07 PM
  • seqadmin
    Multiomics Techniques Advancing Disease Research
    by seqadmin


    New and advanced multiomics tools and technologies have opened new avenues of research and markedly enhanced various disciplines such as disease research and precision medicine1. The practice of merging diverse data from various ‘omes increasingly provides a more holistic understanding of biological systems. As Maddison Masaeli, Co-Founder and CEO at Deepcell, aptly noted, “You can't explain biology in its complex form with one modality.”

    A major leap in the field has
    ...
    02-08-2024, 06:33 AM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, 02-28-2024, 06:12 AM
0 responses
25 views
0 likes
Last Post seqadmin  
Started by seqadmin, 02-23-2024, 04:11 PM
0 responses
72 views
0 likes
Last Post seqadmin  
Started by seqadmin, 02-21-2024, 08:52 AM
0 responses
80 views
0 likes
Last Post seqadmin  
Started by seqadmin, 02-20-2024, 08:57 AM
0 responses
69 views
0 likes
Last Post seqadmin  
Working...
X