Entering edit mode
7.4 years ago
Anny
▴
30
Hi all,
I have a file containing blast hits for each gene, one hit per line, as following
ID E-value
g1 0.0
g1 1e-20
g1 1e-5
g2 1e-5
g2 1e-13
g3 0.01
I would like to grep the lowest e-value for each gene, how to write a script to make it?
Thank you!
Alexie
Why not just do the blast over with
-max_target_seqs 1
?Actually, they are a merged file of results generated by different methods, I want to check which method generally give the best results or more reliable.
my data: blast_scores.txt
code:
output:
datamash is command line tool from GNU