User: graeme.thorn

gravatar for graeme.thorn
graeme.thorn50
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
PCA on exon count data
... 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 ...
R splicing rna-seq written 6 days ago by graeme.thorn50
0
votes
0
answers
100
views
0
answers
Fixing genotypes from split vcf
... 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 . ...
vcf written 4 months ago by graeme.thorn50
0
votes
1
answer
239
views
1
answers
Comment: C: GATK VariantFiltration multiple filters including genotype for multi-sample VCF
... 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.thorn50
0
votes
1
answer
239
views
1
answers
Answer: A: GATK VariantFiltration multiple filters including genotype for multi-sample VCF
... 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.thorn50
0
votes
1
answer
239
views
1
answer
GATK VariantFiltration multiple filters including genotype for multi-sample VCF
... 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 ...
vcf gatk written 4 months ago by graeme.thorn50
0
votes
0
answers
133
views
0
answers
Variant data from covid-19 patients
... 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. ...
vcf covid-19 written 5 months ago by graeme.thorn50
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.thorn50
0
votes
1
answer
155
views
1
answer
RNAseq and PAM50 prediction
... 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 ...
genefu R rna-seq written 5 months ago by graeme.thorn50 • updated 5 months ago by Kevin Blighe66k
0
votes
1
answer
392
views
1
answers
Comment: C: Downloading hg38 centromere and telomere positions from UCSC table browser
... 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.thorn50
1
vote
1
answer
392
views
1
answer
Downloading hg38 centromere and telomere positions from UCSC table browser
... 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 ...
genome R ucsc written 6 months ago by graeme.thorn50 • updated 5 weeks ago by ATpoint40k

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

Help
Access

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