error in intron clustering
0
0
Entering edit mode
3.7 years ago

Hello

I am trying to use leafcutter for differential alternative analysis. As per its requirement, I have generated bam files using STAR.

Now with leafcutter, I have generated bam.juncs and test_juncfiles.txt files successfully using script:

    for bamfile in `ls run/geuvadis/*.bam`
do
    echo Converting $bamfile to $bamfile.junc
    sh ../scripts/bam2junc.sh $bamfile $bamfile.junc
    echo $bamfile.junc >> test_juncfiles.txt
done

Now the next step is intron clustering and differential analysis.

But while doing intron clustering with the script leafcutter_cluster.py, I am not getting numers.counts.gz files with the format as given in the documentation.

According to documentation file should be in the following format:

RNA.NA06986_CEU.chr1.bam RNA.NA06994_CEU.chr1.bam RNA.NA18486_YRI.chr1.bam RNA.NA06985_CEU.chr1.bam RNA.NA18487_YRI.chr1.bam RNA.NA06989_CEU.chr1.bam RNA.NA06984_CEU.chr1.bam RNA.NA18488_YRI.chr1.bam RNA.NA18489_YRI.chr1.bam RNA.NA18498_YRI.chr1.bam
chr1:17055:17233:clu_1 21 13 18 20 17 12 11 8 15 25
chr1:17055:17606:clu_1 4 11 12 7 2 0 5 2 4 4
chr1:17368:17606:clu_1 127 132 128 55 93 90 68 43 112 137

But, Generated counts.gz files contains only my samples name, no other information is present.

Please suggest.

leafcutter intron clustering • 1.1k views
ADD COMMENT
1
Entering edit mode

It would be useful to provide a link for where you got the leafcutter program from (or the package it is part of) and the exact command line you are using.

ADD REPLY
0
Entering edit mode

thanks for the guidance.

link for the leafcutter: http://davidaknowles.github.io/leafcutter/

ADD REPLY

Login before adding your answer.

Traffic: 2645 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