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
      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
    29 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 07-16-2024, 05:49 AM
    0 responses
    42 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 07-15-2024, 06:53 AM
    0 responses
    51 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