Question: Extracting read counts from bam file for CDF calculation
gravatar for phdbiochemist
3.6 years ago by
phdbiochemist0 wrote:

Hi, I need some help if possible. I have nascent RNA-seq data (bam files) and would like to calculate an average CDF (cumulative distribution function) of reads downstream of the TSS, say from +1 to +500. I have no idea how to extract this information from a bam file. Once I get that I think I can use the ecdf function in R.

Thanks, Brian

rna-seq next-gen • 893 views
ADD COMMENTlink modified 3.6 years ago • written 3.6 years ago by phdbiochemist0

I'm sorry. I should add more detail about what I want to extract. I need the total counts for each nucleotide, from +1 to +500, which I should be able to use in R using ecdf.


ADD REPLYlink written 3.6 years ago by phdbiochemist0

Please choose a more descriptive title for your thread. "Principal Investigator and beginning bioinformaticist" doesn't say anything about the content of your question.

ADD REPLYlink written 3.6 years ago by WouterDeCoster45k
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.3.0
Traffic: 1136 users visited in the last hour