Avatar

Juan Pablo Aguilar Cabezas

  • Total activity 19
  • Last activity
  • Member since
  • Following 0 users
  • Followed by 0 users
  • Votes 5
  • Subscriptions 5

Activity overview

Latest activity by Juan Pablo Aguilar Cabezas
  • Avatar

    Juan Pablo Aguilar Cabezas commented,

    Genevieve Brandt (she/her) You are welcome! Yes, I know that it is not an error but since I am doing demographic inference, and those sites will be skipped, right? Then, I will have less SNPs and v...

  • Avatar

    Juan Pablo Aguilar Cabezas commented,

    Sorry the command is: singularity exec /users/PHS0338/jpac1984/appz/gatk_latest.sif gatk --java-options "-Xmx145g" GenotypeGVCFs \   -R myse-hapog.fasta \   -V combine.8-03.vcf.gz \   -O gvcf.8-03....

  • Avatar

    Juan Pablo Aguilar Cabezas commented,

    Hi, my GATK version is v4.2.6.1Everything goes normal at the beginning... 06:32:00.602 INFO  NativeLibraryLoader - Loading libgkl_compression.so from jar:file:/gatk/gatk-package-4.2.6.1-local.jar!/...

  • Avatar

    Juan Pablo Aguilar Cabezas commented,

    Hi, I am in the same position. Not as many samples as Nils, but I have 37 samples that I put together with CombineGVCFs, and got the error message of more than 50 alleles when trying GenotypeGVCFs....

  • Avatar

    Juan Pablo Aguilar Cabezas created a post,

    Variant Filtering problem after merging

    REQUIRED for all errors and issues: a) GATK version used: 4.2.6.1b) Exact command used: java -jar picard.jar MergeVcfs  \          I=PA_113_DP.vcf.gz \          I=No_Pa90-16_DP.vcf.gz \          O=...

  • Avatar

    Juan Pablo Aguilar Cabezas created a post,

    GATK4.2.0.0 docker - Haplotypecaller optimization - how to activate Samtools asynchronous read/write

    Answered

    If you are seeing an error, please provide(REQUIRED) : a) GATK version used: Docker 4.2.0.0.b) Exact command used:singularity exec gatk_latest.sif \gatk --java-options "-DGATK_STACKTRACE_ON_USER_EX...

  • Avatar

    Juan Pablo Aguilar Cabezas commented,

    Genevieve Brandt (she/her) Do you know how to add it for the mapping step?I am trying the AddOrReplaceReadGroups (Picard) but since there is no multi-threading option it for one sample it is taking...

  • Avatar

    Juan Pablo Aguilar Cabezas commented,

    Hi, I am having the same issue, I found that I needed to add the @RG read group information but I am trying and I have found that some of the documents provide incomplete information. For example, ...

  • Avatar

    Juan Pablo Aguilar Cabezas created a post,

    Variant calling/ genotyping with a chromosome-scale genome assembly

    Hi, Hope you are having a great day. I am interested in conservation genomics and estimating recent demographic trends and long-runs-of-homozygosity.For this, I built a genome assembly using a refe...