#!/usr/bin/python#-*-coding:utf-8-*- " Process the Fasta file, place the ID number and sequence on one line " ImportSyswith Open (sys.argv[1]) as F:FW=open ('Out.fasta','W') Line=f.read () Line=line.replace ('\ n',"'). Replace ('>','\n>') forAainchline:fw.write (AA) fw.close ()""">chr1|hos107.1#gene1acactcccgggccccccccccccacctttcaaaaaaaaaaaaaaaaattttccccccaaagggg>chr1|hos107.2# gene2acactcccgggccccccccccccacctttcaaaaaaaaaaaaaaaaattttc>chr1|hos107.4# gene3acactcccgggccccccccccccacctttcaaaaaaaaaaaaaaaaattttc>chr1|hos107.5# GENE4ACACTCCCGGGCCCCCCCCCCCCACCTTTCAAAAAAAAAAAAAAAAATTTTC"""""">chr1|hos107.1#gene1acactcccgggccccccccccccacctttcaaaaaaaaaaaaaaaaattttccccccaaagggg>chr1|hos107.2# gene2acactcccgggccccccccccccacctttcaaaaaaaaaaaaaaaaattttc>chr1|hos107.4# gene3acactcccgggccccccccccccacctttcaaaaaaaaaaaaaaaaattttc>chr1|hos107.5# GENE4ACACTCCCGGGCCCCCCCCCCCCACCTTTCAAAAAAAAAAAAAAAAATTTTC""" # extract target sequence F=open ('./out.fasta','R') FW=open ('Target.fasta','W') forLineinchf.readlines ():ifLine.startswith ('>chr1|hos107.1'): Fw.write (line) F.close () fw.close ()"""You can extract a file of the specified target ID from a single-line file Out.fasta that is processed above and write it to the Target.fasta file"""#Overall idea:#first unify the Fasta file format from Test.fasta----Out.fasta#Take out the target ID sequence: out.fasta----Target.fasta
Python handles Fasta files, IDs and sequences on one line