generate table of counts with strand-specific paired-end RNA-seq data
Entering edit mode
teresati ▴ 10
Last seen 6.7 years ago

Dear All,

I have strand-specific paired-end RNA-seq data, and I am attempting to generate a table of reads counts. I wanted to ask if the following steps are correct:

1. load separately bam files with BamFile function and specify asMates=TRUE for paired-end reads

bamflU1 <- BamFile("PATH/U1.bam", asMates=TRUE)

2. transform the file in a readGAlignmentPairs object and specify strandMode=2 for strand specificity

rGAP_U1 <- readGAlignmentPairs(bamflU1, use.names=TRUE, strandMode=2)

3. use summarizeOverlaps function for each bam file separately after loading the gtf file

countsU1 <- summarizeOverlaps(id_gene, rGAP_U1,mode="Union",singleEnd=FALSE, ignore.strand=FALSE)

4. create an unique object containing all the samples counts

Is there another way to import directly all the bamfiles and define them as paired-end and stranded? (it seems that the BamFilesList function is not compatible with the readGAlignmentPairs object)

Thank you for any help and suggestion,




summarizeoverlaps paired-end reads readGAlignmentPairs • 1.1k views

Login before adding your answer.

Traffic: 505 users visited in the last hour
Help About
Access RSS

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

Powered by the version 2.3.6