In [1]:

    

%matplotlib inline 

import numpy as np
import matplotlib.pyplot as plt


    

In [2]:

    

graph_num = 0


    

In [3]:

    

in_file = open("/mnt/bulk_storage/jackgao/del_reads.fa.counted");

n = 0

while True:
    data = in_file.readline()
    
    if data == '' or data == '\n':
        break
        
    data = data[:-1]
    orig_counts = map(lambda x: int(x), in_file.readline()[:-2].split(" "))
    
    muta_counts = map(lambda x: int(x), in_file.readline()[:-2].split(" "))
     
    refn_counts = map(lambda x: int(x), in_file.readline()[:-2].split(" "))
    
    if n == graph_num:
        plt.figure(figsize=(20,10))
        plt.plot(orig_counts)
        plt.plot(muta_counts)
        plt.plot(refn_counts)
        plt.xticks(range(0, len(orig_counts)), data[:len(orig_counts)])
        
        break
    else:
        n += 1
        
        
#plt.yscale('log', nonposy='clip')
plt.show()


    

In [ ]:

    




    

In [ ]: