Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • akundaje
    started a topic Recommended BWA parameters for SOLiD data ...

    Recommended BWA parameters for SOLiD data ...

    Hi there,

    I am trying to map single-end 33 bp SOLiD data using bwa. I mapped the data using bowtie (-v2 -m1 --best --strata i.e. uniquely mapping reads in the best strata of alignment allowing upto 2 mismatches) and got about 40% of reads mapping. I wanted to try bwa since it supports gapped alignment.

    I have a few questions.

    1. I have a very large number of reads (4,085,198,187) ... it’s nucleosome mapping data and I want the mapping to be conservative. Its ok if I lose reads but I want the mapping to be accurate. It’s the first time I am trying to use BWA since it is able to do gapped alignments, so I was interested in some suggestions regarding parameters. It seems like the default parameters for bwa aln are reasonably conservative. So do you suggest I just go with them?
    The parameters I mainly want feedback on are
    -n
    -k
    -o (max number of gaps)
    -e (max gap extensions)
    -R
    -q

    2. Also, I want to keep only uniquely mapping reads. It wasn’t clear to me at which stage this filtering occurs i.e. bwa aln or samse or postprocessing. Is there some parameter in bwa aln that controls this (similar to –m in bowtie)?
    –R seems to be control the criterion for proceeding with suboptimal alignments so it doesn’t seem to be exactly what I am looking for?
    Or is it the –n parameter in bwa samse?
    Or do I have to postprocess the BAM output from bwa samse using the XT/X0 tag?

    3. For the XT tag, what is the definition of unique. Is it that the number of hits in the best stratum of alignment must be 1 i.e. X0 tag must take value 1? Or does it mean that the total number of hits (including suboptimal ones) based on the bwa aln parameters must be 1?

    4. I’m not clear what the –n INT parameter does in bwa samse : Does that only control what gets written in the XA tag? Or does it also filter out reads that map more than –n INT locations? What is the default value for –n INT? The manual says repetitive hits will be chosen randomly. So for a read that has more than 1 best hit, only one of the locations will be chosen at random and reported?

    I apologize if these are naive questions. I did try to find answers online, by reading the manual and in the forum archives but didn’t get any information to my satisfaction.

    Your input is greatly appreciated.

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...
    04-22-2024, 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
59 views
0 likes
Last Post seqadmin  
Started by seqadmin, 04-10-2024, 10:19 PM
0 responses
57 views
0 likes
Last Post seqadmin  
Started by seqadmin, 04-10-2024, 09:21 AM
0 responses
51 views
0 likes
Last Post seqadmin  
Started by seqadmin, 04-04-2024, 09:00 AM
0 responses
56 views
0 likes
Last Post seqadmin  
Working...
X