python
from Bio import SeqIO
from Bio import pairwise2
sequences = SeqIO.parse("sequences.fasta", "fasta")
for seq in sequences:
print("ID:", seq.id)
print("Sequence:", seq.seq)
print("Length:", len(seq))
seq1 = "ATCG"
seq2 = "AGTC"
alignments = pairwise2.align.globalxx(seq1, seq2)
for alignment in alignments:
print("Score:", alignment.score)
print("Aligned sequence 1:", alignment.seqA)
print("Aligned sequence 2:", alignment.seqB)
