Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
STAR: ultrafast universal RNA-seq aligner alexdobin Bioinformatics 218 04-02-2018 05:59 PM
STAR rna seq Aligner installation arcolombo698 Bioinformatics 10 10-25-2016 09:31 AM
STAR RNA-seq aligner yields different read length than fastqc/zcat? LeonDK RNA Sequencing 2 01-26-2015 01:25 AM
STAR Aligner sbdk82 Bioinformatics 7 09-22-2014 12:29 PM
STAR aligner issue shocker8786 Bioinformatics 3 05-21-2014 01:59 AM

Thread Tools
Old 10-08-2015, 11:36 AM   #1
Location: Philadelphia, PA

Join Date: Aug 2013
Posts: 14
Default STAR Aligner Output BAM Processing

Hi everyone,

I am a new user of STAR Aligner. I have been using GSNAP until now, but my new group likes STAR better. I am very excited about using it and have a couple of questions too. We will be using paired end fastq files as input. I saw there are couple of really cool things with STAR:

1. Output coordinate-sorted/unsorted BAM files using --outSAMtype
2. Output read counts using --quantMode

Question 1: Does the --quantMode perform equally well to htseq-count? Are the read-counts identical in both cases?

When I was using GSNAP, after I obtained the output BAM file, I used to perform the following steps before using htseq-count:
1. samtools fixmate to fill in mate information
2. bamtools filter to keep only "reads in proper pair" and CIGAR string should indicate atleast one match i.e. 'M'

Question 2: So do I have to do these steps after I obtain the bam output from STAR?

Question 3: My command line for htseq-count is as follows with the GSNAP output:
samtools view -f 0x0002 sample.bam | htseq-count - $GTF > sample.counts

Question 3: Do I have to use '-f 0x0002' (i.e. reads mapped in proper pair) while using htseq-count on STAR output bam file as well?

Thanks a lot!
Komal Rathi
Bioinformatics Application Developer
University of Pennsylvania
komalsrathi is offline   Reply With Quote

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 05:29 AM.

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