This is a test version of Biostars. For the public version, visit https://www.biostars.org.
sed or awk command
ENST00000448914.1   13  4.28456     0       0
ENST00000415118.1   8   3.52171     0       0

how to remove the (.*) from column 1 and it looks like 

ENST00000448914 13  4.28456     0       0
ENST00000415118     8   3.52171     0       0

please tell me the sed command or awk command to remove it only .

rna-seq
• 2,136 views • link
updated by AK AK • written by harry • harry

3 answers

Hi harry,

By awk:

awk 'BEGIN{OFS="\t"} {gsub("\\.[0-9]+$", "", $1); print}'

(updated) For sed you can try:

sed -r 's/\.[0-9]+\t/\t/'
• 0 views • link
written by AK AK

Hi harry

Please use the formatting bar (especially the code option) to present your post better. I've done it for you this time.
code_formatting

You could try sed like this

sed 's/\.1//'
• 0 views • link
written by lakhujanivijay lakhujanivijay

This would only address .1s. We should account for .\d+, right?

• 0 views • link
written by Ram Ram

If the gene is always on the first column:

sed 's/\.[0-9]\{1,\}//' yourfile.txt

should work

• 0 views • link
written by vin.darb vin.darb

it will remove other (.) from other column.

• 0 views • link
written by harry • harry

It's weird because I try it and it don't remove the others (.) because I didn't put the 'g' global flag after the last slash

• 0 views • link
written by vin.darb vin.darb

It might if the first . it encounters is not the transcript version. The awk solution, or yours modified to include an anchor and a first-word ensuring regex would be safe.

• 0 views • link
written by Ram Ram

Log in to answer this question.

