New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 891259102 (1.570)5 (0.080) 4/142 7.8 5.2% coding (282/963 nt) ECB_00831 protein similar to RNA‑directed DNA polymerase from retron EC86
?REL606 = 891492 85 (1.390)coding (515/963 nt) ECB_00831 protein similar to RNA‑directed DNA polymerase from retron EC86
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

AGTTTTGACAAAAATATGTTGTTATAAAAATCTGCTACCACAAGGTGCTCCATCATCACCTAAATTAGCTAATCTAATATGTTCTAAACTTGATTATCGTATTCAGGGTTATGCAGGTAGTCGGGGCTTGATATATACGAGATATGCCG  >  REL606/891418‑891566
                                                                          |                                                                          
aGTTTTGACAAAAATATGTTGTTATAAAAATCTGCTACCACAAGGTGCTCCATCATCACCTAAATTAGctaatcta                                                                           <  1:1817224/76‑1 (MQ=255)
aGTTTTGACAAAAATATGTTGTTATAAAAATCTGCTACCACAAGGTGCTCCATCATCACCTAAATTAGctaatcta                                                                           <  2:1465615/76‑1 (MQ=255)
  ttttGACAAAAATATGTTGTTATAAAAATCTGCTACCACAAGGTGCTCCATCATCACCTAAATTAGctaatctaat                                                                         <  2:1198571/76‑1 (MQ=255)
   tttGACAAAAATATGTTGTTATAAAAATCTGCTACCACAAGGTGCTCCATCATCACCTAAATTAGCTAATCTAata                                                                        >  2:1707501/1‑76 (MQ=255)
    ttGACAAACATATGTTGTTATAAAAATCTGCTACCACAAGGTGCTCCATCATCACCTAAATTAGCTaatataatat                                                                       >  2:579019/1‑76 (MQ=255)
     tGACAAAAATATGTTGTTATAAAAATCTGCTACCACAAGGTGCTCCATCATCACCTAAATTAGCTAATCTAATATg                                                                      >  2:1313359/1‑76 (MQ=255)
         aaaaaTATGTTGTTATAAAAATCTGCTACCACAAGGTGCTCCATCATCACCTAAATTAGCTAATCTAATATGTTCt                                                                  >  2:1793003/1‑76 (MQ=255)
         aaaaaTATGTTGTTATAAAAATCTGCTACCACAAGGTGCTCCATCATCACCTAAATTAGCTAATCTAATATGTTCt                                                                  <  1:243497/76‑1 (MQ=255)
         aaaaaTATGTTGTTATAAAAATCTGCTACCACAAGGTGCTCCATCATCACCTAAATTAGCTAATCTAATATGTTCt                                                                  <  2:786384/76‑1 (MQ=255)
         aaaaaTATGTTGTTATAAAAATCTGCTACCACAAGGTGCTCCATCATCACCTAAATTAGCTAATCTAATATGTTCt                                                                  <  1:1459362/76‑1 (MQ=255)
         aaaaaTATGTTGTTATAAAAATCTGCTACCACAAGGTGCTCCATCATCACCTAAATTAGCTAATCTAATATGTTCt                                                                  <  1:1459357/76‑1 (MQ=255)
          aaaaTATGTTGTTATAAAAATCTGCTACCACAAGGTGCTCCATCATCACCTAAATTAGCTAATCTAATATGTTCTa                                                                 <  2:146700/76‑1 (MQ=255)
          aaaaTATGTTGTTATAAAAATCTGCTACCACAAGGTGCTCCATCATCACCTAAATTAGCTAATCTAATATGTTCTa                                                                 <  2:1502197/76‑1 (MQ=255)
          aaaaTATGTTGTTATAAAAATCTGCTACCACAAGGTGCTCCATCATCACCTAAATTAGCTAATCTAATATGTTCTa                                                                 <  2:1665742/76‑1 (MQ=255)
           aaaTATGTTGTTATAAAAATCTGCTACCACAAGGTGCTCCATCATCACCTAAATTAGCTAATCTAATATGTTCTaa                                                                <  1:1012119/76‑1 (MQ=255)
            aaTATGTTGTTATAAAAATCTGCTACCACAAGGTGCTCCATCATCACCTAAATTAGCTAATCTAATATGTTCTaaa                                                               <  1:121085/76‑1 (MQ=255)
               atGTTGTTATAAAAATCTGCTACCACAAGGTGCTCCATCATCACCTAAATTAGCTAATCTAATATGTTCTAAACtt                                                            >  1:2091500/1‑76 (MQ=255)
                tgttgtTATAAAAATCTGCTACCACAAGGTGCTCCATCATCACCTAAATTAGCTAATCTAATATGTTCTAAACTTg                                                           <  1:170159/76‑1 (MQ=255)
                 gttgttATAAAAATCTGCTACCACAAGGTGCTCCATCATCACCTAAATTAGCTAATCTAATATGTTCTAAACTTGa                                                          <  1:1578291/76‑1 (MQ=255)
                 gttgttATAAAAATCTGCTACCACAAGGTGCTCCATCATCACCTAAATTAGCTAATCTAATATGTTCTAAACTTGa                                                          <  1:1098083/76‑1 (MQ=255)
                  ttgttATAAAAATCTGCTACCACAAGGTGCTCCATCATCACCTAAATTAGCTAATCTAATATGTTCTAAACTTGAt                                                         >  1:1854438/1‑76 (MQ=255)
                  ttgttATAAAAATCTGCTACCACAAGGTGCTCCATCATCACCTAAATTAGCTAATCTAATATGTTCTAAACTTGAt                                                         >  2:1952581/1‑76 (MQ=255)
                   tgttATAAAAATCTGCTACCACAAGGTGCTCCATCATCACCTAAATTAGCTAATCTAATATGTTCTAAACTTGAtt                                                        <  1:813216/76‑1 (MQ=255)
                   tgttATAAAAATCTGCTACCACAAGGTGCTCCATCATCACCTAAATTAGCTAATCTAATATGTTCTAAACTTGAtt                                                        >  1:1996631/1‑76 (MQ=255)
                   tgttATAAAAATCTGCTACCACAAGGTGCTCCATCATCACCTAAATTAGCTAATCTAATATGTTCTAAACTTGAtt                                                        >  1:762484/1‑76 (MQ=255)
                   tgttATAAAAATCTGCTACCACAAGGTGCTCCATCATCACCTAAATTAGCTAATCTAATATGTTCTAAACTTGAtt                                                        >  1:1601003/1‑76 (MQ=255)
                     ttATAAAAATCTGCTACCACAAGGTGCTCCATCATCACCTAAATTAGCTAATCTAATATGTTCTAAACTTGATTAt                                                      >  1:1770837/1‑76 (MQ=255)
                      tataAAAATCTGCTACCACAAGGTGCTCCATCATCACCTAAATTAGCTAATCTAATATGTTCTAAACTTGATTATc                                                     <  2:1343717/76‑1 (MQ=255)
                       ataAAAATCTGCTACCACAAGGTGCTCCATCATCACCTAAATTAGCTAATCTAATATGTTCTAAACTTGATTATCg                                                    <  1:379311/76‑1 (MQ=255)
                       ataAAAATCTGCTACCACAAGGTGCTCCATCATCACCTAAATTAGCTAATCTAATATGTTCTAAACTTGATTATCg                                                    >  2:2095692/1‑76 (MQ=255)
                        taAAAATCTGCTACCACAAGGTGCTCCATCATCACCTAAATTAGCTAATCTAATATGTTCTAAACTTGATTATCGt                                                   <  2:291378/76‑1 (MQ=255)
                        taAAAATCTGCTACCACAAGGTGCTCCATCATCACCTAAATTAGCTAATCTAATATGTTCTAAACTTGATTATCGt                                                   >  1:1168330/1‑76 (MQ=255)
                        taAAAATCTGCTACCACAAGGTGCTCCATCATCACCTAAATTAGCTAATCTAATATGTTCTAAACTTGATTATCGt                                                   <  1:200354/76‑1 (MQ=255)
                         aaaaaTCTGCTACCACAAGGTGCTCCATCATCACCTAAATTAGCTAATCTAATATGTTCTAAACTTGATTATCGTa                                                  >  1:303706/1‑76 (MQ=255)
                         aaaaaTCTGCTACCACAAGGTGCTCCATCATCACCTAAATTAGCTAATCTAATATGTTCTAAACTTGATTATCGTa                                                  >  2:1971913/1‑76 (MQ=255)
                           aaaTCTGCTACCACAAGGTGCTCCATCATCACCTAAATTAGCTAATCTAATATGTTCTAAACTTGATTATCGTAtt                                                <  1:1897127/76‑1 (MQ=255)
                            aaTCTGCTACCACAAGGTGCTCCATCATCACCTAAATTAGCTAATCTAATATGTTCTAAACTTGATTATCGTAtta                                               <  2:297919/76‑2 (MQ=255)
                            aaTCTGCTACCACAAGGTGCTCCATCATCACCTAAATTAGCTAATCTAATATGTTCTAAACTTGATTATCGTATTc                                               <  1:1818656/76‑1 (MQ=255)
                               cTGCTACCACAAGGTGCTCCATCATCACCTAAATTAGCTAATCTAATATGTTCTAAACTTGATTATCGTATTCAgg                                            >  2:378465/1‑76 (MQ=255)
                                tGCTACCACAAGGTGCTCCATCATCACCTAAATTAGCTAATCTAATATGTTCTAAACTTGATTATCGTATTCAggg                                           <  1:1107243/76‑1 (MQ=255)
                                  cTACCACAAGGTGCTCCATCATCACCTAAATTAGCTAATCTAATATGTTCTAAACTTGATTATCGTATTCAGGGtt                                         >  1:201494/1‑76 (MQ=255)
                                  cTACCACAAGGTGCTCCATCATCACCTAAATTAGCTAATCTAATATGTTCTAAACTTGATTATCGTATTCAGGGtt                                         <  1:963117/76‑1 (MQ=255)
                                   tACCACAAGGTGCTCCATCATCACCTAAATTAGCTAATCTAATATGTTCTAAACTTGATTATCGTATTCAGGGTTa                                        >  2:1836869/1‑76 (MQ=255)
                                   tACCACAAGGTGCTCCATCATCACCTAAATTAGCTAATCTAATATGTTCTAAACTTGATTATCGTATTCAGGGTTa                                        >  2:1903727/1‑76 (MQ=255)
                                    aCCACAAGGTGCTCCATCATCACCTAAATTAGCTAATCTAATATGTTCTAAACTTGATTATCGTATTCAGGGTTAt                                       <  2:795141/76‑1 (MQ=255)
                                    aCCACAAGGTGCTCCATCATCACCTAAATTAGCTAATCTAATATGTTCTAAACTTGATTATCGTATTCAGGGTTAt                                       >  2:1300126/1‑76 (MQ=255)
                                    aCCACAAGGTGCTCCATCATCACCTAAATTAGCTAATCTAATATGTTCTAAACTTGATTATCGTATTCAGGGTTAt                                       >  2:1559904/1‑76 (MQ=255)
                                    aCCACAAGGTGCTCCATCATCACCTAAATTAGCTAATCTAATATGTTCTAAACTTGATTATCGTATTCAGGGTTAt                                       >  2:1381585/1‑76 (MQ=255)
                                     ccACAAGGTGCTCCATCATCACCTAAATTAGCTAATCTAATATGTTCTAAACTTGATTATCGTATTCAGGGTTATg                                      <  2:1546572/76‑1 (MQ=255)
                                      cacaAGGTGCTCCATTATCGCCTAAATTGGCTAATGTACTATGTTCCAAACTTGGTTATTGTATTCAGGGTTCTGc                                     >  2:719466/1‑76 (MQ=255)
                                      cacaAGGTGCTCCATCATCACCTAAATTAGCTAATCTAATATGTTCTAAACTTGATTATCGTATTCAGGGTTATGc                                     <  2:1181426/76‑1 (MQ=255)
                                      cacaAGGTGCTCCATCATCACCTAAATTAGCTAATCTAATATGTTCTAAACTTGATTATCGTATTCAGGGTTATGc                                     <  1:1644856/76‑1 (MQ=255)
                                      cacaAGGTGCTCCATCATCACCTAAATTAGCTAATCTAATATGTTCTAAACTTGATTATCGTATTCAGGGTTATGc                                     <  1:1140813/76‑1 (MQ=255)
                                       acaAGGTGCTCCATCATCACCTAAATTAGCTAATCTAATATGTTCTAAACTTGATTATCGTATTCAGGGTTATGCa                                    >  2:643503/1‑76 (MQ=255)
                                         aaGGTGCTCCATCATCACCTAAATTAGCTAATCTAATATGTTCTAAACTTGATTATCGTATTCAGGGTTATGCAgg                                  >  2:136789/1‑76 (MQ=255)
                                         aaGGTGCTCCATCATCACCTAAATTAGCTAATCTAATATGTTCTAAACTTGATTATCGTATTCAGGGTTATGCAgg                                  <  1:847616/76‑1 (MQ=255)
                                          aGGTGCTCCATCATCACCTAAATTAGCTAATCTAATATGTTCTAAACTTGATTATCGTATTCAGGGTTATGCAGGt                                 <  2:1922623/76‑1 (MQ=255)
                                          aGGTGCTCCATCATCACCTAAATTAGCTAATCTAATATGTTCTAAACTTGATTATCGTATTCAGGGTTATGCAGGt                                 <  2:1644639/76‑1 (MQ=255)
                                          aGGTGCTCCATCATCACCTAAATTAGCTAATCTAATATGTTCTAAACTTGATTATCGTATTCAGGGTTATGCAGGt                                 <  2:512433/76‑1 (MQ=255)
                                          aGGTGCTCCATCATCACCTAAATTAGCTAATCTAATATGTTCTAAACTTGATTATCGTATTCAGGGTTATGCAGGt                                 <  2:1922610/76‑1 (MQ=255)
                                            gTGCTCCATCATCACCTAAATTAGCTAATCTAATATGTTCTAAACTTGATTATCGTATTCAGGGTTATGCAGGTAg                               <  1:482539/76‑1 (MQ=255)
                                            gTGCTCCATCATCACCTAAATTAGCTAATCTAATATGTTCTAAACTTGATTATCGTATTCAGGGTTATGCAGGTAg                               >  2:689366/1‑76 (MQ=255)
                                            gTGCTCCATCATCACCTAAATTAGCTAATCTAATATGTTCTAAACTTGATTATCGTATTCAGGGCTATGCAGGTAg                               >  2:1292727/1‑76 (MQ=255)
                                             tGCTCCATCATCACCTAAATTAGCTAATCTAATATGTTCTAAACTTGATTATCGTATTCAGGGTTATGCAGGTAGt                              >  1:366383/1‑76 (MQ=255)
                                              gCTCCATCATCACCTAAATTAGCTAATCTAATATGTTCTAAACTTGATTATCGTATTCAGGGTTATGCAGGTAGTc                             <  1:776723/76‑1 (MQ=255)
                                               cTCCATCATCACCTAAATTAGCTAATCTAATATGTTCTAAACTTGATTATCGTATTCAGGGTTATGCAGGTAGTCg                            <  1:1052486/76‑1 (MQ=255)
                                               cTCCATCATCACCTAAATTAGCTAATCTAATATGTTCTAAACTTGATTATCGTATTCAGGGTTATGCAGGTAGTCg                            <  1:1033250/76‑1 (MQ=255)
                                                  catcatCACCTAAATTAGCTAATCTAATATGTTCTAAACTTGATTATCGTATTCAGGGTTATGCAGGTAGTCgggg                         <  1:1924923/76‑1 (MQ=255)
                                                   atcatcACCTAAATTAGCTAATCTAATATGTTCTAAACTTGATTATCGTATTCAGGGTTATGCAGGTAGTCGGGGc                        >  2:1841814/1‑76 (MQ=255)
                                                   atcatcACCTAAATTAGCTAATCTAATATGTTCTAAACTTGATTATCGTATTCAGGGTTATGCAGGTAGTCGGGGc                        >  2:1444626/1‑76 (MQ=255)
                                                    tcatcaCCTAAATTAGCTAATCTAATATGTTCTAAACTTGATTATCGTATTCAGGGTTATGCAGGTAGTCGGGGCt                       >  2:1770773/1‑76 (MQ=255)
                                                     catcaCCTAAATTAGCTAATCTAATATGTTCTAAACTTGATTATCGTATTCAGGGTTATGCAGGTAGTCGGGGCtt                      >  1:1793016/1‑76 (MQ=255)
                                                     catcaCCTAAATTAGCTAATCTAATATGTTCTAAACTTGATTATCGTATTCAGGGTTATGCAGGTAGTCGGGGCtt                      <  1:1787827/76‑1 (MQ=255)
                                                      atcaCCTAAATTAGCTAATCTAATATGTTCTAAACTTGATTATCGTATTCAGGGTTATGCAGGTAGTCGGGGCTTg                     <  2:9204/76‑1 (MQ=255)
                                                      atcaCCTAAATTAGCTAATCTAATATGTTCTAAACTTGATTATCGTATTCAGGGTTATGCAGGTAGTCGGGGCTTg                     <  2:1973000/76‑1 (MQ=255)
                                                      atcaCCTAAATTAGCTAATCTAATATGTTCTAAACTTGATTATCGTATTCAGGGTTATGCAGGTAGTCGGGGCTTg                     <  2:1320510/76‑1 (MQ=255)
                                                        caCCTAAATTAGCTAATCTAATATGTTCTAAACTTGATTATCGTATTCAGGGTTATGCAGGTAGTCGGGGCTTGat                   >  1:1768737/1‑76 (MQ=255)
                                                         aCCTAAATTAGCTAATCTAATATGTTCTAAACTTGATTATCGTATTCAGGGTTATGCAGGTAGTCGGGGCTTGata                  >  2:298268/1‑76 (MQ=255)
                                                          ccTAAATTAGCTAATCTAATATGTTCTAAACTTGATTATCGTATTCAGGGTTATGCAGGTAGTCGGGGCTTGatat                 >  2:561474/1‑76 (MQ=255)
                                                          ccTAAATTAGCTAATCTAATATGTTCTAAACTTGATTATCGTATTCAGGGTTATGCAGGTAGTCGGGGCTTGatat                 >  2:1264856/1‑76 (MQ=255)
                                                           cTAAATTAGCTAATCTAATATGTTCTAAACTTGATTATCGTATTCAGGGTTATGCAGGTAGTCGGGGCTTGatata                >  2:263517/1‑76 (MQ=255)
                                                            tAAATTAGCTAATCTAATATGTTCTAAACTTGATTATCGTATTCAGGGTTATGCAGGTAGTCGGGGCTTGatatat               <  1:303568/76‑1 (MQ=255)
                                                             aaaTTAGCTAATCTAATATGTTCTAAACTTGATTATCGTATTCAGGGTTATGCAGGTAGTCGGGGCTTGAtatata              <  1:1202989/76‑1 (MQ=255)
                                                               aTTAGCTAATCTAATATGTTCTAAACTTGATTATCGTATTCAGGGTTATGCAGGTAGTCGGGGCTTGATATATACg            >  2:472500/1‑76 (MQ=255)
                                                                     taatctaatATGTTCTAAACTTGATTATCGTATTCAGGGTTATGCAGGTAGTCGGGGCTTGATATATACGAGatat      >  2:846283/1‑76 (MQ=255)
                                                                        tctaatATGTTCTAAACTTGATTATCGTATTCAGGGTTATGCAGGTAGTCGGGGCTTGATATATACGAGATATGcc   >  1:990250/1‑76 (MQ=255)
                                                                        tctaatATGTTCTAAACTTGATTATCGTATTCAGGGTTATGCAGGTAGTCGGGGCTTGATATATACGAGATATGcc   >  1:375798/1‑76 (MQ=255)
                                                                         ctaatATGTTCTAAACTTGATTATCGTATTCAGGGTTATGCAGGTAGTCGGGGCTTGATATATACGAGATATGCCg  >  2:769342/1‑76 (MQ=255)
                                                                         ctaatATGTTCTAAACTTGATTATCGTATTCAGGGTTATGCAGGTAGTCGGGGCTTGATATATACGAGATATGCCg  >  1:1549718/1‑76 (MQ=255)
                                                                          |                                                                          
AGTTTTGACAAAAATATGTTGTTATAAAAATCTGCTACCACAAGGTGCTCCATCATCACCTAAATTAGCTAATCTAATATGTTCTAAACTTGATTATCGTATTCAGGGTTATGCAGGTAGTCGGGGCTTGATATATACGAGATATGCCG  >  REL606/891418‑891566

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.