This site is now read-only. You can find our new documentation site and support forum for posting questions here.
Be sure to read our welcome blog!
guidelines for BAM realignment in GATK 4
I have aligned my sequences with BWA and HISAT2, I would like to perform the realignment step. I find confusing reports over the internet on whether performing it or not. Can I ask if it is still a good practice to perform re-alignment and how can I run it?
I created a dictionary with
java -jar picard.jar CreateSequenceDictionary R=<ref.fa> o=<ref.dict>
But I do not have a vcf file. The syntax I have is:
java -jar GenomeAnalysisTK.jar \
-T RealignerTargetCreator \
-R <ref.fa> \ # in the same folder there is the dict file
(the -known option is skipped) followed by
java -Xmx8G -Djava,io,tmpdir=/tmp -jar /
-T IndelRealigner \
-R <ref.fa> \
-I <file_dedup.bam> \
What would be the new syntax for GATK4?