unmapped reads by Tophat can be aligned by Bowtie2
I used Tophat to map my pair-end reads to hg19. I sorted my unmapped.bam file, converted it to sam file and extracted pairs into two seperate fastq files. Finally, I use these two fastq files as input to Bowtie2. Surprisingly, there are 40 percent aligned by Bowti2. I know Tophat uses Bowtie2 to align reads. So why are there reads that can be recognized by Bowtie2 but not Tophat?
• 2,075 views
•
link
updated
by
Alex Reynolds
•
written
by
zhangzhenyusz •
•
written
by
zhangzhenyusz •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Why bamToFastq returns empty fastq files?
written by marongiu.luigi
-
Getting unmapped paired-end reads as fastq from unmapped.bam obtained with tophat2
written by cris.kgs •
-
extract one-side unmapped reads from sam
written by liorglic
-
Bowtie2 output alignments and non-alignments in same file
written by zhangzhenyusz •
-
Problem when having multiple Tophat runs inside one workflow in Galaxy
written by tunl
-
Extracting only paired reads from .bam
written by MS485 •
-
.bam files to fastq including the unmapped reads
written by galozs •
-
Extracting paired unmapped reads into fastq gives empty file
written by vjmorley •
-
Tophat more multiple alignments than mapped reads?
written by bmpbowen •
-
Why # Of Reads From Accepted_Hits.Bam + Unmapped.Bam > # Of Reads In Fastq File?
written by newDNASeqer
Could you post all your commands here? What were the settings for each aligner?
For tophat I used: tophat -G <ref> -o <output> <genome> <read_1> <read_2> Then I got both accepted_hits.bam and unmapped.bam. I sorted the unmapped, converted to sam, extracted pair-end reads to two files to form my new read_1 and read_2. Used Bowtie 2 to align and found almost 40% been mapped.