biostars 105871

disorder.py

with open('input.sequencefile') as fh:
  for line in fh:
    name, seqid, seq, disid, dis = line.split()
    print(' '.join([name, seqid]))
    print(seq)
    print(disid)
    print('Pos R')
    for i, (s, x) in enumerate(zip(seq, dis)):
      if x == 'X':
        print(' '.join([i + 1, s]))