small .bam files to test GATK Variant Discovery pipeline
Could someone refer me to where I could download the smallest possible .bam files to test my GATK Best Practices Variant Discovery pipeline? My pipeline uses the -L option to parallelize over different chromosomes - I would like to test this functionality and so I would like a full .bam file that has data from all chromosomes and will not cause GATK to crash.
Thanks a lot!
• 3,052 views
•
link
written
by
alongalor •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
GATK - Output differences in Calls when Split by chromosome intervals and without splitting
written by Maverick •
-
Seeking Advice on Multisample Variant Calling Pipeline with BWA-MEME and GATK DRAGEN Mode
written by George •
-
Problem with Calling Variants from RNA-Seq data
written by Esraa •
-
no variant remained after GenotypeGVCFs gatk4
written by Mehmet
-
GATK variant calling error
written by endretoth •
-
Does GATK provide datasets to test their pipelines? (Learning purpose)
written by yussab •
-
RNA-Seq Variant Calling
written by user5212 •
-
GATK Variant calling
written by jaqx008
-
How to do a local denovo assembly including unmapped paired reads for many samples to genotype a la…
written by William
-
running GATK or Mutect or a somatic variant caller individual chromosomes
written by Bogdan
Just use any BAM that you have on disk, make a little BED file with one interval per chromosome, e.g. chr1-22 from 6000000-6100000 respectively, and use SAMtools view to get a subset of the whole BAM:
That should be sufficient for testing purposes.
This worked perfectly in terms of splitting the bam file but it caused GATK to crash with a strange error... my pipeline of course still works with the original bam file before it was split.