Body: <p>I'm analyzing mygenome.fa with RepeatScout and RepeatMasker to find transposable elements. I produced a library of Repetitive elements with repeatscout and the I masked TE using Repeatmasker. </p>
  
  <pre><code>./build_lmer_table -l 14 -sequence mygenome.fa -freq mygenome.freq
  ./RepeatScout -sequence mygenome.fa -output mygenome_repeats.fa -freq mygenome.freq -l 14
  cat mygenome_repeats.fa| ./filter-stage-1.prl &gt;mygenome_repeats_filtered1.fasta
  ./RepeatMasker -s -lib mygenome_repeats_filtered1.fa mygenome.fa
  </code></pre>
  
  <p>Generate a masked genome using (non-low-complexity, non-tandem) repeats</p>
  
  <pre><code>cat mygenome_repeats_filtered1.fa | ./filter-stage-2.prl --cat mygenome.fa.out &gt;mygenome_repeats_filtered2.fa
  </code></pre>
  
  <p>Filter out all (non-low-complexity, non-tandem) repeats that have less than 10 repeats</p>
  
  <pre><code>./RepeatMasker -pa 4 -s -lib mygenome_repeats_filtered2.fa -nolow -norna -no_is -gff mygenome.fa
  </code></pre>
  
  <p>After this I ran Repeatmasker: </p>
  
  <pre><code>./RepeatMasker no_is mygenome_repeats_filtered2.fa
  </code></pre>
  
  <p>and I produced:</p>
  
  <pre><code>mygenome_repeats_filt2.fa.masked  mygenome_repeats_filt2.fa.tbl
  mygenome_repeats_filt2.fa.cat  mygenome_repeats_filt2.fa.out
  </code></pre>
  
  <p>I would to know a way to find the location of transposable elements masked by RepeatMasker on my genome. </p>
  
  <p>Supposing that this is a part of my Repetitive elements library produced with RepeatScout as I showed above:</p>
  
  <pre><code>less mygenome_repeats_filt2.fa
  
  &gt;R=3 (RR=4.  TRF=0.000 NSEG=0.000)
  TAAGGCGGCGAGCTGGCAGAATCGTTAGCACGCCGGGCGAAATGCTTAGCGGTATTTCGTCTGTCTTTACGTTCTGAGTT
  CAAATTCCGCCGAGGTCGACTTTGCCTTTCATCCTTTCGGGGTCGATAAAATAAGTACCAGTTGAGCACTGGGGTCGATG
  TAATCGACTTACCCCCTCCCCCAAAATTTCTGGCCTTGTGCCTATATTAGAAACGATTATT
  &gt;R=4 (RR=5.  TRF=0.122 NSEG=0.226)
  ACACACACACACACACACACACACACATATATATATATATACATATATACGACGGGCTTCTTTCAGTTTCCGTCTACCAA
  ATCCACTCACAAGGCTTTGGTCGGCCCGAGGCTATAGTAGAAGACACTTGCCCAAGGTGCCACGCAGTGGGACTGAACCC
  GGAACCATGTGGTTGGTAAGCAAGCTACTTACCACACAGCCACTCCTGCGCCTATATATAT
  &gt;R=6 (RR=7.  TRF=0.134 NSEG=0.247)
  TTGTTTCAGTCATTTGACTGCGGCCATGCTGGAGCACCGCCTTTAGTCGAGCAAATCGACCCCAGGACTTATTCTTTGTA
  AGCCTAGTACTTATTCTATCGGTCTCTTTTGCCGAACCGCTAAGTTACGGGGACGTAAACACACCAGCATCGGTTGTCAA
  GCGATGTTGGGGGGACAAACACAGACACACAAACACACACACACACATACATATATATATATATATATATA
  </code></pre>
  
  <p>from the file mygenome_repeats_filt2.fa.out I can see that in R=3 there is a trasposable element:</p>
  
  <pre><code>   SW   perc perc perc  query     position in query    matching        repeat              position in repeat
  score   div. del. ins.  sequence  begin end   (left)   repeat          class/family      begin  end    (left)  ID
  
  348   20.6  0.0  0.0  R=3          78   140   (81) + AmnSINE2        SINE/tRNA-Deu         67    129  (229) 234  
  </code></pre>
  
  <p>As you can see I have the coordinates to find this element in the library, but I would to find its exactly in the file that contain my assembled genome. </p>
  
