Calculates both scaffold and contig statistics (N50, L50, etc.) from a scaffold FASTA file.
A Python script that takes a FASTA file as input and calculates both scaffold and contig statistics (N50, L50, etc.) from a scaffold FASTA file. It does this by breaking each scaffold wherever there is more than one N and then calculating statistics for both the scaffolds and contigs.
This is a re-write of fasta_metadata_parser to speed up the old implementation, and -- most importantly -- to learn how to install Python scripts onto the Smithsonian HPC.
The only Python library that this script relies on is NumPy, but this is installed and set up in a conda environment on the HPC.
# Load the HPC assembly_stats module, which activates a conda environment # with the appropriate versions of Python and NumPy installed. $ module load bioinformatics/assembly_stats # Now that the module is loaded, assembly_stats.py should be accessible directly # in your PATH. $ assembly_stats.py [genome_file.fasta]
After calculating the statistics for the genome assembly, they will be printed out in JSON format.
- Add ability to save NumPy sequence length arrays for further visualization, since generating these are what takes the most time.
Download the file for your platform. If you're not sure which to choose, learn more about installing packages.
|Filename, size||File type||Python version||Upload date||Hashes|
|Filename, size assembly_stats-0.1-py3-none-any.whl (5.1 kB)||File type Wheel||Python version py3||Upload date||Hashes View|
|Filename, size assembly_stats-0.1.tar.gz (3.2 kB)||File type Source||Python version None||Upload date||Hashes View|
Hashes for assembly_stats-0.1-py3-none-any.whl