User: MMa

gravatar for MMa
MMa230
Reputation:
230
Status:
Trusted
Location:
Last seen:
1 month, 2 weeks ago
Joined:
1 year, 7 months ago
Email:
M**@mdanderson.org

Posts by MMa

<prev • 35 results • page 1 of 4 • next >
0
votes
0
answers
139
views
0
answers
How likely do cDNAs with start-loss variants appear in mRNA-seq?
... I know that many transcripts having frameshift mutations will be degraded by nonsense- or readthrough-mediated degradation, subject to rules that are well-known. How how about start-loss? Most likely they won't be translated if they don't have an alternate start codon, but will they be degraded qui ...
rna-seq written 6 weeks ago by MMa230
0
votes
1
answer
179
views
1
answers
Comment: C: Strategies to call variants from a cancer sample
... Well, the somatic callers I used (MuTect, MuTect2, VarScan, MuSE) do not call somatic reversions, which I hope they could. That aside, using only default parameters and nothing else on a particular pair I'm currently working, MuTect gives less than 2,000 variants and MuSE gives less than 500. I'd b ...
written 7 weeks ago by MMa230
0
votes
1
answer
179
views
1
answers
Comment: C: Strategies to call variants from a cancer sample
... I understand this. The purpose of the variant calling in question, however, is not to study their biological significance. I plan to do some RNA-editing research, so I need to identify DNA-level variations to act blacklist positions. ...
written 7 weeks ago by MMa230
0
votes
1
answer
179
views
1
answers
Comment: C: Strategies to call variants from a cancer sample
... Hi @igor, the intention is to identify all variants regardless of source. ...
written 7 weeks ago by MMa230
1
vote
1
answer
179
views
1
answer
Strategies to call variants from a cancer sample
... Hi all, I am looking to identify DNA-level variations from a matched tumor-normal WES data. Specifically, I just want to know the variations in the tumor sample in relate to the *reference genome*, not the normal sample. I have noticed two possible approaches here: 1. Simply use a germline-varia ...
next-gen snp written 7 weeks ago by MMa230 • updated 7 weeks ago by d-cameron1.0k
1
vote
4
answers
450
views
4
answers
Comment: C: Mapping between genomic coordinates to transcriptomic coordinates
... Eventually I got to try out this one. The important `GenomicFeatures` function here is `mapToTranscripts`: library(GenomicFeatures) library(rtracklayer) gencodeTxDb <- makeTxDbFromGFF (file="gencode.annotation.gtf") gencodeTx <- transcripts (gencodeTxDb) names (gencodeTx) ...
written 8 weeks ago by MMa230
0
votes
1
answer
1.5k
views
1
answers
Answer: A: ssGSEA on RNA-Seq data from TCGA
... I know this is an old thread, but SSGSEA can be calculated using the Bioconductor package GSVA. If you use RPKM, use `ssgsea <- gsva (RPKM, method="ssgsea", kcdf="Gaussian", ...)`; if you use raw counts, use `ssgsea <- gsva (counts, method="ssgsea", kcdf="Poisson", ...)` ...
written 9 weeks ago by MMa230
0
votes
1
answer
352
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 6 months ago by MMa230
0
votes
1
answer
352
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 6 months ago by MMa230 • updated 6 months ago by Pierre Lindenbaum103k
0
votes
2
answers
375
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 6 months ago by MMa230

Latest awards to MMa

Teacher 5 months ago, created an answer with at least 3 up-votes. For A: Different results with edgeR when adding gene symbols
Teacher 8 months ago, created an answer with at least 3 up-votes. For A: Different results with edgeR when adding gene symbols
Scholar 8 months ago, created an answer that has been accepted. For C: I get the same modules but different number of genes WGCNA
Scholar 8 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: 1760 users visited in the last hour