<div dir="ltr"><div><div><div><div>Dear Carson:<br><br></div>Thank you again for all you valuable suggestions. Now I am generating the species specific repeat library. I wonder whether I need to remove the masked the regions by existing repeatMasker library, before I run repeatModeler? I think there may be some redundancy if I run  repeatModeler directly on the genome and then use both existing repeatMasker library and the repeatModeler library to mask the genome. Does it matter, if there is such redundancy?<br><br></div>Thanks<br><br></div>Best<br></div>Quanwei<br></div><div class="gmail_extra"><br><div class="gmail_quote">2017-08-23 14:10 GMT-04:00 Carson Holt <span dir="ltr"><<a href="mailto:carsonhh@gmail.com" target="_blank">carsonhh@gmail.com</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word"><div><span class=""><br><blockquote type="cite"><div><div dir="ltr"><div><div>(1) For the predicted unknown (unclassified) repeat sequences (those in Modelerunknown.lib), it mentioned "Sequences in Modelerunknown.lib were searched against a transposase database (derived from <a rel="nofollow" class="m_-3343096424557117677external m_-3343096424557117677gmail-text" href="http://www.repeatmasker.org/" target="_blank">RepeatMaske</a>r) and sequences matching transposase were considered as transposons belonging to the relevant superfamily".  <br></div><div>I wonder how to do this search. Annotate the "unknown" repeat sequences using the Repeatmaker? Then what to do, if for an "unknown" repeat sequence, only part of the sequence match the known repeat elements.<br></div></div></div></div></blockquote><div><br></div></span><div>You can use RepBase match I guess, but I would not be overly worried about classification. MAKER won’t use any classification info you give it.</div><span class=""><div><br></div><br><blockquote type="cite"><div><div dir="ltr"><div>(2) To exclude gene fragments, I need map the predicted repeat sequences against a protein database, and then run the package "ProExcluder"<b>. </b>Right?<b>  </b>I wonder how to get such protein database.<b> </b>Since I am working on a new rodent species, can I use all the rodent proteins from Uniprot (both Swiss-Prot and TrEMBL)?<br></div></div></div></blockquote><div><br></div></span><div>Try Swiss-Prot. That is a well curated cross species set.</div><span class=""><div><br></div><br><blockquote type="cite"><div><div dir="ltr">(3) After I generate the species specific repeat library, do I still need to select a model organism for RepBase masking (as shown below). <br><div><br></div><div>In the file "maker_opts.ctl"<br>#-----Repeat Masking (leave values blank to skip repeat masking)<br>model_org=Mammalia #select a model organism for RepBase masking in RepeatMasker<br>rmlib=myRepeat.fa #provide an organism specific repeat library in fasta format for RepeatMasker<br></div></div></div></blockquote></span></div><br><div>Yes. Supply both.</div><span class="HOEnZb"><font color="#888888"><div><br></div><div><br></div><div>—Carson</div></font></span></div></blockquote></div><br></div>