Thread: ASEReadCounter
View Single Post
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