Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- #! /usr/bin/env python3
- # trim off first ten bases from the read
- import sys
- filename = sys.argv[1]
- count = 0
- for line in open(filename):
- line = line.rstrip()
- if count == 0:
- name = line
- count += 1
- elif count == 1:
- seq = line
- count += 1
- elif count == 2:
- count += 1
- elif count == 3:
- qual = line
- count = 0
- new_seq = seq[10:]
- new_qual = qual[10:]
- print(name, new_seq, "+", new_qual, sep = "\n")
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement