New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 3375198 =24 (1.200)3 (0.170) 3/148 1.6 12.9% coding (1097/1368 nt) degQ serine endoprotease DegQ
?NZ_CP009273 3375228 = 19 (1.070)coding (1127/1368 nt) degQ serine endoprotease DegQ
Rejected: Frequency below/above cutoff threshold.

CGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCA  >  NZ_CP009273/3375140‑3375315
                                                                                        |                                                                                       
tcgattttaatacctttgccgccatcttttagctgaccatcgctcaacgttgcaccttccagcgctggcgtgatcatttcagcgctggaa                                                                                        <  1:361566‑M2/2‑1 (MQ=255)
cGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGaa                                                                                        <  2:425366/90‑1 (MQ=255)
       gCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCaa                                                                                 <  1:588500/90‑1 (MQ=255)
                     ccatcttttagctgaccatcgctcaacgttgcaccttccagcgctggcgtgatcatttcagcgctggAAGGTGCAACGTTGAGCGATGGt                                                                   <  2:510240‑M2/23‑1 (MQ=255)
                      catcttttagctgaccatcgctcaacgttgcaccttccagcgctggcgtgatcatttcagcgctggAAGGTGCAACGTTGAGCGATGGTc                                                                  <  1:118170‑M2/24‑1 (MQ=255)
                        gAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAg                                                                >  2:309874/1‑90 (MQ=255)
                         aaGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGc                                                               <  1:238357/90‑1 (MQ=255)
                           gTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTa                                                             >  1:278089/1‑90 (MQ=255)
                               cGCTCGATACCAGCCCCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGt                                                                   <  2:564854/80‑1 (MQ=255)
                               cGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGt                                                                   >  1:564854/1‑80 (MQ=255)
                                  gaccatcgctcaacgttgcaccttCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATg                                                      <  1:118170‑M1/66‑1 (MQ=255)
                                   accatcgctcaacgttgcaccttCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATgg                                                     <  2:510240‑M1/67‑1 (MQ=255)
                                           gCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCa                                                                 >  1:332174/1‑70 (MQ=255)
                                           gCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCa                                                                 <  2:332174/70‑1 (MQ=255)
                                           gCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCaaa                                              >  2:549803/1‑89 (MQ=255)
                                           gCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCaaa                                              <  1:549803/89‑1 (MQ=255)
                                           gCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAg                                             <  1:64033/90‑1 (MQ=255)
                                           gCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAg                                             >  1:48034/1‑90 (MQ=255)
                                           gCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAg                                             >  1:375360/1‑90 (MQ=255)
                                               ctctTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTAt                                         >  1:486515/1‑90 (MQ=255)
                                                    cgtcgGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTaaaa                                    <  2:464363/90‑1 (MQ=255)
                                                        ttCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGa                                <  1:361566‑M1/88‑1 (MQ=255)
                                                             GCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGt                                                                   >  1:583367‑M1/1‑50 (MQ=255)
                                                             GCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGt                                                                   <  2:583367‑M1/50‑1 (MQ=255)
                                                              CGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGt                          >  2:278089‑M1/1‑90 (MQ=255)
                                                                  GAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTc                                                                  >  1:469767‑M1/1‑46 (MQ=255)
                                                                  GAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTc                                                                  <  2:469767‑M1/46‑1 (MQ=255)
                                                                       gATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAgg                 >  2:442461/1‑90 (MQ=255)
                                                                            cGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAATAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGcc            >  1:328284/1‑90 (MQ=255)
                                                                              ccAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCa          <  2:118170/90‑1 (MQ=255)
                                                                                agcgctggAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAgc        <  1:242336‑M2/82‑1 (MQ=255)
                                                                                 gcgctggAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATg                                                      >  2:147423‑M2/8‑43 (MQ=255)
                                                                                 gcgctggAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATg                                                      <  1:147423‑M2/36‑1 (MQ=255)
                                                                                 gcgctggAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAgct       >  2:394605‑M2/8‑90 (MQ=255)
                                                                                      ggAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCa  >  2:7977‑M2/3‑90 (MQ=255)
                                                                                        |                                                                                       
CGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCA  >  NZ_CP009273/3375140‑3375315

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