User: graeme.thorn
graeme.thorn • 50
- Reputation:
- 50
- Status:
- Trusted
- Location:
- London, United Kingdom
- Last seen:
- an hour ago
- Joined:
- 5 years ago
- Email:
- g***********@gmail.com
Posts by graeme.thorn
<prev
• 54 results •
page 1 of 6 •
next >
1
vote
0
answers
49
views
0
answers
... I have 67 clinical tissue samples which were sequenced for RNA expression. They are in five main groups, with three of the groups linked by being triple samples from the same people taken from different locations, and the other two being taken from two groups of individuals with slightly different ...
written 6 days ago by
graeme.thorn • 50
0
votes
0
answers
100
views
0
answers
... I have split a multi-sample vcf into each sample, and I was wondering if there was a simple method for fixing the genotype for variants to 0/1 for each single-sample vcf, and changing the alternate allele to match the new genotype. For instance, if the vcf contains a read like
`
chr1 945122 . ...
written 4 months ago by
graeme.thorn • 50
0
votes
1
answer
239
views
1
answers
... Thanks! A quick test has shown that variations of this (see the answer) are what I required. The GATK help pages don't seem that helpful. ...
written 4 months ago by
graeme.thorn • 50
0
votes
1
answer
239
views
1
answers
... As per @Pierre Lindenbaum's comment above, the correct filter for the genotype is
````vc.getGenotype("SAMPLEn").isHomRef()```` to select the "0/0" genotype and
````!vc.getGenotype("SAMPLEn").isHet() && !vc.getGenotype("SAMPLEn").isHomVar()```` to select the "0/0" or "./." genotypes.
This c ...
written 4 months ago by
graeme.thorn • 50
0
votes
1
answer
239
views
1
answer
... I have a multisample VCF and I want to filter it based on:
1. The sequencing depth from the first sample (which is a germline sample)
and
2. The genotype of the first (and/or the third) samples (which are both germline samples)
I've been investigating ````gatk VariantFiltration```` for doing thi ...
written 4 months ago by
graeme.thorn • 50
0
votes
0
answers
133
views
0
answers
... Are there any resources with summarised variant data from covid-19 patients who have had their genomes sequenced? The underlying sequencing data is not required, just the called variants. There are plenty of resources with viral genomes but nothing from the patients suffering from covid-19. ...
written 5 months ago by
graeme.thorn • 50
0
votes
1
answer
155
views
1
answers
Comment:
C: RNAseq and PAM50 prediction
... Thanks Kevin, but this is work in collaboration with a commercial company who will be running PAM50 on the non-deduplicated data (the sequencing included UMIs, which we are taking into account, and they aren't), so I was looking for the most robust way of running PAM50 on the data so we can do a dir ...
written 5 months ago by
graeme.thorn • 50
0
votes
1
answer
155
views
1
answer
... I've a set of RNAseq data from breast cancer tissue samples (counts and post-cqn-normalised log(RPKM) values) and wish to use the PAM50 classifier to classify them.
I've seen the question https://www.biostars.org/p/335640/ and the question https://www.biostars.org/p/102212/, and neither are particu ...
written 5 months ago by
graeme.thorn • 50
• updated
5 months ago by
Kevin Blighe ♦ 66k
0
votes
1
answer
392
views
1
answers
... I will provide feedback to them about this. However, I'm just looking for a table of positions in hg38 that I can bolt on to the existing removed regions to ease the workflow.
EDIT: it does look like this is a frequent question to them, see for instance here: https://groups.google.com/a/soe.ucsc.ed ...
written 6 months ago by
graeme.thorn • 50
1
vote
1
answer
392
views
1
answer
... I'm converting an R code (from here: https://github.com/cancer-genomics/delfi_scripts) from hg19 to hg38 assembly, and it relies on automatically downloading telomeric and centromeric regions from the UCSC table browser:
genome <- "hg19"
mySession <- browserSession()
genome(mySes ...
written 6 months ago by
graeme.thorn • 50
• updated
5 weeks ago by
ATpoint ♦ 40k
Latest awards to graeme.thorn
Popular Question
4 months ago,
created a question with more than 1,000 views.
For Error in GSNAP - unable to find genome in directory
Scholar
12 months ago,
created an answer that has been accepted.
For A: Error in GSNAP - unable to find genome in directory
Popular Question
17 months ago,
created a question with more than 1,000 views.
For Error in GSNAP - unable to find genome in directory
Scholar
19 months ago,
created an answer that has been accepted.
For A: Error in GSNAP - unable to find genome in directory
Popular Question
24 months ago,
created a question with more than 1,000 views.
For Error in GSNAP - unable to find genome in directory
Scholar
2.0 years ago,
created an answer that has been accepted.
For A: Error in GSNAP - unable to find genome in directory
Teacher
2.1 years ago,
created an answer with at least 3 up-votes.
For A: TCGA (legacy) barcodes for identifying tumour/normal samples
Scholar
2.3 years ago,
created an answer that has been accepted.
For A: Error in GSNAP - unable to find genome in directory
Popular Question
3.0 years ago,
created a question with more than 1,000 views.
For Error in GSNAP - unable to find genome in directory
Scholar
4.9 years ago,
created an answer that has been accepted.
For A: Error in GSNAP - unable to find genome in directory
Use of this site constitutes acceptance of our User
Agreement
and Privacy
Policy.
Powered by Biostar
version 2.3.0
Traffic: 1866 users visited in the last hour