SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Bismark - A New Tool for Mapping and Analysis of Bisulfite-Seq Data fkrueger Bioinformatics 649 10-05-2018 01:43 AM
why low mapping rates for RNAseq? NGSfan RNA Sequencing 49 09-09-2014 12:34 PM
mapping tool for resequencing glacerda Bioinformatics 0 04-01-2011 07:22 AM
Is a Gene mapping tool available? ritzriya Bioinformatics 0 02-15-2011 12:45 AM
PerM is an ultra-fast and sensitive SOLiD reads mapping tool KevinLam Bioinformatics 7 06-18-2010 03:03 AM

Reply
 
Thread Tools
Old 03-04-2012, 06:10 PM   #1
sissi
Junior Member
 
Location: Singapore

Join Date: Mar 2012
Posts: 5
Default best RNAseq mapping tool for bacteria

Hi,

I am new to NGS data and am trying to find the best mapping tool for bacterial genome. My data consists of FASTQ files for paired end RNA-seq on Illumina. I have downloaded the reference genome from NCBI website (fasta file).
So far I have identified a few softwares (MAQ, Tophat, BWA, Glimmer, SHRIMP...), but I don't know which one to choose.

Thank you for your help!

Regards,
S
sissi is offline   Reply With Quote
Old 03-04-2012, 11:07 PM   #2
swbarnes2
Senior Member
 
Location: San Diego

Join Date: May 2008
Posts: 912
Default

Bacteria have almost no introns, so you are just mapping to the genome. And since the genome is so small, you can use as slow an algorithm as you like, it'll still be done pretty fast.

As long as it'll output .bam files, you'll probably be fine.
swbarnes2 is offline   Reply With Quote
Old 03-05-2012, 03:33 AM   #3
swaraj
Member
 
Location: Naples, Italy

Join Date: Feb 2012
Posts: 50
Default

For aligning RNA-seq reads consider the following points

1) Paired-end or single-end sequencing?
2) spliced reads?
4) interested in finding novel transcripts?

I suggest subread (http://subread.sourceforge.net/) as it is fast, reported to be accurate and integrates within R environment.
swaraj is offline   Reply With Quote
Old 03-05-2012, 04:55 AM   #4
sissi
Junior Member
 
Location: Singapore

Join Date: Mar 2012
Posts: 5
Default

Hi,
Thx for your answers.

It is for pair end sequencing data with 75bp library.
I will be looking at SNPs and differential expression data.
I don't know exactly if there are many spliced read, but my guess would be very little as it's a bacteria.

I read a few papers on the subject and it seems that the quality of the reference genome is also important. How do I know if that condition is verified?

I will have a look at the R script
Thank you!

Last edited by sissi; 03-05-2012 at 04:58 AM.
sissi is offline   Reply With Quote
Old 03-05-2012, 05:01 AM   #5
Hobbe
Member
 
Location: Uppsala, Sweden

Join Date: Apr 2010
Posts: 29
Default

Shrimp2 is easy to use and runs very well on several cores if you have them, so it is fairly quick. In my experience, it maps well too.
Hobbe is offline   Reply With Quote
Old 03-05-2012, 05:16 AM   #6
swaraj
Member
 
Location: Naples, Italy

Join Date: Feb 2012
Posts: 50
Default

Shrimp2 is a good choice for bacterial genomes. For the genome quality you have to look into the publication which reported the genome.
swaraj is offline   Reply With Quote
Old 03-05-2012, 03:48 PM   #7
sissi
Junior Member
 
Location: Singapore

Join Date: Mar 2012
Posts: 5
Default

Thanks a lot.
sissi is offline   Reply With Quote
Reply

Tags
bacteria, illumina, mapping, rnaseq

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 12:48 AM.


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