Avatar

Shaun Clare

  • Total activity 16
  • Last activity
  • Member since
  • Following 0 users
  • Followed by 0 users
  • Votes 1
  • Subscriptions 4

Activity overview

Latest activity by Shaun Clare
  • Avatar

    Shaun Clare commented,

    I managed to fix it by modifying ls to write as one column using: ls -1 *.sorted.bam | <rest of command> I'm not sure why it worked in the first place, or maybe I accidentally deleted part of my co...

  • Avatar

    Shaun Clare created a post,

    Cannot parallelize GATK

    So I am rather confused what is going on, I am trying to use HaplotyperCaller to genotype ~1100 samples. I ran this is exact code on the entire set, getting through ~800 when it was actively workin...

  • Avatar

    Shaun Clare commented,

    Unfortunately since its 44k separate amplicons, --merge-input-intervals wasn't able to work since it can't detect if they were abutting. I used --merge-contigs-into-num-partitions 50 instead and I ...

  • Avatar

    Shaun Clare commented,

    It seems to be working but I'm guessing this is going to take a long time. It got passed initializing the engine but gave a warning to not use more than 100 intervals, this uses 44,000 amplicons as...

  • Avatar

    Shaun Clare commented,

    Our pipeline was built using CombineGVCFs along time ago so I was familiar with it. Isn't the intervals option required to make the GenomicsDB datastore? I just want all SNP across the whole 'genome'

  • Avatar

    Shaun Clare created a post,

    Issue with combining g.vcf files

    Answered

    I have an issue with combining my g.vcfs. It is 96 samples called against an amplicon based reference so it's approximately 38k (called the 50k) sequences of 250 bp or less. These codes previously ...

  • Avatar

    Shaun Clare commented,

    I do have another issue on another project. This one uses the exact same codes but it's on 96 samples called against a much smaller amplicon based reference.  gatk CombineGVCFs --java-options "-Xmx...

  • Avatar

    Shaun Clare commented,

    Genevieve Brandt (she/her) Thank you, that seems to have worked!

  • Avatar

    Shaun Clare created a post,

    CombineGVCFs

    Answered

    So I noticed I was having trouble combining my g.vcf files, which is saying my index is out of bounds. My HaplotypeCaller command seemed to work fine and all of these codes work fine when I use amp...

  • Avatar

    Shaun Clare commented,

    Running GATK worked fine and that code had worked previously. I couldn't find the error so I rebuilt my Ubuntu from scratch a while ago and it fixed everything so the code works again. Thanks, Shaun