outfile.write(header)
finalfile = open(finalfileName)
+ #TODO: QForAli - the output lines are driven by finalfile. If there are genes in the first 2 that
+ # are not in the finalfile then they will be lost.
for line in finalfile:
fields = line.strip().split()
gene = fields[0]
if __name__ == "__main__":
- main(sys.argv)
+ main(sys.argv)
\ No newline at end of file