Mercurial > repos > drosofff > sam_to_fastq
view sam_to_fastq.py @ 0:4c60ceadc414 draft
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/sam_to_fastq commit 8a467239b4e3aa42c448f5a088c399e53d50c3fd-dirty
author | drosofff |
---|---|
date | Mon, 21 Mar 2016 14:40:06 -0400 |
parents | |
children | c1ab6747fb66 |
line wrap: on
line source
#!/usr/bin/python # import sys import argparse def Parser(): the_parser = argparse.ArgumentParser() the_parser.add_argument( '--input', action="store", type=str, help="input SAM file") the_parser.add_argument( '--output', action="store", type=str, help="output FASTQ file") args = the_parser.parse_args() return args def print_fastq_sequence(samline, file): samfields = samline[:-1].split("\t") print >> file, '@%s\n%s\n+\n%s' % (samfields[0], samfields[9], samfields[10]) def main(input, output): infile = open (input, "r") outfile = open (output, "w") for line in infile: if line[0] == "@": continue if line.split("\t")[1] != "4": print_fastq_sequence (line, outfile) infile.close() outfile.close() if __name__ == "__main__": args = Parser() main (args.input, args.output)