New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 3022554 =41 (0.820)3 (0.070) 3/136 5.8 7.4% coding (131/1140 nt) ECB_02821 hypothetical protein
?REL606 3022717 = 37 (0.810)intergenic (‑33/+541) ECB_02821/ECB_02822 hypothetical protein/bifunctional Ribulose 5‑phosphate reductase/CDP‑ribitol pyrophosphorylase
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TTTTTTTTAATAAACTTAAACGTTATATTATGTTTATTAAATACCTCGACAGCTAAATGTTCTATTTCCTTATAATATTTACAAAAAACACTAATACAATCCTTGTTATTAAATACTTGTGATTGAAAGTAACGAATA  >  REL606/3022482‑3022619
                                                                        |                                                                 
ttttttttAATAAACTTAAACGTTATATTATGTTTATTAAATACCTCGACAGCTAAATGTTCTATTTCCTTataa                                                                 >  1:699427/1‑75 (MQ=255)
ttttttttAATAAACTTAAACGTTATATTATGTTTATTAAATACCTCGACAGCTAAATGTTCTATTTCCTTataa                                                                 >  1:699422/1‑75 (MQ=255)
ttttttttAATAAACTTAAACGTGATATTATGTTTATTAAATACCTCGACAGCTAAATGTTCTATTTCCTTataa                                                                 >  1:699406/1‑75 (MQ=255)
 tttttttAATAAACTTAAACGTTATATTATGTTTATTAAATACCTCGACAGCTAAATGTTCTATTTCCTTataat                                                                <  1:414523/75‑1 (MQ=255)
    ttttAATAAACTTAAACGTTATATTATGTTTATTAAATACCTCGACAGCTAAATGTTCTATTTCCTTATAATAtt                                                             >  2:198661/1‑75 (MQ=255)
     tttAATAAACTTAAACGTTATATTATGTTTATTAAATACCTCGACAGCTAAATGTTCTATTTCCTTATAATAttt                                                            <  1:1251011/75‑1 (MQ=255)
      ttAATAAACTTAAACGTTATATTATGTTTATTAAATACCTCGACAGCTAAATGTTCTATTTCCTTATAATATTTa                                                           <  2:1602905/75‑1 (MQ=255)
       taataaACTTAAACGTTATATTATGTTTATTAAATACCTCGACAGCTAAATGTTCTATTTCCTTATAATATTTAc                                                          <  1:583718/75‑1 (MQ=255)
            aaCTTCAACGTTATATTATGTTTATTAAATACCTCGACAGCTAAATGTTCTATTTCCTTATAATATTTACaaaaa                                                     >  2:1172632/1‑75 (MQ=255)
              cTTAAACGTTATATTATGTTTATTAAATACCTCGACAGCTAAATGTTCTATTTCCTTATAATATTTACAAAAAac                                                   <  1:1521968/75‑1 (MQ=255)
                tAAACGTTATATTATGTTTATTAAATACCTCGACAGCTAAATGTTCTATTTCCTTATAATATTTACAAAAAacac                                                 >  1:317172/1‑75 (MQ=255)
                 aaaCGTTATATTATGTTTATTAAATACCTCGACAGCTAAATGTTCTATTTCCTTATAATATTTACAAAAAACACt                                                >  1:757648/1‑75 (MQ=255)
                  aaCGTTATATTATGTTTATTAAATACCTCGACAGCTAAATGTTCTATTTCCTTATAATATTTACAAAAAACACTa                                               <  2:202104/75‑1 (MQ=255)
                    cGTTATATTATGTTTATTAAATACCTCGACAGCTAAATGTTCTATTTCCTTATAATATTTACAAAAAACACTAAt                                             <  1:338068/75‑1 (MQ=255)
                      ttATATTATGTTTATTAAATACCTCGACAGCTAAATGTTCTATTTCCTTATAATATTTACAAAAAACACTAATAc                                           >  2:232750/1‑75 (MQ=255)
                       tataTTGTGTTTATTAAATACCTCGACAGCTAAATGTTCTATTTCCTTATAATATTTACAAAAAACACTAATACa                                          <  1:190757/75‑1 (MQ=255)
                       tataTTGTGTTTATTAAATACCTCGACAGCTAAATGTTCTATTTCCTTATAATATTTACAAAAAACACTAATACa                                          <  1:190745/75‑1 (MQ=255)
                        atatTATGTTTATTAAATACCTCGACAGCTAAATGTTCTATTTCCTTATAATATTTACAAAAAACACTAATACaa                                         >  1:1132186/1‑75 (MQ=255)
                           ttatGTTTATTAAATACCTCGACAGCTAAATGTTCTATTTCCTTATAATATTTACAAAAAACACTAATACAATcc                                      <  1:852466/75‑1 (MQ=255)
                            tatGTTTATTAAATACCTCGACAGCTAAATGTTCTATTTCCTTATAATATTTACAAAAAACACTAATACAATCCt                                     <  2:1529692/75‑1 (MQ=255)
                                 ttattaAATACCTCGACAGCTAAATGTTCTATTTCCTTATAATATTTACAAAAAACACTAATACAATCCTTGtta                                >  1:395082/1‑75 (MQ=255)
                                  tattaAATACCTCGACAGCTAAATGTTCTATTTCCTTATAATATTTACAAAAAACACTAATACAATCCTTGttat                               <  1:345355/75‑1 (MQ=255)
                                  tattaAATACCTCGACAGCTAAATGTTCTATTTCCTTATAATATTTACAAAAAACACTAATACAATCCTTGttat                               >  2:1547333/1‑75 (MQ=255)
                                   attaAATACCTCGACAGCTAAATGTTCTATTTCCTTATAATATTTACAAAAAACACTAATACAATCCTTGttatt                              >  1:737748/1‑75 (MQ=255)
                                    ttaAATACCTCGACAGCTAAATGTTCTATTTCCTTATAATATTTACAAAAAACACTAATACAATCCTTGttatta                             >  2:722764/1‑75 (MQ=255)
                                      aaaTACCTCGACAGCTAAATGTTCTATTTCCTTATAATATTTACAAAAAACACTAATACAATCCTTGTTATTaaa                           >  1:92129/1‑75 (MQ=255)
                                         tACCTCGACAGCTAAATGTTCTATTTCCTTATAATATTTACAAAAAACACTAATACAATCCTTGTTATTAAATAc                        >  2:1282788/1‑75 (MQ=255)
                                          aCCTCGACAGCTAAATGTTCTATTTCCTTATAATATTTACAAAAAACACTAATACAATCCTTGTTATTAAATACt                       >  2:1590120/1‑75 (MQ=255)
                                             tctaCAGCTAAATGTTCTATTTCCTTATAATATTTACAAAAAACACTAATACAATCCTTGTTATTAAATACTtgt                    <  2:1101671/72‑1 (MQ=255)
                                               gACAGCTAAATGTTCTATTTCCTTATAATATTTACAAAAAACACTAATACAATCCTTGTTATTAAATACTTGTGa                  <  1:52551/75‑1 (MQ=255)
                                                 cAGCTAAATGTTCTATTTCCTTATAATATTTACAAAAAACACTAATACAATCCTTGTTATTAAATACTTGTGAtt                >  2:353100/1‑75 (MQ=255)
                                                     tAAATGTTCTATTTCCTTATAATATTTACAAAAAACACTAATACAATCCTTGTTATTAAATACTTGTGATTGaaa            >  2:632687/1‑75 (MQ=255)
                                                       aaTGTTCTATTTCCTTATAATATTTACAAAAAACACTAATACAATCCTTGTTATTAAATACTTGTGATTGAAAGt          >  2:1135371/1‑75 (MQ=255)
                                                       aaTGTTCTATTTCCTTATAATATTTACAAAAAACACTAATACAATCCTTGTTATTAAATACTTGTGATTGAAAGt          <  1:1337692/75‑1 (MQ=255)
                                                         tGTTCTATTTCCTTATAATATTTACAAAAAACACTAATACAATCCTTGTTATTAAATACTTGTGATTGAAAGTaa        <  2:1588101/75‑1 (MQ=255)
                                                          gTTCTATTTCCTTATAATATTTGAAAAAAACACTAATACAATCCTTGTTATTAAATACTTGTGATTGAAAGTAAc       <  1:433887/75‑1 (MQ=255)
                                                           ttCTATTTCCTTATAATATTTACAAAAAACACTAATACAATCCTTGTTATTAAATACTTGTGATTGAAAGTAACg      >  1:744365/1‑75 (MQ=255)
                                                           ttCTATTTCCTTATAATATTTACAAAAAACACTAATACAATCCTTGTTATTAAATACTTGTGATTGAAAGTAACg      >  1:1476212/1‑75 (MQ=255)
                                                           ttCTATTTCCTTATAATATTTACAAAAAACACTAATACAATCCTTGTTATTAAATACTTGTGATTGAAAGTAACg      >  1:744367/1‑75 (MQ=255)
                                                              tATTTCCTTATAATATTTACAAAAAACACTAATACAATCCTTGTTATTAAATACTTGTGATTGAAAGTAACGAAt   >  1:282077/1‑75 (MQ=255)
                                                               aTTTCCTTATAATATTTACAAAAAACACTAATACAATCCTTGTTATTAAATACTTGTGATTGAAAGTAACGAATa  >  1:1562955/1‑75 (MQ=255)
                                                                        |                                                                 
TTTTTTTTAATAAACTTAAACGTTATATTATGTTTATTAAATACCTCGACAGCTAAATGTTCTATTTCCTTATAATATTTACAAAAAACACTAATACAATCCTTGTTATTAAATACTTGTGATTGAAAGTAACGAATA  >  REL606/3022482‑3022619

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