Skip to content

Instantly share code, notes, and snippets.

@mscook

mscook/q64_q33.py

Created May 12, 2014 06:41
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 mscook/5fafd11794fcc762d9e9 to your computer and use it in GitHub Desktop.
Save mscook/5fafd11794fcc762d9e9 to your computer and use it in GitHub Desktop.
Download ZIP
Use seqret to go q64 --> q33
Raw
q64_q33.py
import glob
import os
inf = glob.glob("*.gz")
for f in inf:
tmp = f.split(".gz")[0]
os.system("cp "+f+" .original/")
os.system("gunzip "+f)
os.system("seqret fastq-illumina::"+tmp+" fastq::"+tmp+".33")
os.system("gzip "+tmp+".33")
os.system("mv "+tmp+".33.gz "+f)
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment