Calculates both scaffold and contig statistics (N50, L50, etc.) from a scaffold FASTA file.
Project description
assembly_stats
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.
Usage
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.
Next steps
- Add ability to save NumPy sequence length arrays for further visualization, since generating these are what takes the most time.
Project details
Download files
Download the file for your platform. If you're not sure which to choose, learn more about installing packages.
Source Distribution
Built Distribution
Hashes for assembly_stats-0.1-py3-none-any.whl
Algorithm | Hash digest | |
---|---|---|
SHA256 | ffcd46e91edeaf3d177b229207a1d3c882019fd2ecc1fdefea83a02835a2bce9 |
|
MD5 | bc66635e690f6d9773a64084d33ab9ae |
|
BLAKE2b-256 | f1dde6f3423112af8a9486388f65e0410c396db7116481043b256321cf15c216 |