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
    Advancing Precision Medicine for Rare Diseases in Children
    by seqadmin




    Many organizations study rare diseases, but few have a mission as impactful as Rady Children’s Institute for Genomic Medicine (RCIGM). “We are all about changing outcomes for children,” explained Dr. Stephen Kingsmore, President and CEO of the group. The institute’s initial goal was to provide rapid diagnoses for critically ill children and shorten their diagnostic odyssey, a term used to describe the long and arduous process it takes patients to obtain an accurate...
    12-16-2024, 07:57 AM
  • seqadmin
    Recent Advances in Sequencing Technologies
    by seqadmin



    Innovations in next-generation sequencing technologies and techniques are driving more precise and comprehensive exploration of complex biological systems. Current advancements include improved accessibility for long-read sequencing and significant progress in single-cell and 3D genomics. This article explores some of the most impactful developments in the field over the past year.

    Long-Read Sequencing
    Long-read sequencing has seen remarkable advancements,...
    12-02-2024, 01:49 PM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, 12-17-2024, 10:28 AM
0 responses
27 views
0 likes
Last Post seqadmin  
Started by seqadmin, 12-13-2024, 08:24 AM
0 responses
43 views
0 likes
Last Post seqadmin  
Started by seqadmin, 12-12-2024, 07:41 AM
0 responses
29 views
0 likes
Last Post seqadmin  
Started by seqadmin, 12-11-2024, 07:45 AM
0 responses
42 views
0 likes
Last Post seqadmin  
Working...
X