Skip to content
Snippets Groups Projects

REPET V4 PASTEC release

Closed BARDET ETIENNE requested to merge dev into main
3 files
+ 201
11
Compare changes
  • Side-by-side
  • Inline
Files
3
+ 20
0
from annot_pipe import NTE50
def get_genome_lenght(infile):
with open(infile) as file:
cumLength = 0
for line in file:
if line[0] != ">":
# cumLength += len(string.rstrip(line))
cumLength += len(line.rstrip())
return cumLength
genome_size=get_genome_lenght(snakemake.params.genome)
iNTE= NTE50(FileName=snakemake.input.FLC, directory=snakemake.params.dir, genomeSize=genome_size)
iNTE.run()
iNTE= NTE50(FileName=snakemake.input.FLF, directory=snakemake.params.dir, genomeSize=genome_size)
iNTE.run()
\ No newline at end of file
Loading