Getbamcounts In Exomedepth Package
when I run the function getBamCounts in ExomeDepth pakage , there is a error as follow. Can anyone tell me what is the problem?
H9.counts <- getBamCounts(bed.frame= exons.hg19,bam.files="H9_exome_b37.bam.pr efix.bam",include.chr=FALSE,referenceFasta="/rd/data/genomes/human/GRCh37/human_g1k_v37.fasta")
Reference fasta file provided so exomeDepth will compute the GC content in each window
Floating point exception
• 3,033 views
•
link
0 answers
No answers yet.
Log in to answer this question.
Is the issue coming from the GC content computation? If you do not specify the referenceFasta option at all, it will skip this step. Do you get the same error then?