Hi,
I have some rnaseq data from xenograft models. I am aligning the reads to human genome and mouse genome separately.
the idea is to throw the reads that are mapped to both human and mouse. How can I achieve this? How can I identify reads that are common to both the human and mouse genome when I run tophat (humanaligned.bam) and mousealigned.bam), then discard these reads and then create a file called "humanminuscommonreads.bam" for the rest of the analysis steps?
Please help.
I have some rnaseq data from xenograft models. I am aligning the reads to human genome and mouse genome separately.
the idea is to throw the reads that are mapped to both human and mouse. How can I achieve this? How can I identify reads that are common to both the human and mouse genome when I run tophat (humanaligned.bam) and mousealigned.bam), then discard these reads and then create a file called "humanminuscommonreads.bam" for the rest of the analysis steps?
Please help.