how to merge multiple cufflink data in one file. what i am expecting is rows will represent gene name and column will represent the FPKM values. i tried cuffdiff too. and all the time i am getting error: Error: atal error: Matched on error [13:22:33] Loading reference annotation. [13:22:36] Inspecting maps and determining fragment length distributions. BAM record error: found spliced alignment without XS attribute BAM record error: found spliced alignment without XS attribute BAM record error: found spliced alignment without XS attribute BAM record error: found spliced alignment without XS attribute i am using mm10 gtf file for annotation from i genome. please let me know what is the fix! Please!
The root problem might be a genome mismatch problem. Was the data mapped within Galaxy against the mm10 natively indexed genome with HISAT? (or optionally, Tophat - which has been deprecated and is not recommended, or RNA STAR - which does not have indexes yet but will soon, by end of the week is the plan).
If the genome mapped against is not exactly the same as that in the other input datasets, this kind of result can occur (with any tool - or it may just manifest as an "exceeding compute resources" error).
CuffMerge or Stringtie Merge are the tools to use with Cufflinks/Stringtie output (gtf) and an optional reference GTF (example: iGenomes) to produce a merged GTF result.
Cuffdiff will give these warnings if the XS attribute is not present in the input BAM datasets (example: if Bowtie was used). Using HISAT will avoid the problem.
Support: https://galaxyproject.org/support/ (see Troubleshooting)
Tutorials: https://galaxyproject.org/learn/ (see RNA-seq)
Thanks! Jen, Galaxy team
Hi Jen Thanks for following up! i used mm10 iGenomes GTF file with the tool along with the natively indexed mm10 genome/and mm10 genome uploaded by me which is actually came from igenome. the problem remains. can i send you the link so that you can see what is happening! thanks
Firs of all Thank you Jennifer! yes i did alignment using in built mm10 genome and then did the cufflink using the igenome gtf file which i uploaded into the system. lets try to upload mm10 genome from i genome and use that as refference genome with HISAT. I am not using UCSC annotation file because that does not produce gene name! Thanks for your suggestion. very helpful! Amalendu