User: fwuffy

gravatar for fwuffy
fwuffy50
Reputation:
50
Status:
Trusted
Location:
Michigan, USA
Website:
http://glyphdataservic...
Last seen:
1 day, 8 hours ago
Joined:
3 years, 9 months ago
Email:
p******@gmail.com

Software developer working in NGS analytics. 

Posts by fwuffy

<prev • 25 results • page 1 of 3 • next >
0
votes
2
answers
112
views
2
answers
Comment: C: Looking for a database containing the least-conserved, most highly polymorphic S
... Ok, thanks. I'll try that... ...
written 4 days ago by fwuffy50
5
votes
2
answers
112
views
2
answers
Looking for a database containing the least-conserved, most highly polymorphic SNPs (or regions) of the human genome (hg38)
... Hi- I need a reliable method to identify SNPs in the human genome with the most variability. Does anyone know of a data file I can download (a VCF or something easy to parse) that contains snps and a conservation score like PhyloP or PhastCons? In lieu of that if you can recommend something to iden ...
conservation human snp written 5 days ago by fwuffy50 • updated 1 day ago by colindaven930
2
votes
0
answers
610
views
0
answers
Filtering VCF with vcf-tools using --maf=0.0001 filters out all homozygous calls ?!
... I know there are some good homozygous calls in this VCF, yet vcf-tools can not see them for whatever reason. Any filtering for MAF or alt read count results in total elimination of all 1/1 homozygous calls. an example VCF row: chr3 69313555 rs1483895 C T 222 . AC1=2;AC=4;AF1=1;AN=4;ANN=T|i ...
vcf vcf-tools filtering variants written 10 months ago by fwuffy50
0
votes
0
answers
362
views
0
answers
Can anyone tell me why BCFTools Call would call a homozygous variant at this position with 1/2 alt variant reads?
... Here's a variant SNP with raw depth of 2, and DP4 field shows 1 alt-forward read. Can someone help me understand why this is called homozygous 1/1 at this position intead of hetero 0/1? (I am filtering out calls with DP < 5 after this, but it still raises questions) Thanks! samtools mpil ...
vcf bcftools variants written 10 months ago by fwuffy50
0
votes
1
answer
479
views
1
answers
Comment: C: Picard MarkDuplicates fatal error
... Updating to 2.17 fixed it Thanks! ...
written 10 months ago by fwuffy50
0
votes
1
answer
479
views
1
answers
Comment: C: Picard MarkDuplicates fatal error
... Using 1.119 Yes Googled but ran into a different error, but this looks like it. Thanks, I will try the upgrade. ...
written 10 months ago by fwuffy50
2
votes
1
answer
479
views
1
answer
Picard MarkDuplicates fatal error
... Hi, I'm running Picard MarkDuplicates on a sorted mapped BAM file. It's mapped with bwa mem, sorted with samtools sort. Pretty standard. I'm running it from a cron job. The process seems to run almost to completion and generates a bam slightly larger than the input. This is genomic data, an exo ...
picard markduplicates written 10 months ago by fwuffy50 • updated 10 months ago by h.mon22k
0
votes
1
answer
473
views
1
answers
Comment: C: VCF PL Values not making sense
... Oh right!!! Duh. It's been a while... so I'd have A/A:72, A/B:25, A/C:45, A/D: 27, B/B: 0, B/C: 47, B/D: 40, C/C: 13, C/D: 15, D/D: 59 The coverage is low cause it's a sub-sampled fastq. I also know there are issues around indels but just for t ...
written 11 months ago by fwuffy50
2
votes
1
answer
473
views
1
answer
VCF PL Values not making sense
... Hi- Can someone help me identify these PL (phred likelihood) genotype values for this multi-allelic indel row in a VCF? The VCF row: chr1 36418926 . caagaaa cAAAAaagaaa,caa,cAAaagaaa 9.95049 . INDEL;IDV=4;IMF=0.666667;DP=6;VDB=0.0340507;SGB=-0.556411;MQSB=1;MQ0F=0;AF1=0.502508;AC1=1;DP4=1,0,3 ...
vcf written 11 months ago by fwuffy50
0
votes
2
answers
2.7k
views
2
answers
Answer: A: Number of Threads for BWA MEM
... Bump. You can not run bwa with 1000 threads on one host and expect it to be 1000x faster than single thread. There will be a theoretical optimum which is a function of memory available per thread and reference genome size. Has anyone actually done tests? ...
written 12 months ago by fwuffy50

Latest awards to fwuffy

Supporter 1 day ago, voted at least 25 times.
Popular Question 15 months ago, created a question with more than 1,000 views. For SnpSift annotate reporting extra alleles not present in the input VCF
Popular Question 15 months ago, created a question with more than 1,000 views. For Convert VCF to BCF
Popular Question 3.0 years ago, created a question with more than 1,000 views. For Convert VCF to BCF
Popular Question 3.0 years ago, created a question with more than 1,000 views. For Samtools mpileup quality scores all zero

Help
Access

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