gagneurlab / FRASER

FRASER - Find RAre Splicing Events in RNA-seq
MIT License
36 stars 20 forks source link

filter out reads with MAPQ=0 & duplicated reads #32

Open yangyangzijiajia opened 2 years ago

yangyangzijiajia commented 2 years ago

for fraser junction count analysis, in default setting: ScanBamParam(mapqFilter=0), but sometime, we can see several duplicated reads with mapQ=0 and reads with multiple locations with mapQ=0, it can also cause the false positive abrrent splicing events. what's your opinion for mapqfilter setting? if we set mapqFilter >0, whether it results other false negative events?