kvn95ss
- Total activity 46
- Last activity
- Member since
- Following 0 users
- Followed by 0 users
- Votes 1
- Subscriptions 17
Activity overview
Latest activity by kvn95ss-
kvn95ss commented,
Any suggestions guys?
-
kvn95ss created a post,
GenomicsDB and recalibration recommended for cohort of mixed exome
REQUIRED for all errors and issues: a) GATK version used: GATK 4.1.6.0b) Exact command used: Genomics DB We have a cohort of around 3k exomes, most of them from different capture kits, such as Agil...
-
kvn95ss commented,
Hay, since the number of samples is high, have you tried ReblockGVCF for the samples? Do look into it, it reduces memory and I/O requirements by quite a lot.
-
kvn95ss created a post,
GATK Haplotypecaller is calling Heterozygous variant, but bamout doesn't reflect the read frequencies in gvcf file
REQUIRED for all errors and issues: a) GATK version used: 4.2.6.1 and 4.3.0.0b) Exact command used: gatk HaplotypeCaller -R "$ref_hg38" -I input.bam -L twist.bed -ERC GVCF -ip 50 -O test_latest.gvc...
-
kvn95ss commented,
Hay, When using hg38, were the alternate contigs retained or removed? If they were retained, you could try removing them (or using a hg38 no alt reference) and then check. Is the data generated via...
-
kvn95ss commented,
Any updates on this? Am I doing something wrong?
-
kvn95ss created a post,
GatherVCFs error - provider not installed
I am trying to create a panel of normal following this guide - https://gatk.broadinstitute.org/hc/en-us/articles/360035531132--How-to-Call-somatic-mutations-using-GATK4-Mutect2 I have split the gen...
-
kvn95ss commented,
You can try writing the bam list into a file, then use parallel -a to use that list as input For Ex. ls *.sorted.bam > input_list # Not the best way to be frankfind $PWD -name "*.sorted.bam" > inpu...
-
kvn95ss commented,
You have got it correct. AFAIK genomicsDB only allows samples to be updated, but doesn't support removal / replacement of samples. A new genomicsDB would have to be created. Do check out the Rebloc...
-
kvn95ss commented,
Just a note, you could run the haplotype caller again with the switch --bam-output to understand which reads were used as evidence for the variant calling.