from qiskit import QuantumCircuit
from qiskit.quantum_info import Kraus
import numpy as np

noise_ops = [np.sqrt(0.9) * np.array([[1, 0],
                                      [0, 1]]),
             np.sqrt(0.1) * np.array([[0, 1],
                                      [1, 0]])]
kraus = Kraus(noise_ops)

qc = QuantumCircuit(2)
qc.append(kraus, [0])
qc.append(kraus, [1])
qc.measure_all()
qc.draw()

from qiskit import Aer, transpile

backend = Aer.get_backend("aer_simulator")
tqc = transpile(qc, backend)
job = backend.run(tqc, shots=1000)
result = job.result()
counts = result.get_counts(tqc)
print(counts)