New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 337520657 (0.960)3 (0.060) 3/146 3.7 5.4% coding (1105/1368 nt) degQ serine endoprotease DegQ
?NZ_CP009273 = 3375218 54 (1.020)coding (1117/1368 nt) degQ serine endoprotease DegQ
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGG  >  NZ_CP009273/3375129‑3375300
                                                                                         |                                                                                  
ttGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGc                                                                                    >  2:95315/1‑90 (MQ=255)
 tGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGcc                                                                                   <  1:1403398/90‑1 (MQ=255)
  ggCCGGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGcct                                                                                  >  2:435700/1‑89 (MQ=255)
   gCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAg                                                                                 >  2:1770877/1‑90 (MQ=255)
       gctgcGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCt                                                                             >  2:22367/1‑90 (MQ=255)
       gctgcGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCt                                                                             >  2:1614265/1‑90 (MQ=255)
           cGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGaa                                                                         <  1:973149/90‑1 (MQ=255)
            gTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAg                                                                        >  2:1063617/1‑90 (MQ=255)
                 ggCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCa                                                                   >  1:259854/1‑90 (MQ=255)
                    aaaCCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACg                                                                >  2:901343/1‑90 (MQ=255)
                    aaaCCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACg                                                                >  2:526829/1‑90 (MQ=255)
                     aaCCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGt                                                               <  2:532933/90‑1 (MQ=255)
                           tGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCg                                                         >  2:415568/1‑90 (MQ=255)
                            ggAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGa                                                        <  1:1164246/90‑1 (MQ=255)
                            ggAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGa                                                        >  2:572338/1‑90 (MQ=255)
                                 tAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTc                                                   <  2:629638/90‑1 (MQ=255)
                                 tAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTc                                                   >  1:51697/1‑90 (MQ=255)
                                 tAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTc                                                   >  2:1076353/1‑90 (MQ=255)
                                 tAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTc                                                   <  1:572338/90‑1 (MQ=255)
                                 tAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTc                                                   >  2:872324/1‑90 (MQ=255)
                                 tAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAACGATGGTc                                                   >  1:802493/1‑90 (MQ=255)
                                  aGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCa                                                  <  1:537534/90‑1 (MQ=255)
                                   gAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAg                                                 >  1:599119/1‑90 (MQ=255)
                                        gACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCACCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTaaa                                            >  2:1637601/1‑90 (MQ=255)
                                          cGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAg                                          >  1:1531326/1‑90 (MQ=255)
                                             tCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATg                                       <  1:95315/90‑1 (MQ=255)
                                             tCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATg                                       <  1:901343/90‑1 (MQ=255)
                                             tCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATg                                       <  1:526829/90‑1 (MQ=255)
                                              cGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATgg                                      <  2:112010/90‑1 (MQ=255)
                                              cGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATgg                                      <  2:590674/90‑1 (MQ=255)
                                                aTACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGt                                                    <  1:280523/74‑1 (MQ=255)
                                                aTACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGt                                                    >  2:280523/1‑74 (MQ=255)
                                                 tACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGATAAAAGATggcgg                                   <  1:1754713/90‑1 (MQ=255)
                                                   ccAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCa                                 <  1:67818/90‑1 (MQ=255)
                                                    cAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCaa                                <  1:1770877/90‑1 (MQ=255)
                                                     aGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCaaa                               >  2:1251402/1‑90 (MQ=255)
                                                     aGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCaaa                               >  2:1081008/1‑90 (MQ=255)
                                                      gCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGGGCGATGGTCAGCTAAAAGATGGCGGCAAAg                              <  2:40186/90‑1 (MQ=255)
                                                      gCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAg                              >  1:1419396/1‑90 (MQ=255)
                                                        aCCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGt                            >  2:1613794/1‑90 (MQ=255)
                                                          ctctTCGTCGGCCAGCGCTGAAGTGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTAt                          >  1:1342920/1‑90 (MQ=255)
                                                          ctctTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTAt                          >  2:1415258/1‑90 (MQ=255)
                                                           tctTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTAtt                         >  2:1319444/1‑90 (MQ=255)
                                                           tctTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTAtt                         >  1:996909/1‑90 (MQ=255)
                                                            ctTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAg                              <  1:1111087/84‑1 (MQ=255)
                                                            ctTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAg                              >  2:1111087/1‑84 (MQ=255)
                                                                   ggCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAATGTATTAAAATCGa                 >  2:1380690/1‑90 (MQ=255)
                                                                        gcgcTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAg            >  2:1186750/1‑90 (MQ=255)
                                                                        gcgcTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAg            >  2:474480/1‑90 (MQ=255)
                                                                         cgcTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTc                                                   >  1:98287/1‑50 (MQ=255)
                                                                         cgcTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTc                                                   <  2:98287/50‑1 (MQ=255)
                                                                         cgcTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGt           <  2:1531326/90‑1 (MQ=255)
                                                                         cgcTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGt           <  2:1282043/90‑1 (MQ=255)
                                                                             gAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTc       >  2:864809/1‑90 (MQ=255)
                                                                                  gATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAg                              <  2:1267642/62‑1 (MQ=255)
                                                                                  gATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAg                              >  1:1267642/1‑62 (MQ=255)
                                                                                  gATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAgg  >  2:116689/1‑90 (MQ=255)
                                                                                  gATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAgg  <  2:952878/90‑1 (MQ=255)
                                                                                       cGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTaaaa                     <  2:1781514/66‑1 (MQ=255)
                                                                                       cGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTaaaa                     >  1:1781514/1‑66 (MQ=255)
                                                                                         |                                                                                  
TTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGG  >  NZ_CP009273/3375129‑3375300

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.