Suppose I have a bam file of paired-end reads and a gff file for exon locations. Is there an easy way to count the number of split reads (which can be identified by CIGAR string) and the number of paired-end reads that overlap any exon junctions? Assume we do not have any isoform information here...
Thanks a lot.
Thanks a lot.
Comment