scripts/last-dotplot
changeset 838 7b972338fe04
parent 837 2c3511fa1734
child 839 bbc6f00e683b
     1.1 --- a/scripts/last-dotplot	Tue Feb 28 18:43:58 2017 +0900
     1.2 +++ b/scripts/last-dotplot	Thu Mar 02 09:25:14 2017 +0900
     1.3 @@ -92,13 +92,13 @@
     1.4      seqLengths1 = {}
     1.5      seqLengths2 = {}
     1.6      lines = fileinput.input(fileName)
     1.7 -    for chr1, seqlen1, chr2, seqlen2, blocks in alignmentInput(lines):
     1.8 -        if not isWantedSequenceName(chr1, opts.seq1): continue
     1.9 -        if not isWantedSequenceName(chr2, opts.seq2): continue
    1.10 -        aln = chr1, chr2, blocks
    1.11 +    for seqName1, seqLen1, seqName2, seqLen2, blocks in alignmentInput(lines):
    1.12 +        if not isWantedSequenceName(seqName1, opts.seq1): continue
    1.13 +        if not isWantedSequenceName(seqName2, opts.seq2): continue
    1.14 +        aln = seqName1, seqName2, blocks
    1.15          alignments.append(aln)
    1.16 -        seqLengths1[chr1] = seqlen1
    1.17 -        seqLengths2[chr2] = seqlen2
    1.18 +        seqLengths1[seqName1] = seqLen1
    1.19 +        seqLengths2[seqName2] = seqLen2
    1.20      return alignments, seqLengths1, seqLengths2
    1.21  
    1.22  def natural_sort_key(my_string):