New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 385156552 (1.780)5 (0.200) 5/144 1.4 10.8% coding (200/1494 nt) yidJ sulfatase‑like hydrolase/transferase
?NZ_CP009273 = 3851588 39 (1.580)coding (177/1494 nt) yidJ sulfatase‑like hydrolase/transferase
Rejected: Frequency below/above cutoff threshold.

TAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTT  >  NZ_CP009273/3851500‑3851675
                                                                                        |                                                                                       
tAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATa                                                                                        >  2:353152‑M2/1‑89 (MQ=255)
   aGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATagtc                                                                                     >  1:727451‑M2/1‑86 (MQ=255)
      aGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGGCTGGTTAGCGTAGATACCGGTAAATagtccgg                                                                                  >  1:257261‑M2/1‑83 (MQ=255)
      aGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATagtccgg                                                                                  <  2:39639‑M2/90‑8 (MQ=255)
       gcgcCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATagtccggc                                                                                 <  1:102378‑M2/90‑9 (MQ=255)
       gcgcCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATagtccggc                                                                                 <  2:450751‑M2/90‑9 (MQ=255)
       gcgcCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATagtccggc                                                                                 <  2:246628‑M2/90‑9 (MQ=255)
       gcgcCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATagtccggc                                                                                 <  2:4805‑M2/90‑9 (MQ=255)
       gcgcCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGCACCACGGGCCGGACTGGTTAGCGCATATACCGGCAAATagtccggc                                                                                 <  1:367924‑M2/90‑9 (MQ=255)
            ccATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATagtccggcccgtg                                                                            <  1:31745‑M2/90‑14 (MQ=255)
             cATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGc                                                                           <  1:830269/90‑1 (MQ=255)
             cATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGc                                                                           <  1:486764/90‑1 (MQ=255)
                aGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAgg                                                                        <  2:637190/90‑1 (MQ=255)
                 gTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAg                                                                         <  1:657569/88‑1 (MQ=255)
                 gTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAg                                                                         >  2:657569/1‑88 (MQ=255)
                   ggAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATagtccggcccgtggaccaac                                                                     >  1:127464‑M2/1‑70 (MQ=255)
                    gagaTGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTa                                                                    >  2:477652/1‑90 (MQ=255)
                      gaTGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACa                                                                  >  2:663439/1‑90 (MQ=255)
                      gaTGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACa                                                                  >  1:741939/1‑90 (MQ=255)
                          tttttGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAAcc                                                              <  2:364817/90‑1 (MQ=255)
                           ttttGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCg                                                             >  2:49663/1‑90 (MQ=255)
                           ttttGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCg                                                             >  2:262789/1‑90 (MQ=255)
                            tttGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCgg                                                            <  1:639146/90‑1 (MQ=255)
                              tGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATagtccggcccgtggaccaacaacgtcgcgcc                                                          >  1:241725‑M2/1‑59 (MQ=255)
                               gCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATagtccggcccgtggaccaacaacgtcgcgcca                                                         >  1:520185‑M2/1‑58 (MQ=255)
                                ccTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGaa                                                        <  2:577214/90‑1 (MQ=255)
                                   ggCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGc                                                                           <  1:220802/68‑1 (MQ=255)
                                   ggCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGc                                                                           >  2:220802/1‑68 (MQ=255)
                                   ggCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAg                                                     >  2:788520/1‑90 (MQ=255)
                                   ggCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAg                                                     >  2:332334/1‑90 (MQ=255)
                                   ggCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAg                                                     <  1:159142/90‑1 (MQ=255)
                                   ggCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAg                                                     >  2:174637/1‑90 (MQ=255)
                                   ggCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAg                                                     >  1:272579/1‑90 (MQ=255)
                                       cgACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGgtgt                                                 <  1:515654/90‑1 (MQ=255)
                                        gACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTa                                                >  2:464839/1‑90 (MQ=255)
                                         aCGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAg                                               >  1:419586/1‑90 (MQ=255)
                                           gttgttGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGc                                             <  1:246628/90‑1 (MQ=255)
                                                        gggCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGa                                >  1:208993/1‑90 (MQ=255)
                                                          gCCGGACTGGTTAGCGTAGATACCGGTAAATagtccggcccgtggaccaacaacgtcgcgccaggcaaaaacatctccactatggggcgc                              >  1:385080‑M2/1‑31 (MQ=255)
                                                           ccGGACTGGTTAGCGTAGATACCGGTAAATagtccggcccgtggaccaacaacgtcgcgccaggcaaaaacatctccactatggggcgct                             >  2:672015‑M2/1‑30 (MQ=255)
                                                                 tGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTcc                       <  1:12732/90‑1 (MQ=255)
                                                                       gCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCgccgcc                 <  2:281081/90‑1 (MQ=255)
                                                                       gCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCgccgcc                 <  1:647132/90‑1 (MQ=255)
                                                                          tAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGa              <  2:208993/90‑1 (MQ=255)
                                                                           aGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGAc             <  1:153374/90‑1 (MQ=255)
                                                                            gATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACt            >  2:178786/1‑90 (MQ=255)
                                                                              tACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTAt          <  2:83703/90‑1 (MQ=255)
                                                                                 cGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAAc                                                       <  1:66976/42‑1 (MQ=255)
                                                                                 cGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAAc                                                       >  2:66976/1‑42 (MQ=255)
                                                                                 cGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCaa       >  2:220511/1‑90 (MQ=255)
                                                                                 cGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCaa       >  1:136109/1‑90 (MQ=255)
                                                                                  ggTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAat      >  2:283926/1‑90 (MQ=255)
                                                                                    tAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAatat    <  1:262789/90‑1 (MQ=255)
                                                                                      aaTAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATAttt  >  2:75333/1‑90 (MQ=255)
                                                                                        |                                                                                       
TAAAGTAGCGCCCCATAGTGGAGATGTTTTTGCCTGGCGCGACGTTGTTGGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTT  >  NZ_CP009273/3851500‑3851675

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