ı didnot do when ı align rnaseq datasets but I did it as soon as I realized I had to do it. ı am following https://gatk.broadinstitute.org/hc/en-us/articles/360035531192-RNAseq-short-variant-discovery-SNPs-Indels but here is not exact information.
Read group addition order in RNA-seq short variant discovery (SNPs + Indels)
Hello, I’m following the GATK best practices for RNA-seq short variant discovery (SNPs + Indels) and wondering about the correct point to add Read Groups (RGs). Does the order (before/after MarkDuplicates or SplitNCigarReads) matter for RNA-seq variant calling with GATK (HaplotypeCaller)? Any official clarification or reference from the GATK team or papers? Pipeline: HISAT2, AddOrReplaceReadGroups, MarkDuplicates, SplitNCigarReads, BaseRecalibrator, HaplotypeCaller Thank you for any kind of help
1 answer
Ideally you would do this at the beginning of the process as you align the data. Read groups would be tagging the samples with relevant info e.g. what lane they were in, if they ran on multiple flowcells etc (all of this may not be applicable in your case).
See Expected Input section: https://gatk.broadinstitute.org/hc/en-us/articles/360035535912-Data-pre-processing-for-variant-discovery
yes this is exactly what ı have done but anway ı couldnot find any clear reference for exact pipeline .so actually ı write here maybe someone can help for reference , because nothing clear when ı search on google. thank you for help, ı also feel ı am right but ı amnot sure anyway without papers.
https://gatk.broadinstitute.org/hc/en-us/articles/360035531192-RNAseq-short-variant-discovery-SNPs-Indels references the https://gatk.broadinstitute.org/hc/en-us/articles/360035535912-Data-pre-processing-for-variant-discovery pipeline linked above. In RNAseq pipeline STAR is replacing default aligner (to account for splicing) so the rest of the steps from there on should be as per the input requirements for DNA pipeline.
We use MergeBamAlignment and MarkDuplicates (similarly to our DNA pre-processing best practices pipeline)
Log in to answer this question.
More posts like this
-
GATK Mutect2
written by sk_24 •
-
Does the order of SplitNCigarReads and MarkDuplicates affect RNA-seq variant calling results?
written by iamsmor •
-
GATK BQSR error — Reference and BAM file chromosome name mismatch (“chr” vs. no “chr”)
written by iamsmor •
-
Seeking Advice on Multisample Variant Calling Pipeline with BWA-MEME and GATK DRAGEN Mode
written by George •
-
HaplotypeCaller: GVCF-Mode or VCF for variant calling of single sample?
written by DavidStreid •
-
no variant remained after GenotypeGVCFs gatk4
written by Mehmet
-
Discussion: Clinical genetics germline variants detection pipeline
written by Jordi •
-
GATK markduplicates / applyBQSR / haplotypecaller tools support unix pipe?
written by roccostorace •
-
GATK GenomicDBImport missing argument --intervals for Illumina HiSeq Exome??
written by yussab •
-
Variant analysis on 1 or 2 samples: are the final steps of the GATK pipeline unneccessary?
written by BioinfGuru
