User: Helene

gravatar for Helene
Helene390
Reputation:
390
Status:
Trusted
Location:
United States
Last seen:
1 month ago
Joined:
3 years ago
Email:
s***********@gmail.com

Posts by Helene

<prev • 93 results • page 1 of 10 • next >
0
votes
2
answers
209
views
2
answers
Comment: C: Mutect 2 on WGS data takes too long to run
... I had a little experience with VarScan2 - don't you need to use samtools pileup? I recall the pileup procedure to be painfully long too... Maybe I did not do it correctly, but please comment on it. Thank you. ...
written 5 weeks ago by Helene390
3
votes
2
answers
209
views
2
answers
Mutect 2 on WGS data takes too long to run
... I used Mutect2 on human exome sequencing (30x) data before and without parallel processing, it took me 5.5 days. That was still doable - However, now I am trying to run Mutect2 on my human WGS data (30x, bqsr bam file size ~120G), the estimation tells me 500+ weeks at the start! I feel like even I ...
wgs mutect written 5 weeks ago by Helene390 • updated 5 weeks ago by donfreed1.1k
0
votes
1
answer
227
views
1
answers
Comment: C: GT field in VCF only has one number?
... By the way, would you please comment on why the software did not take the line3 and line2 as one single indel case? ...
written 7 weeks ago by Helene390
0
votes
2
answers
262
views
2
answers
Comment: C: somatic calls by somatic SNV caller differ a lot, when comparing with cancer min
... Thank you for the reply. By low allele freq, do you mean that the one of the allele is only supported by limited number of reads, but somehow the germline caller captured it? Of curiosity, have you done such comparison using your data? ...
written 7 weeks ago by Helene390
0
votes
1
answer
227
views
1
answers
Comment: C: GT field in VCF only has one number?
... Thank you so much for the reply. Great points! I added some adjacent lines from the original file, and interestingly all the "1" GT cases are next to a indel case. I wonder why IVC did not put line3 with line2 as an indel case... ...
written 7 weeks ago by Helene390
2
votes
2
answers
262
views
2
answers
somatic calls by somatic SNV caller differ a lot, when comparing with cancer minors normal germline calls
... I have a few cancer samples that were analyzed using GATK germline pipeline (call SNVs of each sample, not the cohort study setting). Recently we got the corresponding normal samples sequenced, and I did GATK on them as well. I obtained one sets of somatic calls by subtracting the germline cancer ...
somatic snv snp written 7 weeks ago by Helene390 • updated 7 weeks ago by d-cameron850
2
votes
1
answer
227
views
1
answer
GT field in VCF only has one number?
... I recently got some VCF files from Illumina (IVC output), and I noticed that there are certain lines having GT field only with 1 number. For example, see the 3rd line below (I'm including the adjacent lines.): chr1 21723729 . T C 123.0 PASS SNVSB=-17.2;SNVHPOL=4 ...
vcf written 7 weeks ago by Helene390
0
votes
1
answer
170
views
1
answers
Comment: C: Effectively search for positions in large bedgraph files using python (or other
... Hi @AlexReynolds Thank you so much for the reply and introducing bedops! I will give a short:) ...
written 9 weeks ago by Helene390
2
votes
1
answer
170
views
1
answer
Effectively search for positions in large bedgraph files using python (or other language if needed)
... I guess it is more of an algorithm question. I have a group of large bedgraph files (shown below, contains every single position of the genome), each is about 20-30Gb in size. I need to look into these files and locate the lines having locations match another file (small, < 100kb) (e.g. file test ...
bedgraph python written 10 weeks ago by Helene390 • updated 9 weeks ago by Alex Reynolds20k
0
votes
2
answers
188
views
2
answers
Comment: C: how can I keep only those SNP in file "a" which match to SNPs in file "b"?
... You are right about the line. I edited my answer. As for dictionary, I'm used to it to avoid any duplicated record in filtered.id file. That said, I guess even there are duplicated records, it does not matter in this case. Thank you for your reply. ...
written 10 weeks ago by Helene390

Latest awards to Helene

Appreciated 5 weeks ago, created a post with more than 5 votes. For Develop an R package for CRAN vs BioConductor
Popular Question 7 weeks ago, created a question with more than 1,000 views. For Non-cancer somatic mutation calling
Popular Question 8 weeks ago, created a question with more than 1,000 views. For Remove both pair end reads with low mapping quality using Samtools
Popular Question 3 months ago, created a question with more than 1,000 views. For Remove both pair end reads with low mapping quality using Samtools
Student 6 months ago, asked a question with at least 3 up-votes. For Non-cancer somatic mutation calling
Popular Question 6 months ago, created a question with more than 1,000 views. For Remove both pair end reads with low mapping quality using Samtools
Appreciated 8 months ago, created a post with more than 5 votes. For Develop an R package for CRAN vs BioConductor
Good Question 8 months ago, asked a question that was upvoted at least 5 times. For Develop an R package for CRAN vs BioConductor
Student 8 months ago, asked a question with at least 3 up-votes. For Non-cancer somatic mutation calling
Popular Question 10 months ago, created a question with more than 1,000 views. For Identify point mutations from each read in sam files
Student 10 months ago, asked a question with at least 3 up-votes. For Non-cancer somatic mutation calling
Student 16 months ago, asked a question with at least 3 up-votes. For Non-cancer somatic mutation calling
Scholar 16 months ago, created an answer that has been accepted. For A: Can't locate dbSNP 131 anymore?
Teacher 16 months ago, created an answer with at least 3 up-votes. For A: Illumina's HiSeq strange Phred quality score
Scholar 18 months ago, created an answer that has been accepted. For A: Can't locate dbSNP 131 anymore?
Supporter 18 months ago, voted at least 25 times.
Scholar 19 months ago, created an answer that has been accepted. For A: Can't locate dbSNP 131 anymore?

Help
Access

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