Test-drive the GATK tools and Best Practices pipelines on Terra
Check out this blog post to learn how you can get started with GATK and try out the pipelines in preconfigured workspaces (with a user-friendly interface!) without having to install anything.
Sequence dictionaries are not the same size (613, 0) when trying to merge BAMs with MergeSamFiles
I get the error
Sequence dictionaries are not the same size (613, 0) when trying to merge 4 BAMs using the following expression:
SAMPLE_BAMs_RG=$(ls -1 $TMP/*.bam | grep $SAMPLE | sed 's/\/projekte/I=\/projekte/g') #this adds "I=" before the BAM path/name java -jar $PICARD/picard.jar MergeSamFiles \ $(echo $SAMPLE_BAMs_RG) \ # passing all BAMs at once O=$TMP/$SAMPLE.RG.merged.sorted.bam \ SORT_ORDER=coordinate
I am not sure what I am missing since there is no flag to specify the dictionary path.
Thanks in advance