Reference span vs. width of read in riboseqR
0
0
Entering edit mode
ni2239 • 0
@ni2239-23489
Last seen 4.5 years ago

Hello,

I hope this is the right forum for my question or that someone can help at all. I'm somewhat new to analyzing NGS data in general, but I am trying to analyze some published ribo-seq data. I have aligned it myself and trimmed the reads by quality. When I view these alignments in IGV, for example, I can see that the reads cover a reference span that makes sense for ribo-seq data (between 25 and about 34 bp), but when I load this same BAM file into the readRibodata() function it produces read widths of only 50bp (presumably the pre-trimmed read length). I am unsure if I have done something wrong in my alignment procedure (done with bowtie2), or if there is a way to have readRibodata() utilize the CIGAR string or some other info from the BAM file. Any help is appreciated.

Here is what the code shows:

riboDat <- readRibodata('Ribobt2local1milsort.bam', replicates = 'LB') Reading ribosomal files....done! table(width(riboDat@riboGR$Ribobt2local1milsort.bam))

50

484584

riboSeqR • 612 views
ADD COMMENT

Login before adding your answer.

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