New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 457159079 (1.260)3 (0.050) 3/158 3.9 3.7% coding (1276/1590 nt) hsdM type I restriction‑modification system methyltransferase
?NZ_CP009273 4571594 = 81 (1.330)coding (1272/1590 nt) hsdM type I restriction‑modification system methyltransferase
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACG  >  NZ_CP009273/4571505‑4571667
                                                                                     |                                                                             
cACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCg                                                                           <  1:645186/90‑1 (MQ=255)
 aCTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGGTGGTCGGTGGTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCgg                                                                          >  1:1676884/1‑90 (MQ=255)
 aCTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCgg                                                                          >  2:1113258/1‑90 (MQ=255)
   tCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCg                                                                        >  2:1301077/1‑90 (MQ=255)
     aCGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGtt                                                                      <  2:55803/90‑1 (MQ=255)
      cGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTa                                                                     >  2:958716/1‑90 (MQ=255)
        gCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTaaa                                                                   <  1:1252894/90‑1 (MQ=255)
        gCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTaaa                                                                   <  1:1750435/90‑1 (MQ=255)
         cTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTaaaa                                                                  <  1:1113257/90‑1 (MQ=255)
          tGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAAc                                                                 <  1:1901934/90‑1 (MQ=255)
           gAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGACAACTTCCGTCTCTTCGGCGTTAAAACt                                                                >  2:1016021/1‑90 (MQ=255)
            aaCTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTc                                                               >  1:1899133/1‑90 (MQ=255)
             aCTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTcc                                                              <  2:295044/90‑1 (MQ=255)
               ttGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCAt                                                            <  1:1653524/90‑1 (MQ=255)
               ttGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCAt                                                            <  1:1777295/90‑1 (MQ=255)
                tGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCAtt                                                           >  2:1082600/1‑90 (MQ=255)
                tGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCAtt                                                           >  2:617503/1‑90 (MQ=255)
                tGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCAtt                                                           >  1:716972/1‑90 (MQ=255)
                tGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCAtt                                                           >  2:603625/1‑90 (MQ=255)
                tGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCAtt                                                           >  1:1693158/1‑90 (MQ=255)
                tGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCAtt                                                           >  2:1687220/1‑90 (MQ=255)
                 gcgcCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTc                                                          <  2:1287115/90‑1 (MQ=255)
                 gcgcCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTc                                                          <  1:1692380/90‑1 (MQ=255)
                 gcgcCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTc                                                          <  2:670423/90‑1 (MQ=255)
                 gcgcCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTc                                                          <  1:422167/90‑1 (MQ=255)
                    ccAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCAcc                                                       <  2:627868/90‑1 (MQ=255)
                     cAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCAt                                                            <  1:823464/84‑1 (MQ=255)
                     cAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCAt                                                            >  2:823465/1‑84 (MQ=255)
                       gCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTc                                                    >  1:1067192/1‑90 (MQ=255)
                              gTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTgcgc                                             <  2:404217/90‑1 (MQ=255)
                              gTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTgcgc                                             <  2:1693159/90‑1 (MQ=255)
                              gTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTgcgc                                             <  2:1501214/90‑1 (MQ=255)
                              gTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTgcgc                                             <  2:216472/90‑1 (MQ=255)
                              gTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTgcgc                                             <  2:1301636/90‑1 (MQ=255)
                              gTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTgcgc                                             <  1:1082599/90‑1 (MQ=255)
                              gTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTgcgc                                             <  1:445241/90‑1 (MQ=255)
                              gTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTgcgc                                             <  1:1908608/90‑1 (MQ=255)
                              gTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTgcgc                                             <  1:53285/90‑1 (MQ=255)
                              gTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTgcgc                                             <  2:83772/90‑1 (MQ=255)
                              gTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTgcgc                                             <  1:617503/90‑1 (MQ=255)
                              gTAGCAAGATGCTGGTCGATGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCTCCTTCAGTgcgc                                             <  1:605309/90‑1 (MQ=255)
                              gTAACAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTgcgc                                             <  1:1824755/90‑1 (MQ=255)
                                aGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCgg                                           <  2:1163318/90‑1 (MQ=255)
                                   aaGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCt                                        >  2:1541444/1‑90 (MQ=255)
                                    aGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTc                                                          <  2:1023925/71‑1 (MQ=255)
                                    aGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTc                                                          >  1:1023924/1‑71 (MQ=255)
                                    aGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCG‑‑TCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTaa                                     >  1:855057/1‑90 (MQ=255)
                                       tGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTgcgc                                             >  1:1507207/1‑81 (MQ=255)
                                       tGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTgcgc                                             <  2:1507208/81‑1 (MQ=255)
                                        gCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTgcgc                                             <  1:138388/80‑1 (MQ=255)
                                        gCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTgcgc                                             <  2:220327/80‑1 (MQ=255)
                                        gCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTgcgc                                             >  1:220327/1‑80 (MQ=255)
                                        gCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTgcgc                                             >  2:138388/1‑80 (MQ=255)
                                        gCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTAAAc                                   >  1:842419/1‑90 (MQ=255)
                                                gtgtTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGcgcg                                            >  1:1472851/1‑73 (MQ=255)
                                                gtgtTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGcgcg                                            <  2:1472852/73‑1 (MQ=255)
                                                gtgtTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTC‑‑‑TCGGCGTTAAAACTcc                                                              >  1:1454874/1‑52 (MQ=255)
                                                gtgtTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTC‑‑‑TCGGCGTTAAAACTcc                                                              <  2:1454875/52‑1 (MQ=255)
                                                  gtTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTAAACCGTGCGGGTc                         <  2:1492719/90‑1 (MQ=255)
                                                     tttGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCTCGGGCTTAAACCGTGCGGGTCTTc                      <  1:301698/90‑1 (MQ=255)
                                                     tttGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTAAACCGTGCGGGTCTTc                      >  2:1565271/1‑90 (MQ=255)
                                                     tttGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTAAACCGTGCGGGTCTTc                      <  1:1565270/90‑1 (MQ=255)
                                                          tctcTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTgcgc                                             <  1:1048815/62‑1 (MQ=255)
                                                          tctcTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTgcgc                                             >  2:1048816/1‑62 (MQ=255)
                                                          tctcTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTAAACCGTGCGGGTc                         >  1:519469/1‑82 (MQ=255)
                                                          tctcTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTAAACCGTGCGGGTc                         <  2:519469/82‑1 (MQ=255)
                                                          tctcTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTAAACCGTGCGGGTCTTCGCCAt                 >  1:1864880/1‑90 (MQ=255)
                                                          tctcTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTAAACCGTGCGGGTCTTCGCCAt                 >  2:817458/1‑90 (MQ=255)
                                                          tctcTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTAAACCGTGCGGGTCTTCGCCAt                 <  2:716972/90‑1 (MQ=255)
                                                          tctcTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTAAACCGTGCGGGTCTTCGCCAt                 >  1:575084/1‑90 (MQ=255)
                                                            tctTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTgcgc                                             >  2:1197420/1‑60 (MQ=255)
                                                            tctTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTgcgc                                             <  1:1197419/60‑1 (MQ=255)
                                                            tctTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTAAACCGTTCGGGTCTTCGCCATac               >  2:1488488/1‑90 (MQ=255)
                                                             ctTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTgcgc                                             <  1:917715/59‑1 (MQ=255)
                                                             ctTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTgcgc                                             >  2:469758/1‑59 (MQ=255)
                                                             ctTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTgcgc                                             <  1:469758/59‑1 (MQ=255)
                                                             ctTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTgcgc                                             >  2:917716/1‑59 (MQ=255)
                                                              ttCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCCCCTTCAGTGCGCGGGCTTAAACCGTGCGGGTCTTCGCCATacac             <  1:236743/90‑1 (MQ=255)
                                                              ttCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTAAACCGTGCGGGTCTTCGCCATacac             <  2:1085473/90‑1 (MQ=255)
                                                              ttCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTAAACCGTGCGGGTCTTCGCCATacac             <  2:250453/90‑1 (MQ=255)
                                                              ttCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTAAACCGTGCGGGTCTTCGCCATacac             <  2:1899134/90‑1 (MQ=255)
                                                                         cAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACg  <  1:581075/90‑1 (MQ=255)
                                                                                     |                                                                             
CACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACG  >  NZ_CP009273/4571505‑4571667

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