Skip to content

Latest commit

 

History

History
23 lines (14 loc) · 620 Bytes

README.md

File metadata and controls

23 lines (14 loc) · 620 Bytes

fast_stats.py

Returns N50 (or N of your choice) and other stats of a fasta/fastq file (can be gzipped) of sequences. N50 is calculated as the sequence length above which 50% of the total sequence length lies when the sequences are sorted in order of descending length.

Dependencies

Requires Numpy for calculating the median read length.

Usage

python fast_stats.py -i <infile> -n 50
Arguments

-i The fasta/q file (can be gzipped) for which you want to calculate the stats.

-n The value of n (usually 50) that you want to use.

-h Displays help.