Go Back   SEQanswers > Applications Forums > RNA Sequencing

Thread Tools
Old 02-19-2017, 06:41 AM   #1
Junior Member
Location: Boston

Join Date: May 2012
Posts: 5
Default STAR BAM output processing

Hi all,

I'm running STAR to align RNA-seq reads separately to paternal and maternal genomes for subsequent analysis of allele-specific expression with GATK tools. So as an output I have two BAM files.

I want to keep reads aligned uniquely to only one genome and in cases where reads mapped uniquely to both genomes, use the alignment with the higher alignment quality.

Are there any existing tools for it? I can write a script doing this but I would prefer to use existing tool if there is one.

Thank you!

Last edited by kintany; 02-20-2017 at 09:35 AM. Reason: reformulating better
kintany is offline   Reply With Quote

ase, gatk, rna-seq

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:07 AM.

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