User: MMa
MMa • 240
- Reputation:
- 240
- Status:
- Trusted
- Location:
- Last seen:
- 6 months, 2 weeks ago
- Joined:
- 2 years ago
- Email:
- M**@mdanderson.org
Posts by MMa
<prev
• 35 results •
page 1 of 4 •
next >
0
votes
0
answers
221
views
0
answers
... 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 ...
written 6 months ago by
MMa • 240
0
votes
1
answer
351
views
1
answers
... 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 6 months ago by
MMa • 240
0
votes
1
answer
351
views
1
answers
... 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 6 months ago by
MMa • 240
0
votes
1
answer
351
views
1
answers
... Hi @igor, the intention is to identify all variants regardless of source. ...
written 6 months ago by
MMa • 240
1
vote
1
answer
351
views
1
answer
... 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 ...
1
vote
4
answers
703
views
4
answers
... 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 7 months ago by
MMa • 240
0
votes
2
answers
2.1k
views
2
answers
... 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 7 months ago by
MMa • 240
0
votes
1
answer
516
views
1
answers
... The docs for the current versions of pindel states depth is the same I'd get from `samtools depth`. ...
written 11 months ago by
MMa • 240
0
votes
1
answer
516
views
1
answer
... 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 ...
written 11 months ago by
MMa • 240
• updated
11 months ago by
Pierre Lindenbaum ♦ 108k
0
votes
2
answers
592
views
2
answers
... If you did that successfully on a different file yesterday, then check whether the smallfile.vcf you're using today have been indexed. ...
written 11 months ago by
MMa • 240
Latest awards to MMa
Teacher
10 months ago,
created an answer with at least 3 up-votes.
For A: Different results with edgeR when adding gene symbols
Teacher
13 months ago,
created an answer with at least 3 up-votes.
For A: Different results with edgeR when adding gene symbols
Scholar
13 months ago,
created an answer that has been accepted.
For C: I get the same modules but different number of genes WGCNA
Scholar
13 months ago,
created an answer that has been accepted.
For C: I get the same modules but different number of genes WGCNA
Use of this site constitutes acceptance of our User
Agreement
and Privacy
Policy.
Powered by Biostar
version 2.3.0
Traffic: 600 users visited in the last hour