SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
bowtie2-align dided with the signal 11 (SEGV) (core dumped) with bowtie2 harijay Bioinformatics 0 07-28-2015 02:49 AM
Bowtie2: Uniquely mapped reads hoardin Bioinformatics 19 07-28-2015 02:36 AM
Bowtie2 - getting a bam containing mapped reads only Cantabulous Bioinformatics 3 05-01-2014 01:53 PM
Amazon Web Services: Tophat cannot find Bowtie2 even though Bowtie2 is in path Dr_Funkenstein RNA Sequencing 4 03-06-2014 05:55 AM
IGV shows bowtie2 mapped reads not corresponding to genomic sequence Sarah_Tulin Bioinformatics 2 05-09-2013 02:02 PM

Reply
 
Thread Tools
Old 03-11-2016, 12:27 AM   #1
frymor
Senior Member
 
Location: Germany

Join Date: May 2010
Posts: 148
Question Hisat2 vs. bowtie2 - why does bowtie2 mapped better?

Hi all,

I have 358 pairs pair-end samples from human genome. I am looking for the best way to map them to the reference for differential expression analysis.
I have tested the STAR aligner as well as tophat and Hisat2, but somehow they all fail to map.

Than I have tested the bowtie2 algorithm and found out, that it maps a lot better (AFAIK) than all the others.
As the mapping with bowtie2 takes longer than Hisat2 and or STAR, I would like to try and understand first if the results are correct.
I thought, Hisat2 and STAR work not only faster, but suppose to also better.

These are the commands I am using to map my files

Code:
hisat2 -p 15 -x /home/yeroslaviz/genomes/Mus_musculus/HISAT2/grcm38/genome -1 700APAAXX_1_1.fastq.gz -2 700APAAXX_1_2.fastq.gz -S 700APAAXX_1_2.hisat2.sam

bowtie2 --very-sensitive -p 15 --un tmp.unmapped -x /home/yeroslaviz/genomes/Homo_sapiens/GRCh38.p3/bowtie2indexed/Hsapiens.GRCh38 -1 700APAAXX_1_1.fastq.gz -2 700APAAXX_1_2.fastq.gz -S 700APAAXX_1_2.sam
and the results i am getting from the mapping:

Code:
hisat2

33778368 reads; of these:
  33778368 (100.00%) were paired; of these:
    32551450 (96.37%) aligned concordantly 0 times
    564637 (1.67%) aligned concordantly exactly 1 time
    662281 (1.96%) aligned concordantly >1 times
    ----
    32551450 pairs aligned concordantly 0 times; of these:
      20823 (0.06%) aligned discordantly 1 time
    ----
    32530627 pairs aligned 0 times concordantly or discordantly; of these:
      65061254 mates make up the pairs; of these:
        62207013 (95.61%) aligned 0 times
        1237583 (1.90%) aligned exactly 1 time
        1616658 (2.48%) aligned >1 times
7.92% overall alignment rate

bowtie2


33778368 reads; of these:
  33778368 (100.00%) were paired; of these:
    11345742 (33.59%) aligned concordantly 0 times
    16268991 (48.16%) aligned concordantly exactly 1 time
    6163635 (18.25%) aligned concordantly >1 times
    ----
    11345742 pairs aligned concordantly 0 times; of these:
      2699536 (23.79%) aligned discordantly 1 time
    ----
    8646206 pairs aligned 0 times concordantly or discordantly; of these:
      17292412 mates make up the pairs; of these:
        12165533 (70.35%) aligned 0 times
        3229431 (18.68%) aligned exactly 1 time
        1897448 (10.97%) aligned >1 times
81.99% overall alignment rate
Can someone explain to me, why bowtie2 performs better here?


thanks
frymor is offline   Reply With Quote
Old 03-11-2016, 12:30 AM   #2
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,479
Default

You're aligning to human with bowtie2 and mouse with hisat2...

dpryan is offline   Reply With Quote
Old 03-11-2016, 12:34 AM   #3
frymor
Senior Member
 
Location: Germany

Join Date: May 2010
Posts: 148
Default

yes I have seen that. I am trying to delete the post until I can recheck the results, but it was too late
frymor is offline   Reply With Quote
Reply

Tags
bowtie2, hisat2, mapping

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 08:29 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2018, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO