Hi,
Does anyone here know how to extract lines from a gtf file using a list/subset of gene id obtained from the same gtf file? I basically want a 'recoded' (in vcf terminology) gtf file containing information for only those genes which I am interested in.
I tried awk awk 'FNR==NR {a[$0];next} {for (i in a) if (i~$1) print i}' and grep grep -Fwf but these have not yielded what I want. Thank you for your help.
•
written
by
berge2015
2 answers
try the powerfull CSV/TSV toolkit csvtk, usage of csvtk grep
csvtk grep -H -t -f 1 -r -P gene_id_list.txt ref_CDS.gtf
you may change the column index -f 1 where the gene id locates
After playing with awk for a while, I came up with a solution:
awk -F'"' 'FNR==NR {block[$0];next} $2 in block' gene_id_list.txt ref_CDS.gtf > out.txt
[Note the quote delimeter]
While not the most elegant solution, this does what's asked in the question. Hope it helps anyone else looking for something similar.
Log in to answer this question.
-
Add columns to a file based on matched column and maintain the order of lines using awk
written by AP •
-
sequence extract problem
written by mxlsherry1992 •
-
Output bigger than input. VCFtools
written by victor.agrs •
-
Create multiple files using awk
written by A_Lh •
-
Ignore SAM header section using awk
written by max_19
-
Exacting lines containing Id from one file from ID from other file.
written by 1234anjalianjali1234 •
-
how to extract/grep lines containing specific word?
written by arsilan324
-
search a large dbsnp file using another most efficiently
written by bioguy24
-
Removing reads from fastq file based on position
written by firestar
-
Extracting regions from .bam file using a .gff or .gtf file
written by Tom A •
can you please paste some sample data?