Hi all,
I have build the following data frame
 head(href)
   chr tx_start   tx_end         g_id strand cds_start  cds_end exon_count
1 chr1  8384389  8404227 NM_001080397      +   8384389  8404073          8
2 chr1 16767166 16786584 NM_001145277      +  16767256 16785491          7
3 chr1 16767166 16786584 NM_001145278      +  16767256 16785385          8
4 chr1 16767166 16786584    NM_018090      +  16767256 16785385          8
5 chr1 48998526 50489626    NM_032785      -  48999844 50489468         14
6 chr1 33546713 33585995    NM_052998      +  33547850 33585783         12

Now when I'm trying to export it into a bed file I did:
output_href<- write.table(new_CTTS, file="new_href.bed", quote=FALSE, sep =
"\t", na="NA", row.names=FALSE, col.names=TRUE)

The thing is when I look at the file I only have:
 head(test)
  chr1 X564620 X564649 X564644 X565645  X94 X. X10
1 chr1  565369  565404  565371  566372  217  +   8
2 chr1  565463  565541  565480  566481 1214  +  15
3 chr1  565653  565697  565662  566663 1031  +  28
4 chr1  565861  565922  565883  566884  316  +  12
5 chr1  566537  566573  566564  567565  119  +  11
6 chr1  567535  567579  567562  568563 2085  +  39

but what I want is to include the corresponding gene ID as well in a last
column. How can I do that?

Best,
Nanami

	[[alternative HTML version deleted]]

