search

mcfrith / last-genome-alignments

47 stars 5 forks source link

reference:once; query: no limit #13

Closed AlisaGU closed 2 years ago

AlisaGU commented 2 years ago

Hi, If I want the reference sequence to appear only once but no limit to query, can I do like that?

  1. query as reference and reference as query
  2. after get the many to one maf, use maf-swap to change the order of orginal reference and query and use this as the final result
AlisaGU commented 2 years ago

The reason I proposed this question is I noticed the many to one maf(4.4G, after remove the quality score line) is so big compared to one to one maf(21M).

mcfrith commented 2 years ago

Yes, swapping query and reference should do what you want. (I don't think it's necessary to use maf-swap afterwards, unless you have some reason I don't know about.)