How to extract upstream and downstream of lncRNA in particular range?
1
0
Entering edit mode
10 months ago

Hello All I am trying to extract the upstream and downstream genes of differentially expressed lncRNAs. Initially, I converted the merged gtf file to bed file using bedops tool, then I used hommer Perl script to convert it annotation file using the following command

perl annotatePeaks.pl sample_peak.bed Genome.fasta -gft genome.gft >homer_annotation.txt upstream genes 300 kb downstream genes -300 kb

sample output of annotation.txt https://ibb.co/42s7Thj Then sorted the 300 kb in TSS column, further how do I get the nearby gene names, please check this file and guide me!!!

RNA-Seq lncRNA hommer bed gtf • 410 views
ADD COMMENT
2
Entering edit mode
10 months ago

bedtools flank and bedtools closest might be helpful for your analysis. Might be worth to take a look at them!

ADD COMMENT

Login before adding your answer.

Traffic: 3247 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6