User: Mehulsharma.253

Reputation:
10
Status:
New User
Location:
Last seen:
12 months ago
Joined:
1 year, 2 months ago
Email:
M**************@gmail.com

Posts by Mehulsharma.253

<prev • 22 results • page 1 of 3 • next >
0
votes
0
answers
303
views
0
answers
Mate mismatch errors in 1000genomes CRAM files
... I downloaded some CRAM files for variant calling from the [1000 genomes FTP server][1]. I also downloaded the reference genome and MD5 cache as per the instructions in this [README doc][2]. However, running Picard tools' *ValidateSAMfile* gave errors (refer below) in both the CRAM as well as the s ...
validatesam 1000genomes fixmate matepairs cram written 12 months ago by Mehulsharma.25310
0
votes
1
answer
1.1k
views
1
answers
Comment: C: Convert specific regions in bam file to vcf and annotating the vcf
... This link doesn't work ...
written 13 months ago by Mehulsharma.25310
0
votes
1
answer
707
views
1
answers
Comment: C: How to choose good tools for identifying functional SNPs?
... what tools would you recommend to asses the effect of variants on ligand binding ?Schrodinger maestro seems good but it's a paid software. ...
written 13 months ago by Mehulsharma.25310
2
votes
1
answer
689
views
1
answer
CRAM reference registry and the GRch38 reference genome
... Hi. I'm downloading some CRAM files from 1000genomes for use in variant calling. I will convert them to BAM (since most tools can't call from CRAM files). I'm a bit confused about how to go about this process since I assumed a straightforward Samtools based conversion based on a given reference FAST ...
cram bam samtools variantcalling grch38 written 13 months ago by Mehulsharma.25310 • updated 13 months ago by h.mon29k
0
votes
1
answer
402
views
1
answer
Dealing with a large sample for variant calling
... I have a large number of samples. Like a 100 samples at least (plus these are paired end reads). I aim to call variants on these samples and then predict their effects on protein structure dynamics. The only way that seems possible for now is to align each sample individually, pre-process them ind ...
jointcalling ngs alignment variantcalling written 13 months ago by Mehulsharma.25310 • updated 13 months ago by Pierre Lindenbaum125k
0
votes
2
answers
842
views
2
answers
Comment: C: Bash vs WDL for running GATK
... Thank you ! Pardon me if the question sounds dumb but what would be the difference vis a vis just piping commands together in a bash shell using the standard tools-specific commands ? What would be the primary advantage of a WDL script / Cromwell engine in that aspect ? Does it streamline my process ...
written 13 months ago by Mehulsharma.25310
0
votes
0
answers
332
views
0
answers
Effect of dbSNP build differences in variant calling
... I want to call variants on prealigned CRAM files for a 1000genomes sample, available for download. However, the dbSNP build used to recalibrate base quality scores (the dataset is already recalibrated) for all 1000genomes phase 3 pipelines is dbSNP build 142 which is the 2014 version (current build ...
assembly grch38 dbsnp variantcalling gatk written 13 months ago by Mehulsharma.25310
0
votes
1
answer
410
views
1
answers
Comment: C: GATK best practices for Broad-produced NGS data
... Thank you ! GATK forums seems to be having some issues with enabling posting and commenting because of recent spam reports but I will retry posting there. ...
written 13 months ago by Mehulsharma.25310
4
votes
1
answer
410
views
1
answer
GATK best practices for Broad-produced NGS data
... This is a more generalised question. I wish to discover variants in my raw WGS data which has been produced by studies carried out by Broad Institute. What modifications would you recommend for the standard GATK best practices workflow to better suit Broad-produced data ? For e.g., while GATK/Broa ...
gatk4 ngs broadinstitute variantcalling written 13 months ago by Mehulsharma.25310 • updated 13 months ago by vdauwera960
1
vote
1
answer
727
views
1
answer
Warning in GATK Picard MarkDuplicates
... I ran MarkDuplicates on my sorted BAM file with GATK 4.0.11 picard tools but got the following warning: > WARNING 2018-12-07 > 10:14:23 AbstractOpticalDuplicateFinderCommandLineProgram A field > field parsed out of a read name was expected to contain an integer and > did not. Read name ...
variant calling gatk picard markduplicates written 13 months ago by Mehulsharma.25310 • updated 13 months ago by Pierre Lindenbaum125k

Latest awards to Mehulsharma.253

Scholar 13 months ago, created an answer that has been accepted. For A: Setting temporary directory for GATK4 SortSam

Help
Access

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