comparison get_codon_frequency.py @ 14:344bacf6acb8 draft

Uploaded

13:7c944fd9907e

                    if len(read) == 0:
                        continue
                    len_read = len(read.split('\t')[9])        
                    # if it's read of good length
                    if len_read == kmer and (tag in line or multi_tag not in line):
                        feat = read.split('\t')
                        seq = feat[9]
                        # if it's a reverse read
                        if feat[1] == '16' :
                            if site == "A" :

14:344bacf6acb8

                    if len(read) == 0:
                        continue
                    len_read = len(read.split('\t')[9])        
                    # if it's read of good length
                    if len_read == kmer and (tag in read or multi_tag not in read):
                        feat = read.split('\t')
                        seq = feat[9]
                        # if it's a reverse read
                        if feat[1] == '16' :
                            if site == "A" :