insertion.py 文件源码

def igv_reads(bamfile, virus_contig):
    igv_chimeric_file = os.path.splitext(args.bamfile)[0] + ".chimeric.igv.bam"
    if os.path.exists(igv_chimeric_file):
        return igv_chimeric_file
    with pysam.Samfile(args.bamfile, "rb") as in_handle, \
         pysam.Samfile(igv_chimeric_file, "wb", template=in_handle) as out_handle:
        contig = in_handle.gettid(args.virus_contig)
        for read in in_handle:
            if skip_read(read, contig, args.virus_contig, in_handle, False):
                continue
            chrom = in_handle.getrname(read.tid)
            out_handle.write(read)
    return igv_chimeric_file