SEQanswers

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

Reply
 
Thread Tools
Old 12-10-2009, 08:41 AM   #1
DrD2009
Member
 
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:

Histogram.txt
log file
unpaired.tsv.gz
_simple_paired_samechr.tsv.gz
_simple_paired_spanchr.tsv.gz
_forked_paired.tsv.gz
_complex_paired.tsv.gz

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,
Brandon
DrD2009 is offline   Reply With Quote
Old 12-12-2009, 08:15 PM   #2
frozenlyse
Senior Member
 
Location: Australia

Join Date: Sep 2008
Posts: 136
Default

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

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

Join Date: Oct 2009
Posts: 88
Default

: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
Reply

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