Skip to content

Commit

Permalink
info added to fasta2homozygous
Browse files Browse the repository at this point in the history
  • Loading branch information
Leszek Pryszcz committed Nov 12, 2014
1 parent 98cff09 commit d928197
Showing 1 changed file with 3 additions and 1 deletion.
4 changes: 3 additions & 1 deletion fasta2homozygous.py
Expand Up @@ -236,7 +236,7 @@ def main():
parser = argparse.ArgumentParser(description=desc, epilog=epilog, \
formatter_class=argparse.RawTextHelpFormatter)

parser.add_argument('--version', action='version', version='1.0b')
parser.add_argument('--version', action='version', version='1.01b')
parser.add_argument("-v", "--verbose", default=False, action="store_true",
help="verbose")
parser.add_argument("-i", "-f", "--fasta", nargs="+", type=file,
Expand All @@ -259,6 +259,8 @@ def main():
sys.stderr.write("Options: %s\n"%str(o))

#process fasta
sys.stderr.write("Homozygous assembly/ies will be written with input name + '.homozygous.fa.gz'\n")
sys.stderr.write("#file name\tgenome size\tcontigs\theterozygous size\t[%]\theterozygous contigs\t[%]\tidentity [%]\tpossible joins\thomozygous size\t[%]\thomozygous contigs\t[%]\n")
for fasta in o.fasta:
out = gzip.open(fasta.name+".homozygous.fa.gz", "w")
fasta2homozygous(out, fasta, o.identity, o.overlap, \
Expand Down

0 comments on commit d928197

Please sign in to comment.