- Last Active
An update: I tried this with gatk 188.8.131.52 and it worked!
I got this error, when working with a manually generated fasta file. "ERROR MESSAGE: Bad input:We encountered a non-standard non-IUPAC base in the provided reference: '45'. It took a bit of panicked searching, but then I realised that the referenc…
Hi @bhanuGandham, The file is called twovcfs1.tar.gz. Thank you!
Hi @bhanuGandham, Adding --genomicsdb-segment-size 100000000 still gives the protocol buffer error. These g.vcf files were generated using gatk4's HaplotypeCaller. I'd be happy to share my g.vcf files, how should I go about sharing them? Will yo…
Hi @bhanuGandham Thanks for looking into this! Here are the first 10 lines for my interval.00.list: ctg7180024710462 ctg7180024952859 ctg7180025043953 ctg7180024953445 ctg7180024780449 ctg7180024953454 ctg7180024361837 ctg7180024716175 ctg71…
I have tried with the following: gatk GenomicsDBImport \ --java-options "-Xmx250G -XX:+UseParallelGC -XX:ParallelGCThreads=24" \ -V input.list \ --genomicsdb-workspace-path 5sp_45ind_assmb_00 \ --intervals interval.00.list \ --batch-size 9 w…
Hi Sheila, Thank you for the inspiration! I'll try with --batchSize and splitting up the intervals (which is a list of contigs) instead of trying to combine everything at one go. Will update with my progress Cheers, Quinn
Thank you, Mauro! Will upgrade to 184.108.40.206 and give it a try!
I am currently using GATK 220.127.116.11 and looking to analyse SNPs from a target capture dataset done on a non-model organism without linkage or chromosome information. There are 2900 targeted regions/contigs each across 45 individuals. I have run Haplot…