Skip to content

Repository files navigation

🧬 AMITO

▶️ How to run

1. Clone the repository

git clone https://github.com/pratas/amito.git
cd amito

2. Install dependencies

On Ubuntu/Debian:

sudo apt update
sudo apt install -y ncbi-blast+ samtools python3 python3-pip
python3 -m pip install --user numpy pandas matplotlib

3. Run the NUMT pipeline

chmod +x *.sh
./download_and_extract_numts_T2T.sh

This step downloads or processes the references, searches the human nuclear genome for NUMTs, filters the alignments and extracts the candidate regions.

4. Generate the NUMT plots

python3 plot_numts.py human_NUMTs_coordinates.tsv
python3 plot_numt_substitution_spectrum.py human_NUMTs_substitution_spectrum.tsv

5. Generate the bacterial substitution spectrum

python3 paracoccus_substitution_spectrum.py   human_mitogenome.fasta   paracoccus_denitrificans_genome.fasta

6. Run the competitive mapping benchmark

THREADS=8 READS_PER_GROUP=50000 READ_LENGTHS=25,30,35,50,75 FORCE=1 ./run_competitive_mapping_benchmark.sh

📜 License

This project is licensed under the GNU General Public License v3.0 (GPLv3).

See the LICENSE file for the complete license text.

About

Ancient Mitogenomes Letter

Resources

Stars

Watchers

Forks

Releases

Packages

Contributors

Languages