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.

CCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGA  >  NZ_CP009273/3375110‑3375285                                                                                         |                                                                                        tgagcagctgggcttcctttgacaacttcatcgattttaatacctttgccgccatcttttagctgaccatcgctcaacgttgcaccttcc                                                                                        >  2:7977‑M1/89‑90 (MQ=255)   gggCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAg                                                                                      <  2:226364/90‑1 (MQ=255)      cACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAgcgc                                                                                   >  2:5100/1‑90 (MQ=255)      agctgggcttcctttgacaacttcatcgattttaatacctttgccgccatcttttagctgaccatcgctcaacgttgcaccttCCAgcgc                                                                                   >  2:394605‑M1/84‑90 (MQ=255)       gctgggcttcctttgacaacttcatcgattttaatacctttgccgccatcttttagctgaccatcgctcaacgttgcaccttCCAGCGCt                                                                                  <  1:242336‑M1/8‑1 (MQ=255)               gAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGAt                                                                          >  1:266121/1‑90 (MQ=255)                     tGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGcc                                                                    <  2:425298/90‑1 (MQ=255)                               cGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTggcc                                                          <  2:521858/90‑3 (MQ=255)                               cGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGaa                                                          <  2:425366/90‑1 (MQ=255)                                 tAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTggccga                                                        <  1:5100/90‑5 (MQ=255)                                      gCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCaa                                                   <  1:588500/90‑1 (MQ=255)                                                     tAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGcc                                                                    <  1:103452/58‑1 (MQ=255)                                                     tAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGcc                                                                    >  2:103452/1‑58 (MQ=255)                                                     catcttttagctgaccatcgctcaacgttgcaccttCCAgcgc                                                                                   <  1:147423‑M1/7‑1 (MQ=255)                                                     catcttttagctgaccatcgctcaacgttgcaccttCCAgcgc                                                                                   >  2:147423‑M1/37‑43 (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                                   <  2:332174/70‑1 (MQ=255)                                                                          gCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCa                                   >  1:332174/1‑70 (MQ=255)                                                                          gCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCaaa                <  1:549803/89‑1 (MQ=255)                                                                          gCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCaaa                >  2:549803/1‑89 (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)                                                                                         |                                                                                        CCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGA  >  NZ_CP009273/3375110‑3375285

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.