INNER CODE UNIT · Python

alignment_length

MariaNattestad/assemblytics · assemblytics/index.py:137

            alignment_length = abs(int(fields[2]) - int(fields[1]))
        else:
            alignment_length = abs(int(fields[3]) - int(fields[2]))
        fields.append(alignment_length)
        if fields[8] == "unique":
            uniques.append(fields)
        else:
            repetitives.append(fields)
    f.close()

    uniques.sort(key=lambda x: x[alignment_length_column],reverse=True)
    repetitives.sort(key=lambda x: x[alignment_length_column],reverse=True)
    
    fout_info = open(os.path.join(output_dir, "assemblytics_info.csv"),'w')
    fout_info.write("key,value\n")
    fout_info.write("unique alignments,%d\n" % len(uniques))
    fout_info.write("repetitive alignments,%d\n" % len(repetitives))

View source record →