SEQanswers

SEQanswers (http://seqanswers.com/forums/index.php)
-   RNA Sequencing (http://seqanswers.com/forums/forumdisplay.php?f=26)
-   -   ASEReadCounter (http://seqanswers.com/forums/showthread.php?t=74353)

kintany 02-19-2017 06:41 AM

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!


All times are GMT -8. The time now is 11:55 PM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.