¿Cómo puedo realizar un dotplot en python, que no visualice cuadrados sino líneas? Que permita comparar de forma adecuada las dos secuencias fasta. He conseguido obtener un dotplot con el siguiente código, pero es para mí imposible interpretar algo.
import numpy as np
import os
import matplotlib.pyplot as plt
def dotplot (ruta1,ruta2):
fichero = open (ruta1,'r')
secuencia = ''
for linea in fichero:
if '>' not in linea:
lineasinreto = linea.strip()
secuencia += lineasinreto
fichero2 = open(ruta2,'r')
secuencia2 = ''
for linea in fichero2:
if '>' not in linea:
lineasinreto = linea.strip()
secuencia2+=lineasinreto
matriz = np.zeros ((len(secuencia),len(secuencia2)))
for i in range(len(secuencia)):
for x in range (len(secuencia2)):
if secuencia[i] ==secuencia2[x]:
matriz[i,x]=1
return matriz
a,=dotplot(os.getcwd()+'/fichero.fa',os.getcwd()+'/fichero1.fa')
plt.matshow(a)
Me gustaría obtener algo así: http://biopython.org/DIST/docs/tutorial/images/dot_plot_scatter.png