User: jtwalker
jtwalker • 20
- 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
... 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
jtwalker • 20
0
votes
1
answer
416
views
1
answer
... 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 ...
written 7 months ago by
jtwalker • 20
0
votes
1
answer
553
views
1
answers
... 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
jtwalker • 20
1
vote
1
answer
455
views
1
answer
... 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 ...
written 9 months ago by
jtwalker • 20
• updated
9 months ago by
Pierre Lindenbaum ♦ 108k
0
votes
1
answer
553
views
1
answers
... They were all created with the same method, namely GATK's HaplotypeCaller ...
written 9 months ago by
jtwalker • 20
0
votes
1
answer
553
views
1
answers
... 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
jtwalker • 20
0
votes
1
answer
553
views
1
answers
... 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
jtwalker • 20
0
votes
1
answer
553
views
1
answer
... 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 ...
written 9 months ago by
jtwalker • 20
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
jtwalker • 20
0
votes
1
answer
558
views
1
answers
... It didn't give me an error, but I'll try that now. Thanks! ...
written 11 months ago by
jtwalker • 20
Latest awards to jtwalker
No awards yet. Soon to come :-)
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