Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Bowtie output question

    Hi ,
    i'm running the following code to align two primer pair files to my database.

    Code:
    bowtie  -f -t -p 10 -a --best --pairtries 1000 -n0  --fr  --minins 600 --maxins 2500 Bact_archea  -1 fortest.fa -2 revtest.fa
    The output looks like this...
    Code:
    f_1/1   +       GQ246403.1.1498 0       AGAGTTTGATCATGGC        IIIIIIIIIIIIIIII        79309   
    r_2/2   -       GQ246403.1.1498 1479    AAGTCGTAACAAGGTA        IIIIIIIIIIIIIIII        79309   
    f_1/1   +       JF775628.1.1474 0       AGAGTTTGATCATGGC        IIIIIIIIIIIIIIII        79309   
    r_2/2   -       JF775628.1.1474 1455    AAGTCGTAACAAGGTA        IIIIIIIIIIIIIIII        79309
    What I look at this output i know that i have 79309 other hits in the database. This very helpful information that i want.

    Now , i have tried to use the sam output by using the same command adding the "-S" tag . The sam output looks like this.

    Code:
    @HD     VN:1.0  SO:unsorted
    @PG     ID:Bowtie       VN:1.1.2        CL:"bowtie --wrapper basic-0 -f -t -p 1 -a --best --pairtries 1000 -n0 --fr --minins 600 --maxins 2500 Bact_archea -1 fortest.fa -2 revtest.fa -S"
    f_1     99      GQ246403.1.1498 1       255     16M     =       1480    1495    AGAGTTTGATCATGGC        IIIIIIIIIIIIIIII        XA:i:0  MD:Z:16 NM:i:0
    r_2     147     GQ246403.1.1498 1480    255     16M     =       1       -1495   AAGTCGTAACAAGGTA        IIIIIIIIIIIIIIII        XA:i:0  MD:Z:16 NM:i:0
    f_1     99      JF775628.1.1474 1       255     16M     =       1456    1471    AGAGTTTGATCATGGC        IIIIIIIIIIIIIIII        XA:i:0  MD:Z:16 NM:i:0
    r_2     147     JF775628.1.1474 1456    255     16M     =       1       -1471   AAGTCGTAACAAGGTA        IIIIIIIIIIIIIIII        XA:i:0  MD:Z:16 NM:i:0
    f_1     99      EU265948.1.1520 1       255     16M     =       1502    1517    AGAGTTTGATCATGGC        IIIIIIIIIIIIIIII        XA:i:0  MD:Z:16 NM:i:0
    I don't have the total number of alignments as i had before (79309). Is there a way to get the number of alignments found.

    Since my forward and reverse primers contains several primer pairs i want to get the total number of alignments per pair.

    Of course I can do it from my first output file but since it seems that sam format is a standard i was wondering how to get the same information ??

    Thanks,
    david

Latest Articles

Collapse

  • seqadmin
    Strategies for Sequencing Challenging Samples
    by seqadmin


    Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
    03-22-2024, 06:39 AM
  • seqadmin
    Techniques and Challenges in Conservation Genomics
    by seqadmin



    The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

    Avian Conservation
    Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
    03-08-2024, 10:41 AM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, 03-27-2024, 06:37 PM
0 responses
12 views
0 likes
Last Post seqadmin  
Started by seqadmin, 03-27-2024, 06:07 PM
0 responses
11 views
0 likes
Last Post seqadmin  
Started by seqadmin, 03-22-2024, 10:03 AM
0 responses
53 views
0 likes
Last Post seqadmin  
Started by seqadmin, 03-21-2024, 07:32 AM
0 responses
69 views
0 likes
Last Post seqadmin  
Working...
X