Breaking News: Grepper is joining You.com. Read the official announcement!

double digest fasta files

Purple Team answered on February 25, 2022 Popularity 2/10 Helpfulness 1/10

answer double digest fasta files

double digest fasta files

Comment

import re,sys
from Bio import SeqIO

def print_ReSites(id,seqence):
    pattern=r"GATC"
    seq_len=len(seqence)
    sites = [str(m.start()) for m in re.finditer(pattern,seqence)]
    sites.append(str(seq_len))
    for start,end in zip(sites,sites[1:]):
        print id+"\t"+start+"\t"+end

for seq in SeqIO.parse(sys.argv[1],"fasta"):
    print_ReSites( str( seq.id),str(seq.seq))

Popularity 2/10 Helpfulness 1/10 Language python

Source: www.biostars.org

Tags: digest double fasta python

Link to this answer
Share Copy Link

Contributed on Feb 25 2022

Purple Team

0 Answers Avg Quality 2/10

Continue with Google

By continuing, I agree that I have read and agree to Greppers's Terms of Service and Privacy Policy.

double digest fasta files

Contents

More Related Answers

double digest fasta files

Grepper

Documentation

Social

Legal

Contact

Oops, You will need to install Grepper and log-in to perform this action.