Question: RNA Star problems
0
gravatar for fraheto
17 days ago by
fraheto • 10
fraheto • 10 wrote:

I am pretty new in Galaxy. I have been following some tutorials in order to perform some DEG with my data( https://galaxyproject.github.io/training-material/topics/transcriptomics/tutorials/ref-based/tutorial.html). Everything seems to be easy, however I am getting problems. I am trying to alinng single-ends reads from Illumina as follow:

Single-end or paired-end reads:
    Single-end reads in my history.
Custom or built-in reference genome:
    Use a built-in index.
Reference genome with or without an annotation:
    Use genome reference without builtin gene-model.
Select reference genome:
    Mouse (Mus musculus):mm10.
Gene model (gff3,gtf) file for splice junctions:
 ftp://ftp.ensembl.org/pub/release94/gff3/mus_musculus/Mus_musculus.GRCm38.94.gff3.gz.

I launch the program but sooafter I get this message:
This job was terminated because it used more memory than it was allocated.

I have tried to download the mouse genome from ENSMBL ftp://ftp.ensembl.org/pub/release-94/fasta/mus_musculus/dna/Mus_musculus.GRCm38.dna.primary_assembly.fa.gz, and select it to perform the alinment instead of Mouse (Mus musculus):mm10 in your database, but this time i get the mesasage: Fatal error: Matched on FATAL ERROR Fatal INPUT FILE error, no valid exon lines in the GTF file: /galaxy-repl/main/files/028/060

I have tried to transform Mus_musculus.GRCm38.94.gff3 file to gtf with gffread tool, since gft3 file from ENSEMBL is not recognized by RNA star. I have also tried to use mmu.gff3 file from miRbase (my reads come from small RNA libreries) and I get the same problem.

Could you be so kind to hellp me?

Thanks

ADD COMMENT • link • modified 16 days ago by Jennifer Hillman Jackson ♦ 25k • written 17 days ago by fraheto • 10
0
gravatar for Jennifer Hillman Jackson
16 days ago by
United States
Jennifer Hillman Jackson ♦ 25k wrote:

Hello,

The annotation needs to have chromosome identifiers that are a match for the genome used for mapping.

Try the GeneCode Genes GTF version of the annotation.

Full details are in this prior Q&A, including where to source the GTF and FAQs about how to detect and correct format/chromosome mismatch problems: https://biostar.usegalaxy.org/p/28878/

Thanks! Jen, Galaxy team

ADD COMMENT • link written 16 days ago by Jennifer Hillman Jackson ♦ 25k

Thank you so much. I will try it again taking into account your advices

ADD REPLY • link written 16 days ago by fraheto • 10

It worked perfectly. Thank you so much

ADD REPLY • link written 15 days ago by fraheto • 10
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 16.09
Traffic: 182 users visited in the last hour