cufflinks doesn't process the loci, all the output files are empty
0
0
Entering edit mode
9.0 years ago
aleka ▴ 110

Hi,

I am running paired alignments with tophat and when I run cufflinks

./cufflinks -o emrsa15.cuff/ -g annotation.gtf accepted_hits.bam

cufflinks doesn't process the loci, so all the output files are empty (only header is appeared)

Cufflinks output

 Inspecting reads and determining fragment length distribution.
> Processed 1 loci.                            [*************************] 100%
> Map Properties:
>       Normalized Map Mass: 19021202.87
>       Raw Map Mass: 19021202.87
>       Fragment Length Distribution: Empirical (learned)
>                     Estimated Mean: 213.28
>                  Estimated Std Dev: 91.26

Here is the summary of the alignment, which seems fine.

Left reads:
          Input     :  18922184
           Mapped   :  18759527 (99.1% of input)
            of these:     41369 ( 0.2%) have multiple alignments (1463 have >20)
Right reads:
          Input     :  18922184
           Mapped   :  18494767 (97.7% of input)
            of these:     41039 ( 0.2%) have multiple alignments (1463 have >20)
98.4% overall read mapping rate.

Aligned pairs:  18412669
     of these:     40472 ( 0.2%) have multiple alignments
                  143680 ( 0.8%) are discordant alignments
96.5% concordant pair alignment rate.

When I run tophat with one end, cufflinks works fine. But not with the paired ends.

Do you have any suggestion why is that happening or how can I merge the results of cufflinks from the two ends?

I would appreciate any help.

SNP next-gen-sequencing alignment • 2.0k views
ADD COMMENT

Login before adding your answer.

Traffic: 2723 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6