Hi --

I would like to suggest an addition to the PSL format example given here:
http://genome.ucsc.edu/FAQ/FAQformat#format2

The current example is

0         1         2         3 tens position in query
0123456789012345678901234567890 ones position in query
            ++++          +++++ plus strand alignment on query
    --------    ----------      minus strand alignment on query

Plus strand:
     qStart=12
     qEnd=31
     blockSizes=4,5
     qStarts=12,26

Minus strand:
     qStart=4
     qEnd=26
     blockSizes=10,8
     qStarts=5,19

To an ignoramus (me) trying to puzzle this out, the Minus strand
qStart and qEnd values seem ambiguous. The stand is 30 bases long, so
the coordinates 4 and 26 are not informative about the correct
directional relationship. By adding 10 bases to the query length, the
ambiguity is resolved:

0         1         2         3         4 tens position in query
01234567890123456789012345678901234567890 ones position in query
            ++++          +++++           plus strand alignment on query
    --------    ----------                minus strand alignment on query

Plus strand:
     qStart=12
     qEnd=31
     blockSizes=4,5
     qStarts=12,26

Minus strand:
     qStart=4
     qEnd=26
     blockSizes=10,8
     qStarts=15,29

Is my extension of the example correct?

Thanks
_______________________________________________
Genome maillist  -  [email protected]
https://lists.soe.ucsc.edu/mailman/listinfo/genome

Reply via email to