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
    Advanced Methods for the Detection of Infectious Disease
    by seqadmin




    The recent pandemic caused worldwide health, economic, and social disruptions with its reverberations still felt today. A key takeaway from this event is the need for accurate and accessible tools for detecting and tracking infectious diseases. Timely identification is essential for early intervention, managing outbreaks, and preventing their spread. This article reviews several valuable tools employed in the detection and surveillance of infectious diseases.
    ...
    11-27-2023, 01:15 PM
  • seqadmin
    Strategies for Investigating the Microbiome
    by seqadmin




    Microbiome research has led to the discovery of important connections to human and environmental health. Sequencing has become a core investigational tool in microbiome research, a subject that we covered during a recent webinar. Our expert speakers shared a number of advancements including improved experimental workflows, research involving transmission dynamics, and invaluable analysis resources. This article recaps their informative presentations, offering insights...
    11-09-2023, 07:02 AM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, 12-01-2023, 09:55 AM
0 responses
21 views
0 likes
Last Post seqadmin  
Started by seqadmin, 11-30-2023, 10:48 AM
0 responses
20 views
0 likes
Last Post seqadmin  
Started by seqadmin, 11-29-2023, 08:26 AM
0 responses
15 views
0 likes
Last Post seqadmin  
Started by seqadmin, 11-29-2023, 08:12 AM
0 responses
18 views
0 likes
Last Post seqadmin  
Working...
X