#!/usr/bin/env python
##coding=utf-8

import re
import os
import sys
import argparse
import gzip

parser = argparse.ArgumentParser(description="")
parser.add_argument('-i', '--input', help = 'the pathway of the input fastq file ', required = True)
parser.add_argument('-o', '--output', help = 'the pathway of the output vcf file,filter vcf', required = True)
argv = vars(parser.parse_args())
ifile = os.path.abspath(argv['input'].strip())
ofile = os.path.abspath(argv['output'].strip())

def lenght(inf,ouf):
    ouc=open(ouf,'w')
    with open(inf,'rb') as v:
        for vi in v:
            vic = vi.strip().split('/')
            lv=str(vic[-1]).strip().split('_')
            lvs=int(lv[0])
            lve=int(lv[1])
            lenlv=lve-lvs
            zwid='/'.join(vic)
            ouc.write(str(zwid)+"\t"+str(lenlv)+"\n")

默默加班2天
众人之中，独自摇晃的红酒杯
空无一人，独自干活的身影
略微懂得
follow you want to do
世界很大
--2023年1月15日