|
use lib "$RealBin/../lib" ; use_ok 'Bio::Kmer' ;
my @correctCounts =(
0,
16087,
17621,
12868,
6857,
3070,
1096,
380,
105,
17,
);
my %query =(
TTGGAGCA => 3,
TTGGAGCT => 6,
TTGGAGCTA => -1,
AAAAAAAA => 0,
);
my $kmer =Bio::Kmer->new(dirname($0). "/../data/rand.fastq.gz" ,{ kmerlength =>8});
my $hist = $kmer ->histogram();
for ( my $i =0; $i < @correctCounts ; $i ++){
is $$hist [ $i ], $correctCounts [ $i ], "Freq of $i checks out" ;
}
for my $query ( keys ( %query )){
is $query { $query }, $kmer ->query( $query );
}
my $kmerJf =Bio::Kmer->new(dirname($0). "/../data/rand.fastq.gz" ,{ kmerlength =>8, kmerCounter => "jellyfish" });
my $histJf = $kmerJf ->histogram();
for ( my $i =0; $i < @correctCounts ; $i ++){
is $$histJf [ $i ], $correctCounts [ $i ], "Freq of $i checks out" ;
}
for my $query ( keys ( %query )){
is $query { $query }, $kmerJf ->query( $query );
}
|