Latest
Open
RNA-Seq
ChIP-Seq
SNP
Assembly
Tutorials
Tools
Jobs
Forum
Planet
All »
Home
Log In

User: MMa

gravatar for MMa
MMa • 280
Reputation:
280
Status:
Trusted
Location:
Last seen:
4 months, 3 weeks ago
Joined:
4 years, 2 months ago
Email:
M**@mdanderson.org

Posts by MMa

<prev • 37 results • page 1 of 4 • next >
1
vote
2
answers
1.7k
views
2
answers
Answer: A: Narrow & Broad peaks from 'macs2 bdgdiff'
... I'll be a bit more technical here. `bdgdiff` first generate three tracks from the 4 input tracks, which contains the log-likelihoods of t1 > t2, t1=t2, and t1 < t2. The three output files are in fact peaks called by `bdgpeakcall`--the subfunction for calling narrow peaks--from these three tra ...
written 18 months ago by MMa • 280
1
vote
0
answers
633
views
0
answers
Does a no-normal pipeline exist for MuSiC?
... Hi, We have a large set of tumor-only WES data that we want to run `music smg`. However, it's unclear from the documentations on whether: 1. Matched tumor-normal bam files are required for `music smg` *per se*, 1. Running `music bmr calc-bmr` is necessary before running `music smg`, and 1. I can u ...
music written 2.1 years ago by MMa • 280
0
votes
0
answers
515
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 2.7 years ago by MMa • 280
0
votes
1
answer
1.1k
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 2.7 years ago by MMa • 280
0
votes
1
answer
1.1k
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 2.7 years ago by MMa • 280
0
votes
1
answer
1.1k
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 2.7 years ago by MMa • 280
2
votes
1
answer
1.1k
views
6 follow
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 2.7 years ago by MMa • 280 • updated 2.7 years ago by d-cameron • 2.2k
1
vote
4
answers
2.1k
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 2.7 years ago by MMa • 280
1
vote
1
answer
7.2k
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 2.7 years ago by MMa • 280
0
votes
1
answer
1.3k
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 3.1 years ago by MMa • 280

Latest awards to MMa

Popular Question 5 months ago, created a question with more than 1,000 views. For Converting MuTect v1 text outputs to VCF
Popular Question 18 months ago, created a question with more than 1,000 views. For Converting MuTect v1 text outputs to VCF
Popular Question 18 months ago, created a question with more than 1,000 views. For Mapping between genomic coordinates to transcriptomic coordinates
Popular Question 2.0 years ago, created a question with more than 1,000 views. For Adding Coverage Stats into MuTect2 output VCFs
Popular Question 2.0 years ago, created a question with more than 1,000 views. For Is log-transformation needed to calculate Euclidean distance metric from FPKM/TPM?
Popular Question 2.0 years ago, created a question with more than 1,000 views. For Mapping between genomic coordinates to transcriptomic coordinates
Student 2.7 years ago, asked a question with at least 3 up-votes. For Adding Coverage Stats into MuTect2 output VCFs
Teacher 3.0 years ago, created an answer with at least 3 up-votes. For A: Different results with edgeR when adding gene symbols
Teacher 3.2 years ago, created an answer with at least 3 up-votes. For A: Different results with edgeR when adding gene symbols
Scholar 3.2 years ago, created an answer that has been accepted. For C: I get the same modules but different number of genes WGCNA
Scholar 3.3 years ago, created an answer that has been accepted. For C: I get the same modules but different number of genes WGCNA
Content
Help
Access
Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1855 users visited in the last hour