User: jtwalker

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

Posts by jtwalker

<prev • 18 results • page 1 of 2 • next >
0
votes
1
answer
808
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 15 months ago by jtwalker20
0
votes
1
answer
808
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 15 months ago by jtwalker20
0
votes
1
answer
1000
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 16 months ago by jtwalker20
2
votes
1
answer
1.1k
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 17 months ago by jtwalker20 • updated 17 months ago by Pierre Lindenbaum116k
0
votes
1
answer
1000
views
1
answers
Comment: C: Badly formed genome location? GATK GenotypeGVCFs
... They were all created with the same method, namely GATK's HaplotypeCaller ...
written 17 months ago by jtwalker20
0
votes
1
answer
1000
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 17 months ago by jtwalker20
0
votes
1
answer
1000
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 17 months ago by jtwalker20
0
votes
1
answer
1000
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 17 months ago by jtwalker20
0
votes
1
answer
746
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 19 months ago by jtwalker20
0
votes
1
answer
884
views
1
answers
Comment: C: GATK variantFiltration unexpected behaviour?
... It didn't give me an error, but I'll try that now. Thanks! ...
written 19 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: 2404 users visited in the last hour