Incorporate results from GEO into jBrowse

Today maybe very specific hint, but some of you may find it useful. I was looking for a way of incorporating the data from GEO in jBrowse. First, I’ve generated BED file out of results in proprietary format, then I have converted this to bedGraph and bigWig following the hint from BioStars.
transcript2genomic.py is available through github.

# convert GEO to BED
zcat GSE53693_allinfo_final_ORFs.txt.gz | cut -f3,4,9 | transcript2genomic.py ../../ref/DANRE.gtf > GSE53693.bed
 
# BED 2 BedGraph
sort -k1,1 -k2,2n GSE53693.bed | bedtools merge -scores max | awk '{printf "%s\t%d\t%d\t%2.3f\n" , $1,$2,$3,$4}' > GSE53693.bg
# using bedtools2
#sort -k1,1 -k2,2n GSE53693.bed | bedtools merge -c 5 -o max | awk '{printf "%s\t%d\t%d\t%2.3f\n" , $1,$2,$3,$4}' > GSE53693.bg
 
# bedGraph 2 BigWig
bedGraphToBigWig GSE53693.bg ../../ref/DANRE.fa.fai GSE53693.bw

Leave a Reply

Your email address will not be published. Required fields are marked *