Para selecionar a partir de blast2lca
selReadsExcPlants.pl
use por exemplo: selReadsExPlant.pl -i nrxUame-Purple_samples.results -k Bacteria | more
para somente bacterias
ou não use o -k e pegue tudo que não for planta
se quiser apenas o nome das reads (este deve ser o caso)
use o "cut -f 1" no lugar do "more" (comando acima)
[]s
Comentários
Postar um comentário