User: ekwame
ekwame • 60
- Reputation:
- 60
- Status:
- Trusted
- Location:
- NYC
- Last seen:
- 1 week, 3 days ago
- Joined:
- 2 years, 8 months ago
- Email:
- e********@gmail.com
Posts by ekwame
<prev
• 26 results •
page 1 of 3 •
next >
0
votes
1
answer
147
views
1
answers
... Thanks! this was helpful!
...
written 6 weeks ago by
ekwame • 60
2
votes
1
answer
147
views
1
answer
... Hi all,
I'm I found overlaps between couple of genes and TCGA copy number calls but I've been unsuccessful in appendding metadata (hgnc_symbol) to subjectHits
data look like this;
1. sample of genes of interest
> my.cordinates.hg19.gr
GRanges object with 3 ranges and 3 metadata column ...
0
votes
0
answers
262
views
0
answers
... thanks @genomax for the explanation. yes it seems these genes are also in the patches.
what is usually done in the community? these are list of DNA repair genes which i want to the intersect with my `.VCFs`
any suggestions on how to go about this correctly? ...
written 5 months ago by
ekwame • 60
0
votes
0
answers
262
views
0
answers
... Hi all,
i download few ensensemble genes with cordinates (Human genes (GRCh38.p13)) from biomart and want to convert to ucsc hg38 cordinates.
there are some contigs not in hg38 ref dictionary
see this example (full file below);
CHR_HSCHR5_1_CTG1_1 70376388 70411404 ENSG00000276910 GTF2H2
i did ...
written 5 months ago by
ekwame • 60
2
votes
1
answer
394
views
1
answers
... Specifying an interval list restricts your variant calling to only regions in the interval list. by not specifying an interval list, you are calling variants across the whole genome hence the big size of your output file. interval list are particularly useful when for whole exomes sequence data whe ...
written 5 months ago by
ekwame • 60
0
votes
0
answers
156
views
0
answers
... Hi
I'm using annovar to filter germline variants by `--score_threshold` but nothing is either filtered or dropped.
please I'm I missing something in the code?
`perl annotate_variation.pl --filter /sc/hydra/projects/canan/Breast_Cancer/results/germline/new_germline07302020/Germline_108_07302020.vc ...
written 5 months ago by
ekwame • 60
1
vote
1
answer
342
views
1
answers
... Thanks @Kevin Blighe,
it turns out the 'strand' column was the source of error. i think because of `ignore.strand=TRUE` flag in `GenomicRanges::makeGRangesFromDataFrame`
set `ignore.strand=TRUE` inistially because `GenomicRanges::makeGRangesFromDataFrame` wasn't recognizing values in `strand` colu ...
written 7 months ago by
ekwame • 60
0
votes
1
answer
342
views
1
answers
... class(merged.biomart.hg19.cna)
[1] "data.frame"
typeof(merged.biomart.hg19.cna)
[1] "list"
str(merged.biomart.hg19.cna)
'data.frame': 57812 obs. of 24 variables:
$ seqnames : Factor w/ 24 levels "1","10","11",..: 1 1 1 1 1 1 1 1 1 1 ...
$ start ...
written 7 months ago by
ekwame • 60
1
vote
1
answer
342
views
1
answer
... Hello Biostars,
converting dataframe to Granges object using `GenomicRanges::makeGRangesFromDataFrame` but keeeps failing due to presumably reserved names in metadata column
but my metadata columns do not contain these
"seqnames", "ranges", "strand", "seqlevels", "seqlengths", "isCircular", "s ...
written 7 months ago by
ekwame • 60
• updated
7 months ago by
Kevin Blighe ♦ 69k
1
vote
1
answer
368
views
1
answers
... ouch! I forget to post my solution!
just in case someone else encounters similar solution.
`bcftools reheader --samples` requires a file with sample names so I specified correctly file in local directory by inserting quotes `"${sampleName}.txt"` to make it work.
thanks everyone for the support
S
...
written 7 months ago by
ekwame • 60
Latest awards to ekwame
Popular Question
4 months ago,
created a question with more than 1,000 views.
For .gstmp extension on .bam files downloaded from google bucket; and EOF marker is absent in bam file; is this a problem with downloading or my bam file is actually corrupted?
Use of this site constitutes acceptance of our User
Agreement
and Privacy
Policy.
Powered by Biostar
version 2.3.0
Traffic: 884 users visited in the last hour