MetaMaps/convertMetaMapsToKraken.pl at master · DiltheyLab/MetaMaps · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
#!/usr/bin/perl

use strict;
use warnings;
use Data::Dumper;
use FindBin;
use lib "$FindBin::Bin/perlLib";
use Cwd qw/getcwd abs_path/;
use File::Copy;

$| = 1;

use SimulationsKraken;
use Util;

my $kraken_binPrefix = SimulationsKraken::getKrakenBinPrefix();
my $kraken2_binPrefix = SimulationsKraken::getKraken2BinPrefix();

my $Bracken_dir = SimulationsKraken::getBrackenDir();

my $krakenDBTemplate = SimulationsKraken::getKrakenDBTemplate();
my $kraken2DBTemplate = SimulationsKraken::getKraken2DBTemplate();

unless(scalar(@ARGV) == 1)
{
	print_help();
}

my $database = $ARGV[0];

# test that database is valid
# todo reinstate
#my %taxonID_2_contigs;
#my %contigLength;
#Util::read_taxonIDs_and_contigs($database, \%taxonID_2_contigs, \%contigLength);

my $outputDir_kraken = $database . '/kraken';
my $outputDir_kraken2 = $database . '/kraken2';

warn "Kraken 1 conversion currently deactivated.";

# SimulationsKraken::translateMetaMapToKraken (
	# $outputDir_kraken,
	# $database,
	# $krakenDBTemplate,
	# $kraken_binPrefix,
	# $Bracken_dir
# );

SimulationsKraken::translateMetaMapToKraken2 (
	$outputDir_kraken2,
	$database,
	$kraken2DBTemplate,
	$kraken2_binPrefix,
);


sub print_help
{
	print qq(
convertMetaMapsToKraken.pl

  Convert a MetaMaps DB to Kraken/Bracken.

Usage:

  perl convertMetaMapsToKraken.pl dbNAME

Example:

  perl convertMetaMapsToKraken.pl databases/miniSeq


);
exit;
}