Bowtie2: Single End Reads Reporting 100% or 0% Alignment

I am having trouble using single end reads with Bowtie2. This pipeline has previously been used for paired end reads and works just fine. When I use the following command:

bowtie2 -x HP_bacteria -f -U singleendread.fasta -S singleendread.sam

The files are reported as aligned either 100% or 0%.

1 reads; of these:
  1 (100.00%) were unpaired; of these:
    0 (0.00%) aligned 0 times
    1 (100.00%) aligned exactly 1 time
    0 (0.00%) aligned >1 times
100.00% overall alignment rate

1 reads; of these:
  1 (100.00%) were unpaired; of these:
    1 (100.00%) aligned 0 times
    0 (0.00%) aligned exactly 1 time
    0 (0.00%) aligned >1 times
 0.00% overall alignment rate

Based on the previous paired reads, I would expect the files to align between 80 and 90%. Has anyone encountered this problem and know why I am getting this error? Thank you.

sequencing alignment

• 2,517 views • link

updated 6.3 years ago by GenoMax

• written 6.3 years ago by akm126 •

The software is reporting that your input file is one read long.

• 0 views • link

written 6.3 years ago by swbarnes2

Are these two reads the mates of a pair? Please be more clear. What are these two reads and what do they have to do with this 80% alignment?

• 0 views • link

written 6.3 years ago by ATpoint

0 answers

No answers yet.

More posts like this

Discrepancy in HISAT2 Alignment Rates for Paired-End vs. Single-End Reads
written 2.0 years ago by Ujjawal •

Hi everyone, I'm analyzing RNA-seq data using HISAT2 and encountered an interesting discrepancy between paired-end and single-end alignment rates. Here's a summary of the alignment …
unusual alignment rates for IgG and H3K27me3 in cuttag
written 2.5 years ago by david •

I am using Bowtie2 to align my cut&tag data with the reference genome mm10, and the spike-in species is Ecoli. However, I obtained abnormal alignment …
Zero Mapping rate
written 3.0 years ago by waqaskhokhar999

Hi, I have 150 bp long paired-end virus sequencing (fastq) data, I have performed the QC analysis (Fastqc), trimmed out the bad quality reads (Trim_galore) …
miRNA-seq alignment bowtie2
written 4.3 years ago by imaparna27 •

I am trying to align miRNA-seq reads against miRBase reference using Bowtie2, but I am getting a very low overall alignment rate. using Bowtie2: > …
bowtie2 results with NGS data
written 4.9 years ago by adnan.sabar7865 •

Hello, I made index from my reference file and run command to align my metagenomic data by bowtie2. command is `bowtie2 -x <index_referance> -1 <paired_end_read_path_1.fastq> …
Differences in results of Bowtie2
written 5.5 years ago by snehu.sambare •

Bowtie 2 yields different results when **"-a"** option is used. bowtie2 **-a** -x ../../masked_rDNA ../hTERT-TP2-R1_p1.fq --no-unal > temp.sam 16914184 reads; of these: 16914184 (100.00%) were …
Overall alignment 0.00%
written 7.7 years ago by heeyairenkhaba •

hi everyone, I have done alignment using HiSAT2 the result comes likes this 19887884 (100.00%) were unpaired; of these: 19887828 (100.00%) aligned 0 times 36 …
bowtie2 paired end insertion length
written 8.1 years ago by jesselee516

Hi, I am using bowtie2 to do alignment to reference genome. Actually, instead of getting real data, I simulated the paired-end reads(100bp) with insertion length …
alignment rate in bowtie2
written 10.5 years ago by zizigolu

Hi, I used bowtie2 to map GSE69414 data sets but rate of alignment is 0 or max 0.02 How is that possible? Am I doing …
unexpected bowtie2 unpaired alignment behavior
written 12.0 years ago by mark.rose •

Hello Seemingly according to the manual, bowtie2 will attempt to align paired-end reads first as pairs but then as unpaired reads if necessary. If this …