def test_write_trunc_fasta(self):
""" writes fasta seq to output filepath correctly """
seq_order = ['seq1', 'seq2', 'seq3']
seqs = {'seq1': 'ATCG', 'seq3': 'ACCC', 'seq2': 'GGACC'}
output_dir = '/tmp/truncate_fasta_qual_test/'
create_dir(output_dir)
fasta_out_fp = output_dir + 'seqs_filtered.fna'
write_trunc_fasta(seqs, fasta_out_fp, seq_order)
expected_seqs = ['>seq1', 'ATCG', '>seq2', 'GGACC', '>seq3', 'ACCC']
actual_fasta = open(fasta_out_fp, "U")
actual_fasta = [line.strip() for line in actual_fasta]
self.assertEqual(actual_fasta, expected_seqs)
def test_write_trunc_fasta(self):
""" writes fasta seq to output filepath correctly """
seq_order = ['seq1', 'seq2', 'seq3']
seqs = {'seq1':'ATCG', 'seq3':'ACCC', 'seq2':'GGACC'}
output_dir = '/tmp/truncate_fasta_qual_test/'
create_dir(output_dir)
fasta_out_fp = output_dir + 'seqs_filtered.fna'
write_trunc_fasta(seqs, fasta_out_fp, seq_order)
expected_seqs = ['>seq1', 'ATCG', '>seq2', 'GGACC', '>seq3', 'ACCC']
actual_fasta = open(fasta_out_fp, "U")
actual_fasta = [line.strip() for line in actual_fasta]
self.assertEqual(actual_fasta, expected_seqs)
def test_write_trunc_fasta(self):
""" writes fasta seq to output filepath correctly """
seq_order = ["seq1", "seq2", "seq3"]
seqs = {"seq1": "ATCG", "seq3": "ACCC", "seq2": "GGACC"}
output_dir = "/tmp/truncate_fasta_qual_test/"
create_dir(output_dir)
fasta_out_fp = output_dir + "seqs_filtered.fna"
write_trunc_fasta(seqs, fasta_out_fp, seq_order)
expected_seqs = [">seq1", "ATCG", ">seq2", "GGACC", ">seq3", "ACCC"]
actual_fasta = open(fasta_out_fp, "U")
actual_fasta = [line.strip() for line in actual_fasta]
self.assertEqual(actual_fasta, expected_seqs)
