def run(self):
self.logger.info("Starting analysis")
fasta = Fasta(self.logger,
self.plasmid_data,
self.kmer,
self.homopolyer_compression,
max_kmer_count=self.max_kmer_count)
fastq = Fastq(self.logger,
self.input_fastq,
self.kmer,
fasta.all_kmers_in_file(),
self.min_fasta_hits,
self.print_interval,
self.output_file,
self.filtered_reads_file,
fasta,
self.homopolyer_compression,
max_gap=self.max_gap,
min_block_size=self.min_block_size,
margin=self.margin,
start_time=self.start_time,
min_kmers_for_onex_pass=self.min_kmers_for_onex_pass,
min_perc_coverage=self.min_perc_coverage,
max_kmer_count=self.max_kmer_count,
no_gene_filter=self.no_gene_filter)
fastq.read_filter_and_map()
def test_writting_to_output_file(self):
logger = logging.getLogger(__name__)
logger.setLevel(logging.ERROR)
fasta = Fasta(logger, os.path.join(data_dir,'plasmid_data.fa'),4, True)
fastq = Fastq(logger, os.path.join(data_dir,'query_gz.fastq.gz'), 4 , fasta.all_kmers_in_file(), 1, 50, 'outputfile', None, fasta, True)
fastq.read_filter_and_map()
self.assertTrue(os.path.exists('outputfile'))
self.assertTrue(filecmp.cmp(os.path.join(data_dir, 'expected_outputfile'), 'outputfile'))
os.remove('outputfile')
def test_gzipped_input(self): logger = logging.getLogger(__name__) logger.setLevel(logging.ERROR) fasta = Fasta(logger, os.path.join(data_dir,'plasmid_data.fa'),4, True) fastq = Fastq(logger, os.path.join(data_dir,'query_gz.fastq.gz'), 4 , fasta.all_kmers_in_file(), 1, 50, None, None, fasta, True) self.assertTrue(fastq.read_filter_and_map())