Hello everybody, 

Here i my issue, I have a Trinity.fasta file like this 

    >TRINITY_DN5631_c0_g2_i1 len=947 path=[0:0-946]
    TACAACTTGAACATCAACAATGGTTGCGCAGCTATTGCCATCCGCGACGTTCGAGGACTGCGTGCGAA
    >TRINITY_DN62279_c1_g1_i1 len=298 path=[0:0-297]
    TATTACCATTATTATTATTATCATATTTATGTTCATTATTATCATTATCATAATCATTATCATCTTGATA
     ...
And I also have a list of id in a id.txt file :
     
    TRINITY_DN16359_c0_g1_i4
    TRINITY_DN62279_c1_g1_i1
    ...

I am trying to extract for my fasta file the sequences that have their id in my txt.file.
I am using seqkit for that, but with no success :

     seqkit grep -n -f id.txt Trinity.fasta -o result.fa

Does anyone know how to fix it ?
     


Referring to the [usage](https://bioinf.shenwei.me/seqkit/usage/#grep), you should not switch on `-n` , just use `seqkit grep -f id.txt seqs.fa`.

    -n, --by-name               match by full name instead of just id
    -i, --ignore-case           ignore case
    -r, --use-regexp            patterns are regular expression

For `-nrif`,  it's **partly matching** full FASTA header  by regular expression and case-ignored, with patterns from file. 
*This may produce some unwanted results*. For example, `seq_1` matches `seq_10` with `-nri`.


