User: jtwalker

gravatar for jtwalker
jtwalker20
Reputation:
20
Status:
New User
Location:
Last seen:
7 months, 2 weeks ago
Joined:
1 year, 1 month ago
Email:
j*******@iastate.edu

Posts by jtwalker

<prev • 18 results • page 1 of 2 • next >
0
votes
1
answer
416
views
1
answers
Comment: C: Speed up GATK joint genotype calling?
... I had considered this, but the gatk doc says that we'd need to have the same amount of ram available for each thread. As our ram usage is already a concern, would this still be a viable option? ...
written 7 months ago by jtwalker20
0
votes
1
answer
416
views
1
answer
Speed up GATK joint genotype calling?
... I'm using GATK's GenotypeGVCFs tool to jointly genotype ~1000 samples. It's my understanding that because of the genome wide annotations that are calculated, I can't speed things up by using CombineVCFs on smaller jointly called groups. Am I correct? Is there some way to speed up my joint genotypin ...
gatk joint genotyping written 7 months ago by jtwalker20
0
votes
1
answer
553
views
1
answers
Answer: A: Badly formed genome location? GATK GenotypeGVCFs
... Just to update this question in case anyone is having a similar problem - for some reason when I used BWA mem the output for a few of my files was malformed, which led to problems down the line. I simply realigned the files and my joint call worked. ...
written 8 months ago by jtwalker20
1
vote
1
answer
455
views
1
answer
VCF format: what does AD 1,0 denote?
... I'm looking at one of my VCF files and I noticed that allele depth (AD) isn't just one number, but rather two values seperated by a comma e.g.: GT:AD:DP:GQ:PL 0/0:**1,0**:1:3:0,3,27 The [GATK doc](https://gatkforums.broadinstitute.org/gatk/discussion/1268/what-is-a-vcf-and-how-should-i-i ...
vcf written 9 months ago by jtwalker20 • updated 9 months ago by Pierre Lindenbaum108k
0
votes
1
answer
553
views
1
answers
Comment: C: Badly formed genome location? GATK GenotypeGVCFs
... They were all created with the same method, namely GATK's HaplotypeCaller ...
written 9 months ago by jtwalker20
0
votes
1
answer
553
views
1
answers
Comment: C: Badly formed genome location? GATK GenotypeGVCFs
... cut -f 1 S1097.g.vcf | grep -F 024218.1 | uniq ##contig= NC_024218.1 I have over a thousand g.vcfs. Could it be that one of them contains a malformed contig name causing the problem? ...
written 9 months ago by jtwalker20
0
votes
1
answer
553
views
1
answers
Comment: C: Badly formed genome location? GATK GenotypeGVCFs
... GCF.dict:@SQ SN:NC_024218.1 LN:77392008 M5:0ad5ac74565c0b48329eec6020994b16 UR:file:/Users/lukehoekstra/GenomeAnalysisTK-3.7/GCF.fna GCF.fna.fai:NC_024218.1 77392008 125 80 81 ...
written 9 months ago by jtwalker20
0
votes
1
answer
553
views
1
answer
Badly formed genome location? GATK GenotypeGVCFs
... I'm trying to genotype >1000 samples using the GATK pipeline. I've already created gVCFs for my samples, but when attempting to use GATK's GenotypeGVCFs tool, I get the following error BiasBySample values to calculate strand bias annotation values. If no sample has the SB genotype annotatio ...
gatk joint genotyping written 9 months ago by jtwalker20
0
votes
1
answer
495
views
1
answers
Comment: C: Regarding vcf output
... You might also consider filtering on QD rather than QUAL and DP: https://software.broadinstitute.org/gatk/documentation/article.php?id=6925 ...
written 11 months ago by jtwalker20
0
votes
1
answer
558
views
1
answers
Comment: C: GATK variantFiltration unexpected behaviour?
... It didn't give me an error, but I'll try that now. Thanks! ...
written 11 months ago by jtwalker20

Latest awards to jtwalker

No awards yet. Soon to come :-)

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 690 users visited in the last hour