faith-hplc / huasmalltools

Some bioinformic domestic scripts.

HuaSmallTools

This repository contains some home-made bioinformic scripts. NOTICE: All python scripts in this repository were writen in python3!

Why I write this repository?

I am a student working on non-model plants. Unlike model organism, reams of annotation need to be done on non-model plants mannually. So I put my home-made scripts in this repository.

What's in this repository

For now, I assort the scripts into four functions:convert, parse, spider and some tools for sequence.

convert

1. convert_aspera.py
   Convert SRP/SRX/SRS accession to SRR url which could be used for aspera directly.
   Notes: I sincerely recommend use the fastq mode.

parse

These scripts mainly used for parse results for multiple Databases. eg.KEGG, GO, NCBI

1. parse_eggNOG.py
   Parse KEGG and GO result from eggNOG 5.0 web server results to a table format.
2. parse_KEGGkolist.py
   Parse KOALA web server results with pathway information. The result file could be used in clusterProfiler.
3. parse_IPSGO.py
   Parse GO terms in InterPro Scan results.
4. parse_BioProject.py
   Parse Biosample information from the directly downloaded BioProject txt file.
5. parse_go_obofile.py
   Parse GO database OBO file to a more user-friendly format. Some scripts in this repository may use its result.
6. parse_go_obofile.py
   Parse Uniprot database seq headers. You could use this interpret your BLAST result against Uniprot database

spider

1. spider_chinese_name.py
   Convert plant latin name to chinese name based on 植物智
   Notes: I used the pyquery module in this scripts, which is not a default module in anaconda.

sequence

These scripts was used for simple sequence processing (including MSA, GFF and other related files)

1. calculate_psy_chemi.py
   ExPASy - ProtParam tool is recommended for calculating physicochemical properties of DNA sequence. However, it could not be applied to batch processing.
2. trim_align.py
   A simple scirpt to trim alignment before construct the tree.
3. extract_cds_prokka
   Prokka results exclude the cds sequence file. This script could do this job.
   Notes: This script could only be used for prokaryote. And gff file from different software,in my experience were different more or less. So, if you use it in places other than Prokka results, Please be careful.
4. calculate_GC123.py
   The in-house scirpt for calculating GC123, which are necessary in codon usage bias analysis.

Contact me

If you have any question in using these scripts or any advices, please feel free to contact me by [email protected].