Mercurial > repos > padge > clipkit
view clipkit_repo/tests/integration/test_gappy_mode.py @ 0:49b058e85902 draft
"planemo upload for repository https://github.com/jlsteenwyk/clipkit commit cbe1e8577ecb1a46709034a40dff36052e876e7a-dirty"
author | padge |
---|---|
date | Fri, 25 Mar 2022 13:04:31 +0000 |
parents | |
children |
line wrap: on
line source
import pytest from pathlib import Path from clipkit.clipkit import execute from clipkit.files import FileFormat from clipkit.modes import TrimmingMode here = Path(__file__) @pytest.mark.integration class TestGappyMode(object): def test_simple_no_change(self): """ test gappy where no changes are expected in the resulting output alignment. usage: clipkit simple.fa """ input_file = f"{here.parent}/samples/simple.fa" output_file = "output/simpla.fa.clipkit" kwargs = dict( input_file=input_file, output_file=output_file, input_file_format='fasta', output_file_format='fasta', complement=False, gaps=0.9, mode=TrimmingMode.gappy, use_log=False, ) execute(**kwargs) with open(input_file, "r") as expected: expected_content = expected.read() with open(output_file, "r") as out_file: output_content = out_file.read() assert expected_content == output_content def test_simple_no_change_long_description(self): """ test gappy where no changes are expected in the resulting output alignment. usage: clipkit simple.fa """ input_file = f"{here.parent}/samples/simple_long_description.fa" output_file = "output/simple_long_description.fa.clipkit" kwargs = dict( input_file=input_file, output_file=output_file, input_file_format='fasta', output_file_format='fasta', complement=False, gaps=0.9, mode=TrimmingMode.gappy, use_log=False, ) execute(**kwargs) with open(input_file, "r") as expected: expected_content = expected.read() with open(output_file, "r") as out_file: output_content = out_file.read() assert expected_content == output_content def test_12_YIL115C_Anc_2_253_codon_aln(self): """ test gappy with codon alignment of yeast sequences usage: clipkit 12_YIL115C_Anc_2.253_codon_aln.fasta """ input_file = f"{here.parent}/samples/12_YIL115C_Anc_2.253_codon_aln.fasta" output_file = "output/12_YIL115C_Anc_2.253_codon_aln.fasta.clipkit" kwargs = dict( input_file=input_file, output_file=output_file, input_file_format='fasta', output_file_format='fasta', complement=False, gaps=0.9, mode=TrimmingMode.gappy, use_log=False, ) execute(**kwargs) with open( f"{here.parent}/expected/12_YIL115C_Anc_2.253_codon_aln.fasta_gappy", "r" ) as expected: expected_content = expected.read() with open(output_file, "r") as out_file: output_content = out_file.read() assert expected_content == output_content def test_12_YIL115C_Anc_2_253_aa_aln(self): """ test gappy with amino acid alignment of yeast sequences usage: clipkit 12_YIL115C_Anc_2.253_aa_aln.fasta """ input_file = f"{here.parent}/samples/12_YIL115C_Anc_2.253_aa_aln.fasta" output_file = "output/12_YIL115C_Anc_2.253_aa_aln.fasta.clipkit" kwargs = dict( input_file=input_file, output_file=output_file, input_file_format='fasta', output_file_format='fasta', complement=False, gaps=0.9, mode=TrimmingMode.gappy, use_log=False, ) execute(**kwargs) with open( f"{here.parent}/expected/12_YIL115C_Anc_2.253_aa_aln.fasta_gappy", "r" ) as expected: expected_content = expected.read() with open(output_file, "r") as out_file: output_content = out_file.read() assert expected_content == output_content def test_24_ENSG00000163519_aa_aln(self): """ test gappy with amino acid alignment of mammalian sequences usage: clipkit 24_ENSG00000163519_aa_aln.fasta """ input_file = f"{here.parent}/samples/24_ENSG00000163519_aa_aln.fasta" output_file = "output/24_ENSG00000163519_aa_aln.fasta.clipkit" kwargs = dict( input_file=input_file, output_file=output_file, input_file_format='fasta', output_file_format='fasta', complement=False, gaps=0.9, mode=TrimmingMode.gappy, use_log=False, ) execute(**kwargs) with open( f"{here.parent}/expected/24_ENSG00000163519_aa_aln.fasta_gappy", "r" ) as expected: expected_content = expected.read() with open(output_file, "r") as out_file: output_content = out_file.read() assert expected_content == output_content def test_24_ENSG00000163519_codon_aln(self): """ test gappy with codon alignment of mammalian sequences usage: clipkit 24_ENSG00000163519_codon_aln.fasta """ input_file = f"{here.parent}/samples/24_ENSG00000163519_codon_aln.fasta" output_file = "output/24_ENSG00000163519_codon_aln.fasta.clipkit" kwargs = dict( input_file=input_file, output_file=output_file, input_file_format='fasta', output_file_format='fasta', complement=False, gaps=0.9, mode=TrimmingMode.gappy, use_log=False, ) execute(**kwargs) with open( f"{here.parent}/expected/24_ENSG00000163519_codon_aln.fasta_gappy", "r" ) as expected: expected_content = expected.read() with open(output_file, "r") as out_file: output_content = out_file.read() assert expected_content == output_content def test_EOG091N44M8_aa(self): """ test gappy with amino acid alignment of Penicillium sequences usage: clipkit EOG091N44M8_aa.fa """ input_file = f"{here.parent}/samples/EOG091N44M8_aa.fa" output_file = "output/EOG091N44M8_aa.fa.clipkit" kwargs = dict( input_file=input_file, output_file=output_file, input_file_format='fasta', output_file_format='fasta', complement=False, gaps=0.9, mode=TrimmingMode.gappy, use_log=False, ) execute(**kwargs) with open(f"{here.parent}/expected/EOG091N44M8_aa.fa_gappy", "r") as expected: expected_content = expected.read() with open(output_file, "r") as out_file: output_content = out_file.read() assert expected_content == output_content def test_EOG091N44M8_nt(self): """ test gappy with nucleotide alignment of Penicillium sequences usage: clipkit EOG091N44M8_nt.fa """ input_file = f"{here.parent}/samples/EOG091N44M8_nt.fa" output_file = "output/EOG091N44M8_nt.fa.clipkit" kwargs = dict( input_file=input_file, output_file=output_file, input_file_format='fasta', output_file_format='fasta', complement=False, gaps=0.9, mode=TrimmingMode.gappy, use_log=False, ) execute(**kwargs) with open(f"{here.parent}/expected/EOG091N44M8_nt.fa_gappy", "r") as expected: expected_content = expected.read() with open(output_file, "r") as out_file: output_content = out_file.read() assert expected_content == output_content @pytest.mark.slow def test_EOG092C4VOX_aa(self): """ test gappy with amino alignment of fungal sequences usage: clipkit EOG092C4VOX_aa_aln.fasta """ input_file = f"{here.parent}/samples/EOG092C4VOX_aa_aln.fasta" output_file = "output/EOG092C4VOX_aa_aln.fasta.clipkit" kwargs = dict( input_file=input_file, output_file=output_file, input_file_format='fasta', output_file_format='fasta', complement=False, gaps=0.9, mode=TrimmingMode.gappy, use_log=False, ) execute(**kwargs) with open( f"{here.parent}/expected/EOG092C4VOX_aa_aln.fasta_gappy", "r" ) as expected: expected_content = expected.read() with open(output_file, "r") as out_file: output_content = out_file.read() assert expected_content == output_content @pytest.mark.integration class TestGappyModeCustomGapsParameter(object): def test_simple(self): """ test gappy with a custom gaps parameter usage: clipkit simple.fa -g 0.2 """ input_file = f"{here.parent}/samples/simple.fa" output_file = "output/simpla.fa.clipkit" kwargs = dict( input_file=input_file, output_file=output_file, input_file_format='fasta', output_file_format='fasta', complement=False, gaps=0.2, mode=TrimmingMode.gappy, use_log=False, ) execute(**kwargs) with open( f"{here.parent}/expected/simple.fa_gappy_gaps_set_to_0.2", "r" ) as expected: expected_content = expected.read() with open(output_file, "r") as out_file: output_content = out_file.read() assert expected_content == output_content def test_12_YIL115C_Anc_2_253_codon_aln(self): """ test gappy with codon alignment of yeast sequences usage: clipkit 12_YIL115C_Anc_2.253_codon_aln.fasta -g 0.3 """ input_file = f"{here.parent}/samples/12_YIL115C_Anc_2.253_codon_aln.fasta" output_file = "output/12_YIL115C_Anc_2.253_codon_aln.fasta.clipkit" kwargs = dict( input_file=input_file, output_file=output_file, input_file_format='fasta', output_file_format='fasta', complement=False, gaps=0.3, mode=TrimmingMode.gappy, use_log=False, ) execute(**kwargs) with open( f"{here.parent}/expected/12_YIL115C_Anc_2.253_codon_aln.fasta_gappy_custom_gaps", "r", ) as expected: expected_content = expected.read() with open(output_file, "r") as out_file: output_content = out_file.read() assert expected_content == output_content def test_24_ENSG00000163519_codon_aln(self): """ test gappy with codon alignment of mammalian sequences usage: clipkit 24_ENSG00000163519_codon_aln.fasta -g .4 """ input_file = f"{here.parent}/samples/24_ENSG00000163519_codon_aln.fasta" output_file = "output/24_ENSG00000163519_codon_aln.fasta.clipkit" kwargs = dict( input_file=input_file, output_file=output_file, input_file_format='fasta', output_file_format='fasta', complement=False, gaps=0.4, mode=TrimmingMode.gappy, use_log=False, ) execute(**kwargs) with open( f"{here.parent}/expected/24_ENSG00000163519_codon_aln.fasta_gappy_custom_gaps", "r", ) as expected: expected_content = expected.read() with open(output_file, "r") as out_file: output_content = out_file.read() assert expected_content == output_content def test_EOG091N44M8_nt(self): """ test gappy with nucleotide alignment of Penicillium sequences usage: clipkit EOG091N44M8_nt.fa -g .1 """ input_file = f"{here.parent}/samples/EOG091N44M8_nt.fa" output_file = "output/EOG091N44M8_nt.fa.clipkit" kwargs = dict( input_file=input_file, output_file=output_file, input_file_format='fasta', output_file_format='fasta', complement=False, gaps=0.1, mode=TrimmingMode.gappy, use_log=False, ) execute(**kwargs) with open( f"{here.parent}/expected/EOG091N44M8_nt.fa_gappy_custom_gaps", "r" ) as expected: expected_content = expected.read() with open(output_file, "r") as out_file: output_content = out_file.read() assert expected_content == output_content @pytest.mark.slow def test_EOG092C0CZK_aa(self): """ test gappy with amino alignment of fungal sequences usage: clipkit EOG092C0CZK_aa_aln.fasta -g .5 """ input_file = f"{here.parent}/samples/EOG092C0CZK_aa_aln.fasta" output_file = "output/EOG092C0CZK_aa_aln.fasta.clipkit" kwargs = dict( input_file=input_file, output_file=output_file, input_file_format='fasta', output_file_format='fasta', complement=False, gaps=0.5, mode=TrimmingMode.gappy, use_log=False, ) execute(**kwargs) with open( f"{here.parent}/expected/EOG092C0CZK_aa_aln.fasta_gappy_custom_gaps", "r" ) as expected: expected_content = expected.read() with open(output_file, "r") as out_file: output_content = out_file.read() assert expected_content == output_content @pytest.mark.slow def test_EOG092C4VOX_aa(self): """ test gappy with amino alignment of fungal sequences usage: clipkit EOG092C4VOX_aa_aln.fasta -g .25 """ input_file = f"{here.parent}/samples/EOG092C4VOX_aa_aln.fasta" output_file = "output/EOG092C4VOX_aa_aln.fasta.clipkit" kwargs = dict( input_file=input_file, output_file=output_file, input_file_format='fasta', output_file_format='fasta', complement=False, gaps=0.25, mode=TrimmingMode.gappy, use_log=False, ) execute(**kwargs) with open( f"{here.parent}/expected/EOG092C4VOX_aa_aln.fasta_gappy_custom_gaps", "r" ) as expected: expected_content = expected.read() with open(output_file, "r") as out_file: output_content = out_file.read() assert expected_content == output_content