Seqanswers Leaderboard Ad

Collapse

Announcement

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

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

  • #2
    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).

    Comment

    Latest Articles

    Collapse

    • seqadmin
      The Impact of AI in Genomic Medicine
      by seqadmin



      Article Coming Soon......
      Today, 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-23-2024, 04:11 PM
    0 responses
    31 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 02-21-2024, 08:52 AM
    0 responses
    46 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 02-20-2024, 08:57 AM
    0 responses
    37 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 02-14-2024, 09:19 AM
    0 responses
    63 views
    0 likes
    Last Post seqadmin  
    Working...
    X