Coverage for biobb_dna / utils / transform.py: 83%

1def inverse_complement(sequence, dna=True): 

2 """compute inverse complement sequence.""" 

3 if dna: 

4 A_complement = "T" 

5 else: 

6 A_complement = "U" 

7 complement = { 

8 "A": A_complement, 

9 A_complement: "A", 

10 "G": "C", 

11 "C": "G" 

12 } 

13 return "".join([complement[b] for b in sequence])[::-1]