Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
RNA-Seq: Characterizing short read sequencing for gene discovery and RNA-Seq analysis Newsbot! Literature Watch 0 01-17-2012 06:50 AM
RNA-seq analysis-feedback on a commercial package rnaseq RNA Sequencing 7 04-05-2011 11:03 AM
ShortRead R package Aligned Read Output?? nazanin Bioinformatics 0 01-31-2011 12:34 AM
PubMed: Pash 3.0: A versatile software package for read mapping and integrative analy Newsbot! Literature Watch 0 11-26-2010 03:00 AM
Hello from Vancouver rwarren Introductions 0 10-23-2008 08:26 AM

Thread Tools
Old 12-10-2009, 08:41 AM   #1
Location: Kansas City

Join Date: Oct 2009
Posts: 88
Default BowtieToBedFormat.jar from Vancouver Short Read Analysis Package

Hello everyone,

I've been having problems using the BowtieToBedFormat tool that comes with the Vancouver Short Read Analysis Package. I have bowtie files in the default output and when I go to process them with the conversion tool it produces the output files:

log file

Great right? Well, when I unzip the files they contain nothing or txt files that are blank.

Any suggestions? All I want to do is create a GFF or 3-column BED file from the bowtie output.

Thanks in advance,
DrD2009 is offline   Reply With Quote
Old 12-12-2009, 08:15 PM   #2
Senior Member
Location: Australia

Join Date: Sep 2008
Posts: 136

you can use awk to turn a bowtie text output file into a 3 column bed file

awk -v bp=36 '{print $3"\t"$4"\t"($4+bp)}' > outfile.bed
set bp to however long your reads are, in this case 36bp
frozenlyse is offline   Reply With Quote
Old 12-14-2009, 08:39 AM   #3
Location: Kansas City

Join Date: Oct 2009
Posts: 88

:O It works!! Haha.

I had to learn what awk was, but after that I installed it on my Debian box and it worked perfectly.

Is there a way for it to work for sequence length ranges? I do have a couple bowtie files where all of the sequence lengths are 17, but there are a few where the lengths range 15-26bps, is there a setting for that using awk?

Thank you so much.
DrD2009 is offline   Reply With Quote

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 02:08 PM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO