      SUBROUTINE WTRAXY
C++
C
C WTRAXY.F V2.41 Convert X,Y pixel position to drizzle output position
C
C This task uses coefficients held in a text file to define
C geometrical distortions.
C
C This code uses the F77/VOS for IRAF i/o.
C It may be compiled and linked using a command like:
C
C xc -p stsdas tranback.f drutil.f x_tranback.x <dir>libiraf77.a
C
C History:
C
C First quick version, 17th April 1997
C
C V0.2 modified to have the centre position as parameters, 4th Sept 97
C
C V0.3 added a rotation around the centre, 19th August 1998, STScI
C      also specified "align" and chip size.
C
C V0.31 possibility to have coeffs="" for simple linear case, 25th August 98
C
C V0.4 with full set of drizzle parameter for compatibility with
C      TRANBACK and release of STECF IRAF package. 17th December 1999
C
C V0.5 modified to use new calling sequences for some DRUTIL routines
C      and hence fix bugs, November 2001
C
C V1.0 using the full geometrical possibilities available in Drizzle V2.6.
C      December 2001
C
C V1.1 return to 3 decimal places and small code tinyup
C      January 2002
C
C V1.2 using new DRIVAL call
C      Richard Hook, ST-ECF, 27th February 2002
C
C V1.3 using packaged coefficient file access and hence fixing bug
C      Richard Hook, ST-ECF, 11th March 2002
C
C V1.4 write out results to IRAF parameters and added silent flag
C      for more convenient scripting application
C      Richard Hook, ST-ECF/STScI, 4th September 2002
C
C V1.5 WTRAXY version for WCS tests,
C      Richard Hook, ST-ECF/STScI, 19th March 2003
C
C V1.6 include distortion image support,,
C      Richard Hook, ST-ECF/STScI, 26th August 2003
C
C V1.7 using double precision for almost everything
C      Richard Hook, ST-ECF/STScI, 21st October 2003
C
C V1.8 added PREFIX parameter
C       Richard Hook, ST-ECF/STScI, 23rd October 2003
C
C V1.9 added refpix support
C       Richard Hook, ST-ECF/STScI, 5th December 2003
C
C V2.0 added pixel-area output
C       Richard Hook, ST-ECF/STScI, 12th March 2004
C
C V2.1 allowed the input image header to be passed to
C      GETCOE to read coefficients, if appropriate. This is
C      part of Spitzer support.
C       Richard Hook, ST-ECF/STScI, 26th March 2004
C
C V2.2 add xylist input option.
C       Richard Hook, ST-ECF/STScI, April 2004
C
C V2.3 add output X,Y image support.
C       Richard Hook, ST-ECF/STScI, April 12th 2004.
C
C V2.4 modified to work with common callable code
C       Warren Hack, STScI, June 25th 2004.
C
C V2.5 bug fix in SIP support (was trying to read SIP header from wrong file)
C       Richard Hook, ST-ECF, October 2008.
C
C--
      IMPLICIT NONE

C Local variables
      INTEGER ISTAT,NXIN,NYIN,NXOUT,NYOUT
      INTEGER COTY,CONUM,COMAX
      PARAMETER (COMAX=100)

      DOUBLE PRECISION XCO(COMAX),YCO(COMAX)
      DOUBLE PRECISION XIN,YIN,XOUT,YOUT,X(4),Y(4),YO(4),XO(4)
      DOUBLE PRECISION RA,DEC
      INTEGER IID,IRD,DATTYP,NDIMS,DDIMS(7)
      CHARACTER*80 INIM,OUTIM
      CHARACTER*1 PREFIX
      CHARACTER*132 CHARS
      CHARACTER*45 VERS
      CHARACTER*8 CTYPE1,CTYPE2

      CHARACTER*80 OUTXIM,OUTYIM
      INTEGER MAXBUF,I,J,IDX,IDY
      PARAMETER (MAXBUF=100000)
      DOUBLE PRECISION XIA(2),YIA(2),XOA(MAXBUF),YOA(MAXBUF)
      REAL XOAS(MAXBUF),YOAS(MAXBUF)

C Geometrical parameters, the standard set
      DOUBLE PRECISION SCALE,ROT,XSH,YSH,LAM
      DOUBLE PRECISION XSH2,YSH2,ROT2,XSCALE,YSCALE
      DOUBLE PRECISION RACEN,DECCEN,XREFP,YREFP,WCSIN(8),WCSOUT(8)
      DOUBLE PRECISION ORIENT,OUTSCL,PAM
      CHARACTER*80 COEFFS,SHFTUN,XGEOIM,YGEOIM,XYLIST,LINE
      CHARACTER*8 SHFR2,SHFTFR,ALIGN,GEOMOD
      INTEGER PXG,PYG,LUN,XGDIM,YGDIM
      LOGICAL DISIM

      REAL MEMR(1)
      COMMON /MEM/MEMR

      DOUBLE PRECISION PI
      PARAMETER (PI=3.141592653589793D0)

      LOGICAL ROTFIR
      LOGICAL USEWCS
      LOGICAL ROTF2
      LOGICAL SECPAR
      LOGICAL SILENT

C-- Start of executable code
      VERS='WTRAXY Version 2.5.1 (1st May 2009)'

C Get the noise flag
      CALL UCLGSB('silent',SILENT,ISTAT)

