join snp list
1
0
Entering edit mode
7.0 years ago
mra8187 ▴ 20

Dear all

i am trying to join 2 list of SNPs to find out common SNP's ... whats wrong with my script ?

join‬‬ ‫‪-1‬‬ ‫‪1‬‬ ‫‪-2‬‬ ‫‪1‬‬ ‫‪<( sort‬‬ ‫‪-k1,1‬‬ SNPList1.filtered-D-1.txt) ‫‪<( sort‬‬ ‫‪-k1,1‬‬ SNPList1.filtered-D-2.txt) | sort -V > '/home/mra/Desktop/snps/D/ ‫‪CommonSNPlistD1-2.txt‬‬'

error :

sort‬‬: command not found
join‬‬: command not found
sort‬‬: command not found
rna-seq SNP • 1.1k views
ADD COMMENT
1
Entering edit mode

Which operating system are you using?

ADD REPLY
0
Entering edit mode

linux mint cinnamon 17

ADD REPLY
0
Entering edit mode

It's rather odd that even sort and join are not found on your system. There is nothing wrong with your command, but your operating system can't find the tools it needs. Which should be installed by default if I'm not mistaken.

Can you check your path using echo $PATH?

ADD REPLY
0
Entering edit mode

thank you for help

echo $PATH /home/mra/miniconda3/bin:/home/mra/miniconda3/bin:/usr/local/sbin:/usr/local/bin:/usr/sbin:/usr/bin:/sbin:/bin:/usr/games:/usr/local/games

ADD REPLY
0
Entering edit mode

Does it work if you do separately sort‬‬ ‫‪-k1,1‬‬ SNPList1.filtered-D-2.txt without the remaining pipes and parts?

ADD REPLY
0
Entering edit mode

no didnt work sort‬‬ ‫‪-k1,1‬‬ SNPList1.filtered-D-2.txt sort‬‬: command not found

ADD REPLY
0
Entering edit mode
7.0 years ago
Bioaln ▴ 360

Have you tried the comm command?

It is trivial with that:

comm -12 snplist1.txt snplist2.txt

input files must be sorted, nothing a e.g. <(sort snplist1.txt) couldn't handle.

ADD COMMENT

Login before adding your answer.

Traffic: 2344 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6