Advertisement
Guest User

Untitled

a guest
Mar 30th, 2017
69
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.50 KB | None | 0 0
  1. #! /usr/bin/env python3
  2.  
  3. # trim off first ten bases from the read
  4. import sys
  5.  
  6. filename = sys.argv[1]
  7.  
  8. count = 0
  9. for line in open(filename):
  10. line = line.rstrip()
  11.  
  12. if count == 0:
  13. name = line
  14. count += 1
  15.  
  16. elif count == 1:
  17. seq = line
  18. count += 1
  19.  
  20. elif count == 2:
  21. count += 1
  22.  
  23. elif count == 3:
  24. qual = line
  25. count = 0
  26.  
  27. new_seq = seq[10:]
  28. new_qual = qual[10:]
  29.  
  30. print(name, new_seq, "+", new_qual, sep = "\n")
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement