This is a test version of Biostars. For the public version, visit https://www.biostars.org.
which file I can use for cufflinks from the aligner STAR

Hello,

I have three files generated from the software of STAR after mapping: 2-pass_Aligned.out.bam, 2-pass_Aligned.sortedByCoord.out.bam, 2-pass_Aligned.toTramscriptome.out.bam.

Which one I can use for cufflilnks?

Thanks in advance!

rna-seq
• 2,006 views • link
updated by zx8754 zx8754 • written by yueli7 yueli7

Please post your code for generating these bams with aligner.

• 1 views • link
written by cpad0112 cpad0112

Hi yueli7,

Please follow up on all your previous questions and provide feedback.

If an answer was helpful you should upvote it, if the answer resolved your question you should mark it as accepted.
Upvote|Bookmark|Accept

Note that future questions might be closed if you do not provide feedback on older threads.

Cheers,
Wouter

• 1 views • link
written by WouterDeCoster WouterDeCoster

1 answer

You'll want to use the 2-pass_Aligned.sortedByCoord.out.bam file. The transcriptome file is, as you guessed by the name alignments to a transcriptome, which is completely useless if you're trying to do de novo discovery. The 2-pass_Aligned.out.bam file is the same as the first one I mentioned, just not sorted. You can delete this and the transcriptome BAM file.

Also, use stringTie rather than cufflinks. Cufflinks should never be used for new projects.

• 1 views • link
written by Devon Ryan Devon Ryan

Thank you so much for your great help!

• 0 views • link
written by yueli7 yueli7

Log in to answer this question.

More posts like this
  • software or R package to change protein *.msf file to any readable file
    written by yueli7 yueli7
    Hello I have some protein *.msf file. I can not open it. Which software or R package that can change this kind of *.msf to …
  • SILAC protein data analysis
    written by yueli7 yueli7
    Hello, I have SILAC protein data from three patients. How can I normalize the A10: Light/Heavy value over these patients? Any available software or tools? …
  • SILAC protein data analysis
    written by yueli7 yueli7
    Hello I am new in analysis SILAC protein data. I have three patients samples which collected from three different points of time, and no replicates. …
  • which gtf file I can use?
    written by yueli7 yueli7
    Hello, I use star to make human index. I have hg38.fa file which downloaded from: http://hgdownload.soe.ucsc.edu/goldenPath/hg38/bigZips/ There are four gtf files:`hg38.refGene.gtf.gz, hg38.ncbiRefSeq.gtf.gz, hg38.knownGene.gtf.gz and hg38.ensGene.gtf.gz …
  • software of single-cell RNA-seq from fastq or fasta
    written by yueli7 yueli7
    Hello, Anyone know which software or pipeline can treat the single-cell RNA-seq data from the first beginning? That means from the form of single or …
  • zUMIs software usage
    written by yueli7 yueli7
    Hello, Everyone, I tried to use zUMIs. Thanks in advance for any help! Yue EXITING because of FATAL ERROR: Genome version: 20201 is INCOMPATIBLE with …
  • software of duplicates in RNA-seq data
    written by yueli7 yueli7
    Hello, I have already used star to mapping all the reads to the human genome. Which software has the best perform in detecting duplicates in …
  • differential express of RNA seq
    written by yueli7 yueli7
    Hello, This is my first time analysis RNA seq data. I have already used STAR, featureCounts. Then, I will try to find the differential expressed …
  • software after human STAR mapping
    written by yueli7 yueli7
    Hello, I have already finished mapping reads to the huamn genome by STAR. Which software I can use for transcript abundance estiamtes or Alignment/Junction? Cufflinkd, …
  • could not open genomeParameters.txt
    written by yueli7 yueli7
    Hello, By using STAR, I have already successfully generated genome indexes. Apr 26 05:20:22 ..... finished successfully Then I try to running mapping. It said …