mmiladi / MutaRNA

1 stars 3 forks source link

new option '--idxPos0' (default=1) to guide indexing #17

Open martin-raden opened 4 years ago

martin-raden commented 4 years ago

(similar to IntaRNA --t/qIdxPos0)

mmiladi commented 4 years ago

I didn't get the idea..

martin-raden commented 4 years ago

right now, indexing always starts with 1. but mRNAs are usually indexed by the first position of the start codon. all indices upstream of the start codon are numbered negatively (skipping 0).

all you need for proper indexing supportign this scenario is to know what is the index of the first sequence position (idxPos0). thus, you need to map all input/output indices to internal 0-based indexing, best via two respective functions. that way, you can change sequence context length of an RNA around a mutation without changing the mutation index (by setting respective idxPos0 values).

That's quite useful if you want to provide figures/results for mRNAs that can be directly used afterwards without reindexing..

better explained? ;)

mmiladi commented 4 years ago

Thanks :-)