C First announce the version
      IF(.NOT.SILENT) CALL UMSPUT('+ '//VERS,1,0,ISTAT)

C Allow WCS mode
      CALL UCLGST('geomode',GEOMOD,ISTAT)

C Get all the geometrical parameters
      CALL GTGEPA(GEOMOD,COEFFS,LAM,XGEOIM,YGEOIM,
     :             SCALE,ROT,XSH,YSH,SHFTFR,SHFTUN,ALIGN,
     :             XSCALE,YSCALE,XSH2,YSH2,ROT2,SHFR2,
     :             OUTSCL,RACEN,DECCEN,XREFP,YREFP,ORIENT)

C If so, check for image names (output is optional)
      CALL UCLGST('inimage',INIM,ISTAT)
      CALL UCLGST('refim',OUTIM,ISTAT)

      IF(GEOMOD.EQ.'wcs') THEN
         USEWCS=.TRUE.
      ELSE
         USEWCS=.FALSE.
      ENDIF

C Open the files and read the WCS from the headers
      IF(USEWCS .OR. COEFFS.EQ.'wcs') THEN
         IF(INIM.EQ.' ') THEN
            CALL UMSPUT('! No input image specified',
     :                  1,0,ISTAT)
            GO TO 99
         ELSE
            CALL UIMOPN(INIM,1,IID,ISTAT)
            IF(ISTAT.NE.0) THEN
               CALL UMSPUT('! Unable to open input image',
     :                     1,0,ISTAT)
               GO TO 99
            ELSE
               CALL UMSPUT('-Reading WCS from '//INIM,1,0,ISTAT)
            ENDIF

            CALL GETWCS(IID,WCSIN,CTYPE1,CTYPE2,ISTAT)
            IF(ISTAT.NE.0) THEN
               CALL UMSPUT(
     : '! Unable to get WCS from reference image',
     :                     1,0,ISTAT)
               GO TO 99
            ENDIF

C Get the dimensions and use them
            CALL UIMGID(IID,DATTYP,NDIMS,DDIMS,ISTAT)
            IF(NDIMS.NE.2) THEN
               CALL UMSPUT(
     :   '! Input image is not two dimensional',
     :                          1,0,ISTAT)
               GO TO 99
            ELSE
               NXIN=DDIMS(1)
               NYIN=DDIMS(2)
            ENDIF
         ENDIF
      ENDIF

      IF(USEWCS) THEN

C If we have a reference image we need to try to get its WCS for
C use as the output
C
C Note that this can have a prefix (CCCD1_1 etc)
         IF(OUTIM.NE.' ') THEN

C Get the prefix character
            CALL UCLGST('refpref',PREFIX,ISTAT)

            CALL UIMOPN(OUTIM,1,IRD,ISTAT)
            IF(ISTAT.NE.0) THEN
               CALL UMSPUT('! Unable to open output image',
     :                     1,0,ISTAT)
               GO TO 99
            ELSE
               CALL UMSPUT('-Reading WCS from '//OUTIM,1,0,ISTAT)
               IF(PREFIX.NE.' ') THEN
                  CALL UMSPUT('-Using prefix: '//PREFIX,1,0,ISTAT)
               ENDIF
            ENDIF

            IF(PREFIX.NE.' ') THEN
               CALL GTPWCS(IRD,PREFIX,WCSOUT,CTYPE1,CTYPE2,ISTAT)
            ELSE
               CALL GETWCS(IRD,WCSOUT,CTYPE1,CTYPE2,ISTAT)
            ENDIF

            IF(ISTAT.NE.0) THEN
               CALL UMSPUT(
     : '! Unable to get WCS from reference image',
     :                     1,0,ISTAT)
               GO TO 99
            ENDIF

C Get the dimensions and use them
            CALL UIMGID(IRD,DATTYP,NDIMS,DDIMS,ISTAT)
            IF(NDIMS.NE.2) THEN
               CALL UMSPUT(
     :   '! Reference image is not two dimensional',
     :                          1,0,ISTAT)
               GO TO 99
            ELSE
               NXOUT=DDIMS(1)
               NYOUT=DDIMS(2)
            ENDIF
         ELSE

C If there is no output image then we read the output WCS from the
C user supplied values and convert to a WCS representation
            CALL SETWCS(OUTSCL,ORIENT,XREFP,RACEN,
     :                  YREFP,DECCEN,WCSOUT)
            CALL UMSPUT(
     : '-Defining output WCS from user-specified parameters',
     :      1,0,ISTAT)

         ENDIF
      ENDIF

C Convert the secondary parameters into suitable units
      ROT2=ROT2*PI/180.0

      IF(SHFR2.eq.'input') THEN
         ROTF2=.FALSE.
      ELSE
         ROTF2=.TRUE.
      ENDIF

C Give a warning message if secondary parameters will have an effect
      IF(XSCALE.NE.1.0D0 .OR. YSCALE.NE.1.0D0 .OR.
     :   XSH2.NE.0.0D0 .OR. YSH2.NE.0.0D0 .OR. ROT2.NE.0.0D0) THEN
       IF(.NOT.SILENT) CALL UMSPUT(
     : '! Warning, secondary geometric transform is being used',
     :               1,0,ISTAT)
        SECPAR=.TRUE.
      ELSE
        SECPAR=.FALSE.
      ENDIF

      IF(SHFTFR.EQ.'input') THEN
         ROTFIR=.FALSE.
      ELSE
         ROTFIR=.TRUE.
      ENDIF

      IF(SHFTUN.EQ.'output') THEN
         XSH=XSH*SCALE
         YSH=YSH*SCALE
      ENDIF

C Get the geometric distortion coefficient information
C Note that an input image is not always available so
C coeffs="header"  or coeffs="wcs" may fail
C
C A bug occured here in earlier versions when coeffs="wcs" - the SIP coefficients were read from the
C output rather than the input.
C
C Richard Hook, ST-ECF, October 2008
      CALL GETGEO(COEFFS,IID,LAM,
     :            COTY,COMAX,CONUM,XCO,YCO,ISTAT)
      IF(ISTAT.NE.0) THEN
         CALL UMSPUT(
     :  '! WTRAXY: Error, failed to get distortion coefficients',
     :  1,0,ISTAT)
          GO TO 99
      ENDIF

C  Read in any distortion image corrections, if specified
C  PXG,PYG will default to 2x2 arrays when DISIM returns as FALSE.
C
      CALL GGEOIM(XGEOIM,YGEOIM,PXG,PYG,XGDIM,YGDIM,DISIM)

C Also get the chip size - in and out
C Not needed for GEOMOD='wcs'
      IF(.NOT.USEWCS) THEN
         CALL UCLGSI('nxin',NXIN,ISTAT)
         CALL UCLGSI('nyin',NYIN,ISTAT)
      ENDIF

      IF(.NOT.USEWCS .AND. OUTIM.EQ.' ') THEN
         CALL UCLGSI('nxout',NXOUT,ISTAT)
         CALL UCLGSI('nyout',NYOUT,ISTAT)
      ENDIF

C Convert to radians
      ROT=ROT*PI/180.0D0

C Either we have a list of X,Y positions or a single one
      CALL UCLGST('xylist',XYLIST,ISTAT)
      IF(XYLIST.NE.' ') THEN
         CALL UFOPEN(XYLIST,1,LUN,ISTAT)
         IF(ISTAT.NE.0) THEN
            CALL UMSPUT('! Unable to open xylist file',1,0,
     :                  ISTAT)
            ISTAT=1
            RETURN
         ELSE
            CALL UMSPUT('-Opening xylist file: '//XYLIST,
     :                      1,0,ISTAT)
         ENDIF
      ENDIF

C Get the X,Y pixel position to be transformed
      IF(XYLIST.EQ.' ') THEN
         CALL UCLGSD('xin',XIN,ISTAT)
         CALL UCLGSD('yin',YIN,ISTAT)
      ENDIF

C Transform the position
C Note that this is exactly the same as that in Drizzle
C Double precision version
      DO WHILE(.TRUE.)

         IF(XYLIST.NE.' ') THEN
            CALL UFGLIN(LUN,LINE,ISTAT)
            IF(ISTAT.NE.0) GO TO 88
            READ(LINE,*) XIN,YIN
         ENDIF

C         IF(NINT(XIN).GT.1 .AND. NINT(XIN).LT.NXIN-1 .AND.
C     :         NINT(YIN).GT.1 .AND. NINT(YIN).LT.NYIN-1) THEN

            CALL DRIVAL(XIN,YIN,1,NXIN,NYIN,NXOUT,NYOUT,.FALSE.,
     :            XSH,YSH,ROT,SCALE,ALIGN,ROTFIR,
     :            SECPAR,XSH2,YSH2,ROT2,XSCALE,YSCALE,SHFR2,ROTF2,
     :            USEWCS,WCSIN,WCSOUT,
     :            COTY,CONUM,XCO,YCO,
     :            DISIM,MEMR(PXG),MEMR(PYG),XGDIM,YGDIM,XOUT,YOUT)
C         ELSE
C            CALL DRIVAL(XIN,YIN,1,NXIN,NYIN,NXOUT,NYOUT,.FALSE.,
C     :            XSH,YSH,ROT,SCALE,ALIGN,ROTFIR,
C     :            SECPAR,XSH2,YSH2,ROT2,XSCALE,YSCALE,SHFR2,ROTF2,
C     :            USEWCS,WCSIN,WCSOUT,
C     :            COTY,CONUM,XCO,YCO,
C     :            .FALSE.,MEMR(PXG),MEMR(PYG),XGDIM,YGDIM,XOUT,YOUT)
C         ENDIF

C Transform a unit square to calculate the pixel area
         X(1)=XIN-0.5D0
         X(2)=XIN+0.5D0
         X(3)=XIN+0.5D0
         X(4)=XIN-0.5D0
         Y(1)=YIN-0.5D0
         Y(2)=YIN-0.5D0
         Y(3)=YIN+0.5D0
         Y(4)=YIN+0.5D0

         CALL DRIVAL(X,Y,4,NXIN,NYIN,NXOUT,NYOUT,.FALSE.,
     :            0.0D0,0.0D0,0.0D0,1.0D0,ALIGN,ROTFIR,
     :            .FALSE.,0.0D0,0.0D0,0.0D0,1.0D0,1.0D0,SHFR2,ROTF2,
     :            .FALSE.,WCSIN,WCSOUT,
     :            COTY,CONUM,XCO,YCO,
     :            .FALSE.,MEMR(PXG),MEMR(PYG),XGDIM,YGDIM,XO,YO)

         PAM=DABS(0.5*((XO(2)-XO(4))*(YO(1)-YO(3)) -
     :              (XO(1)-XO(3))*(YO(2)-YO(4))))

C Write out the result
         IF(.NOT.SILENT) THEN
            WRITE(CHARS,
     :  '('' Xin,Yin: '',2F13.6'' Xout,Yout: '',2F13.6,
     :    '' PAM: '',F12.5)')
     : XIN,YIN,XOUT,YOUT,PAM
            CALL UMSPUT(CHARS,1,0,ISTAT)

C If in WCS mode also report the sky position
            IF(USEWCS) THEN
               CALL XY2RD(XOUT,YOUT,RA,DEC,WCSOUT)
               WRITE(CHARS,
     : '('' RA: '',F12.7,'' Dec: '',F12.7,'' (deg)'')')
     :    RA,DEC
               CALL UMSPUT(CHARS,1,0,ISTAT)
            ENDIF
         ENDIF

         IF(XYLIST.EQ.' ') GO TO 88
      ENDDO

 88   CONTINUE

C Also write the result to the CL
      IF(XYLIST.EQ.' ') THEN
         CALL UCLPSD('xout',XOUT,ISTAT)
         CALL UCLPSD('yout',YOUT,ISTAT)
      ELSE
         CALL UFCLOS(LUN,ISTAT)
      ENDIF

C Image creation support
      CALL UCLGST('outxim',OUTXIM,ISTAT)
      CALL UCLGST('outyim',OUTYIM,ISTAT)

C Skip this if not image option
      IF(OUTXIM.EQ.' ' .AND. OUTYIM.EQ.' ') GO TO 99

      DDIMS(1)=NXIN
      DDIMS(2)=NYIN

C Check for huge images
      IF(NXIN.GT.MAXBUF .OR. NYIN.GT.MAXBUF) THEN
         CALL UMSPUT('! Image too large for output image option',
     :               1,0,ISTAT)
         GO TO 99
      ENDIF

C Create images if not null
      IF(OUTXIM.NE.' ') THEN
         CALL UIMCRE(OUTXIM,6,2,DDIMS,IDX,ISTAT)
         IF(ISTAT.NE.0) THEN
            CALL UMSPUT('! Failed to created output X image',
     :                  1,0,ISTAT)
            GO TO 99
         ELSE
            CALL UMSPUT('-Created output X position image: '//
     :                  OUTXIM,1,0,ISTAT)
         ENDIF
      ENDIF

      IF(OUTYIM.NE.' ') THEN
         CALL UIMCRE(OUTYIM,6,2,DDIMS,IDY,ISTAT)
         IF(ISTAT.NE.0) THEN
            CALL UMSPUT('! Failed to created output Y image',
     :                  1,0,ISTAT)
            GO TO 99
         ELSE
            CALL UMSPUT('-Created output Y position image: '//
     :                  OUTYIM,1,0,ISTAT)
         ENDIF
      ENDIF

C Loop over the input image grid
C We are going to use the "regular" feature so need to
C initialise some things
      XIA(1)=1.0D0
      XIA(2)=0.0D0
      YIA(2)=0.0D0

      DO J=1,NYIN
         YIA(1)=DBLE(J)

         CALL DRIVAL(XIA,YIA,NXIN,NXIN,NYIN,NXOUT,NYOUT,.TRUE.,
     :            XSH,YSH,ROT,SCALE,ALIGN,ROTFIR,
     :            SECPAR,XSH2,YSH2,ROT2,XSCALE,YSCALE,SHFR2,ROTF2,
     :            USEWCS,WCSIN,WCSOUT,
     :            COTY,CONUM,XCO,YCO,
     :            DISIM,MEMR(PXG),MEMR(PYG),XGDIM,YGDIM,XOA,YOA)

C Copy to single precision
         DO I=1,NXIN
            XOAS(I)=SNGL(XOA(I))
            YOAS(I)=SNGL(YOA(I))
         ENDDO

C Write out the lines
         IF(OUTXIM.NE.' ') CALL UIPL2R(IDX,J,XOAS,ISTAT)
         IF(OUTYIM.NE.' ') CALL UIPL2R(IDY,J,YOAS,ISTAT)
      ENDDO

C Close the files
      IF(OUTXIM.NE.' ') CALL UIMCLO(IDX,ISTAT)
      IF(OUTYIM.NE.' ') CALL UIMCLO(IDY,ISTAT)

C End of main TRAXY module
 99   CONTINUE

C If allocated, free the distortion image memory
      IF(PXG.NE.0) CALL UDMFRE(PXG,6,ISTAT)
      IF(PYG.NE.0) CALL UDMFRE(PYG,6,ISTAT)

      RETURN
      END


      PRINT 25 / 5
      END


c               Copyright 1993 Colorado State University
c                       All Rights Reserved

C ... READBLK.F

      subroutine readblk

      implicit none
      include 'chrvar.inc'
      include 'const.inc'
      include 'schvar.inc'
      include 'seq.inc'
      include 'site.inc'
      include 't0par.inc'
      include 'timvar.inc'
      include 'zztim.inc'

c ... Reads the next block of events from the schedule file

c ... Function declarations
      integer   readint
      double precision readreal
      character readstring*20

c ... Local variables
      integer   blknum, ctemp, ii, modt, ttemp, yrnum
      integer   prevsys
      double precision pltmo

c ... Saved variables
      save      prevsys

c ... Local variables for new input block
      integer   ststrt, stend, lch, istat
      logical   dbgp
      character buffr*90

      parameter (dbgp = .FALSE.)
c      parameter (dbgp = .TRUE.)

c ... Set the starting year of the block
c ... Using ANINT which rounds to the nearest whole number.  -rm
      strtyr = anint(time)

c ... Read the block number, ending year, and # of years set up
      read(15,*) blknum
      read(15,*) blktnd
      read(15,*) rptyrs
      read(15,*) strplt
      read(15,*) pltmo
      read(15,10) dtpl
10    format(f6.3)

c ... Add 1 more to blktnd for January of the ending year
      blktnd = blktnd + 1

c ... Reset the dtpl value; the units are years, but it should
c ... have a value representing an integral number of months as
c ... accurately as possible.
      if (dtpl .lt. .01) then

c ..... Reset to monthly if lower than lowest limit
        dtpl = dt
      else

c ..... Rounding
        modt = int(12.*dtpl+.05)
        dtpl = dble(modt)/12.
      endif

c ... Set tplt, the next time to spit out output
c      pltmo = pltmo - 1
      tplt = strplt + dble(pltmo)/12.
c ... Round tplt to 2 digits of precision, cak - 03/22/02
      tplt =  anint(tplt * 100.0) / 100.0
c      if (tplt-int(tplt) .eq. 0.0) then
c        tplt = tplt + dtpl
c      endif

c ... Determine the weather type and call wthini
c ... to initialize the weather for this block
      read(15,20) wthr
20    format(a1)
      if (wthr .eq. 'F') then
        read(15,30) wthnam
30      format(a20)
c ..... This call just opens the weather file for RAMS/Century linkage,
c ..... mdh 12/96
        call wthini
      endif

      if (dbgp) write(*,'(/a,i5)') ' reading block ', blknum

c=======================================================================
c
c   The following block reads the events from the schedule file.
c   Uses subroutines that allow free form number and string input.
c   It also simplifies the input format by relaxing the year input requirement
c   (it will assume the same year as the previous command)
c
c  History
c    Modified : 04/03/95 : Kendrick Killian
c              added error checking on Century events and dates
c    Written : 1/25/95 : Kendrick Killian
c
c  New External References
c    character function READSTRING (UNIT, BUFFR, LCH)
c    double precision function      READREAL   (UNIT, BUFFR, LCH, ISTAT)
c    integer function   READINT    (UNIT, BUFFR, LCH, ISTAT)
c
c  Additional comments
c    1) input commands are converted to upper case
c    2) input fields can be separated by commas, white space or endline
c    3) UNIX like comments can be inserted in the input stream
c       - # is the comment character
c       - comments extend to the end of the line
c    4) Century events can be entered without the year.  The year is
c        pulled from the previous event.
c    5) Century commands with NO date arguments generate a fatal error
c
c========================== BEGIN MODIFIED CODE ========================

c ... Read and save the events for this block
      ttlind = 0
c ... set the buffer pointer (lch) to beyond the buffer length to
c ... force the "get" routines to read an input line
      lch = len(buffr)+1
40    continue
      yrnum = readint(15,buffr,lch,istat)
c ... check for the existance of the date arguments
      if (istat .le. 0) then
        write(*,*) ' fatal error: missing date information :',buffr
        stop 'Abnormal Termination'
      endif

c ... check for end of block
      if (yrnum .ne. -999) then
        ttlind = ttlind + 1
c ..... Fill in the date parameters
        timary(ttlind, 1) = yrnum
        timary(ttlind, 2) = readint(15,buffr,lch,istat)

