sort the homo sapiens genes by chromosome with Unix command lines
Hi,
I want to sort homo sapiens gene by chromosome. I tried with sort -t _ -k 3 -g genes.txt > genes_by_chr.txt, while the result is not what I want. In my genes.txt file, the chromosome column includes numbers like 1, 2, 3 ... and characters like CHR_HG1_PATCH, CHR_HG26_PATCH, and GL000009.2. What is the easiest way to sort my genes by chromosome, if I only use Unix commands?
Thank you in advance!
• 1,951 views
•
link
written
by
Huichen03 •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
How to use a variable to install an R package
written by nonaddldy •
-
Converting U to T in my miRNA sequences
written by aranyak111 •
-
Filtering pileup from site lists
written by Eliran Turgeman •
-
How to tell a sequence found in NCBI ("Nucleotide") is from Human Reference Genome or not.
written by freeheart.steven •
-
Proteolysis Of A Protein Sequence
written by karunaworah •
-
How to find out read counts from bedtools coverage?
written by ganeshram •
-
make a different color for a given amino acid
written by Learner
-
What cause the differences between genes annotations from different databases?
written by darklings
-
data file format
written by ahmedakhokhar
-
Number of SNPs from data from BioMart
written by jamesxli2007 •
How about
sort -kz,z genes.txt > genes_by_chr.txt, where z is the column of interest. In case you have positional information for your genes, you may add-kx,xn, where x is the column of the genomic coordinate and "n" indicates to sort numerically.Paste add some example data, in additon to data description. Data description is good, but it is more helpful if you can post some example data.
try
sort -V -k1if chromsomes are in first column.