Guest User

Untitled

a guest
Nov 22nd, 2017
71
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.54 KB | None | 0 0
  1. >line1
  2. ACCGAGCTACTAGXXXXX
  3. >line2
  4. ACGTAX
  5.  
  6. print """Input file must be named FILE.fasta"""
  7. fasta_file = raw_input('Input file name:') # Input fasta file
  8. char = raw_input('Which sequence should be stripped?:')
  9. OutFileName = fasta_file.strip('.fasta') + '_stripped.fasta'
  10. OutFile = open(OutFileName, 'w')
  11. WriteOutFile = True
  12. data = open(fasta_file, "r")
  13. for line in data:
  14. if line.startswith('>'):
  15. OutPut = line
  16. else:
  17. OutPut = line.strip(char)
  18. print OutPut
  19. OutFile.write(OutPut)
  20. print(char)
  21. OutFile.close()
  22. quit()
Add Comment
Please, Sign In to add comment