[Bioperl-l] LocatableSeq NSE format - suggested changes
Chris Fields
cjfields at illinois.edu
Tue Dec 15 01:23:05 UTC 2009
All,
The current output for NSE format (Name/Start-End) via Bio::LocatableSeq::get_nse() currently doesn't allow for strandedness. I have seen two variations of NSE that incorporate strandedness:
1) Stockholm Rfam reverses start and end if the strand == -1
chrY/598-1
2) Sheldon McKay's Gbrowse_syn uses Name(strand)/start-end
rice-3(+)/16598648-16600199
The former breaks fewer things within BioPerl, but the latter seems more explicit. Any preferences? Do we want a new method that creates this, and deprecate out simple non-stranded NSE?
chris
More information about the Bioperl-l
mailing list