c ..... check for the explicit month argument
        if (istat.eq.0) then
          if (ttlind.gt.1) then
            timary(ttlind, 1) = timary(ttlind-1, 1)
            timary(ttlind, 2) = yrnum
c ....... bad input
          else
            write(*,*) ' missing year designator :',buffr
            stop 'Abnormal Termination'
          endif
        endif

c ..... check date input consistancy
c ..... use blknum to mark an error  (blknum < 0  bad input)
        if (timary(ttlind, 1) .le.0) then
          write(*,'(a,i5/a)') ' fatal error: year less than zero ',
     &          timary(ttlind, 1), buffr
          timary(ttlind, 1) = 1
          blknum = -abs(blknum)
        elseif (timary(ttlind, 1) .gt.rptyrs) then
          write(*,'(a,i5/a)') ' fatal error: year greater than '//
     &          ' rotation length ',timary(ttlind, 1), buffr
          timary(ttlind, 1) = rptyrs
          blknum = -abs(blknum)
        elseif ((timary(ttlind,2).lt.0).or.
     &          (timary(ttlind,2).gt.365)) then
          write(*,'(a,i5/a)') ' fatal error: illegal day of year ',
     &          timary(ttlind, 2), buffr
          timary(ttlind, 2) = 0
          blknum = -abs(blknum)
        endif
        if (ttlind.gt.1) then
          if ((timary(ttlind-1,1).gt.timary(ttlind,1)  ) .or.
     &       ((timary(ttlind,1)  .eq.timary(ttlind-1,1)) .and.
     &       (timary(ttlind-1,2).gt.timary(ttlind,2)  )))     then
            write(*,'(a/a)') ' fatal error: dates out of sequence ',
     &            buffr
            write(*,'(3(i3,a),i3)') timary(ttlind-1,1),'/',
     &            timary(ttlind-1,2)," comes before ",
     &            timary(ttlind,1),'/',timary(ttlind,2)
            timary(ttlind, 1) = timary(ttlind-1, 1)
            timary(ttlind, 2) = timary(ttlind-1, 2)
            blknum = -abs(blknum)
          endif
        endif

c ..... Read and record the command
        cmdary(ttlind) = readstring(15,buffr,lch,ststrt,stend)

c ..... Fill in additional information (if required)
        if(INDEX('CROP CULT FERT FIRE GRAZ HARV IRRI OMAD TREE TREM ',
     &           cmdary(ttlind)(1:4)) .ne. 0) then

c ....... get string options (types)
          typary(ttlind) = readstring(15,buffr,lch,ststrt,stend)
        elseif ('EROD' .eq. cmdary(ttlind)) then

c ....... get real number for erosion
          fltary(ttlind, 1) = readreal(15,buffr,lch,istat)
          if (istat .eq. 0)  write(*,*)
     &        "bad errosion  ",fltary(ttlind, 1),buffr,lch,istat
        else

c ....... check for valid commands
          if(INDEX('PLTM FRST LAST SENM TFST TLST ',
     &             cmdary(ttlind)(1:4)) .eq. 0) then
            write(*,'(2a/a)')' fatal error: unknown CENTURY command : ',
     &            cmdary(ttlind)(1:4), buffr
            cmdary(ttlind)(1:4) = 'err'
            blknum = -abs(blknum)
          endif
        endif
c ..... ***************************************************
        if (dbgp) write(*,'(2(a,i5),i3,a,a4,3a,g12.5)') ' Read event #',
     &    ttlind,'   ->', timary(ttlind, 1),timary(ttlind, 2),
     &    ' "',cmdary(ttlind),'" "', typary(ttlind),'"',
     &    fltary(ttlind, 1)
c ..... ***************************************************
        goto 40
      endif

c ... abort on illegal event input
      if(blknum .lt. 0) stop 'Abnormal Termination'
c ... ============================ End Modified Code ========================

c ... Reset evtptr, the array index into cmdary, timary, typary
      evtptr = 1

c ... Set up cursys, the current system(s) in use if
c ... this is the first block read
c ... cursys: CRPSYS = crop/grass; FORSYS = forest; SAVSYS = savanna
c ... Reset system for each block to handle changing systems.

c ... Store the value of the previous system
      if (blknum .gt. 1) then
        prevsys = cursys
      endif

c ... Initialize the temporary variables
      ctemp = 0
      ttemp = 0
      do 80 ii = 1, ttlind
        if (cmdary(ii) .eq. 'CROP') then
          ctemp = CRPSYS
        else if (cmdary(ii) .eq. 'TREE') then
          ttemp = FORSYS
        endif
80    continue
      cursys = ctemp + ttemp

c ... If no crop or forest designation was given in the
c ... first block, use the initial system designated
      if (cursys .eq. 0) then
        if (blknum .eq. 1) then
          if (initcp .ne. ' ') then
            ctemp = CRPSYS
          endif
          if (initre .ne. ' ') then
            ttemp = FORSYS
          endif
          cursys = ctemp + ttemp
        else

c ....... If no crop or forest designation was given in this
c ....... block and it is not the first, use the previous
c ....... block's system designation.
          cursys = prevsys
        endif
      endif

c ... Check if decsys needs to be updated.
      if (decsys .eq. CRPSYS .and. ttemp .eq. FORSYS) then
        decsys = FORSYS
      endif

      return
      end


