This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to find shortest lenth or longest length from fasta file

Hi, I am trying to sort shortest length and longest length fasta seq of amino acids from fasta file. If anyone know anything about it please let me know.

GNNVVIGAGAKVLGSFKVGDNVKIGAGSVVNKVVPSDTTVVGIPGRIVLHNGVPIKDPDL
RHDELPDPVNEMLKCLMQRVEYLEERLNDEESRYNVFEIIQHHDTSKRGVPPKGERESCD
VRVRSNDL
>tr|A0A3D3WF91|A0A3D3WF91_9PROT
MRGFIDELRSIRERDPAAGGLLGLMFLTPSVHVMLAYRLGHRLWRWRMRFLARFTMQLAR
WFTGIEIHPAARIGKRFFIDHGMGVVIGETAVIGDDVTFYHGVTLGGLLPSVDSDGQRQK
KRHPTIRNNVIVGAGAQILGAITVNECARVGANSVVVKDVPQGTTVTGIPARPVSARRVT
DETTFQPYGTPSDLSSDARDKAIKGLLREVDILHNKLKMLEDLQGEAKQSESSFQPEATS
LSASDRN
>tr|Q2JWT3|Q2JWT3_SYNJA
MLTESSRVSGTKSLPSAVESPETGGSDSPSSAKPEPKLGFWQQFWEDIDCVFERDPAARN
RWEVLLTYPGVHALFLHRIAHWLWKRRCFFLARLLSFISRSFTLIEIHPAARIGRRFFID
HGCGVVIGETAEIGDDVTLYHGVTLGGTSWTKGKRHPTLEDGVIVGTGAKILGPVRIGAR
ARIGANAVVIQDVAPGMTVVGIPGRAVIPPHQRRIPAHGIDLDHHLMPDPVGRAIEQLLH
RIQELEAQVARLNQEREQDRQDQVRCE
>sp|P0A9D5|CYSE_ECOL6
MSCEELEIVWNNIKAEARTLADCEPMLASFYHATLLKHENLGSALSYMLANKLSSPIMPA
IAIREVVEEAYAADPEMIASAACDIQAVRTRDPAVDKYSTPLLYLKGFHALQAYRIGHWL
WNQGRRALAIFLQNQVSVTFQVDIHPAAKIGRGIMLDHATGIVVGETAVIENDVSILQSV
TLGGTGKSGGDRHPKIREGVMIGAGAKILGNIEVGRGAKIGAGSVVLQPVPPHTTAAGVP
ARIVGKPDSDKPSMDMDQHFNGINHTFEYGDGI

Thank you!

fasta
• 3,073 views • link
written by Neel • Neel

use seqkit sort sequences by length

$ seqkit -w 0 --quiet sort -l test.fa
• 0 views • link
written by cpad0112 cpad0112

this looks suspiciously like homework ...

• 0 views • link
written by JC JC

Hi, i used this script, which result is given below but i want to find the seq which is shortest and longest in all seq of fasta file.

"awk '/^>/{if (l!="") print l; print; l=0; next}{l+=length($0)}END{print l}' file.fasta"

>tr|A0A0B4WZN6|A0A0B4WZN6_9HYPH Serine acetyltransferase OS=Rhizobium gallicum bv. gallicum R602sp OX=1041138 GN=cysE PE=3 SV=1
276
>tr|A0A2A5GE89|A0A2A5GE89_9GAMM Serine acetyltransferase OS=Porticoccaceae bacterium OX=2026782 GN=cysE PE=3 SV=1
280
>tr|A0A5F0DFU6|A0A5F0DFU6_9MICO Serine acetyltransferase OS=Cryobacterium sp. HLT2-28 OX=1259146 GN=cysE PE=3 SV=1
192

Thank you|

• 0 views • link
updated by GenoMax GenoMax • written by Neel • Neel

Thank you so much to all for your reply

• 0 views • link
written by Neel • Neel

1 answer

seqmagick has a feature for summarizing sequence files that gives minimum and maximum length:

$ seqmagick info file.fasta 
name       alignment    min_len   max_len   avg_len  num_seqs
file.fasta FALSE            128       273    228.75         4
• 0 views • link
written by Jesse • Jesse

Log in to answer this question.

More posts like this

No similar posts found.