I'm trying to use my
.sam file from
GTF file with
bedtools complement to get the following the number of intergenic reads, the number of genic reads, and the number reads that didn't map at all.
How can I do this? I saw this post which suggested some other tools but I need to do a lot of preprocessing for RNA-SeqQC and would like to use bedtools (and any downstream tool that would be necessary for this process): Does STAR give the reads that were mapped at intergenic regions? How can I get the regions and count reads?
Does anyone know how I would do this?
I converted my sam file to a bed file and ran the following:
bedtools -i mapped.bed -g contig.lengths.tsv
but it doesn't know where the intergenic regions are...