Reads with ambiguous annotation
1
0
Entering edit mode
raya.fai ▴ 60
@rayafai-9396
Last seen 13 months ago
Israel

Hi,

We have a question about the proper handing of ambiguous annotation reads when using DESeq2. In our bacterial RNAseq data there is a fraction of reads (0.5-2%) that map inconclusively to two adjacent genes (gene A and gene B). We count such reads twice, one for gene A and one for gene B. Can we use DESeq2 for such counting? We prefer not to exclude these reads because they may have biological meaning.

Thank you very much, Raya

deseq2 • 613 views
ADD COMMENT
0
Entering edit mode
@mikelove
Last seen 1 day ago
United States

You can use a method such as Salmon to perform quantification followed by tximport and DESeqDataSetFromTximport. This is our official pipeline for taking in the proportioned reads, and has been thoroughly benchmarked.

Additionally, we have a method Swish in the fishpond Bioconductor package which focuses on when there is substantial uncertainty in read assignment, as occurs in transcript level analysis or in gene level for scRNA-Seq.

ADD COMMENT
0
Entering edit mode

Thank you very much! Raya

ADD REPLY

Login before adding your answer.

Traffic: 668 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