Skip to content
Snippets Groups Projects
Commit eaad78a2 authored by Amandine  PERRIN's avatar Amandine PERRIN
Browse files

Adapt test for sup gene in ffn to cove all cases

parent 0e33e554
No related branches found
No related tags found
No related merge requests found
Pipeline #38961 passed
>JGIKIPIJ_00001
aacgcgcatccagcgctagagacgctctcgacgagc
cagatctgca
cacaacaggtcgcgctcg
>JGIKIPIJ_00005
cagagatgcccc
ccgcgttt
>JGIKIPIJ_00008
aaccggcctcgcgcatcgcatcagcag
>JGIKIPIJ_00009
acgctagagagctcgcgctaagagatc
>JGIKIPIJ_00010
CCGAATAGCGCGCTCAGAGAGAGAGGA
CGATAGCTCTCGCA
CCGCATAGC
>JGIKIPIJ_00011
AAGCTCGCTAGACNNCTCGATACNNCTCGGATAGC
>JGIKIPIJ_00011
CAGATAGCNNCGCTAGAGGAGAGCTCTCGGAGA
>JGIKIPIJ_00011
CCGAGATCGCGCGCGCTCTTCTCGAGA
>JGIKIPIJ_00013
AGCTCTCGAGAGGAGCGCTCGA
CGCTCGCGCA
CCTATAGGAACCACCGGGG
>JGIKIPIJ_00014
CAGGATAGCGCGCTCAGAG
>JGIKIPIJ_00016
AAGGGATATGGCCACACAGTT
>CRISPR_toto
CCGAGGTGTGCG
CCGGTGTGCGCG
>PROKKA_02927
CCAGAGATCGCGCGC
>PROKKA_02928
AGATATGGCGCGCTCTCTC
>JGIKIPIJ_03015
AGATCCGCGCGCTATAGAGC
>JGIKIPIJ_03016
AAGGATCTCTCGCGAGAGGA
......@@ -25,7 +25,7 @@ CGCTCGCGCA
CCTATAGGAACCACCGGGG
>JGIKIPIJ_00014
CAGGATAGCGCGCTCAGAG
>JGIKIPIJ_00015
>JGIKIPIJ_03050
CAGGATAGCGCGCTCAGAG
>JGIKIPIJ_00016
AAGGGATATGGCCACACAGTT
......
......@@ -350,7 +350,7 @@ def test_create_gen_supgen(caplog):
lstfile = os.path.join(EXP_ANNOTE, "res_create_lst-prokka.lst")
res_gen_file = os.path.join(GENEPATH, "prodigal_res.gen")
assert not prokkafunc.create_gen(ffnfile, lstfile, res_gen_file)
assert ("Missing info for gene >JGIKIPIJ_00015 (from test/data/annotate/test_files/prokka_out_for_test-supGene.ffn) in test/data/annotate/exp_files/res_create_lst-prokka.lst. If it is actually present in the lst file, check that genes are ordered by increasing number in both lst and ffn files.") in caplog.text
assert ("Missing info for gene >JGIKIPIJ_03050 (from test/data/annotate/test_files/prokka_out_for_test-supGene.ffn) in test/data/annotate/exp_files/res_create_lst-prokka.lst. If it is actually present in the lst file, check that genes are ordered by increasing number in both lst and ffn files.") in caplog.text
def test_create_gen_missingSeq(caplog):
......@@ -368,21 +368,6 @@ def test_create_gen_missingSeq(caplog):
assert tutil.compare_order_content(exp_gen, res_gen_file)
def test_create_gen_missing_lastSeq(caplog):
"""
Check create gen file. The last gene in lst does not have a sequence in ffn.
Just skip it, and go to next sequence for gen file.
"""
caplog.set_level(logging.DEBUG)
logger = logging.getLogger("test_prodigal")
ffnfile = os.path.join(TEST_ANNOTE, "prokka_out_for_test-noLastSeq.ffn")
lstfile = os.path.join(EXP_ANNOTE, "res_create_lst-prokka.lst")
res_gen_file = os.path.join(GENEPATH, "prodigal_res.gen")
assert prokkafunc.create_gen(ffnfile, lstfile, res_gen_file)
exp_gen = os.path.join(EXP_ANNOTE, "res_create_gene_prokka-miss-lastGene.gen")
assert tutil.compare_order_content(exp_gen, res_gen_file)
def test_create_prt(caplog):
"""
Check that prt file is generated as expected
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment