This is a test version of Biostars. For the public version, visit https://www.biostars.org.
filter blast output

I did blastx against swissprot database in linux and the e value was set to 1e-3. The output is like

>TCONS_00006082    sp|P43298|TMK1_ARATH    57.89    57    23    2    459    292    751    806    8e-11    67.4
>TCONS_00006082    sp|P0DH62|Y4407_SELML    69.05    42    12    1    435    313    373    414    1e-08    60.1
>TCONS_00006082    sp|Q8RWZ5|SD25_ARATH    56.52    46    19    1    435    301    659    704    1e-08    60.1
>TCONS_00006082    sp|Q8VYA3|WAKLJ_ARATH    41.43    70    40    1    477    271    586    655    5e-08    58.2
>TCONS_00006082    sp|Q9C9L5|WAKLH_ARATH    40.00    70    41    1    477    271    600    669    5e-08    58.2
>TCONS_00006082    sp|P93749|Y2197_ARATH    40.54    74    43    2    501    283    225    297    6e-08    57.8
>TCONS_00006082    sp|Q9M342|WAKLP_ARATH    34.52    84    54    2    519    271    504    580    6e-08    57.8

Now I want to filter those id with e value > 0.001. How it can be done?

sequence blast

1 answer

export LANG=C; awk '$11 > 0.001' yourFile

Log in to answer this question.