Question: why the bed produced by bowtie2 and tophat2 are different?
gravatar for A
5.2 years ago by
A3.9k wrote:

hi friends,

my accepted_hits.bam produced by tophat2 is different from bed produced by bowtie2 because is has not gene name in column one anymore and when i tried to extract the average of reads by biophysconnector package it told:

> source('/usr/data/nfs6/izadi/dropoff_matrix_modified_Subramaniam_yesleu.R')
Error in `colnames<-`(`*tmp*`, value = c("Gene ID", "Gene length", "Gene ID",  :
  length of 'dimnames' [2] not equal to array extent

what can i do to have a bed but hasing the gene name in column one please?

myposts tophat2 bowtie2 • 1.1k views
ADD COMMENTlink modified 5.2 years ago by Devon Ryan97k • written 5.2 years ago by A3.9k
gravatar for Devon Ryan
5.2 years ago by
Devon Ryan97k
Freiburg, Germany
Devon Ryan97k wrote:

Bowtie2 doesn't produce BED files, presumably you meant a BAM file. In any case, I wouldn't expect a splicing-aware aligner to produce the same output as an unspliced aligner (there is no reason to use bowtie2 when aligning RNAseq reads unless you're aligning directly against the transcriptome).

SInce we can't see the contents of the .R file, there's no possible way we could help at all. In any case, try to figure out what that error means and run the .R file line by line so you know exactly what line causes it (if you don't already know).

ADD COMMENTlink written 5.2 years ago by Devon Ryan97k

thank you Devon

ADD REPLYlink written 5.2 years ago by A3.9k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 710 users visited in the last hour