Question: Cufflinks On Chromosome 2
0
gravatar for komal.rathi
5.2 years ago by
komal.rathi3.4k
Children's Hospital of Philadelphia, Philadelphia, PA
komal.rathi3.4k wrote:

Hello everyone,

I have a sorted and index BAM file. I want to run cufflinks on just Chromosome2 from it to save time. How can I do it?

I was using:

cufflinks -p 2 -g <chr2 filtered GTF file> -o outputdir input_filter_sort.bam 2> input.out

I dont know whether this is the best way to do it. I also want to assemble novel transcripts along with the known ones on chromosome2 so I am using the -g option instead of -G. But this command assembles novel transcripts on all the chromosomes not just on chromosome 2. Is there a way to assemble known and novel transcripts on just one chromosome?

Thanks

cufflinks • 1.6k views
ADD COMMENTlink modified 5.2 years ago by Vivek2.2k • written 5.2 years ago by komal.rathi3.4k
2
gravatar for Vivek
5.2 years ago by
Vivek2.2k
Denmark
Vivek2.2k wrote:

You can use samtools to just extract chromosome 2 from your bam and use cufflinks on it.

samtools view -bh allChr.bam chr2 > onlyChr2.bam
ADD COMMENTlink written 5.2 years ago by Vivek2.2k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 988 users visited in the last hour