I have a problem in TopHat-Cufflinks pipeline. I want to get splicing information of a gene family, but the annotation is bad so far, that is, nothing is annotated or huge exons are annotated.
This family has a number of processed pseudogenes with almost same sequence. It seems that reads that are supposed to be mapped onto "real" genes are mapped onto these pseudogenes.
Now, I have a gff file defining the location of those pseudogenes in the genome. Is there a way to mask these pseudogenes in TopHat mapping?
I do not think -M option in Cufflinks works for me, because the reads are already mapped in .bam file from TopHat.
I really need help. If someone knows, please tell me
This family has a number of processed pseudogenes with almost same sequence. It seems that reads that are supposed to be mapped onto "real" genes are mapped onto these pseudogenes.
Now, I have a gff file defining the location of those pseudogenes in the genome. Is there a way to mask these pseudogenes in TopHat mapping?
I do not think -M option in Cufflinks works for me, because the reads are already mapped in .bam file from TopHat.
I really need help. If someone knows, please tell me
Comment