New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 570046 =36 (0.790)4 (0.090) 3/142 6.1 7.3% intergenic (‑110/+569) yhhI/ECB_00526 predicted transposase/hypothetical protein
?REL606 570481 = 67 (1.550)intergenic (‑545/+134) yhhI/ECB_00526 predicted transposase/hypothetical protein
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

AATAAAAAACAGCTTCAGGTAATAAGGAATATCTCAATTTTTAAACATAAAATGCGAATTATTTAGTACAAAAAGTGAACTTCTGATGATCTTTTCCTGGTTAATATCCAGGCATCATGTTACGTGACATGGTGTAGCTA  >  REL606/569982‑570121
                                                                |                                                                           
aaTAAAAAACAGCTTCAGGTAATAAGGAATATCTCAATTTTTAAACATAAAATGCGAATTATTTAGTACAAAAAGt                                                                  <  1:889967/76‑1 (MQ=31)
aaTAAAAAACAGCTTCAGGTAATAAGGAATATCTCAATTTTTAAACATAAAATGCGAATTATTTAGTACAAAAAGt                                                                  <  1:889955/76‑1 (MQ=31)
aaTAAAAAACAGCTTCAGGTAATAAGGAATATCTCAATTTTTAAACATAAAATGCGAATTATTTAGTACAAAAAGt                                                                  >  1:300996/1‑76 (MQ=31)
  tAAAAAACAGCTTCAGGTAATAAGGAATATCTCAATTTTTAAACATAAAATGCGAATTATTTAGTACAAAAAGTGa                                                                <  2:1100343/76‑1 (MQ=34)
     aaatCAGCTTCAGGTAATAAGGAATATCTCAATTTTTAGACATAAAATGCGAATTATTTAGTACAAAACGTGCACt                                                             >  2:1445534/1‑76 (MQ=11)
      aaaCAGCTTCAGGTAATAAGGAATATCTCAATTTTTAAACATAAAATGCGAATTATTTAGTACAAAAAGTGAACtt                                                            <  1:243906/76‑1 (MQ=34)
          aGCTTCAGGTAATAAGGAATATCTCAATTTTTAAACATAAAATGCGAATTATTTAGTACAAAAAGTGAACTTCtga                                                        >  1:1345347/1‑76 (MQ=34)
             ttCAGGTAATAAGGAATATCTCAATTTTTAAACATAAAATGCGAATTATTTAGTACAAAAAGTGAACTTCtgatga                                                     <  2:402351/76‑1 (MQ=34)
              tCAGGTAATAAGGAATATCTCAATTTTTAAACATAAAATGCGAATTATTTAGTACAAAAAGTGAACTTCTgatgat                                                    >  2:815564/1‑76 (MQ=34)
              tCAGGTAATAAGGAATATCTCAATTTTTAAACATAAAATGCGAATTATTTAGTACAAAAAGTGAACTTCTgatgat                                                    <  2:1323387/76‑1 (MQ=34)
              tCAGGTAATAAGGAATATCTCAATTTTTAAACATAAAATGCGAATTATTTAGTACAAAAAGTGAACTTCTgatgat                                                    <  2:515534/76‑1 (MQ=34)
               cAGGTAATACGGAATATCTCAATTTTTAAACATAAAATGCGAATTATTTAGTACAAAAAGTGAACTTCTGATGATc                                                   >  1:349035/1‑76 (MQ=21)
               cAGGTAATAAGGAATATCTCAATTTTTAAACATAAAATGCGAATTATTTAGTACAAAAAGTGAACTTCTGATGATc                                                   >  2:1006407/1‑76 (MQ=34)
               cAGGTAATAAGGAATATCTCAATTTTTAAACATAAAATGCGAATTATTTAGTACAAAAAGTGAACTTCTGATGATc                                                   >  1:349038/1‑76 (MQ=34)
                aGGTAATAAGGAATATCTCAATTTTTAAACATAAAATGCGAATTATTTAGTACAAAAAGTGAACTTCTGATGATCt                                                  >  1:994540/1‑76 (MQ=34)
                   taataaGGAATATCTCAATTTTTAAACATAAAATGCGAATTATTTAGTACAAAAAGTGAACTTCTGATGATCtttt                                               >  2:332347/1‑76 (MQ=35)
                   taataaGGAATATCTCAATTTTTAAACATAAAATGCGAATTATTTAGTACAAAAAGTGAACTTCTGATGATCtttt                                               >  2:332346/1‑76 (MQ=35)
                    aataaGGAATATCTCAATTTTTAAACATAAAATGCGAATTATTTAGTACAAAAAGTGAACTTCTGATGATCTTTTc                                              >  2:1296048/1‑76 (MQ=38)
                            atatCTCAATTTTTAAACATAAAATGCGAATTATTTAGTACAAAAAGTGAACTTCTGATGATCTTTTCCTGGTTaa                                      >  2:1262408/1‑76 (MQ=40)
                            atatCTCAATTTTTAAACATAAAATGCGAATTATTTAGTACAAAAAGTGAACTTCTGATGATCTTTTCCTGGTTaa                                      >  2:1351195/1‑76 (MQ=40)
                                 tcAATTTTTAAACATAAAATGCGAATTATTTAGTACAAAAAGTGAACTTCTGATGATCTTTTCCTGGTTAATATcc                                 <  2:863340/76‑1 (MQ=255)
                                  cAATTTTTAAACATAAAATGCGAATTATTTAGTACAAAAAGTGAACTTCTGATGATCTTTTCCTGGTNAATATCCa                                <  1:1417953/76‑1 (MQ=255)
                                        ttAAACATAAAATGCGAATGATTTAGTACAAAAAGTGAACTTCTGATGATCTTTTCCTGGTTAATATCCAGGcatc                          >  1:1258061/1‑76 (MQ=255)
                                         tAAACATAAAATGCGAATTATTTAGTACAAAAAGTGAACTTCTGATGATCTTTTCCTGGTTAATATCCAGGcatca                         <  1:193955/76‑1 (MQ=255)
                                          aaaCATAAAATGCGAATTATTTAGTACAAAAAGTGAACTTCTGATGATCTTTTCCTGGTTAATATCCAGGcatcat                        <  2:597998/76‑1 (MQ=255)
                                             cATAAAATGCGAATTATTTAGTACAAAAAGTGAACTTCTGATGATCTTTTCCTGGTTAATATCCAGGCATCATGtt                     <  2:1333686/76‑1 (MQ=255)
                                              aTAAAATGCGAATTATTTAGTACAAAAAGTGAACTTCTGATGATCTTTTCCTGGTTAATATCCAGGCATCATGTTa                    >  1:397141/1‑76 (MQ=255)
                                                    tGCGAATTATTTAGTACAAAAAGTGAACTTCTGACGATCTTNTCCTGGTTAATATCCAGGNATNATGTTACGTGAc              >  1:1202934/1‑76 (MQ=255)
                                                        aaTTATTTAGTACAAAAAGTGAACTTCTGATGATCTTTTCCTGGTTAATATCCAGGCATCATGTTACGTGACATgg          >  2:1107123/1‑76 (MQ=255)
                                                         attattTAGTACAAAAAGTGAACTTCTGATGATCTTTTCCTGGTTAATATCCAGGCATCATGTTACGTGACATGgt         <  2:1312688/76‑1 (MQ=255)
                                                         attattTAGTACAAAAAGTGAACTTCTGATGATCTTTTCCTGGTTAATATCCAGGCATCATGTTACGTGACATGgt         <  1:1006409/76‑1 (MQ=255)
                                                           tattTAGTACAAAAAGTGAACTTCTGATGATCTTTTCCTGGTTAATATCCAGGCATCATGTTACGTGACATGgtgt       >  2:1395394/1‑76 (MQ=255)
                                                           tattTAGTACAAAAAGTGAACTTCTGATGATCTTTTCCTGGTTAATATCCAGGCATCATGTTACGTGACATGgtgt       >  1:808290/1‑76 (MQ=255)
                                                           tattTAGTACAAAAAGTGAACTTCTGATGATCTTTTCCTGGTTAATATCCAGGCATCATGTTACGTGACATGgtgt       >  2:725289/1‑76 (MQ=255)
                                                           tattTAGTACAAAAAGTGAACTTCGGATNTTCTTTTCCTTGTTAATATACNGGCATCATGTTACGTGACATGgtgt       >  1:221896/1‑76 (MQ=25)
                                                            attTAGTACAAAAAGTGAACTTCTGATGATCTTTTCCTGGTTAATATCCAGGCATCATGTTACGTGACATGGTGTa      <  2:1371974/76‑1 (MQ=255)
                                                                aGTACAAAAAGTGAACTTCTGATGATCTTTTCCTGGTTAATATCCAGGCATCATGTTACGTGACATGGTGTAGCTa  <  2:881607/76‑1 (MQ=255)
                                                                aGTACAAAAAGTGAACTTCNGNTGATCTTNTCCTGGTTAATATCCAGGCATNATNTTCCGTGACATGGTGTAGCTa  <  2:214840/76‑1 (MQ=37)
                                                                |                                                                           
AATAAAAAACAGCTTCAGGTAATAAGGAATATCTCAATTTTTAAACATAAAATGCGAATTATTTAGTACAAAAAGTGAACTTCTGATGATCTTTTCCTGGTTAATATCCAGGCATCATGTTACGTGACATGGTGTAGCTA  >  REL606/569982‑570121

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 34 ≤ 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.