I have two files with genes
File one (with 40000 genes)
Gene 1
Gene 2
Gene 3
Gene b
Gene f
Gene c
Gene r
Gene z
File two (with 39000 genes)
Gene 1
Gene 3
Gene 2
Gene b
I would like to know if there is a command line (with awk or bash) to extract that lines that exist in the one file and not in the two file
•
written
by
biostarsb •
I would like to know if there is a command line (with awk or bash) to extract that lines that exist in the one file and not in the two file
use comm : http://man7.org/linux/man-pages/man1/comm.1.html
comm -3 <(sort file1.txt) <(sort file2.txt )
i tested this but i have all genes not only those in the file1 me i need to extract only genes present in file1 and not in file2
only in 1st:
comm -23 <(sort file1.txt) <(sort file2.txt )
only in 2nd
comm -13 <(sort file1.txt) <(sort file2.txt )
use uniq -u instead of uniq -c |awk '$1==1'
What not simply use grep?
grep -f file2.txt -v file1.txt
You're right! Thanks
Log in to answer this question.
