Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • kmcarr
    replied
    Originally posted by sshen View Post
    Hi group,

    Happy new year to your all! Does anyone happen to know why this happens?

    --strata has no effect unless combined with -k, -m or -a
    Command: ../bowtie-0.12.7/bowtie -S -k 1 --best --strata -v 2 -p 12

    Thanks in advance,
    sshen,

    These options can be difficult to fully understand at first. If you haven't yet, carefully read the Bowtie manual (http://bowtie-bio.sourceforge.net/manual.shtml) paying special attention to all of the examples under the Reporting Modes section.

    There are two primary groups of options in Bowtie which impact the alignments which are reported. The first group (-n, -v, -e, -l, -I, -X, --fr, --rf, --ff) relate to defining what is or is not a valid alignment. The second group, and the one you are asking about (-k, -a, -m, -M, --best, --strata) affect which and how many of these valid alignments are reported in the output.

    -k simply limits the maximum number of alignments reported. The argument to -k is any positive integer (1..infinity); -a is simply a special case of -k meaning "all". Bowtie must use one or the other of these options, but if you don't explicitly include it on the command line Bowtie will use its default of k=1.

    If by the setting of -k (>1) or -a the valid, reportable alignments fall into more than one stratum (please read the manual for the definition of stratum), the adding the --strata option will further limit the reported alignments, only including those with the same number of mismatches in the seed (if -n mode was used) or in the whole alignments (if -v was used).

    Note that if you use --strata you MUST also use --best.

    -m is not required to use --strata but the two options can combine to further limit the reported alignments.

    It should also be pointed out that in the example command line you provided --strata is meaningless. You have already limited the output to one best alignment (-k 1 --best). --strata will only come into play if you are first allowing more than one reportable alignment (-k > 1 or -a).

    Leave a comment:


  • sshen
    started a topic bowtie related questions

    bowtie related questions

    Hi group,

    Happy new year to your all! Does anyone happen to know why this happens?

    --strata has no effect unless combined with -k, -m or -a
    Command: ../bowtie-0.12.7/bowtie -S -k 1 --best --strata -v 2 -p 12

    Thanks in advance,

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
53 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