https://github.com/fjruizruano/ngs-protocols
Raw File
Tip revision: 39a091d1fa569a7fc717ac73c4b3de07f0a1204d authored by fjruizruano on 03 August 2023, 11:48:27 UTC
adding gfa2fas.py and extract_gfa.py
Tip revision: 39a091d
rm_join_tbl.py
#! /usr/bin/python

files = open("lista_tbl.txt")
outout = open("test.all.tbl", "w")

t = 0
m = 0

for file in files:
        file = open(file[:-1]).readlines()
        total = file[3]
	total = total.split()
	total = int(total[2])

	masked = file[5]
	masked = masked.split()
	masked = int(masked[2])

	t += total
	m += masked

	print str(t) + "\t" + str(m) + "\t" + str(100.0*m/t) + "%"

outout.write("total bases\tmasked bases\tproportion of masked bases\n")
outout.write(str(t) + "\t" + str(m) + "\t" + str(100.0*m/t) + "%")
		
outout.close()

back to top