Question: generate table of counts with strand-specific paired-end RNA-seq data
gravatar for teresati
8 months ago by
teresati10 wrote:

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,




ADD COMMENTlink modified 8 months ago • written 8 months ago by teresati10
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.2.0
Traffic: 295 users visited in the last hour