Biostar Beta. Not for public use.
cufflinks large bundles warning
0
Entering edit mode
18 months ago
Sam • 110

Dear All

I used star to align the PE lib with this command and then I converted sam file to bam and sorted it with samtools but after running the cufflinks I got "Warning: Skipping large bundle" how I can fix it ?

Thanks

STAR --runThreadN 12 --genomeDir  /genome_index_star -sjdbGTFfile /genome_index/Egene_exons.gtf --readFilesIn R1.fastq,R2.fastq --outFileNamePrefix test  --outSAMstrandField intronMotif

cufflinks -o cufflink_out_gtf -g ./Egene_exons.gtf test_Aligned.out_sorted.bam
cufflinks star • 346 views
ADD COMMENTlink
0
Entering edit mode

Try to increase the --max-bundle-length in cufflinks command, try 10000000

similar error

cufflinks options

ADD REPLYlink
0
Entering edit mode

but I tried also(the same lib ) with Tophat2 and Hisat2 and cufflinks worked well without any warning and error. should I change any flag in the star to be compatible with cufflinks?

ADD REPLYlink
1
Entering edit mode

In STAR manual

In addition, it is recommended to remove the non-canonical junctions for Cufflinks runs using --outFilterIntronMotifs RemoveNoncanonical.

jM jI attributes require samtools 0.1.18 or later, and were reported to be incompatible with some downstream tools such as Cufflinks.

I would say also that what you got is just a warning not an error

ADD REPLYlink
0
Entering edit mode

Thanks for your help, but with this warning, the processing time of cufflinks increased at least 10 times.

ADD REPLYlink
0
Entering edit mode

Did you try the --max-bundle-length and the --outFilterIntronMotifs RemoveNoncanonical ?

ADD REPLYlink

Login before adding your answer.

Similar Posts
Loading Similar Posts
Powered by the version 2.3.1