Bowtie output question

danova

Member

Join Date: Sep 2010
Posts: 27

Bowtie output question

07-22-2015, 02:34 AM

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

Tags: bowtie, samtools

Previous template Next

From Collection to Sequencing: Why Sample Preparation and Preservation Define Sequencing Data

by SEQadmin2

Data variability is still an issue in sequencing technologies despite the advances in reproducibility and accuracy of these platforms. But the problem does not originate in the sequencing itself, but in the previous steps, before the sample reaches the sequencer.

The first step is collection, followed by preservation and sample preparation for analysis. Most scientists overlook those steps, but not being careful might just be skewing the experiment’s results.
...
- Channel: Articles
06-02-2026, 10:05 AM
Single-Cell Sequencing at an Inflection Point: Early Impacts of New Platforms and Emerging Trends

by SEQadmin2

With the launch of new single-cell sequencing platforms in 2026, the field stands at an exciting inflection point. This article surveys the most impactful advances in the field and discusses how they’re reshaping research in cancer, immunology, and beyond.

Introduction

Single-cell sequencing technologies have undergone remarkable advances over the past decade, transitioning from low-throughput experimental approaches to highly scalable platforms capable of...
- Channel: Articles
05-22-2026, 06:42 AM

Topics	Statistics	Last Post
A New Method Makes Hantavirus Genome Analysis Faster and More Accessible by SEQadmin2 Started by SEQadmin2, 06-05-2026, 10:09 AM	0 responses 12 views 0 reactions	Last Post by SEQadmin2 06-05-2026, 10:09 AM
A New Single-Cell Method Maps DNA-Protein Interactions by SEQadmin2 Started by SEQadmin2, 06-04-2026, 08:59 AM	0 responses 23 views 0 reactions	Last Post by SEQadmin2 06-04-2026, 08:59 AM
Long-Read RNA Sequencing Uncovers a Hidden Layer of Immune Cell Regulation by SEQadmin2 Started by SEQadmin2, 06-02-2026, 12:03 PM	0 responses 28 views 0 reactions	Last Post by SEQadmin2 06-02-2026, 12:03 PM
DNA Methylation Study Reveals How Epigenetic Changes Pass Between Generations by SEQadmin2 Started by SEQadmin2, 06-02-2026, 11:40 AM	0 responses 22 views 0 reactions	Last Post by SEQadmin2 06-02-2026, 11:40 AM

Unconfigured Ad

Bowtie output question

Latest Articles

ad_right_rmr

News