Entering edit mode
12.0 years ago
siyu
▴
150
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
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?