Change chromosome names in bed file

I have a bed file like this:

chr1    11869   12227   +       processed_transcript
chr2    12613   12721   +       processed_transcript
chr3    13221   14409   +       processed_transcript
chr4    12010   12057   +       transcribed_unprocessed_pseudogene
chr5    12179   12227   +       transcribed_unprocessed_pseudogene
`

I want to change the chromosomes to NCBI format which is like this:

NC_000001.11    11869   12227   +       processed_transcript
NC_000002.12    12613   12721   +       processed_transcript
NC_000003.12    13221   14409   +       processed_transcript
NC_000004.12    12010   12057   +       transcribed_unprocessed_pseudogene
NC_000005.10    12179   12227   +       transcribed_unprocessed_pseudogene

How to do that?

rna-seq next-gen

• 2,543 views • link

written 4.5 years ago by rk.khayami94 •

1 answer

find a file mapping the two notations e.g: https://github.com/dpryan79/ChromosomeMappings/blob/master/GRCh37_NCBI2UCSC.txt

sort both files on chromosome

join both files on the chromosome using join https://linux.die.net/man/1/join

• 0 views • link

written 4.5 years ago by Pierre Lindenbaum

Thank you very much.

• 0 views • link

written 4.5 years ago by rk.khayami94 •

More posts like this

Is there a way to safely remove all underscore from gene id for hg19 from genecode
written 4.6 years ago by simplitia •

Hi, I realize that https://www.gencodegenes.org/human/release_33lift37.html gene codes hg19 gtf files has a strange annotation with an underscore append to the end of each of the …
bedtools getfasta command
written 5.4 years ago by harry •

Hi, I want to extract sequences from the coordinates with header include the name with transcript id and exon rank. Till now I used bedtools …
Error with Rseqc infer_experiment - unknown datatype
written 6.0 years ago by newbie

I'm using rseqc on my samples bam file to check the strand specific. For this I used gencode gtf file. Initially I converted gencode gtf …
Counting how many transcripts per gene in gtf file
written 6.4 years ago by sallyey2 •

Hi I want to count the number of transcripts per gene in GTF file, using python 2 on Ubuntu 18.04. Could you help me how …
featureCount vs GENCODE gene lengths
written 8.1 years ago by user31888

I used featureCount for counting reads in a bam file. featureCount notably outputs the length of each gene. However, I noticed huge differences in number …
Choice of Gene Length in TPM Calculations
written 8.5 years ago by dec986

Hello, I am attempting to convert read count output from STAR into TPM. However, there is some ambiguity in what gene length to use. from …
Comprehensive gene annotation for drosophila melanogaster
written 9.1 years ago by newscient •

Is there a trustworthy and so-called comprehensive annotation for drosophila melanogaster of all genomic regions in order to construct a TSS (transcription start site) database? …
How to get the chromosome numbers from RefSeq accession IDs?
written 9.9 years ago by birdguy.saikat •

I have an array of RefSeq accession IDs, which looks like the following: NC_000001.11 NC_000002.12 NC_000003.12 NC_000004.12 NC_000005.10 NC_000006.12 NC_000007.14 NC_000008.11 NC_000009.12 NC_000010.11 NC_000011.10 . …
Parsing GENCODE GTF to simpler BED files: am I reinventing the wheel ?
written 10.6 years ago by Charles Plessy

Dear Biostars, I am using [a shell script][1] to transform a GENCODE GTF file into smaller BED files that I use to annotate transcriptome (CAGE) …
What's the number at the begining of some lines of transcriptome GTF file means ?
written 11.8 years ago by jack

Hi, I need to parse human GTF file for my work. I downloaded it from Ensembl. Basically I don't know what is the number "1" …