Hello, I am having a problem whereby Cufflinks sometimes doesn't generate models in regions with abundant reads (while performing de novo gene model prediction), but does in other scenarios in which reads are even less abundant.
I am using cufflinks v2.0.2 and running the following command:
cufflinks -j 0.3 --library-type ff-firststrand -o TestRegion testregion.bam
As you can see from the attached image, cufflinks generates a model for several lowly abundant regions (blue and red regions in the "TestRegion" track), but completely skips a region with high read coverage on the + strand (the orange region in the "TestFwd+" track).
Does anybody have an idea of why regions can be skipped like this?
Any help would be greatly helpful, thanks!
I am using cufflinks v2.0.2 and running the following command:
cufflinks -j 0.3 --library-type ff-firststrand -o TestRegion testregion.bam
As you can see from the attached image, cufflinks generates a model for several lowly abundant regions (blue and red regions in the "TestRegion" track), but completely skips a region with high read coverage on the + strand (the orange region in the "TestFwd+" track).
Does anybody have an idea of why regions can be skipped like this?
Any help would be greatly helpful, thanks!
Comment