More posts like this
-
Issue when trying to extract sequences from a genome using bed coordinates and samtools faidx
written by epi_vincent •Hi, I am trying to convert ChIP peak coordinates (.bed output from Pepr or sicer) to fasta sequences for motif analysis using RSAT. To that …
-
Knowing the number of exons in a range
written by zizigoluHi I have copy number segment from exome seq like > head(CN) # A tibble: 6 x 6 file Chromosome Start End Total_CN Minor_CN <chr> …
-
Plot raw coverage for copy number data
written by zizigoluI have raw copy number like > head(cndata[,-1]) Chromosome Start End Total_CN Minor_CN 1 1 10583 3457311 2 0 2 1 3458681 143542402 1 0 …
-
Calculating human genome covering each segment
written by zizigoluHi I have some genomic ranges like Chr start end 1 10583 863511 1 12841835 12854479 I want to know for instance; The range of …
-
Creating a marker file for copy number analysis
written by zizigoluHi I have a markers file contains snp position across the genome like below Marker chr position rs62635286 1 13116 rs75454623 1 14930 rs806731 1 …
-
Error in Gviz BiomartGeneRegionTrack
written by codezy •Dear friends, I would like to run the "BiomartGeneRegionTrack" in package Gviz to get the annotation file. But it seems I can not use the …
-
Creating marker file for GICTIC2
written by zizigoluHi I have generated a segmentation file (copy number) like this > head(segg) Illumina_Barcode_ful Chromosome Start End Segments_mean number_of_snps 1 t_005 1 13116 2063094 -0.722466 …
-
Getting the number of SNPs in some ranges
written by zizigoluHi, I have called copy number and I have 2 files (I have shared link of my files ); One contains some ranges > head(cndata[,c(2,3,4)]) …
-
Extracting this data frame from a .vcf file
written by zizigoluHi, I have one .vcf file of whole genome sequencing of tumour Vs normal samples of 21 patients. I need a data from like this …
-
Finding overlapping ranges in R
written by EVRHi, I have a set of intervals in a data frame and a query interval range. All I want to find the interval ranges that …
If it is from human/ a model organism with reference annotation available, then use UCSC Table Browser or Ensembl BioMart to download exon coordinates and then use bedtools to check overlap.
Alternatively, if your query segments are few, you could use the online interfaces (UCSC Table Browser) to provide query segment coordinates and retrieve relevant exons only.
Hello A!
We believe that this post does not fit the main topic of this site.
Cross-posted and zero effort. Users in multiple communities keep asking you to show effort. At least what you've tried, even if it failed. But you keep ignoring it, I personally find this utterly disrespectful, even though you conststantly keep apologizing for it. I do not buy it. To me this reads like "do my work for me". You are not an inexperienced user who is entering the field. New users might be overwhelmed with new information so it can be difficult to express yourself, difficult to find the right words to put together a good question. That excuse does not hold true for you.
https://bioinformatics.stackexchange.com/questions/12806/calculating-the-number-of-probes-for-a-given-genomic-range
For this reason we have closed your question. This allows us to keep the site focused on the topics that the community can help with.
If you disagree please tell us why in a reply below, we'll be happy to talk about it.
Cheers!