BOL

#Intersect your bed data with coverage $ bedtools genomecov -ibam aln.bam -bga | awk '$4==0' | bedtools intersect -a cds.bed -b - -f 1.0