Entering edit mode
3.1 years ago
Vaishnavi
•
0
Hi, I would like to differentiate between the Agilent SureSelect .bed file version v6 & v8. The release notes is unclear and couldn't arrive at any inference. I would like to differentiate and get the genomic regions and gene names as a separate .bed file comparing both the versions. Kindly help.
If it is BED files then use any of the
bedtools
operations to find overlaps and differences, e.g.intersect
.I tried doing it, but 6k+ rows in my bed file returned "-" in the place of gene regions. I've created those regions as a separate .bed file. I would like to annotate those regions alone in the file with gene names. Any lead would be helpful.
Without code or data help is close to impossible, please add something to work with.
I annotated the Agilent SureSelect bed file v6 with refFlat.txt which returned a bed file with gene names annotated, which looked something like this
I intersected v6.bed with v8.bed to pull out the common region. using the code
$ bedtools intersect -a v6.bed -b v8.bed -f 1 > results.bed
but i could notice that there are some regions in the bed file that had no gene names and the bed file looked like this after filtering out those regions alone
So now, i would like to annotate these regions alone with gene names. Help!!