Swissprot annotation

Swiss-Prot protein database (release 2024-07-24) for 25 grassland plant species.Species list:Species_code Species_nameGerPra Geranium pratenseLeuVul Leucanthemum vulgare agg.TriFla Trisetum flavescensLotCor Lotus corniculatusCreBie Crepis biennisBelPer Bellis perennisDauCar Daucus carotaPruVul Prunella vulgarisTriPra Trifolium pratensePriVer Primula verisHolLan Holcus lanatusAloPra Alopecurus pratensisGalMol Galium mollugo agg.KnaArv Knautia arvensisMedVar Medicago x variaPlaLan Plantago lanceolataLuzCam Luzula campestrisVerCha Veronica chamaedrysArrEla Arrhenatherum elatiusTriDub Trifolium dubiumLatPra Lathyrus pratensisVicCra Vicia craccaRanAcr Ranunculus acrisPlaMed Plantago mediaAjuRep Ajuga reptansAnalyse pseudo code:### blast annotationdatabase="path/to/swissprot.fasta"for i in {01..25}; do# the number corresponding to a particular speciesblastx -query path/to/BMK230426-BJ202-01P000$i.flnc.clustered.isoforms.fa -db $database -out BMK230426-BJ202-01P000${i}.blasx -evalue 1e-5 -outfmt "6 qseqid sseqid pident length mismatch gapopen qstart qend sstart send evalue bitscore qlen slen" -num_threads 16 -max_target_seqs 1done