SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
HELP!! XS tag to filter unique reads Do NOT match bowtie2 report momo Bioinformatics 24 10-11-2019 08:11 AM
Filter unique hits in Bowtie2 for stranded RNA-seq MattB Bioinformatics 2 05-05-2015 05:37 AM
STAR vs TopHat viewing junctions biocomputer Bioinformatics 1 11-27-2014 11:45 PM
Generating genome with splice junctions - STAR shivshiv Bioinformatics 8 09-27-2014 02:21 AM
Filter reads by mappability, unique reads, not more than 1 mismatches hari_iyer16 Bioinformatics 2 02-28-2014 10:47 AM

Reply
 
Thread Tools
Old 10-04-2015, 09:20 PM   #1
frymor
Senior Member
 
Location: Germany

Join Date: May 2010
Posts: 149
Default RNA-STAR filter out unique Junctions

Hi all,

I am trying to filter un-annotated splice junction from my STAR mapping run, but with no success.
I would like to exclude all junction from both my *.bam and *.SJ.out.tab files which have less than 5 reads overlapping the junction motifs.
I have trid it several times with different parameters until I came these construction:
Code:
--alignSJoverhangMin 8
--alignSJDBoverhangMin 5
--outFilterType BySJout
--outSJfilterCountTotalMin 8 5 5 5
--outSJfilterCountUniqueMin -1 -1 -1 -1
With this settings I get a lot less junctions than with the default values. This is fine by me, as I hope only to get those with high coverage, but even here I still get this kind of junctions in my SJ.tab.out file:

Code:
awk 'FNR==5442 {print}' STAR.SJ.out.tab
 2L      11438179        11487284        1       1       0       1       0       22
This is an unannotated splice junction with only one read overlapping the intron motifs.

What do I need to change to get rid of these junctions as well?

Thanks in advance for the help.

Assa
frymor is offline   Reply With Quote
Reply

Tags
splice junctions, 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 08:58 AM.


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