Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Bowtie behavior with --best

    ok, i apologize it may be just something very obvious but despite reading manual several times, i can not understand the following behavior. A sequence searched against build in e_coli index:

    bowtie e_coli -c TTCTTTTTTACCTTCCCC -n 3 -e 90 -a
    0 + gi|110640213|ref|NC_008253.1| 2745489 TTCTTTTTTACCTTCCCC IIIIIIIIIIIIIIIIII 0
    0 - gi|110640213|ref|NC_008253.1| 3948415 GGGGAAGGTAAAAAAGAA IIIIIIIIIIIIIIIIII 0 9:C>T,10:T>G,12:G>A
    0 + gi|110640213|ref|NC_008253.1| 2345717 TTCTTTTTTACCTTCCCC IIIIIIIIIIIIIIIIII 0 13:G>T,16:T>C,17:G>C
    0 + gi|110640213|ref|NC_008253.1| 1032847 TTCTTTTTTACCTTCCCC IIIIIIIIIIIIIIIIII 0 11:G>C,13:A>T,16:G>C
    0 + gi|110640213|ref|NC_008253.1| 2899814 TTCTTTTTTACCTTCCCC IIIIIIIIIIIIIIIIII 0 11:T>C,15:A>C
    0 + gi|110640213|ref|NC_008253.1| 2939803 TTCTTTTTTACCTTCCCC IIIIIIIIIIIIIIIIII 0 9:C>A,10:A>C,15:G>C
    0 + gi|110640213|ref|NC_008253.1| 550112 TTCTTTTTTACCTTCCCC IIIIIIIIIIIIIIIIII 0 0:G>T,1:C>T
    0 + gi|110640213|ref|NC_008253.1| 4879006 TTCTTTTTTACCTTCCCC IIIIIIIIIIIIIIIIII 0 8:C>T,13:G>T,15:T>C
    0 + gi|110640213|ref|NC_008253.1| 4630520 TTCTTTTTTACCTTCCCC IIIIIIIIIIIIIIIIII 0 8:G>T,12:A>T,14:G>C
    # reads processed: 1
    # reads with at least one reported alignment: 1 (100.00%)
    # reads that failed to align: 0 (0.00%)
    Reported 9 alignments to 1 output stream(s)

    =========================================================

    bowtie e_coli -c TTCTTTTTTACCTTCCCC -n 3 -e 90 -a --best
    0 + gi|110640213|ref|NC_008253.1| 2745489 TTCTTTTTTACCTTCCCC IIIIIIIIIIIIIIIIII 0
    0 + gi|110640213|ref|NC_008253.1| 550112 TTCTTTTTTACCTTCCCC IIIIIIIIIIIIIIIIII 0 0:G>T,1:C>T
    0 + gi|110640213|ref|NC_008253.1| 2899814 TTCTTTTTTACCTTCCCC IIIIIIIIIIIIIIIIII 0 11:T>C,15:A>C
    0 + gi|110640213|ref|NC_008253.1| 4760194 TTCTTTTTTACCTTCCCC IIIIIIIIIIIIIIIIII 0 0:A>T,2:A>C,14:A>C
    0 + gi|110640213|ref|NC_008253.1| 2345717 TTCTTTTTTACCTTCCCC IIIIIIIIIIIIIIIIII 0 13:G>T,16:T>C,17:G>C
    # reads processed: 1
    # reads with at least one reported alignment: 1 (100.00%)
    # reads that failed to align: 0 (0.00%)
    Reported 5 alignments to 1 output stream(s)

    =========================================================

    Q.1 what happened to other 4 alignments with best (no quality is considered as you see)
    Q.2 Where is the alignment 4760194 coming from when --best is used

    Thanks for your response
    Last edited by analyst; 05-05-2011, 06:10 AM.

Latest Articles

Collapse

  • seqadmin
    Non-Coding RNA Research and Technologies
    by seqadmin




    Non-coding RNAs (ncRNAs) do not code for proteins but play important roles in numerous cellular processes including gene silencing, developmental pathways, and more. There are numerous types including microRNA (miRNA), long ncRNA (lncRNA), circular RNA (circRNA), and more. In this article, we discuss innovative ncRNA research and explore recent technological advancements that improve the study of ncRNAs.

    Nobel Prize for MicroRNA Discovery
    This week,...
    10-07-2024, 08:07 AM
  • seqadmin
    Recent Developments in Metagenomics
    by seqadmin





    Metagenomics has improved the way researchers study microorganisms across diverse environments. Historically, studying microorganisms relied on culturing them in the lab, a method that limits the investigation of many species since most are unculturable1. Metagenomics overcomes these issues by allowing the study of microorganisms regardless of their ability to be cultured or the environments they inhabit. Over time, the field has evolved, especially with the advent...
    09-23-2024, 06:35 AM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, 10-02-2024, 04:51 AM
0 responses
103 views
0 likes
Last Post seqadmin  
Started by seqadmin, 10-01-2024, 07:10 AM
0 responses
111 views
0 likes
Last Post seqadmin  
Started by seqadmin, 09-30-2024, 08:33 AM
1 response
114 views
0 likes
Last Post EmiTom
by EmiTom
 
Started by seqadmin, 09-26-2024, 12:57 PM
0 responses
20 views
0 likes
Last Post seqadmin  
Working...
X