<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>

CC BY-SA 4.0
2019-05-09T13:12:22.397Z
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.
llrs
CC BY-SA 4.0
2019-05-09T13:30:08.840Z
ok, I modified the answer :)
jonny
CC BY-SA 4.0
2019-05-09T13:47:45.930Z
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).
llrs
CC BY-SA 4.0
2019-05-09T13:48:44.843Z
sorry,I'll modified again the answer and I'll try to specify it better
jonny
CC BY-SA 4.0
2019-05-09T13:58:06.120Z
I modified the file with some other informations
jonny
<repeatmasker>
find repetitive elements on genome
# 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

