SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
STAR mapping error Him26 Bioinformatics 2 03-09-2016 03:17 PM
STAR vs tophat2 - mapping reads to long exons frymor Bioinformatics 7 08-18-2015 02:17 PM
STAR Mapping taking too long mtiwaridros RNA Sequencing 1 05-02-2014 08:26 AM
STAR multi-mapping id0 Bioinformatics 1 03-28-2014 07:41 AM
Using Star/ bowtie on cluster babi2305 Bioinformatics 7 02-06-2013 11:11 AM

Reply
 
Thread Tools
Old 04-13-2018, 08:00 AM   #1
boro2014
Junior Member
 
Location: Belgium

Join Date: Jan 2014
Posts: 3
Default Mapping using STAR and Bowtie

Hi all,

I want to perform an alignment of short reads against a reference genome (human genome from Ensembl). For that I've used STAR, and Bowtie (not Bowtie2). I've some questions about that.

1- I've used STAR by setting the more important parameters as follows:

--outFilterScoreMinOverLread 0.4 --outFilterMatchNminOverLread 0.4 --outFilterMatchNmin 15

Here is the final log file content generated by STAR:

Number of input reads | 5845995
Average input read length | 27
UNIQUE READS:
Uniquely mapped reads number | 446299
Uniquely mapped reads % | 7.63%
Average mapped length | 25.10
Number of splices: Total | 4043
Number of splices: Annotated (sjdb) | 622
Number of splices: GT/AG | 3930
Number of splices: GC/AG | 108
Number of splices: AT/AC | 5
Number of splices: Non-canonical | 0
Mismatch rate per base, % | 1.74%
Deletion rate per base | 0.08%
Deletion average length | 1.36
Insertion rate per base | 0.04%
Insertion average length | 1.07
MULTI-MAPPING READS:
Number of reads mapped to multiple loci | 4346149
% of reads mapped to multiple loci | 74.34%
Number of reads mapped to too many loci | 963965
% of reads mapped to too many loci | 16.49%
UNMAPPED READS:
% of reads unmapped: too many mismatches | 0.00%
% of reads unmapped: too short | 0.31%
% of reads unmapped: other | 1.22%
CHIMERIC READS:
Number of chimeric reads | 1503
% of chimeric reads | 0.03%

According to this, I got a very low rate of uniquely mapped reads (I used the parameter --outFilterMultimapNmax 10).

How can I improve the uniquely mapped reads rate ?

2 - I've also used Bowtie to map the same sample to the same reference genome. Bowtie didn't output the number of uniquely mapped reads and the number of multimapped reads as STAR did. I checked the log file generated by Bowtie:

# reads processed: 5845995
# reads with at least one reported alignment: 4206928 (71.96%)
# reads that failed to align: 1639067 (28.04%)
Reported 98304317 alignments to 1 output stream(s)

I want to output the number of uniquely mapped reads and multimapped reads from the sam file generated by bowtie. So I have done the following:

samtools view -Sb myfile.sam > myfile.bam
samtools view -F 4 myfile.bam | grep -v "XS:" | wc -l

But it didn't work for me, and I got a wrong number of uniquely mapped reads. I know that those commands work for files generated by bowtie2.

Do you have any solution able to retrieve the number of uniquely mapped, and multimapped reads from the sam file generated by bowtie ?

Thanks a lot
boro2014 is offline   Reply With Quote
Old 04-13-2018, 08:10 AM   #2
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 6,761
Default

Cross-posted: https://www.biostars.org/p/309209/
GenoMax is offline   Reply With Quote
Reply

Tags
bowtie, rnaseq alignment, short read alignment, star aligner

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 01:16 AM.


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