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