User: ekwame

gravatar for ekwame
ekwame60
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
Comment: C: [GRanges] How to append query metadata to subjectHits
... Thanks! this was helpful! ...
written 6 weeks ago by ekwame60
2
votes
1
answer
147
views
1
answer
[GRanges] How to append query metadata to subjectHits
... 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 ...
R bioconductor findoverlaps granges written 6 weeks ago by ekwame60 • updated 6 weeks ago by Papyrus620
0
votes
0
answers
262
views
0
answers
Comment: A: Convert list of ensembl Human genes -GRCh38.p13- from biomart to hg38 coordinate
... 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 ekwame60
0
votes
0
answers
262
views
0
answers
Convert list of ensembl Human genes -GRCh38.p13- from biomart to hg38 coordinates
... 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 ...
genome ensembl hg38 liftover assembly written 5 months ago by ekwame60
2
votes
1
answer
394
views
1
answers
Answer: A: what is the difference between calling variants with and without interval list?
... 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 ekwame60
0
votes
0
answers
156
views
0
answers
filtered and dropped vcfs empty when filtering 1000 gemomes variants with Annova
... 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 ...
annovar germline annotation 1000genomes written 5 months ago by ekwame60
1
vote
1
answer
342
views
1
answers
Answer: C: [R] makeGRangesFromDataFrame errors that metadata columns (in dframe) has reserv
... 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 ekwame60
0
votes
1
answer
342
views
1
answers
Comment: C: [R] makeGRangesFromDataFrame errors that metadata columns (in dframe) has reserv
... 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 ekwame60
1
vote
1
answer
342
views
1
answer
[R] makeGRangesFromDataFrame errors that metadata columns (in dframe) has reserved words. this is not the case
... 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 ...
R bioconductor genomicranges written 7 months ago by ekwame60 • updated 7 months ago by Kevin Blighe69k
1
vote
1
answer
368
views
1
answers
Answer: A: Rename sample names in g.vcf files fails [Error Segmentation fault (core dumped)
... 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 ekwame60

Latest awards to ekwame


Help
Access

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
_