Ways to make a variant calling for RNA Seq (paired-end)

Hello everyone! I’m a newer in bioinformatic and galaxy - sorry :mask:
I trying to use RNA STAR tool. I have full human transcriptom (RNAseq) in 8files (with R1 and R2 reads). So I have to get names of genes and their characteristics (pathogenic etc) in the end for clinicians.
I want to create vcf file for 1-snpsift variant type and 2-snpsift annotate SNPs from dbsnp. And after this I wanted download this vcf file (for ClinVar annotation etc).

I dont need the Trimmomatic so I was trying to do next:
Fastq–>RNA-STAR(hg38)–>Samtoolsmerge(for mapped.bam)—>FreeBayes

But FreeBayes show me the error:

  • { “code_desc”: “”, “desc”: “Fatal error: Exit code 101 ()”, “error_level”: 3, “exit_code”: 101, “type”: “exit_code” }
    Job API ID: bbd44e69cb8906b5339b2fad2d9562d3

Please explain me what I do wrong and which tools I need to use for my purpose?

Hi @Sofi

You should at least run FastQCand probably Fastq Info on your read data to make sure there is not a content/format problem.

These prior topics may also be helpful: Search results for 'FreeBayes order:likes' - Galaxy Community Help

Does the Dataset Details view contain any other informative messages? If not, check the prior topics for common help with this protocol, along with the tutorials below.

If you need more help, please post back a share link to the history containing this error. Please leave all inputs and outputs undeleted, along with the outputs from the QA checks. Sharing your History


Tutorials