More posts like this
  • grep command for fasta header
    written by harry • harry
    I used this command:-- grep -Fw -A 1 -f header.txt test.fa >test_result.fa But it extracts only 1 header, not the whole which are present in …
  • Fasta header contain multiple string
    written by harry • harry
    I have large fasta file. As you see below there are > sign present in some fasta header like >exon2_ENST00000218032|>exon2_ENST00000218032 >exon17_ENST00000253024|>exon17_ENST00000253024 I want to remove …
  • Fasta and csv file reformatting
    written by MSRS MSRS
    Hi, 1. My fasta looks like >AAM01497_1 Glutamate-1-semialdehyde aminotransferase [Methanopyrus kandleri AV19] MGYEDEFPESLELFKRAERVMPGGVSSPVRRFDPYPFYVERAEGSRLYTVDGHVLIDYCLAFGPLILGHAHPEVVEAVVERVR >AAM04025_1 glutamate-1-semialdehyde 2,1-aminomutase [Methanosarcina acetivorans C2A] MVSEVTLDKSRQMYEKAKTLIPGGVSSPVRAIKPYPFYTASADGSKIRDLDGNEYIDYCLAYGPAVLGHNHPVIKAAIKEQLD I want the output like …
  • Script for making exon file
    written by harry • harry
    Hi everyone, I have an exon.fasta file which looks like --- >ENSE00002234944 TTAACTTGCCGTCAGCCTTTTCTTTGACCTCTTCTTTCTGTTCATGTGTATTTGCTGTCTCTTAGCCCAGACTTCCCGTGTCCTTTCCACCGGGCCTTTGAGAGGTCACAGGGTCTTGATGCTGTGGTCTTCATCTGCAGGTGTCTGACTTCCAGCAACTGC TGGCCTGTGCCAGGGTGCAAGCTGAGCACTGGAGTGGAGTTTTCCTGTGGAGAGGAGCCATGCCTAGAGTGGGATGGGCCATTGTTCATCTTCTGGCCCCTGTTGTCTGCATGTAACTTAATACCACAACCAGGCATAGGGGAAAGATTGGAGGAAAGATGA GTGAGAGCATCAACTTCTCTCACAACCTAGGCCA >ENSE00001746346 CAGGGCCATCAGGCACCAAAGGGATTCTGCCAGCATAGTGCTCCTGGACCAGTGATACACCCGGCACCCTGTCCTGGACACGCTGTTGGCCTGGATCTGAGCCCTGGTGGAGGTCAAAGCCACCTTTGGTTCTGCCATTGCTGCTGTGTGGAA >ENSE00002312635 CAGGGCCATCAGGCACCAAAGGGATTCTGCCAGCATAGTGCTCCTGGACCAGTGATACACCCGGCACCCTGTCCTGGACACGCTGTTGGCCTGGATCTGAGCCCTGGTGGAGGTCAAAGCCACCTTTGGTTCTGCCATTGCTGCTGTGTGGAAGTTCACTCC TGCCTTTTCCTTTCCCTAGAGCCTCCACCACCCCGAGATCACATTTCTCACTGCCTTTTGTCTGCCCAGTTTCACCAGAAGTAGGCCTCTTCCTGACAGGCAGCTGCACCACTGCCTGGCGCTGTGCCCTTCCTTTGCTCTGCCCGCTGGAGACGGTGTTTG TCATGGGCCTGGTCTGCAGGGATCCTGCTACAAAGGTGAAACCCAGGAGAGTGTGGAGTCCAGAGTGTTGCCAGGACCCAGGCACAGGCATTAGTGCCCGTTGGAGAAAACAGGGGAATCCCGAAGAAATGGTGGGTCCTGGCCATCCGTGAGATCTTCCCA GGGCAGCTCCCCTCTGTGGAATCCAATCTGTCTTCCATCCTGCGTGGCCGAGGGCCAGGCTTCTCACTGGGCCTCTGCAGGAGGCTGCCATTTGTCCTGCCCACCTTCTTAGAAGCGAGACGGAGCAGACCCATCTGCTACTGCCCTTTCTATAATAACTAA AGTTAGCTGCCCTGGACTATTCACCCCCTAGTCTCAATTTAAGAAGATCCCCATGGCCACAGGGCCCCTGCCTGGGGGCTTGTCACCTCCCCCACCTTCTTCCTGAGTCATTCCTGCAGCCTTGCTCCCTAACCTGCCCCACAGCCTTGCCTGGATTTCTAT CTCCCTGGCTTGGTGCCAGTTCCTCCAAGTCGATGGCACCTCCCTCCCTCTCAACCACTTGAGCAAACTCCAAGACATCTTCTACCCCAACACCAGCAATTGTGCCAAGGGCCATTAGGCTCTCAGCATGACTATTTTTAGAGACCCCGTGTCTGTCACTGA AACCTTTTTTGTGGGAGACTATTCCTCCCATCTGCAACAGCTGCCCCTGCTGACTGCCCTTCTCTCCTCCCTCTCATCCCAGAGAAACAGGTCAGCTGGGAGCTTCTGCCCCCACTGCCTAGGGACCAACAGGGGCAGGAGGCAGTCACTGACCCCGAGACG …
  • assigning the values in matrix in bash
    written by smrutimayipanda • smrutimayipanda
    I have a matrix consisting of gene name column and log fold change column. I want to change the log fold change values in -1,1 …
  • Place wordcount next to each word in a paragraph of text
    written by zack.saud • zack.saud
    Hi, Is there any way to place the count of each word next to the word in GREP SED or AWK (or is there any …
  • bustools transcript to gene file
    written by Barry Digby Barry Digby
    Hi, I want to ask is it appropriate to generate a `tx2gene` file from an ENSEMBL cDNA reference file for bustools input? I want to …
  • How to merge Kallisto TSV files for use with Tximport
    written by Colari19 Colari19
    I have pseudo-aligned some RNASeq reads using Kallisto where each sample was split over two flow cells. I pseudo-aligned the fastq files for each flow …
  • using tximport to prepare salmon quantification files for DEG analysis
    written by tarek.mohamed tarek.mohamed
    Hi All, I am using tximport to prepare quant.sf files generated from salmon for Deseq2 DEG analysis. My quant.sf and txgenes files looks fine. However …
  • bed file run macs error
    written by catherine catherine
    I have wig file downloaded from GEO, and since MACS doesn't accept wig file, so I convert it to bed file using wig2bed command. And …