using only uniquely-mapping reads for alternative splicing using SGSeq
1
0
Entering edit mode
116100016 • 0
@116100016-22156
Last seen 4.9 years ago

I am using SGSeq for alternative splicing analysis of the androgen receptor and was wondering if there was a way to tell the program only to use uniquely mapped reads from the BAM file in its analyses. As there seems to be discrepancy with IGV counting which I'd like to try to tease out.

rna seq SGSeq Alternative splicing • 952 views
ADD COMMENT
0
Entering edit mode
@leonard-goldstein-6845
Last seen 5 months ago
Australia

You can influence which reads are analyzed when setting up the data frame with sample information. To do this column file_bam needs to be a BamFileList and sample information must be stored in a DataFrame object. Please see the manual page for Rsamtools::BamFileList for more details.

ADD COMMENT
0
Entering edit mode

Actually looking at the code in more detail, BamFile details from the sample information currently do not propagate. The hard-coded behavior is that only primary alignments in the BAM file are considered for analysis. In other words, each read should be analyzed only once, but the analysis may include reads that do not align uniquely. For your question it may be easiest to pre-filter BAMs to only contain uniquely mapped reads.

ADD REPLY

Login before adding your answer.

Traffic: 671 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6