https://github.com/fjruizruano/ngs-protocols
Tip revision: 39a091d1fa569a7fc717ac73c4b3de07f0a1204d authored by fjruizruano on 03 August 2023, 11:48:27 UTC
adding gfa2fas.py and extract_gfa.py
adding gfa2fas.py and extract_gfa.py
Tip revision: 39a091d
remove_ns.py
#!/usr/bin/python
import sys
from Bio import SeqIO
print "remove_ns.py MaskedFastaFile"
try:
fasta = sys.argv[1]
except:
fasta = raw_input("Introduce masked FASTA file: ")
data = SeqIO.parse(open(fasta), "fasta")
out = open(fasta+".withoutn", "w")
for s in data:
secu = str(s.seq)
if secu.find("N") == -1:
out.write(">%s\n%s\n" % (str(s.id), secu))
out.close()