def test_read_outfile(self):
"""
check if it can "read" .fna
TODO: have check what happen in the file format is invalid,
assuming its the correct fasta now
"""
model_infile_var = "ErrorModelSolexa36bp.mconf"
taxon_infile_var = "MetaSim_bint.mprf"
outfile_var = "MetaSim_bint-Empirical.fna"
metasim = RunMetaSim(model_file=model_infile_var, no_reads=100,
taxon_infile=taxon_infile_var, pdir=self.data_dir, wdir=self.working_dir,
filename=outfile_var, check_exist=True)
metasim.run()
result = metasim.read_outfile()
self.assertEqual(len(result), 2)
self.assertEqual(result.keys(), ["1", "2"])
expected = [170, 60]
for i, key in enumerate(result):
# print key, i, type(result[key]), result[key]
self.assertEqual(len(result[key]), expected[i])
