User: MMa

gravatar for MMa
MMa210
Reputation:
210
Status:
Trusted
Location:
Last seen:
1 month, 1 week ago
Joined:
1 year, 4 months ago
Email:
M**@mdanderson.org

Posts by MMa

<prev • 28 results • page 1 of 3 • next >
0
votes
1
answer
227
views
1
answers
Comment: C: How to get sequencing depths from VCF with Rsamtools
... The docs for the current versions of pindel states depth is the same I'd get from `samtools depth`. ...
written 9 weeks ago by MMa210
0
votes
1
answer
227
views
1
answer
How to get sequencing depths from VCF with Rsamtools
... I got from our collaborator some pindel output and the BAM files that is used as pindel input. Unfortunately, they used an older of version of pindel (before approximately 0.2.4u) that doesn't record total depth at breakpoint (i.e. the DP field in VCF). While our collaborator's standard practice is ...
R rsamtools pindel written 9 weeks ago by MMa210 • updated 9 weeks ago by Pierre Lindenbaum98k
0
votes
2
answers
220
views
2
answers
Comment: C: bcftools: how to annotate a file with a larger file using bcftools annotate?
... If you did that successfully on a different file yesterday, then check whether the smallfile.vcf you're using today have been indexed. ...
written 10 weeks ago by MMa210
1
vote
2
answers
220
views
2
answers
Comment: C: bcftools: how to annotate a file with a larger file using bcftools annotate?
... I mean the `bcftools annotate` command line. ...
written 11 weeks ago by MMa210
1
vote
2
answers
220
views
2
answers
Comment: C: bcftools: how to annotate a file with a larger file using bcftools annotate?
... What exactly was the command line you used? bcftools output gzipped vcf by default, so try renaming OCT1.genome_min.vcf to OCT1.genome_min.vcf.gz and retry. ...
written 11 weeks ago by MMa210
1
vote
2
answers
220
views
2
answers
Answer: A: bcftools: how to annotate a file with a larger file using bcftools annotate?
... It's actually more of `bcftools annotate -a largefile.vcf -c [comma separated of INFO tags you want to transfer in the form of DST_TAG:=SRC_TAG] smallfile.vcf` ...
written 11 weeks ago by MMa210
0
votes
0
answers
164
views
0
answers
Search VCF by ID in Python
... I've been trying to find a python driver for VCF that allows me to search the VCF by one or more IDs without me writing loops, but I haven't found any that doesn't use databases. Does parsers like that exist, or do I really need to use a database for that? ...
python snp written 3 months ago by MMa210
0
votes
1
answer
174
views
1
answers
Comment: C: Converting MuTect v1 text outputs to VCF
... library (vcfR) MT12VCF <- function (inputFile, outFile){ inputDF <- read.delim (inputFile, stringsAsFactors=F) if ((length(unique(inputDF$normal_name)) > 1) || (length(unique(inputDF$tumor_name)) > 1)){ stop ("This function assumes ...
written 3 months ago by MMa210
0
votes
1
answer
174
views
1
answers
Answer: A: Converting MuTect v1 text outputs to VCF
... Since the comment above implied there isn't, eventually I wrote one in R that should provide all the fields that `-vcf` emits; but since the raw MuTect text out doesn't give out the raw DP, the DP is derived as the sum of the two ADs. This script is also based on the assumption that MuTect1 only em ...
written 3 months ago by MMa210
0
votes
1
answer
174
views
1
answer
Converting MuTect v1 text outputs to VCF
... I have received some MuTect v1 text outputs from our collaborator. While I can probably write scripts that can convert the output into VCF, are there software out there that does this? R, Perl or Python are fine. ...
vcf snp mutect written 3 months ago by MMa210

Latest awards to MMa

Teacher 4 months ago, created an answer with at least 3 up-votes. For A: Different results with edgeR when adding gene symbols
Scholar 4 months ago, created an answer that has been accepted. For C: I get the same modules but different number of genes WGCNA
Scholar 4 months ago, created an answer that has been accepted. For C: I get the same modules but different number of genes WGCNA

Help
Access

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