Skip to content

Instantly share code, notes, and snippets.

@johnsolk
Created June 14, 2017 23:07
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save johnsolk/6946be6b0eb78613c012066f20dc3d03 to your computer and use it in GitHub Desktop.
Save johnsolk/6946be6b0eb78613c012066f20dc3d03 to your computer and use it in GitHub Desktop.
# checks for redundancies in record headers
# converts fastq to fasta format
import screed
f = open('porecamp_killifish.fasta','wb')
for n,r in enumerate(screed.open('porecamp_killifish.fastq')):
if r.name in s:
continue
else:
s.update([r.name])
del r["quality"]
screed.screedRecord.write_fastx(r , f)
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment