Problem bam files 0 counts
Hello, I am trying to get the counts of my experiment from the .bam files by using summarized experiment function in R.
My reads are single ends, not strand specific. Here is some R code:
gtffile <- file.path(dir2,"Homo_sapiens.GRCh38.91.gtf")
(txdb <- makeTxDbFromGFF(gtffile, format="gtf"))
(genes <- exonsBy(txdb, by="gene"))
library("GenomicAlignments")
se <- summarizeOverlaps(features=genes, reads=bamfiles,
mode="Union",
singleEnd=FALSE,
ignore.strand=FALSE,
fragments=FALSE )
I am getting 0 counts as the output for every sample and every gene. Anyone knows if it could be the .gtf file or if my .bam files are not correct? Ive done it many times. Any alternatives?
• 2,281 views
•
link
written
by
Pin.Bioinf
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
I'm not sure what is going on with my DEXSeq object
written by cftattersfield •
-
How to use NCBI annotation instead of Ensembl in analysing Transcriptome from RNAseq data
written by Vojtěch •
-
Help with summarizeOverlaps function in RNASeq analysis using R
written by FJCF •
-
warning when using GenomicAlignments
written by donny.dw •
-
Counting reads over genes with summarizeOverlaps
written by concetta •
-
GTF file problem
written by johnsonn573 •
-
HTSeq count versus summarizeOverlaps, mismathc of exon counts
written by Gama313
-
chromosome name in summarizeOverlaps for RNASeq data
written by BioProg •
-
summarizeOverlap in systemPipeR not working!
written by Parham
-
Error when counting reads in genes with summarizeOverlaps (Genomic-Alignments package)
written by alejandro.colaneri •
Please use the formatting bar (especially the
codeoption) to present your post better. I've done it for you this time.Didn't you get any error? What is the output of
samtools idxstats file.bam?I get :
Please use the formatting bar!
Then you should index
24_control.bamby doingsamtools index 24_control.bam.Top thing to check when zero counts manifest is to ensure that chromosome names match in your alignments and annotation (i.e.
chr1 =\= 1).Thank you a lot, i was helping another person who had stumbled upon this problem, and i realized she was using different genome anotations for mapping and for quantification.