Comments.ContentLicense: CC BY-SA 4.0
Comments.CreationDate: 2019-05-09T13:12:22.397Z
Comments.Text: Hi. In order to help you better could you explain what have you tried to do this? This way our answers can explain how to find the solution (not only the solution itself) or how to improve your skills.
Comments.UserDisplayName: llrs
Comments.ContentLicense: CC BY-SA 4.0
Comments.CreationDate: 2019-05-09T13:30:08.840Z
Comments.Text: ok, I modified the answer :)
Comments.UserDisplayName: jonny
Comments.ContentLicense: CC BY-SA 4.0
Comments.CreationDate: 2019-05-09T13:47:45.930Z
Comments.Text: Great edit!! But this code does not show what you tried to "find the location of transposable elements masked by RepeatMasker" (or at least is not clear enough for me that I am not familiar with RepeatMasker).
Comments.UserDisplayName: llrs
Comments.ContentLicense: CC BY-SA 4.0
Comments.CreationDate: 2019-05-09T13:48:44.843Z
Comments.Text: sorry,I'll modified again the answer and I'll try to specify it better
Comments.UserDisplayName: jonny
Comments.ContentLicense: CC BY-SA 4.0
Comments.CreationDate: 2019-05-09T13:58:06.120Z
Comments.Text: I modified the file with some other informations
Comments.UserDisplayName: jonny
Tags: <repeatmasker>
Title: find repetitive elements on genome
ThreadText: # find repetitive elements on genome
  
  (Asked by: username_0 on 2019-05-09)
  
  I'm analyzing mygenome.fa with RepeatScout and RepeatMasker to find transposable elements. I produced a library of Repetitive elements with repeatscout and the I masked TE using Repeatmasker.
  
  ```
  ./build_lmer_table -l 14 -sequence mygenome.fa -freq mygenome.freq
  ./RepeatScout -sequence mygenome.fa -output mygenome_repeats.fa -freq mygenome.freq -l 14
  cat mygenome_repeats.fa| ./filter-stage-1.prl >mygenome_repeats_filtered1.fasta
  ./RepeatMasker -s -lib mygenome_repeats_filtered1.fa mygenome.fa
  
  ```
  
  Generate a masked genome using (non-low-complexity, non-tandem) repeats
  
  ```
  cat mygenome_repeats_filtered1.fa | ./filter-stage-2.prl --cat mygenome.fa.out >mygenome_repeats_filtered2.fa
  
  ```
  
  Filter out all (non-low-complexity, non-tandem) repeats that have less than 10 repeats
  
  ```
  ./RepeatMasker -pa 4 -s -lib mygenome_repeats_filtered2.fa -nolow -norna -no_is -gff mygenome.fa
  
  ```
  
  After this I ran Repeatmasker:
  
  ```
  ./RepeatMasker no_is mygenome_repeats_filtered2.fa
  
  ```
  
  and I produced:
  
  ```
  mygenome_repeats_filt2.fa.masked  mygenome_repeats_filt2.fa.tbl
  mygenome_repeats_filt2.fa.cat  mygenome_repeats_filt2.fa.out
  
  ```
  
  I would to know a way to find the location of transposable elements masked by RepeatMasker on my genome.
  
  Supposing that this is a part of my Repetitive elements library produced with RepeatScout as I showed above:
  
  ```
  less mygenome_repeats_filt2.fa
  
  >R=3 (RR=4.  TRF=0.000 NSEG=0.000)
  TAAGGCGGCGAGCTGGCAGAATCGTTAGCACGCCGGGCGAAATGCTTAGCGGTATTTCGTCTGTCTTTACGTTCTGAGTT
  CAAATTCCGCCGAGGTCGACTTTGCCTTTCATCCTTTCGGGGTCGATAAAATAAGTACCAGTTGAGCACTGGGGTCGATG
  TAATCGACTTACCCCCTCCCCCAAAATTTCTGGCCTTGTGCCTATATTAGAAACGATTATT
  >R=4 (RR=5.  TRF=0.122 NSEG=0.226)
  ACACACACACACACACACACACACACATATATATATATATACATATATACGACGGGCTTCTTTCAGTTTCCGTCTACCAA
  ATCCACTCACAAGGCTTTGGTCGGCCCGAGGCTATAGTAGAAGACACTTGCCCAAGGTGCCACGCAGTGGGACTGAACCC
  GGAACCATGTGGTTGGTAAGCAAGCTACTTACCACACAGCCACTCCTGCGCCTATATATAT
  >R=6 (RR=7.  TRF=0.134 NSEG=0.247)
  TTGTTTCAGTCATTTGACTGCGGCCATGCTGGAGCACCGCCTTTAGTCGAGCAAATCGACCCCAGGACTTATTCTTTGTA
  AGCCTAGTACTTATTCTATCGGTCTCTTTTGCCGAACCGCTAAGTTACGGGGACGTAAACACACCAGCATCGGTTGTCAA
  GCGATGTTGGGGGGACAAACACAGACACACAAACACACACACACACATACATATATATATATATATATATA
  
  ```
  
  from the file mygenome\_repeats\_filt2.fa.out I can see that in R=3 there is a trasposable element:
  
  ```
     SW   perc perc perc  query     position in query    matching        repeat              position in repeat
  score   div. del. ins.  sequence  begin end   (left)   repeat          class/family      begin  end    (left)  ID
  
  348   20.6  0.0  0.0  R=3          78   140   (81) + AmnSINE2        SINE/tRNA-Deu         67    129  (229) 234  
  
  ```
  
  As you can see I have the coordinates to find this element in the library, but I would to find its exactly in the file that contain my assembled genome.
  
  **Tags:** repeatmasker
  
  > **Comments:**
  > **username_1**: Hi. In order to help you better could you explain what have you tried to do this? This way our answers can explain how to find the solution (not only the solution itself) or how to improve your skills.  
  > **username_0**: ok, I modified the answer :)  
  > **username_1**: Great edit!! But this code does not show what you tried to "find the location of transposable elements masked by RepeatMasker" (or at least is not clear enough for me that I am not familiar with RepeatMasker).  
  > **username_0**: sorry,I'll modified again the answer and I'll try to specify it better  
  > **username_0**: I modified the file with some other informations  
  
  ---
  
  ## 0 answers
  
  