User: Robert Sicko

gravatar for Robert Sicko
Robert Sicko570
Reputation:
570
Status:
Trusted
Location:
United States
Last seen:
4 hours ago
Joined:
6 years, 4 months ago
Email:
r********@gmail.com

about me

Posts by Robert Sicko

<prev • 83 results • page 1 of 9 • next >
0
votes
6
answers
313
views
6
answers
Comment: C: How to create a tab delimited file?
... I haven't seen this command before. Awesome! ...
written 11 weeks ago by Robert Sicko570
0
votes
2
answers
285
views
2
answers
Comment: C: Calling CFTR Poly TG/Poly T Genotype
... Charles, I looked into GET-evidence further. I see their about page does at least states it is not clinical quality: > "GET-Evidence is not a “clinical quality” tool for diagnosis or > treatment of any disease or medical condition." They do allow anyone with a login to [edit][1] their annot ...
written 12 weeks ago by Robert Sicko570
0
votes
2
answers
285
views
2
answers
Comment: C: Calling CFTR Poly TG/Poly T Genotype
... Charles, thank you for the thoughtful response. I completely agree with your statement on the need to throughly validate variant calling, especially in a targeted panel, especially when reporting clinically. We are validating the assay for clinical use (screening). The script I posted above is only ...
written 12 weeks ago by Robert Sicko570
1
vote
2
answers
285
views
2
answers
Answer: A: Calling CFTR Poly TG/Poly T Genotype
... Just to wrap this up. I found this [paper][1]. I contacted Dr. Pagin and he explained their script used grep on fastqs and counted different alleles to make a call per sample. I used that idea in this [script][2]. If anyone uses this script in the future, the HET_CUTOFF (line 554) will likely vary d ...
written 12 weeks ago by Robert Sicko570
3
votes
2
answers
285
views
2
answers
Calling CFTR Poly TG/Poly T Genotype
... Background: The CFTR gene has a Poly TG/Poly T region in intron 8 that has [varying clinical][1] consequences. I'm using targeted NGS for this region (150bp, paried end) and I'm trying to call the PolyTG/PolyT genotypes from fastq files. My first attempt was bbmap run with default parameters and a ...
ngs alignment written 3 months ago by Robert Sicko570 • updated 12 weeks ago by Charles Warden6.6k
0
votes
2
answers
715
views
2
answers
Comment: C: Estimating cross contamination in a set of BAMS
... **verifyBamID** does need a vcf, but it is a population reference VCF (1000genomes) I've used it for detecting contamination in a targeted panel with alright results. see my [question][1] on their user group page. reference_vcf=/media/sf_BigShare/SCID/180213_TSCA_r1_sop_test/work/reference/180 ...
written 11 months ago by Robert Sicko570
1
vote
1
answer
1.0k
views
1
answers
Comment: C: BWA mem running extremely slow
... Any chance you can allocate more memory? I'm guessing your 14Gbp reference will require just over 24gigs. I [had 5gig allocated for human genome][1] and alignment was really slow. After I followed Heng's advice and ran with 6gig it was normal speed. 3.1 human reference/5.5 gig minimum allocation = ...
written 17 months ago by Robert Sicko570
1
vote
1
answer
706
views
1
answers
Answer: A: Finding Gaps in Targeted NGS Coverage Incorporating Replicates
... Thanks Kevin for pointing me in right direction. I was able to accomplish what I need doing exactly what you say. For anyone who may run across a similar problem. I used ````bedtools unionbedg```` to output a combined per-base file with depth of each replicate in columns 4 and 5. Then I could ````a ...
written 18 months ago by Robert Sicko570
0
votes
1
answer
706
views
1
answers
Comment: C: Finding Gaps in Targeted NGS Coverage Incorporating Replicates
... thanks, I figured I was over complicating things. Comparing the per-base output to find the common regions under 20x does seem like a better approach. I'll have to guess and check on finding the regions in common in the per-base output that are under 20x. My initial thought is combine the replicate ...
written 18 months ago by Robert Sicko570
2
votes
1
answer
706
views
1
answer
Finding Gaps in Targeted NGS Coverage Incorporating Replicates
... I'm trying to identify regions of a targeted NGS panel that lack 20x coverage. I've figured out how to get my desired output (described below) per sequencing run. I'm running into trouble figuring out how to use sample replicates and only output regions that are missing in both replicates. Example: ...
bedtools python ngs pandas awk written 18 months ago by Robert Sicko570

Latest awards to Robert Sicko

Scholar 12 weeks ago, created an answer that has been accepted. For A: Question in Database of Genomic variants
Popular Question 4 months ago, created a question with more than 1,000 views. For Functional annotation of network of gene-sets
Popular Question 5 months ago, created a question with more than 1,000 views. For Functional annotation of network of gene-sets
Popular Question 8 months ago, created a question with more than 1,000 views. For Functional annotation of network of gene-sets
Popular Question 11 months ago, created a question with more than 1,000 views. For BWA mem alignment taking a while in a virtual machine
Scholar 18 months ago, created an answer that has been accepted. For A: Question in Database of Genomic variants
Popular Question 21 months ago, created a question with more than 1,000 views. For BWA mem alignment taking a while in a virtual machine
Popular Question 21 months ago, created a question with more than 1,000 views. For BWA mem alignment taking a while in a virtual machine
Appreciated 2.4 years ago, created a post with more than 5 votes. For A: Do You Want Your Exome/Genome Sequenced?
Teacher 2.4 years ago, created an answer with at least 3 up-votes. For A: Question in Database of Genomic variants
Popular Question 2.7 years ago, created a question with more than 1,000 views. For Calculating Depth For Each Amplicon From A Mapped Bam
Great Question 2.9 years ago, created a question with more than 5,000 views. For Renaming All Files In A Directory With A Perl Script
Teacher 3.1 years ago, created an answer with at least 3 up-votes. For A: Question in Database of Genomic variants
Popular Question 3.1 years ago, created a question with more than 1,000 views. For Performing Pathway Analysis On Cnv Data
Popular Question 3.9 years ago, created a question with more than 1,000 views. For Calculating Depth For Each Amplicon From A Mapped Bam
Popular Question 4.3 years ago, created a question with more than 1,000 views. For Calculating Depth For Each Amplicon From A Mapped Bam
Teacher 4.4 years ago, created an answer with at least 3 up-votes. For A: Do You Want Your Exome/Genome Sequenced?
Scholar 4.4 years ago, created an answer that has been accepted. For A: Question in Database of Genomic variants
Popular Question 4.8 years ago, created a question with more than 1,000 views. For Renaming All Files In A Directory With A Perl Script
Voter 5.1 years ago, voted more than 100 times.
Supporter 5.1 years ago, voted at least 25 times.
Teacher 5.7 years ago, created an answer with at least 3 up-votes. For A: Do You Want Your Exome/Genome Sequenced?
Teacher 5.7 years ago, created an answer with at least 3 up-votes. For A: Do You Want Your Exome/Genome Sequenced?
Appreciated 5.7 years ago, created a post with more than 5 votes. For A: Do You Want Your Exome/Genome Sequenced?

Help
Access

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