Calculating Fpkm/Rpkm For Each Read Group Given A Gtf File
I have a couple of BAM files, each containing several read groups, and an associated GTF file. I'd like to calculate the FPKM/RPKM's for each transcript for each read group. One of the files is from single-end reads, the other from paired-end reads.
Is there an existing tool that can calculate the FPKM/RPKMs for each individual read group in a BAM file?
• 4,584 views
•
link
updated
by
Biostar
•
written
by
Nathan S. Watson-Haigh
•
written
by
Nathan S. Watson-Haigh
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Paired-End Bam Files
written by sarahmanderni
-
Removing Phix Data From Ngs Data For Analysis
written by crivenster
-
Samtools Merge/Reheader Of Bam Files With Different @Sq Lines
written by Nathan S. Watson-Haigh
-
Cufflinks Paired End Read Names
written by Frenkiboy
-
Merging Bam Files
written by rickyflintoff •
-
How Do I Pull Out One Read (Given Its Id) From A Bam?
written by zam.iqbal.genome
-
How To Fix A Bam-Readcount Sm Error: "Couldn'T Grab Single-End Mapping Quality For Read"
written by amnonim •
-
Should Read Names Be Unique Over Multiple Readgroups In A Bam?
written by William
-
Problem With Sorting Bam File: Samtools
written by lyfsa •
-
How To Read Bam.Ebd, Bam.Ebi Format Dataset
written by Lerong •
I don't know of any tool that can do it at the read group level, if your bam file contains several read groups. You will have to separate the big bam into several small ones (one for each read group) and process them individually.
You need a way of identifying the read groups, are they annotated in the BAM file (should use
@RGtag). Samtools has the-R/rparameters:maybe that can be used to produce a set of BAM files for which to calculate FPKM separately?