ACCEPTABLE LEFT PRIMERS
                                  0-based     #               self  self  qual-
   # sequence                       start ln  N   GC%     Tm   any   end   lity
   0 CCCTACATACAGCCAAGGGA              21 20  0 55.00 59.948  4.00  0.00  0.052
   1 TCCCTACATACAGCCAAGGG              20 20  0 55.00 59.948  4.00  1.00  0.052
   2 CTAAGAATGTGGCCAAGGGA             122 20  0 50.00 60.066  6.00  0.00  0.066
   3 TCCCTAAGAATGTGGCCAAG             119 20  0 50.00 60.066  6.00  2.00  0.066
   4 GTTCCCTAAGAATGTGGCCA             117 20  0 50.00 59.933  6.00  6.00  0.067
   5 TGTTCCCTAAGAATGTGGCC             116 20  0 50.00 59.933  4.00  4.00  0.067
   6 GACTGCCACCGTGTACACAC              39 20  0 60.00 60.080  8.00  6.00  0.080
   7 ACATACAGCCAAGGGACTGC              25 20  0 55.00 60.142  3.00  2.00  0.142
   8 CCGCCATGTACACACTGAAG              95 20  0 55.00 60.175  6.00  2.00  0.175
   9 TAAGAATGTGGCCAAGGGAA             123 20  0 45.00 60.439  6.00  0.00  0.439
  10 TTCCCTAAGAATGTGGCCAA             118 20  0 45.00 60.439  6.00  4.00  0.439
  11 GCCACCGTGTACACACTGAA              43 20  0 55.00 60.639  8.00  4.00  0.639
  12 AAGAATGTGGCCAAGGGAAT             124 20  0 45.00 60.693  6.00  2.00  0.693
  13 GAAATGTGGTTAAGGGACCG              75 20  0 50.00 59.288  4.00  2.00  0.712
  14 TGGGTGTACCACACACTTCC             183 20  0 55.00 59.283  6.00  0.00  0.717
  15 AATGCACTGAAGACCTTCCC             141 20  0 50.00 59.141  6.00  2.00  0.859
  16 GGGAATGCACTGAAGACCTT             138 20  0 50.00 59.141  4.00  4.00  0.859
  17 AAGGGAATGCACTGAAGACC             136 20  0 50.00 59.141  4.00  3.00  0.859
  18 ACTGCCACCGTGTACACACT              40 20  0 55.00 59.101  8.00  6.00  0.899
  19 GGACTGCCACCGTGTACAC               38 19  0 63.16 60.021  8.00  8.00  1.021
  20 ATACAGCCAAGGGACTGCC               27 19  0 57.89 60.088  3.00  3.00  1.088
  21 AAATGTGGTTAAGGGACCGC              76 20  0 50.00 61.099  4.00  2.00  1.099
  22 TCCTACAGCTGTGGGAAAATG            200 21  0 47.62 60.118 12.00  0.00  1.118
  23 AATGCACTGAAGACCTTCCCT            141 21  0 47.62 60.125  6.00  4.00  1.125
  24 AGGGAATGCACTGAAGACCTT            137 21  0 47.62 60.125  4.00  4.00  1.125
  25 AAGGGAATGCACTGAAGACCT            136 21  0 47.62 60.125  4.00  2.00  1.125
  26 TGCACTGAAGACCTTCCCTAA            143 21  0 47.62 59.861  6.00  1.00  1.139
  27 TACATACAGCCAAGGGACTGC             24 21  0 52.38 60.147  3.00  2.00  1.147
  28 CCAAGGGAATGCACTGAAGA             134 20  0 50.00 61.176  4.00  0.00  1.176
  29 GTGGGTGTACCACACACTTCC            182 21  0 57.14 60.189  8.00  0.00  1.189
  30 AATGTGGTTAAGGGACCGC               77 19  0 52.63 59.805  4.00  2.00  1.195
  31 ATGTTCCCTAAGAATGTGGCC            115 21  0 47.62 60.201  4.00  4.00  1.201
  32 TGCCACCGTGTACACACTG               42 19  0 57.89 60.217  8.00  8.00  1.217
  33 AGAAATGTGGTTAAGGGACCG             74 21  0 47.62 60.228  4.00  2.00  1.228
  34 ATGCACTGAAGACCTTCCCT             142 20  0 50.00 58.744  6.00  4.00  1.256
  35 AGGGAATGCACTGAAGACCT             137 20  0 50.00 58.744  4.00  2.00  1.256
  36 CTGCCACCGTGTACACACTG              41 20  0 60.00 61.274  8.00  8.00  1.274
  37 CATTGTGGGTGTACCACACA             178 20  0 50.00 58.674  8.00  8.00  1.326
  38 TGGGTGTACCACACACTTCCT            183 21  0 52.38 60.328  6.00  0.00  1.328
  39 CATTGTGGGTGTACCACACAC            178 21  0 52.38 59.632  8.00  8.00  1.368
  40 ACATTGTGGGTGTACCACACA            177 21  0 47.62 59.632  8.00  8.00  1.368
  41 CCGTGTACACACTGAACACGT             47 21  0 52.38 59.573  9.00  7.00  1.427
  42 AAGAATGTGGCCAAGGGAA              124 19  0 47.37 60.445  6.00  0.00  1.445
  43 CACGTTCCCTAAGAAATGTGG             63 21  0 47.62 59.480  8.00  6.00  1.520
  44 CCACCGTGTACACACTGAACA             44 21  0 52.38 60.530  8.00  5.00  1.530
  45 TGCACTGAAGACCTTCCCTA             143 20  0 50.00 58.442  6.00  2.00  1.558
  46 AGAATGTGGCCAAGGGAAT              125 19  0 47.37 59.351  6.00  2.00  1.649
  47 TGTGGGTGTACCACACACTTC            181 21  0 52.38 59.347  8.00  2.00  1.653
  48 TAAGAATGTGGCCAAGGGAAT            123 21  0 42.86 60.678  6.00  2.00  1.678
  49 CAAGGGAATGCACTGAAGAC             135 20  0 50.00 58.287  4.00  1.00  1.713
  50 CGCCATGTACACACTGAAGA              96 20  0 50.00 58.276  6.00  1.00  1.724
  51 CCTAAGAATGTGGCCAAGGG             121 20  0 55.00 61.732  6.00  1.00  1.732
  52 CCCTAAGAATGTGGCCAAGG             120 20  0 55.00 61.732  6.00  2.00  1.732
  53 CATACAGCCAAGGGACTGC               26 19  0 57.89 59.256  3.00  2.00  1.744
  54 CCCTACATACAGCCAAGGGAC             21 21  0 57.14 60.748  4.00  1.00  1.748
  55 CCAAGGGAATGCACTGAAG              134 19  0 52.63 59.226  4.00  2.00  1.774
  56 GCCACCGTGTACACACTGA               43 19  0 57.89 59.114  8.00  6.00  1.886
  57 GACTGCCACCGTGTACACA               39 19  0 57.89 59.114  8.00  6.00  1.886
  58 CCGCCATGTACACACTGAA               95 19  0 52.63 59.097  6.00  1.00  1.903
  59 TAAGAATGTGGCCAAGGGA              123 19  0 47.37 59.061  6.00  0.00  1.939
  60 TCCCTAAGAATGTGGCCAA              119 19  0 47.37 59.061  6.00  4.00  1.939
  61 TTCCCTAAGAATGTGGCCA              118 19  0 47.37 59.061  6.00  6.00  1.939
  62 ACACACTGAACACGTTCCCT              53 20  0 50.00 58.039  6.00  4.00  1.961
  63 AGGGACTGCCACCGTGTAC               36 19  0 63.16 61.004  5.00  4.00  2.004
  64 GAATGCACTGAAGACCTTCCC            140 21  0 52.38 61.014  6.00  2.00  2.014
  65 GGAATGCACTGAAGACCTTCC            139 21  0 52.38 61.014  6.00  4.00  2.014
  66 GGGAATGCACTGAAGACCTTC            138 21  0 52.38 61.014  6.00  6.00  2.014
  67 GGGACTGCCACCGTGTAC                37 18  0 66.67 59.962  5.00  4.00  2.038
  68 GACTGCCACCGTGTACACACT             39 21  0 57.14 61.098  8.00  6.00  2.098
  69 AGGGACTGCCACCGTGTA                36 18  0 61.11 60.113  5.00  2.00  2.113
  70 TGCACTGAAGACCTTCCCTAAT           143 22  0 45.45 60.130  6.00  2.00  2.130
  71 ATGCACTGAAGACCTTCCCTAA           142 22  0 45.45 60.130  6.00  1.00  2.130
  72 AATGCACTGAAGACCTTCCCTA           141 22  0 45.45 60.130  6.00  2.00  2.130
  73 AGAATGTGGCCAAGGGAATG             125 20  0 50.00 62.168  6.00  0.00  2.168
  74 TTGTGGGTGTACCACACACTT            180 21  0 47.62 58.817  8.00  5.00  2.183
  75 AAGACATTGTGGGTGTACCACA           174 22  0 45.45 60.193  6.00  6.00  2.193
  76 ATGCACTGAAGACCTTCCCTA            142 21  0 47.62 58.803  6.00  2.00  2.197
  77 TACAGCCAAGGGACTGCC                28 18  0 61.11 59.787  3.00  3.00  2.213
  78 GTTCCCTAAGAATGTGGCCAA            117 21  0 47.62 61.214  6.00  4.00  2.214
  79 TTCCCTACATACAGCCAAGGG             19 21  0 52.38 61.221  4.00  1.00  2.221
  80 AGACATTGTGGGTGTACCACA            175 21  0 47.62 58.765  6.00  6.00  2.235
  81 GAATGTGGCCAAGGGAATG              126 19  0 52.63 61.265  6.00  0.00  2.265
  82 CTTCCCTACATACAGCCAAGG             18 21  0 52.38 58.730  4.00  2.00  2.270
  83 CACGTTCCCTAAGAAATGTGGT            63 22  0 45.45 60.275  8.00  1.00  2.275
  84 ACACGTTCCCTAAGAAATGTGG            62 22  0 45.45 60.275  8.00  6.00  2.275
  85 TTCCCTACATACAGCCAAGG              19 20  0 50.00 57.700  3.00  2.00  2.300
  86 CACACTGAACACGTTCCCTAAG            54 22  0 50.00 59.696  6.00  3.00  2.304
  87 TGCCACCGTGTACACACTGA              42 20  0 55.00 62.305  8.00  6.00  2.305
  88 CACACTGAACACGTTCCCTAA             54 21  0 47.62 58.687  6.00  1.00  2.313
  89 AAGGGACTGCCACCGTGTAC              35 20  0 60.00 62.314  5.00  4.00  2.314
  90 GCCGCCATGTACACACTG                94 18  0 61.11 59.685  6.00  4.00  2.315
  91 TGGGTGTACCACACACTTCCTA           183 22  0 50.00 60.326  6.00  2.00  2.326
  92 AGACATTGTGGGTGTACCACAC           175 22  0 50.00 59.670  6.00  6.00  2.330
  93 CTAAGAATGTGGCCAAGGGAA            122 21  0 47.62 61.331  6.00  0.00  2.331
  94 TTCCCTAAGAATGTGGCCAAG            118 21  0 47.62 61.331  6.00  2.00  2.331
  95 CGCCATGTACACACTGAAGAT             96 21  0 47.62 58.664  6.00  3.00  2.336
  96 GACATTGTGGGTGTACCACAC            176 21  0 52.38 58.609  6.00  6.00  2.391
  97 GACATTGTGGGTGTACCACA             176 20  0 50.00 57.608  6.00  6.00  2.392
  98 ACACACTGAACACGTTCCCTAA            53 22  0 45.45 59.562  6.00  1.00  2.438
  99 ACCGTGTACACACTGAACACGT            46 22  0 50.00 60.445  9.00  7.00  2.445
 100 ACTTCCCTACATACAGCCAAGG            17 22  0 50.00 59.543  4.00  2.00  2.457
 101 GATGTTCCCTAAGAATGTGGC            114 21  0 47.62 58.517  3.00  2.00  2.483
 102 GCCACCGTGTACACACTGAAC             43 21  0 57.14 61.486  8.00  4.00  2.486
 103 ACATTGTGGGTGTACCACACAC           177 22  0 50.00 60.510  8.00  8.00  2.510
 104 AGATGTTCCCTAAGAATGTGGC           113 22  0 45.45 59.487  3.00  2.00  2.513
 105 AAGGGACTGCCACCGTGTA               35 19  0 57.89 61.530  5.00  2.00  2.530
 106 CCACCGTGTACACACTGAAC              44 20  0 55.00 57.400  8.00  3.00  2.600
 107 TGCACTGAAGACCTTCCCT              143 19  0 52.63 58.364  6.00  4.00  2.636
 108 CATTGTGGGTGTACCACACACT           178 22  0 50.00 60.644  8.00  8.00  2.644
 109 TCACACTCACTGAAGACTTCCC             2 22  0 50.00 59.350  7.00  3.00  2.650
 110 GGGTGTACCACACACTTCCT             184 20  0 55.00 57.315  5.00  0.00  2.685
 111 AATGTGGCCAAGGGAATG               127 18  0 50.00 59.296  6.00  0.00  2.704
 112 GAATGCACTGAAGACCTTCC             140 20  0 50.00 57.288  6.00  4.00  2.712
 113 GGAATGCACTGAAGACCTTC             139 20  0 50.00 57.288  6.00  6.00  2.712
 114 CACTGAACACGTTCCCTAAGAA            56 22  0 45.45 59.281  6.00  2.00  2.719
 115 TCCTACAGCTGTGGGAAAAT             200 20  0 45.00 57.268 12.00  2.00  2.732
 116 TCCCTACATACAGCCAAGGGA             20 21  0 52.38 61.748  5.00  0.00  2.748
 117 TTGTGGGTGTACCACACACTTC           180 22  0 50.00 60.765  8.00  3.00  2.765
 118 TGTGGGTGTACCACACACTT             181 20  0 50.00 57.229  8.00  5.00  2.771
 119 TTGTGGGTGTACCACACACT             180 20  0 50.00 57.229  8.00  6.00  2.771
 120 GCCGCCATGTACACACTGA               94 19  0 57.89 61.775  6.00  2.00  2.775
 121 GCCAAGGGAATGCACTGA               133 18  0 55.56 60.784  4.00  1.00  2.784
 122 CCTACATACAGCCAAGGGACT             22 21  0 52.38 58.214  3.00  1.00  2.786
 123 CACACACTTCCTACAGCTGTGG           192 22  0 54.55 60.796  8.00  8.00  2.796
 124 CCACACACTTCCTACAGCTGTG           191 22  0 54.55 60.796  8.00  6.00  2.796
 125 CGTGTACACACTGAACACGTTC            48 22  0 50.00 59.186  9.00  6.00  2.814
 126 ATTGTGGGTGTACCACACACTT           179 22  0 45.45 59.161  8.00  6.00  2.839
 127 CTAAGAATGTGGCCAAGGG              122 19  0 52.63 58.159  6.00  1.00  2.841
 128 CCTAAGAATGTGGCCAAGG              121 19  0 52.63 58.159  6.00  2.00  2.841
 129 CCCTAAGAATGTGGCCAAG              120 19  0 52.63 58.159  6.00  2.00  2.841
 130 CCTACATACAGCCAAGGGAC              22 20  0 55.00 57.143  3.00  1.00  2.857
 131 CACACTGAACACGTTCCCTA              54 20  0 50.00 57.141  6.00  2.00  2.859
 132 TGAACACGTTCCCTAAGAAATG            59 22  0 40.91 59.131  6.00  4.00  2.869
 133 ACACACTGAACACGTTCCCTA             53 21  0 47.62 58.130  6.00  2.00  2.870
 134 TACACACTGAACACGTTCCCT             52 21  0 47.62 58.130  6.00  4.00  2.870
 135 CGTTCCCTAAGAAATGTGGTT             65 21  0 42.86 58.111  4.00  0.00  2.889
 136 CTGCCACCGTGTACACACT               41 19  0 57.89 58.109  8.00  6.00  2.891
 137 TGTACACACTGAACACGTTCCC            50 22  0 50.00 60.906  6.00  2.00  2.906
 138 ATACAGCCAAGGGACTGCCA              27 20  0 55.00 62.908  3.00  3.00  2.908
 139 CATACAGCCAAGGGACTGCC              26 20  0 60.00 62.908  3.00  3.00  2.908
 140 CCTACATACAGCCAAGGGACTG            22 22  0 54.55 60.913  3.00  3.00  2.913
 141 CCGTGTACACACTGAACACGTT            47 22  0 50.00 60.956  9.00  5.00  2.956
 142 GGACTGCCACCGTGTACA                38 18  0 61.11 59.043  6.00  6.00  2.957
 143 CAAGGGAATGCACTGAAGACC            135 21  0 52.38 61.957  4.00  3.00  2.957
 144 CCAAGGGAATGCACTGAAGAC            134 21  0 52.38 61.957  4.00  1.00  2.957
 145 GTACACACTGAACACGTTCCCT            51 22  0 50.00 59.041  6.00  4.00  2.959
 146 CCCTACATACAGCCAAGGG               21 19  0 57.89 58.034  4.00  1.00  2.966
 147 AGAATGTGGCCAAGGGAA               125 18  0 50.00 59.009  6.00  0.00  2.991
 148 AAGAATGTGGCCAAGGGA               124 18  0 50.00 59.009  6.00  0.00  2.991
 149 GTACACACTGAACACGTTCCC             51 21  0 52.38 57.973  6.00  2.00  3.027
 150 CTGAACACGTTCCCTAAGAAATG           58 23  0 43.48 60.042  6.00  4.00  3.042
 151 CTACATACAGCCAAGGGACTGC            23 22  0 54.55 61.052  3.00  2.00  3.052
 152 ATGTGGTTAAGGGACCGCC               78 19  0 57.89 62.052  4.00  1.00  3.052
 153 ACTGCCACCGTGTACACAC               40 19  0 57.89 57.937  8.00  6.00  3.063
 154 ACGTTCCCTAAGAAATGTGGTT            64 22  0 40.91 58.937  6.00  0.00  3.063
 155 CCCTAAGAAATGTGGTTAAGGG            69 22  0 45.45 58.926  4.00  4.00  3.074
 156 GAACACGTTCCCTAAGAAATGTG           60 23  0 43.48 59.920  8.00  8.00  3.080
 157 TGAACACGTTCCCTAAGAAATGT           59 23  0 39.13 59.920  6.00  6.00  3.080
 158 ACACTGAACACGTTCCCTAAGAA           55 23  0 43.48 60.082  6.00  2.00  3.082
 159 CCGCCATGTACACACTGAAGA             95 21  0 52.38 62.092  6.00  1.00  3.092
 160 ACTGCCACCGTGTACACACTG             40 21  0 57.14 62.102  8.00  8.00  3.102
 161 ACACACTTCCTACAGCTGTGG            193 21  0 52.38 57.875  8.00  8.00  3.125
 162 CCACACACTTCCTACAGCTGT            191 21  0 52.38 57.875  8.00  8.00  3.125
 163 ACCACACACTTCCTACAGCTG            190 21  0 52.38 57.875  6.00  6.00  3.125
 164 CACTGAACACGTTCCCTAAGA             56 21  0 47.62 57.871  6.00  1.00  3.129
 165 TGCACTGAAGACCTTCCCTAATA          143 23  0 43.48 60.136  6.00  3.00  3.136
 166 GTGGGTGTACCACACACTTCCT           182 22  0 54.55 61.149  8.00  0.00  3.149
 167 GCCAAGGGAATGCACTGAA              133 19  0 52.63 62.155  4.00  3.00  3.155
 168 TAAGACATTGTGGGTGTACCACA          173 23  0 43.48 60.196  6.00  6.00  3.196
 169 ACCACACACTTCCTACAGCTGT           190 22  0 50.00 58.799  8.00  8.00  3.201
 170 ACACTGAACACGTTCCCTAAGA            55 22  0 45.45 58.769  6.00  1.00  3.231
 171 ATTGTGGGTGTACCACACACT            179 21  0 47.62 57.686  8.00  8.00  3.314
 172 GAAGACTTCCCTACATACAGCCA           13 23  0 47.83 59.666  7.00  0.00  3.334
 173 TGAAGACTTCCCTACATACAGCC           12 23  0 47.83 59.666  7.00  1.00  3.334
 174 TCACACTCACTGAAGACTTCCCT            2 23  0 47.83 60.337  7.00  5.00  3.337
 175 ATCACACTCACTGAAGACTTCCC            1 23  0 47.83 59.659  7.00  3.00  3.341
 176 AAGGGAATGCACTGAAGACCTT           136 22  0 45.45 61.364  4.00  4.00  3.364
 177 ATGCACTGAAGACCTTCCC              142 19  0 52.63 57.630  6.00  2.00  3.370
 178 GGGAATGCACTGAAGACCT              138 19  0 52.63 57.630  4.00  2.00  3.370
 179 AGGGAATGCACTGAAGACC              137 19  0 52.63 57.630  4.00  3.00  3.370
 180 ATGCACTGAAGACCTTCCCTAAT          142 23  0 43.48 60.373  6.00  2.00  3.373
 181 CCACCGTGTACACACTGAACAC            44 22  0 54.55 61.374  8.00  7.00  3.374
 182 ACGTTCCCTAAGAAATGTGGT             64 21  0 42.86 57.593  6.00  0.00  3.407
 183 TACACACTGAACACGTTCCCTAA           52 23  0 43.48 59.588  6.00  1.00  3.412
 184 AAGAAATGTGGTTAAGGGACCG            73 22  0 45.45 61.417  4.00  2.00  3.417
 185 GCTAAGACATTGTGGGTGTACC           171 22  0 50.00 58.513  5.00  4.00  3.487
 186 ACACACTGAACACGTTCCCTAAG           53 23  0 47.83 60.494  6.00  3.00  3.494
 187 CTAAGAATGTGGCCAAGGGAAT           122 22  0 45.45 61.515  6.00  2.00  3.515
 188 CACCGTGTACACACTGAACAC             45 21  0 52.38 57.482  8.00  7.00  3.518
 189 AGCTAAGACATTGTGGGTGTACC          170 23  0 47.83 59.468  5.00  4.00  3.532
 190 CGCCATGTACACACTGAAGATG            96 22  0 50.00 61.537  6.00  5.00  3.537
 191 CACTGAACACGTTCCCTAAGAAA           56 23  0 43.48 60.550  6.00  4.00  3.550
 192 CACACTCACTGAAGACTTCCCT             3 22  0 50.00 58.449  7.00  5.00  3.551
 193 CCCTAAGAAATGTGGTTAAGGGA           69 23  0 43.48 60.553  4.00  0.00  3.553
 194 TCCCTAAGAAATGTGGTTAAGGG           68 23  0 43.48 60.553  4.00  4.00  3.553
 195 CAAAAGCTAAGACATTGTGGG            166 21  0 42.86 57.442  5.00  0.00  3.558
 196 CGTTCCCTAAGAAATGTGGTTAA           65 23  0 39.13 59.441  4.00  4.00  3.559
 197 GGGTGTACCACACACTTCCTA            184 21  0 52.38 57.434  5.00  2.00  3.566
 198 AAGGGACTGCCACCGTGT                35 18  0 61.11 61.583  5.00  3.00  3.583
 199 GCACTGAAGACCTTCCCTAAT            144 21  0 47.62 57.417  6.00  2.00  3.583
 200 GACATTGTGGGTGTACCACACA           176 22  0 50.00 61.591  8.00  8.00  3.591
 201 CCCTACATACAGCCAAGGGACT            21 22  0 54.55 61.601  4.00  1.00  3.601
 202 AAATGTGGTTAAGGGACCG               76 19  0 47.37 57.380  4.00  2.00  3.620
 203 CTACATACAGCCAAGGGACTG             23 21  0 52.38 57.376  3.00  3.00  3.624
 204 GGGTGTACCACACACTTCCTAC           184 22  0 54.55 58.361  5.00  0.00  3.639
 205 ATGTGGTTAAGGGACCGC                78 18  0 55.56 58.359  4.00  2.00  3.641
 206 CACACTCACTGAAGACTTCCC              3 21  0 52.38 57.332  7.00  3.00  3.668
 207 TGTTCCCTAAGAATGTGGCCA            116 21  0 47.62 62.668  6.00  6.00  3.668
 208 AAAAGCTAAGACATTGTGGGTG           167 22  0 40.91 58.314  5.00  0.00  3.686
 209 CAAAAGCTAAGACATTGTGGGT           166 22  0 40.91 58.314  5.00  0.00  3.686
 210 ACAAAAGCTAAGACATTGTGGG           165 22  0 40.91 58.314  5.00  1.00  3.686
 211 AAGATGTTCCCTAAGAATGTGGC          112 23  0 43.48 60.707  4.00  2.00  3.707
 212 CAAGGGAATGCACTGAAGA              135 19  0 47.37 57.292  4.00  0.00  3.708
 213 GAATGTGGCCAAGGGAAT               126 18  0 50.00 58.274  6.00  2.00  3.726
 214 AAGACATTGTGGGTGTACCAC            174 21  0 47.62 57.273  5.00  4.00  3.727
 215 CGCCGCCATGTACACACT                93 18  0 61.11 61.777  6.00  2.00  3.777
 216 TACAGCCAAGGGACTGCCA               28 19  0 57.89 62.781  3.00  3.00  3.781
 217 TACACACTGAACACGTTCCCTA            52 22  0 45.45 58.214  6.00  2.00  3.786
 218 AAGACTTCCCTACATACAGCCAA           14 23  0 43.48 59.209  5.00  2.00  3.791
 219 GAAATGTGGTTAAGGGACCGC             75 21  0 52.38 62.792  4.00  2.00  3.792
 220 CGTTCCCTAAGAAATGTGGTTA            65 22  0 40.91 58.188  4.00  2.00  3.812
 221 AACACGTTCCCTAAGAAATGTG            61 22  0 40.91 58.175  8.00  8.00  3.825
 222 CGTGTACACACTGAACACGTT             48 21  0 47.62 57.155  9.00  4.00  3.845
 223 ATGTGGCCAAGGGAATGC               128 18  0 55.56 61.869  6.00  2.00  3.869
 224 TGTGGTTAAGGGACCGCC                79 18  0 61.11 61.876  4.00  1.00  3.876
 225 AAAAGCTAAGACATTGTGGGTGT          167 23  0 39.13 59.123  5.00  0.00  3.877
 226 ACAAAAGCTAAGACATTGTGGGT          165 23  0 39.13 59.123  5.00  1.00  3.877
 227 GAATGCACTGAAGACCTTCCCT           140 22  0 50.00 61.884  6.00  4.00  3.884
 228 AGGGAATGCACTGAAGACCTTC           137 22  0 50.00 61.884  6.00  6.00  3.884
 229 TCCTACAGCTGTGGGAAAATGA           200 22  0 45.45 61.895 12.00  1.00  3.895
 230 CCAAGGGAATGCACTGAA               134 18  0 50.00 58.095  4.00  3.00  3.905
 231 CAAAAGCTAAGACATTGTGGGTG          166 23  0 43.48 60.909  5.00  0.00  3.909
 232 GTACACACTGAACACGTTCCCTA           51 23  0 47.83 59.087  6.00  2.00  3.913
 233 GATGTTCCCTAAGAATGTGGCC           114 22  0 50.00 61.913  4.00  4.00  3.913
 234 AAGAAATGTGGTTAAGGGACC             73 21  0 42.86 57.076  4.00  3.00  3.924
 235 TGTACACACTGAACACGTTCC             50 21  0 47.62 57.053  6.00  4.00  3.947
 236 GTGTACACACTGAACACGTTCC            49 22  0 50.00 58.046  8.00  4.00  3.954
 237 GCCATGTACACACTGAAGATGTT           97 23  0 43.48 59.046  6.00  4.00  3.954
 238 GTTCCCTAAGAATGTGGCC              117 19  0 52.63 57.030  4.00  4.00  3.970
 239 CTGAACACGTTCCCTAAGAAA             58 21  0 42.86 57.025  6.00  4.00  3.975
 240 GAAGATGTTCCCTAAGAATGTGG          111 23  0 43.48 59.014  6.00  0.00  3.986
 241 ATTGTGGGTGTACCACACACTTC          179 23  0 47.83 60.995  8.00  5.00  3.995
 242 AAGACATTGTGGGTGTACCACAC          174 23  0 47.83 60.995  6.00  6.00  3.995
 243 CACACTGAACACGTTCCCT               54 19  0 52.63 57.003  6.00  4.00  3.997
 244 ACACGTTCCCTAAGAAATGTGGT           62 23  0 43.48 61.013  8.00  6.00  4.013
 245 ACGTTCCCTAAGAAATGTGGTTA           64 23  0 39.13 58.983  6.00  2.00  4.017
 246 TACCACACACTTCCTACAGCTG           189 22  0 50.00 57.969  6.00  6.00  4.031
 247 GTTCCCTAAGAATGTGGCCAAG           117 22  0 50.00 62.039  6.00  2.00  4.039
 248 CTTCCCTACATACAGCCAAGGG            18 22  0 54.55 62.041  4.00  1.00  4.041
 249 AAGACTTCCCTACATACAGCCAAG          14 24  0 45.83 60.076  5.00  4.00  4.076
 250 ACTGAACACGTTCCCTAAGAAA            57 22  0 40.91 57.921  6.00  4.00  4.079
 251 GACTTCCCTACATACAGCCAAG            16 22  0 50.00 57.893  4.00  4.00  4.107
 252 AGACTTCCCTACATACAGCCAA            15 22  0 45.45 57.893  3.00  2.00  4.107
 253 AAGACTTCCCTACATACAGCCA            14 22  0 45.45 57.893  5.00  0.00  4.107
 254 GGGTGTACCACACACTTCCTACA          184 23  0 52.17 61.117  5.00  1.00  4.117
 255 TGGGTGTACCACACACTTCCTAC          183 23  0 52.17 61.117  6.00  0.00  4.117
 256 GTGGGTGTACCACACACTTCCTA          182 23  0 52.17 61.117  8.00  2.00  4.117
 257 AGACTTCCCTACATACAGCCAAG           15 23  0 47.83 58.861  4.00  4.00  4.139
 258 TACCACACACTTCCTACAGCTGT          189 23  0 47.83 58.855  8.00  8.00  4.145
 259 GTACCACACACTTCCTACAGCTG          188 23  0 52.17 58.855  6.00  6.00  4.145
 260 TGTACCACACACTTCCTACAGCT          187 23  0 47.83 58.855  4.00  4.00  4.145
 261 ACGTTCCCTAAGAAATGTGGTTAA          64 24  0 37.50 60.157  6.00  4.00  4.157
 262 TGTACCACACACTTCCTACAGC           187 22  0 50.00 57.817  4.00  2.00  4.183
 263 GCTAAGACATTGTGGGTGTACCA          171 23  0 47.83 61.184  5.00  2.00  4.184
 264 AAAGCTAAGACATTGTGGGTGT           168 22  0 40.91 57.804  5.00  0.00  4.196
 265 GACTTCCCTACATACAGCCAAGG           16 23  0 52.17 61.222  4.00  2.00  4.222
 266 CCGCCATGTACACACTGAAGAT            95 22  0 50.00 62.253  6.00  3.00  4.253
 267 CGTTCCCTAAGAAATGTGGTTAAG          65 24  0 41.67 60.268  4.00  2.00  4.268
 268 GTGTACCACACACTTCCTACAGC          186 23  0 52.17 58.715  5.00  2.00  4.285
 269 CACACTGAAGATGTTCCCTAAGAA         105 24  0 41.67 59.690  6.00  2.00  4.310
 270 GTACCACACACTTCCTACAGCTGT         188 24  0 50.00 59.674  8.00  8.00  4.326
 271 GTGTACCACACACTTCCTACAGCT         186 24  0 50.00 59.674  5.00  4.00  4.326
 272 AATGCACTGAAGACCTTCCCTAA          141 23  0 43.48 61.327  6.00  1.00  4.327
 273 TCACACTCACTGAAGACTTCCCTA           2 24  0 45.83 60.335  7.00  3.00  4.335
 274 TGTACACACTGAAGATGTTCCCT          101 23  0 43.48 58.639  6.00  4.00  4.361
 275 GTACACACTGAACACGTTCCCTAA          51 24  0 45.83 60.367  6.00  1.00  4.367
 276 ATGCACTGAAGACCTTCCCTAATA         142 24  0 41.67 60.370  6.00  3.00  4.370
 277 TAAGAAATGTGGTTAAGGGACCG           72 23  0 43.48 61.379  4.00  2.00  4.379
 278 TTCCCTAAGAAATGTGGTTAAGG           67 23  0 39.13 58.619  4.00  4.00  4.381
 279 GCCATGTACACACTGAAGATGT            97 22  0 45.45 57.601  6.00  4.00  4.399
 280 TGTACACACTGAAGATGTTCCC           101 22  0 45.45 57.574  6.00  2.00  4.426
 281 GGTGTACCACACACTTCCTACAG          185 23  0 52.17 58.570  5.00  2.00  4.430
 282 ACATTGTGGGTGTACCACACACT          177 23  0 47.83 61.441  8.00  8.00  4.441
 283 CACGTTCCCTAAGAAATGTGGTT           63 23  0 43.48 61.453  8.00  0.00  4.453
 284 AACACGTTCCCTAAGAAATGTGG           61 23  0 43.48 61.453  8.00  6.00  4.453
 285 CTAAGACATTGTGGGTGTACCA           172 22  0 45.45 57.543  5.00  2.00  4.457
 286 TCCCTACATACAGCCAAGGGAC            20 22  0 54.55 62.458  5.00  1.00  4.458
 287 CACACTGAACACGTTCCCTAAGA           54 23  0 47.83 61.476  6.00  1.00  4.476
 288 GCACTGAAGACCTTCCCTAATA           144 22  0 45.45 57.520  6.00  3.00  4.480
 289 CACTCACTGAAGACTTCCCTACA            5 23  0 47.83 58.519  7.00  1.00  4.481
 290 CACACTCACTGAAGACTTCCCTA            3 23  0 47.83 58.519  7.00  3.00  4.481
 291 TCCCTAAGAATGTGGCCA               119 18  0 50.00 57.514  6.00  6.00  4.486
 292 TACACACTGAACACGTTCCCTAAG          52 24  0 45.83 60.486  6.00  3.00  4.486
 293 GGTGTACCACACACTTCCTACA           185 22  0 50.00 57.509  5.00  1.00  4.491
 294 GTGGCCAAGGGAATGCAC               130 18  0 61.11 62.504  6.00  2.00  4.504
 295 CTGAAGACTTCCCTACATACAGCC          11 24  0 50.00 60.533  7.00  1.00  4.533
 296 CTGAACACGTTCCCTAAGAAAT            58 22  0 40.91 57.436  6.00  2.00  4.564
 297 CTAAGACATTGTGGGTGTACCAC          172 23  0 47.83 58.433  5.00  4.00  4.567
 298 CCGCCATGTACACACTGA                95 18  0 55.56 57.425  6.00  2.00  4.575
 299 ACCACACACTTCCTACAGCTGTG          190 23  0 52.17 61.580  8.00  6.00  4.580
 300 ATCACACTCACTGAAGACTTCCCT           1 24  0 45.83 60.582  7.00  5.00  4.582
 301 TAAGACATTGTGGGTGTACCAC           173 22  0 45.45 57.389  5.00  4.00  4.611
 302 TACAAAAGCTAAGACATTGTGGG          164 23  0 39.13 58.382  5.00  1.00  4.618
 303 CACACTGAAGATGTTCCCTAAGA          105 23  0 43.48 58.382  6.00  1.00  4.618
 304 GGCCAAGGGAATGCACTG               132 18  0 61.11 62.628  4.00  3.00  4.628
 305 TGGCCAAGGGAATGCACT               131 18  0 55.56 62.628  6.00  4.00  4.628
 306 ACACTCACTGAAGACTTCCCTACA           4 24  0 45.83 59.361  7.00  1.00  4.639
 307 CACACTCACTGAAGACTTCCCTAC           3 24  0 50.00 59.361  7.00  1.00  4.639
 308 GTTCCCTAAGAAATGTGGTTAAGG          66 24  0 41.67 59.360  4.00  4.00  4.640
 309 GCACTGAAGACCTTCCCTAATAC          144 23  0 47.83 58.358  6.00  1.00  4.642
 310 CCTAAGAAATGTGGTTAAGGGA            70 22  0 40.91 57.332  4.00  0.00  4.668
 311 TCCCTAAGAAATGTGGTTAAGG            68 22  0 40.91 57.332  4.00  4.00  4.668
 312 AAGCTAAGACATTGTGGGTGTACC         169 24  0 45.83 60.673  5.00  4.00  4.673
 313 GTGTACACACTGAACACGTTCCC           49 23  0 52.17 61.688  8.00  2.00  4.688
 314 CATGTACACACTGAAGATGTTCCC          99 24  0 45.83 60.700  6.00  2.00  4.700
 315 CCATGTACACACTGAAGATGTTCC          98 24  0 45.83 60.700  6.00  4.00  4.700
 316 GAACACGTTCCCTAAGAAATGT            60 22  0 40.91 57.291  6.00  6.00  4.709
 317 ACTGAACACGTTCCCTAAGAAAT           57 23  0 39.13 58.273  6.00  2.00  4.727
 318 TGAAGATGTTCCCTAAGAATGTG          110 23  0 39.13 58.251  6.00  0.00  4.749
 319 AAGATGTTCCCTAAGAATGTGG           112 22  0 40.91 57.236  4.00  0.00  4.764
 320 CTGAACACGTTCCCTAAGAAATGT          58 24  0 41.67 60.767  6.00  6.00  4.767
 321 ACTGAACACGTTCCCTAAGAAATG          57 24  0 41.67 60.767  6.00  4.00  4.767
 322 CACTGAACACGTTCCCTAAGAAAT          56 24  0 41.67 60.767  6.00  2.00  4.767
 323 CAAGGGAATGCACTGAAGACCT           135 22  0 50.00 62.786  4.00  2.00  4.786
 324 ACACACTGAAGATGTTCCCTAAGA         104 24  0 41.67 59.212  6.00  1.00  4.788
 325 ATGTTCCCTAAGAATGTGGCCA           115 22  0 45.45 62.791  6.00  6.00  4.791
 326 CTGAAGATGTTCCCTAAGAATGTG         109 24  0 41.67 59.196  6.00  0.00  4.804
 327 CACTGAAGATGTTCCCTAAGAATG         107 24  0 41.67 59.196  6.00  0.00  4.804
 328 TGTACACACTGAACACGTTCCCT           50 23  0 47.83 61.806  6.00  4.00  4.806
 329 TAAGAAATGTGGTTAAGGGACC            72 22  0 40.91 57.190  4.00  3.00  4.810
 330 GAATGCACTGAAGACCTTCCCTA          140 23  0 47.83 61.828  6.00  2.00  4.828
 331 AAAAGCTAAGACATTGTGGGTGTA         167 24  0 37.50 59.162  5.00  2.00  4.838
 332 TACAAAAGCTAAGACATTGTGGGT         164 24  0 37.50 59.162  5.00  0.00  4.838
 333 GCCATGTACACACTGAAGATGTTC          97 24  0 45.83 60.841  6.00  6.00  4.841
 334 TGCACTGAAGACCTTCCC               143 18  0 55.56 57.158  6.00  2.00  4.842
 335 CTAAGAAATGTGGTTAAGGGACC           71 23  0 43.48 58.153  4.00  3.00  4.847
 336 CCTAAGAAATGTGGTTAAGGGAC           70 23  0 43.48 58.153  4.00  1.00  4.847
 337 GCACTGAAGACCTTCCCTAATACA         144 24  0 45.83 60.860  6.00  0.00  4.860
 338 TGCACTGAAGACCTTCCCTAATAC         143 24  0 45.83 60.860  6.00  1.00  4.860
 339 GAAGACTTCCCTACATACAGCCAA          13 24  0 45.83 60.860  7.00  2.00  4.860
 340 CAGCCAAGGGACTGCCAC                30 18  0 66.67 62.888  3.00  1.00  4.888
 341 ACAGCCAAGGGACTGCCA                29 18  0 61.11 62.888  3.00  3.00  4.888
 342 TTCCCTACATACAGCCAAGGGA            19 22  0 50.00 62.893  5.00  0.00  4.893
 343 CATTGTGGGTGTACCACACACTT          178 23  0 47.83 61.918  8.00  6.00  4.918
 344 TAAGACATTGTGGGTGTACCACAC         173 24  0 45.83 60.969  6.00  6.00  4.969
 345 TAAGAATGTGGCCAAGGG               123 18  0 50.00 57.022  6.00  1.00  4.978
 346 CCCTAAGAATGTGGCCAA               120 18  0 50.00 57.022  6.00  4.00  4.978
 347 ATGTACACACTGAAGATGTTCCC          100 23  0 43.48 57.976  6.00  2.00  5.024
 348 TGTACACACTGAAGATGTTCCCTAA        101 25  0 40.00 59.975  6.00  1.00  5.025
 349 ATGTACACACTGAAGATGTTCCCT         100 24  0 41.67 58.969  6.00  4.00  5.031
 350 ACTGAAGATGTTCCCTAAGAATGTG        108 25  0 40.00 59.947  6.00  0.00  5.053
 351 CACTGAAGATGTTCCCTAAGAATGT        107 25  0 40.00 59.947  6.00  1.00  5.053
 352 ACACTGAAGATGTTCCCTAAGAATG        106 25  0 40.00 59.947  6.00  0.00  5.053
 353 CACACTGAAGATGTTCCCTAAGAAT        105 25  0 40.00 59.947  6.00  2.00  5.053
 354 CACTGAAGACCTTCCCTAATACAAA        145 25  0 40.00 59.946  6.00  0.00  5.054
 355 CTAAGACATTGTGGGTGTACCACA         172 24  0 45.83 61.086  6.00  6.00  5.086
 356 GATCACACTCACTGAAGACTTCC            0 23  0 47.83 57.900  7.00  5.00  5.100
 357 AAAGCTAAGACATTGTGGGTGTA          168 23  0 39.13 57.892  5.00  2.00  5.108
 358 AAAAGCTAAGACATTGTGGGTGTAC        167 25  0 40.00 59.890  5.00  4.00  5.110
 359 AATACAAAAGCTAAGACATTGTGGG        162 25  0 36.00 59.858  5.00  0.00  5.142
 360 CACTCACTGAAGACTTCCCTACAT           5 24  0 45.83 58.854  7.00  2.00  5.146
 361 GAAGACCTTCCCTAATACAAAAGC         149 24  0 41.67 58.814  6.00  2.00  5.186
 362 CCTAAGAAATGTGGTTAAGGGACC          70 24  0 45.83 61.226  4.00  3.00  5.226
 363 CCCTAAGAAATGTGGTTAAGGGAC          69 24  0 45.83 61.226  4.00  1.00  5.226
 364 CACTGAAGACCTTCCCTAATACAA         145 24  0 41.67 58.759  6.00  0.00  5.241
 365 ACACTGAACACGTTCCCTAAGAAA          55 24  0 41.67 61.267  6.00  4.00  5.267
 366 CGCCATGTACACACTGAAGATGT           96 23  0 47.83 62.286  6.00  4.00  5.286
 367 AAAGCTAAGACATTGTGGGTGTAC         168 24  0 41.67 58.699  5.00  4.00  5.301
 368 TGTACACACTGAAGATGTTCCCTA         101 24  0 41.67 58.699  6.00  2.00  5.301
 369 ATACAAAAGCTAAGACATTGTGGG         163 24  0 37.50 58.698  5.00  0.00  5.302
 370 ACACTCACTGAAGACTTCCCTACAT          4 25  0 44.00 59.643  7.00  2.00  5.357
 371 GACATTGTGGGTGTACCACACAC          176 23  0 52.17 62.364  8.00  8.00  5.364
 372 GAAGACCTTCCCTAATACAAAAGCT        149 25  0 40.00 59.629  6.00  4.00  5.371
 373 CACGTTCCCTAAGAAATGTGGTTA          63 24  0 41.67 61.414  8.00  2.00  5.414
 374 ACACACTGAAGATGTTCCCTAAGAA        104 25  0 40.00 60.434  6.00  2.00  5.434
 375 GATCACACTCACTGAAGACTTCCC           0 24  0 50.00 61.444  7.00  3.00  5.444
 376 AGACATTGTGGGTGTACCACACA          175 23  0 47.83 62.479  8.00  8.00  5.479
 377 AATGCACTGAAGACCTTCCCTAAT         141 24  0 41.67 61.502  6.00  2.00  5.502
 378 TACCACACACTTCCTACAGCTGTG         189 24  0 50.00 61.532  8.00  6.00  5.532
 379 TGTACCACACACTTCCTACAGCTG         187 24  0 50.00 61.532  6.00  6.00  5.532
 380 TGAAGATGTTCCCTAAGAATGTGG         110 24  0 41.67 61.537  6.00  0.00  5.537
 381 CACTGAAGACCTTCCCTAATACA          145 23  0 43.48 57.450  6.00  0.00  5.550
 382 ATACAAAAGCTAAGACATTGTGGGT        163 25  0 36.00 59.432  5.00  1.00  5.568
 383 ATCACACTCACTGAAGACTTCCCTA          1 25  0 44.00 60.571  7.00  3.00  5.571
 384 CACTGAAGACTTCCCTACATACAGC          9 25  0 48.00 60.575  7.00  2.00  5.575
 385 ACACACTGAAGATGTTCCCTAAG          104 23  0 43.48 57.416  6.00  3.00  5.584
 386 CAAAAGCTAAGACATTGTGGGTGT         166 24  0 41.67 61.607  5.00  0.00  5.607
 387 ACAAAAGCTAAGACATTGTGGGTG         165 24  0 41.67 61.607  5.00  0.00  5.607
 388 AAGCTAAGACATTGTGGGTGTAC          169 23  0 43.48 57.388  5.00  4.00  5.612
 389 TTCCCTAAGAAATGTGGTTAAGGG          67 24  0 41.67 61.631  4.00  4.00  5.631
 390 AGATGTTCCCTAAGAATGTGGCC          113 23  0 47.83 62.687  4.00  4.00  5.687
 391 TCCTACAGCTGTGGGAAAATGAG          200 23  0 47.83 62.697 12.00  0.00  5.697
 392 ACTTCCCTACATACAGCCAAGGG           17 23  0 52.17 62.697  4.00  1.00  5.697
 393 CGTGTACACACTGAACACGTTCC           48 23  0 52.17 62.750  9.00  4.00  5.750
 394 CCGTGTACACACTGAACACGTTC           47 23  0 52.17 62.750  9.00  6.00  5.750
 395 TGTACACACTGAACACGTTCCCTA          50 24  0 45.83 61.750  6.00  2.00  5.750
 396 TACACACTGAAGATGTTCCCTAAGA        103 25  0 40.00 59.249  6.00  1.00  5.751
 397 CTGAAGACTTCCCTACATACAGC           11 23  0 47.83 57.219  7.00  2.00  5.781
 398 AAGACCTTCCCTAATACAAAAGC          150 23  0 39.13 57.202  4.00  2.00  5.798
 399 GTGGGTGTACCACACACTTCCTAC         182 24  0 54.17 61.849  8.00  0.00  5.849
 400 CATGTACACACTGAAGATGTTCC           99 23  0 43.48 57.126  6.00  4.00  5.874
 401 CCATGTACACACTGAAGATGTTC           98 23  0 43.48 57.126  6.00  6.00  5.874
 402 GCTAAGACATTGTGGGTGTACCAC         171 24  0 50.00 61.894  5.00  4.00  5.894
 403 AAGACCTTCCCTAATACAAAAGCT         150 24  0 37.50 58.100  4.00  4.00  5.900
 404 ACGTTCCCTAAGAAATGTGGTTAAG         64 25  0 40.00 60.930  6.00  2.00  5.930
 405 ACTGAAGACTTCCCTACATACAGC          10 24  0 45.83 58.067  7.00  2.00  5.933
 406 TGAAGACCTTCCCTAATACAAAA          148 23  0 34.78 57.052  6.00  1.00  5.948
 407 ACACTGAAGATGTTCCCTAAGAA          106 23  0 39.13 57.049  6.00  2.00  5.951
 408 GGGTGTACCACACACTTCCTACAG         184 24  0 54.17 61.959  5.00  2.00  5.959
 409 ATGTACACACTGAAGATGTTCCCTA        100 25  0 40.00 59.015  6.00  2.00  5.985
 410 AGCTAAGACATTGTGGGTGTACCA         170 24  0 45.83 62.000  5.00  2.00  6.000
 411 AGACTTCCCTACATACAGCCAAGG          15 24  0 50.00 62.003  4.00  2.00  6.003
 412 TGAAGACCTTCCCTAATACAAAAG         148 24  0 37.50 57.993  6.00  3.00  6.007
 413 CTGAAGACCTTCCCTAATACAAAA         147 24  0 37.50 57.993  6.00  1.00  6.007
 414 GTACACACTGAAGATGTTCCCTAAGA       102 26  0 42.31 59.993  6.00  1.00  6.007
 415 TCACTGAAGACTTCCCTACATACA           8 24  0 41.67 57.947  7.00  0.00  6.053
 416 ATGCACTGAAGACCTTCCCTAATAC        142 25  0 44.00 61.054  6.00  1.00  6.054
 417 TCACACTCACTGAAGACTTCCCTAC          2 25  0 48.00 61.088  7.00  1.00  6.088
 418 GGTGTACCACACACTTCCTACAGC         185 24  0 54.17 62.095  5.00  2.00  6.095
 419 TCCCTAAGAAATGTGGTTAAGGGA          68 24  0 41.67 62.095  6.00  6.00  6.095
 420 TCACTGAAGACTTCCCTACATACAG          8 25  0 44.00 58.904  7.00  2.00  6.096
 421 CTCACTGAAGACTTCCCTACATACA          7 25  0 44.00 58.904  7.00  0.00  6.096
 422 CACTCACTGAAGACTTCCCTACATA          5 25  0 44.00 58.904  7.00  2.00  6.096
 423 TGAAGACCTTCCCTAATACAAAAGC        148 25  0 40.00 61.112  6.00  2.00  6.112
 424 ACACGTTCCCTAAGAAATGTGGTT          62 24  0 41.67 62.116  8.00  4.00  6.116
 425 AACACGTTCCCTAAGAAATGTGGT          61 24  0 41.67 62.116  8.00  8.00  6.116
 426 CTAAGAAATGTGGTTAAGGGACCG          71 24  0 45.83 62.125  4.00  2.00  6.125
 427 TAATACAAAAGCTAAGACATTGTGGG       161 26  0 34.62 59.872  5.00  0.00  6.128
 428 CTGAAGACCTTCCCTAATACAAAAG        147 25  0 40.00 58.861  6.00  3.00  6.139
 429 ACACACTGAACACGTTCCCTAAGA          53 24  0 45.83 62.186  6.00  1.00  6.186
 430 TGAAGACTTCCCTACATACAGCCA          12 24  0 45.83 62.191  7.00  0.00  6.191
 431 GTACACACTGAACACGTTCCCTAAG         51 25  0 48.00 61.201  6.00  3.00  6.201
 432 ATGTACACACTGAAGATGTTCCCTAA       100 26  0 38.46 60.216  6.00  1.00  6.216
 433 CTCACTGAAGACTTCCCTACATACAG         7 26  0 46.15 59.782  7.00  2.00  6.218
 434 ACTGAAGACTTCCCTACATACAGCC         10 25  0 48.00 61.224  7.00  1.00  6.224
 435 ACTGAAGACCTTCCCTAATACAAAA        146 25  0 36.00 58.743  6.00  1.00  6.257
 436 GAAGATGTTCCCTAAGAATGTGGC         111 24  0 45.83 62.298  6.00  2.00  6.298
 437 ACTCACTGAAGACTTCCCTACATACA         6 26  0 42.31 59.665  7.00  0.00  6.335
 438 CACTCACTGAAGACTTCCCTACATAC         5 26  0 46.15 59.665  7.00  0.00  6.335
 439 ACACTCACTGAAGACTTCCCTACATA         4 26  0 42.31 59.665  7.00  2.00  6.335
 440 GACCTTCCCTAATACAAAAGCTAAGA       152 26  0 38.46 59.649  4.00  1.00  6.351
 441 GAAGACCTTCCCTAATACAAAAGCTA       149 26  0 38.46 59.649  6.00  2.00  6.351
 442 CCTTCCCTAATACAAAAGCTAAGACA       154 26  0 38.46 60.381  4.00  0.00  6.381
 443 TGAACACGTTCCCTAAGAAATGTG          59 24  0 41.67 62.411  8.00  8.00  6.411
 444 TACACACTGAAGATGTTCCCTAAGAA       103 26  0 38.46 60.429  6.00  2.00  6.429
 445 ACTGAACACGTTCCCTAAGAAATGT         57 25  0 40.00 61.443  6.00  6.00  6.443
 446 ACACTGAACACGTTCCCTAAGAAAT         55 25  0 40.00 61.443  6.00  2.00  6.443
 447 ACTGAAGACCTTCCCTAATACAAAAG       146 26  0 38.46 59.556  6.00  3.00  6.444
 448 AATGCACTGAAGACCTTCCCTAATA        141 25  0 40.00 61.463  6.00  3.00  6.463
 449 ACTGAAGACCTTCCCTAATACAAA         146 24  0 37.50 57.525  6.00  0.00  6.475
 450 TACACACTGAAGATGTTCCCTAAG         103 24  0 41.67 57.518  6.00  3.00  6.482
 451 GTGTACACACTGAACACGTTCCCT          49 24  0 50.00 62.520  8.00  4.00  6.520
 452 AATACAAAAGCTAAGACATTGTGGGT       162 26  0 34.62 60.524  5.00  1.00  6.524
 453 CATGTACACACTGAAGATGTTCCCT         99 25  0 44.00 61.548  6.00  4.00  6.548
 454 ACTGAAGATGTTCCCTAAGAATGT         108 24  0 37.50 57.442  6.00  1.00  6.558
 455 ACACTGAAGATGTTCCCTAAGAAT         106 24  0 37.50 57.442  6.00  2.00  6.558
 456 CAAAAGCTAAGACATTGTGGGTGTA        166 25  0 40.00 61.563  5.00  2.00  6.563
 457 TACAAAAGCTAAGACATTGTGGGTG        164 25  0 40.00 61.563  5.00  0.00  6.563
 458 CCTTCCCTAATACAAAAGCTAAGA         154 24  0 37.50 57.427  4.00  1.00  6.573
 459 CACACTGAACACGTTCCCTAAGAA          54 24  0 45.83 62.619  6.00  2.00  6.619
 460 GTACACACTGAAGATGTTCCCTAA         102 24  0 41.67 57.376  6.00  1.00  6.624
 461 ACATTGTGGGTGTACCACACACTT         177 24  0 45.83 62.627  8.00  6.00  6.627
 462 GAAGACTTCCCTACATACAGCCAAG         13 25  0 48.00 61.639  7.00  4.00  6.639
 463 ACACTGAAGATGTTCCCTAAGAATGT       106 26  0 38.46 60.648  6.00  1.00  6.648
 464 ACACACTGAAGATGTTCCCTAAGAAT       104 26  0 38.46 60.648  6.00  2.00  6.648
 465 GTACACACTGAAGATGTTCCCTAAG        102 25  0 44.00 58.344  6.00  3.00  6.656
 466 GACCTTCCCTAATACAAAAGCTAA         152 24  0 37.50 57.300  4.00  1.00  6.700
 467 AAGACCTTCCCTAATACAAAAGCTAA       150 26  0 34.62 59.264  4.00  1.00  6.736
 468 AAAGCTAAGACATTGTGGGTGTACC        168 25  0 44.00 61.766  5.00  4.00  6.766
 469 TCCCTAATACAAAAGCTAAGACATTG       157 26  0 34.62 59.231  4.00  3.00  6.769
 470 CTAAGACATTGTGGGTGTACCACAC        172 25  0 48.00 61.794  6.00  6.00  6.794
 471 TGTACACACTGAAGATGTTCCCTAAG       101 26  0 42.31 60.811  6.00  3.00  6.811
 472 CCTTCCCTAATACAAAAGCTAAGAC        154 25  0 40.00 58.167  4.00  1.00  6.833
 473 ACCTTCCCTAATACAAAAGCTAAGA        153 25  0 36.00 58.167  4.00  1.00  6.833
 474 GACCTTCCCTAATACAAAAGCTAAG        152 25  0 40.00 58.167  4.00  3.00  6.833
 475 AGACCTTCCCTAATACAAAAGCTAA        151 25  0 36.00 58.167  4.00  1.00  6.833
 476 AAGACCTTCCCTAATACAAAAGCTA        150 25  0 36.00 58.167  4.00  2.00  6.833
 477 CACTGAAGACTTCCCTACATACAG           9 24  0 45.83 57.151  7.00  2.00  6.849
 478 GAATGCACTGAAGACCTTCCCTAA         140 24  0 45.83 62.910  6.00  1.00  6.910
 479 GCACTGAAGACCTTCCCTAATACAA        144 25  0 44.00 61.945  6.00  0.00  6.945
 480 CACACTCACTGAAGACTTCCCTACA          3 25  0 48.00 61.960  7.00  1.00  6.960
 481 AAGACCTTCCCTAATACAAAAGCTAAG      150 27  0 37.04 59.997  4.00  3.00  7.003
 482 CACTGAAGACCTTCCCTAATACAAAA       145 26  0 38.46 61.028  6.00  1.00  7.028
 483 AGACCTTCCCTAATACAAAAGCTAAG       151 26  0 38.46 58.971  4.00  3.00  7.029
 484 ACACGTTCCCTAAGAAATGTGGTTA         62 25  0 40.00 62.057  8.00  2.00  7.057
 485 TCCCTAATACAAAAGCTAAGACATTGT      157 27  0 33.33 59.891  5.00  5.00  7.109
 486 TACACACTGAACACGTTCCCTAAGA         52 25  0 44.00 62.121  6.00  1.00  7.121
 487 ACCTTCCCTAATACAAAAGCTAAGAC       153 26  0 38.46 58.862  4.00  1.00  7.138
 488 GTACCACACACTTCCTACAGCTGTG        188 25  0 52.00 62.233  8.00  6.00  7.233
 489 TGTACCACACACTTCCTACAGCTGT        187 25  0 48.00 62.233  8.00  8.00  7.233
 490 GTGTACCACACACTTCCTACAGCTG        186 25  0 52.00 62.233  6.00  6.00  7.233
 491 GTTCCCTAAGAAATGTGGTTAAGGG         66 25  0 44.00 62.241  4.00  4.00  7.241
 492 ACAAAAGCTAAGACATTGTGGGTGT        165 25  0 40.00 62.258  5.00  0.00  7.258
 493 GATCACACTCACTGAAGACTTCCCT          0 25  0 48.00 62.262  7.00  5.00  7.262
 494 TTCCCTAATACAAAAGCTAAGACATTG      156 27  0 33.33 60.279  4.00  3.00  7.279
 495 ATCACACTCACTGAAGACTTCCCTAC         1 26  0 46.15 61.279  7.00  1.00  7.279
 496 GACCTTCCCTAATACAAAAGCTAAGAC      152 27  0 40.74 60.283  4.00  1.00  7.283
 497 CCCTAATACAAAAGCTAAGACATTG        158 25  0 36.00 57.709  4.00  3.00  7.291
 498 CTGAAGATGTTCCCTAAGAATGTGG        109 25  0 44.00 62.299  6.00  0.00  7.299
 499 ACACTCACTGAAGACTTCCCTACATAC        4 27  0 44.44 60.375  7.00  0.00  7.375
 500 AGACCTTCCCTAATACAAAAGCTAAGA      151 27  0 37.04 60.382  4.00  1.00  7.382
 501 GTGTACACACTGAACACGTTCCCTA         49 25  0 48.00 62.441  8.00  2.00  7.441
 502 CACGTTCCCTAAGAAATGTGGTTAA         63 25  0 40.00 62.450  8.00  4.00  7.450
 503 ACTCACTGAAGACTTCCCTACATACAG        6 27  0 44.44 60.484  7.00  2.00  7.484
 504 CATGTACACACTGAAGATGTTCCCTA        99 26  0 42.31 61.505  6.00  2.00  7.505
 505 CTTCCCTAATACAAAAGCTAAGACA        155 25  0 36.00 57.487  4.00  0.00  7.513
 506 TAATACAAAAGCTAAGACATTGTGGGT      161 27  0 33.33 60.517  5.00  1.00  7.517
 507 CCCTAATACAAAAGCTAAGACATTGT       158 26  0 34.62 58.432  5.00  5.00  7.568
 508 CCTTCCCTAATACAAAAGCTAAGACAT      154 27  0 37.04 60.573  4.00  2.00  7.573
 509 CTAATACAAAAGCTAAGACATTGTGGG      160 27  0 37.04 60.616  5.00  0.00  7.616
 510 CCTAATACAAAAGCTAAGACATTGTGG      159 27  0 37.04 60.616  5.00  1.00  7.616
 511 CCCTAATACAAAAGCTAAGACATTGTG      158 27  0 37.04 60.616  5.00  3.00  7.616
 512 TACACACTGAAGATGTTCCCTAAGAAT      103 27  0 37.04 60.637  6.00  2.00  7.637
 513 AGCTAAGACATTGTGGGTGTACCAC        170 25  0 48.00 62.652  5.00  4.00  7.652
 514 GAAGACCTTCCCTAATACAAAAGCTAA      149 27  0 37.04 60.660  6.00  1.00  7.660
 515 TCCCTAAGAAATGTGGTTAAGGGAC         68 25  0 44.00 62.702  6.00  2.00  7.702
 516 ATACAAAAGCTAAGACATTGTGGGTG       163 26  0 38.46 61.720  5.00  0.00  7.720
 517 TGAAGACCTTCCCTAATACAAAAGCT       148 26  0 38.46 61.823  6.00  4.00  7.823
 518 CTGAAGACCTTCCCTAATACAAAAGC       147 26  0 42.31 61.823  6.00  2.00  7.823
 519 ACTCACTGAAGACTTCCCTACATAC          6 25  0 44.00 57.154  7.00  0.00  7.846
 520 GGTGTACCACACACTTCCTACAGCT        185 25  0 52.00 62.860  5.00  4.00  7.860
 521 TGTACACACTGAACACGTTCCCTAA         50 25  0 44.00 62.865  6.00  1.00  7.865
 522 TTCCCTAATACAAAAGCTAAGACATT       156 26  0 30.77 58.116  4.00  3.00  7.884
 523 CTGAAGACTTCCCTACATACAGCCA         11 25  0 48.00 62.932  7.00  0.00  7.932
 524 ACCTTCCCTAATACAAAAGCTAAGACA      153 27  0 37.04 60.995  4.00  0.00  7.995
 525 ATGTACACACTGAAGATGTTCCCTAAG      100 27  0 40.74 61.005  6.00  3.00  8.005
 526 AATGCACTGAAGACCTTCCCTAATAC       141 26  0 42.31 62.082  6.00  1.00  8.082
 527 CTTCCCTAATACAAAAGCTAAGACATT      155 27  0 33.33 58.903  4.00  3.00  8.097
 528 GTACACACTGAAGATGTTCCCTAAGAA      102 27  0 40.74 61.101  6.00  2.00  8.101
 529 CACACTCACTGAAGACTTCCCTACAT         3 26  0 46.15 62.110  7.00  2.00  8.110
 530 TCACTGAAGACTTCCCTACATACAGC         8 26  0 46.15 62.145  7.00  2.00  8.145
 531 CTTCCCTAATACAAAAGCTAAGACAT       155 26  0 34.62 57.811  4.00  2.00  8.189
 532 GATCACACTCACTGAAGACTTCCCTA         0 26  0 46.15 62.195  7.00  3.00  8.195
 533 CAAAAGCTAAGACATTGTGGGTGTAC       166 26  0 42.31 62.195  5.00  4.00  8.195
 534 ACAAAAGCTAAGACATTGTGGGTGTA       165 26  0 38.46 62.195  5.00  2.00  8.195
 535 TACAAAAGCTAAGACATTGTGGGTGT       164 26  0 38.46 62.195  5.00  0.00  8.195
 536 CTAATACAAAAGCTAAGACATTGTGG       160 26  0 34.62 57.766  5.00  1.00  8.234
 537 CCTAATACAAAAGCTAAGACATTGTG       159 26  0 34.62 57.766  5.00  3.00  8.234
 538 CACTGAAGATGTTCCCTAAGAATGTG       107 26  0 42.31 62.327  6.00  0.00  8.327
 539 CACACTGAAGATGTTCCCTAAGAATG       105 26  0 42.31 62.327  6.00  0.00  8.327
 540 CACTGAAGACCTTCCCTAATACAAAAG      145 27  0 40.74 61.738  6.00  3.00  8.738
 541 GTACACACTGAACACGTTCCCTAAGA        51 26  0 46.15 62.761  6.00  1.00  8.761
 542 AAAAGCTAAGACATTGTGGGTGTACC       167 26  0 42.31 62.764  5.00  4.00  8.764
 543 TGAAGACCTTCCCTAATACAAAAGCTA      148 27  0 37.04 61.778  6.00  2.00  8.778
 544 GTGTACCACACACTTCCTACAGCTGT       186 26  0 50.00 62.885  8.00  8.00  8.885
 545 ACTGAAGATGTTCCCTAAGAATGTGG       108 26  0 42.31 62.911  6.00  0.00  8.911
 546 GAATGCACTGAAGACCTTCCCTAATA       140 26  0 42.31 62.924  6.00  3.00  8.924
 547 GCACTGAAGACCTTCCCTAATACAAA       144 26  0 42.31 62.934  6.00  0.00  8.934
 548 CACTCACTGAAGACTTCCCTACATACA        5 27  0 44.44 62.050  7.00  0.00  9.050
 549 CACACTCACTGAAGACTTCCCTACATA        3 27  0 44.44 62.050  7.00  2.00  9.050
 550 TACAAAAGCTAAGACATTGTGGGTGTA      164 27  0 37.04 62.136  5.00  2.00  9.136
 551 ATACAAAAGCTAAGACATTGTGGGTGT      163 27  0 37.04 62.320  5.00  0.00  9.320
 552 TGTACACACTGAAGATGTTCCCTAAGA      101 27  0 40.74 62.355  6.00  1.00  9.355
 553 ACTGAAGACCTTCCCTAATACAAAAGC      146 27  0 40.74 62.397  6.00  2.00  9.397
 554 CTGAAGACCTTCCCTAATACAAAAGCT      147 27  0 40.74 62.483  6.00  4.00  9.483
 555 CATGTACACACTGAAGATGTTCCCTAA       99 27  0 40.74 62.561  6.00  1.00  9.561
 556 AATACAAAAGCTAAGACATTGTGGGTG      162 27  0 37.04 62.689  5.00  0.00  9.689
 557 ACAAAAGCTAAGACATTGTGGGTGTAC      165 27  0 40.74 62.787  5.00  4.00  9.787
 558 GATCACACTCACTGAAGACTTCCCTAC        0 27  0 48.15 62.829  7.00  1.00  9.829
 559 CTCACTGAAGACTTCCCTACATACAGC        7 27  0 48.15 62.853  7.00  2.00  9.853
 560 ACACTGAAGATGTTCCCTAAGAATGTG      106 27  0 40.74 62.934  6.00  0.00  9.934
 561 CACACTGAAGATGTTCCCTAAGAATGT      105 27  0 40.74 62.934  6.00  1.00  9.934
 562 ACACACTGAAGATGTTCCCTAAGAATG      104 27  0 40.74 62.934  6.00  0.00  9.934