Which fasta file for Variant Calling using RNASeq

Hey everyone,

I want to perform variant calling using transcriptomic samples. I will use the recommended pipeline from GATK (https://gatk.broadinstitute.org/hc/en-us/articles/360035531192-RNAseq-short-variant-discovery-SNPs-Indels).

For the 2-pass alignment step with STAR, I was considering to use the same genome that I used for the gene counts step (dna.primary.assembly), but I got the idea this may not be ideal - Fasta file and GTF file for STAR alignment.

I wanted to confirm if I'm using the right file or not.

Thanks in advance.

rna-seq variant-calling

• 537 views • link

updated 2.4 years ago by Ram

• written 2.4 years ago by andrebolerbarros •

0 answers

No answers yet.

More posts like this

Can MarkDuplicates of Picard be used for RNA reads?
written 4.0 years ago by Nemo •

Hello, I have bam files from RNA sequence data. I am following the pipeline of gatk in [Variant calling in RNA sequences][1]. In the second …
Variant calling for matched samples
written 4.6 years ago by shivangi.agarwal800

Hi Guys, I want to call variants for 36 samples having matched normal. Can you suggest me how to run (what are the steps) GATK …
Why GATK4 doesn't carry out local realignment around indels ?
written 4.6 years ago by ManuelDB

I am a trainee working with GATK. In the clinical lab where I am working, we use GATK3 and I am learning creating my own …
What is GenotypeGVCFs?
written 4.7 years ago by wormball •

Hello! This article https://gatk.broadinstitute.org/hc/en-us/articles/360035535932-Germline-short-variant-discovery-SNPs-Indels- says i should use HaplotypeCaller in GVCF mode and GenotypeGVCFs then, and this article https://gatk.broadinstitute.org/hc/en-us/articles/360035531192-RNAseq-short-variant-discovery-SNPs-Indels- advises to use HaplotypeCaller without GenotypeGVCFs. …
no variant remained after GenotypeGVCFs gatk4
written 5.0 years ago by Mehmet

Dear all, I am analyzing a DNA-seq data to identify variants in two genes of human using a single-sample based gatk4 pipeline (https://gatk.broadinstitute.org/hc/en-us/articles/360035535932-Germline-short-variant-discovery-SNPs-Indels-). Haplotypecaller identified …
GATK GenomicDBImport missing argument --intervals for Illumina HiSeq Exome??
written 5.8 years ago by yussab •

Hi everyone, I'm using the GATK germline short variant discovery (SNPs + Indels) pipeline. https://gatk.broadinstitute.org/hc/en-us/articles/360035535932-Germline-short-variant-discovery-SNPs-Indels- After using HaplotypeCaller to call the GVCF, I'm not consolidating …
Does GATK provide datasets to test their pipelines? (Learning purpose)
written 5.9 years ago by yussab •

Hi Everyone, I am working in the oncology field and I would like to test this pipeline . **GATK -Germline short variant discovery (SNPs + …
GATK for RNA seq variant calling
written 6.0 years ago by danson.wooi •

I am very new to this but I was just wondering what commands do people normally use with HaplotypeCaller for variant calling with GATK on …
How can I call both somatic and germline variants from bam file?
written 6.1 years ago by svp

I have 'analysis ready bam file created using 'Best Practice Workflows'. My next plan is to get both somatic and germline mutation from bam file. …
Mutect2 somatic hemizygote variant on chrX
written 6.2 years ago by cocchi.e89

I am working on calling some somatic variants from WES samples, I followed step-by-step the [GATK guidelines][1]. This is not a tumor sample so I …