Nucleotide statistics for fasta sequences
2 answers
awk '/^>/ {N++;next;} {B+=length($0);} END {printf("%f\n",B/N);}' input.fasta
• 0 views
•
link
written
by
Pierre Lindenbaum
seqkit stats input.fasta
Look at avg.len field and link to seqkit: http://bioinf.shenwei.me/seqkit/download/
Log in to answer this question.
More posts like this
-
How to get the average length of sequences in a fasta file
written by joaovmio •
-
How to Return Non-Aligned Proteins - BLAST
written by marcelofolhas •
-
How to extract the last 1000 nt from a group of sequences in a FASTA file?
written by Andrei •
-
Kmer enrichment analysis
written by chipolino •
-
How to divide FASTA file?
written by l.souza
-
How to submit bacterial genomes to the database?
written by olp123 •
-
Mutation analysis with amino acid sequences
written by l.souza
-
Merging FASTA files
written by l.souza
-
how to get read density for the first 1000 nt of each transcript
written by Sara
-
How to extract all non-seqenced positions from a genome (Fasta file)?
written by TheBabyBlueEyes •
I edited the title to be consistent with the question.