New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3381839 =	28 (0.810)	5 (0.150)	5/278	2.4	16.1%	coding (1097/1368 nt)	degQ	serine endoprotease, periplasmic
	?	NC_000913	3381869 =	26 (0.800)					coding (1127/1368 nt)	degQ	serine endoprotease, periplasmic
Rejected: Frequency below cutoff threshold.

CCTCAACGGCAAACCGCTGAATAGCTTTGCTGAGTTGCGCTCTCGTATCGCGACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATG  >  NC_000913/3381690‑3381979                                                                                                                                                      |                                                                                                                                             ccTCAACGGCACACCGCTGATTAGCTTGGCTGAGTTGCGCACTCGTATCGCGACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGcc                                                                                                                                             <  1:78752/151‑1 (MQ=255)        ttcacccgatcgcggttgacgccaatgatcacatcgtctttttgcaagccagcctgagcagctgggcttcctttgacaacttcatcgattttaatacctttgccgccatcttttagctgaccatcgctcaacgttgcaccttCCAgcgc                                                                                                                                        >  1:136542‑M1/143‑149 (MQ=255)           aaaCCGCTGAATAGCTTTGCTGAGTTGCGCTCTCGTATCGCGACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGaaa                                                                                                                                   <  2:55252/151‑1 (MQ=255)                   gAATAGCTTTGCTGAGTTGCGCTCTCGTATCGCGACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACg                                                                                                                           >  1:453296/1‑151 (MQ=255)                                        ctctCGTATCGCGACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACg                                                                                                      >  2:204396/1‑151 (MQ=255)                                         tctcGTATCGCGACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGt                                                                                                     <  1:240998/151‑1 (MQ=255)                                         tctcGTATCGCGACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGt                                                                                                     <  1:437683/151‑1 (MQ=255)                                          ctcGTATCGCGACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACTACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGtt                                                                                                    >  2:180843/1‑151 (MQ=255)                                          ctcGTATCGCGACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGtt                                                                                                    >  2:485090/1‑151 (MQ=255)                                          ctcGTATCGCGACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGtt                                                                                                    >  2:143580/1‑151 (MQ=255)                                              tttttgcaagccagcctgagcagctgggcttcctttgacaacttcatcgattttaatacctttgccgccatcttttagctgaccatcgctcaacgttgcaccttCCAGCGCTGGCGTGATCATTTCAGCGCTGGAAGGTGCAACGTTGAGc                                                                                                >  2:451218‑M1/105‑151 (MQ=255)                                              tATCGCGACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGc                                                                                                <  1:132473/151‑1 (MQ=255)                                                tttgcaagccagcctgagcagctgggcttcctttgacaacttcatcgattttaatacctttgccgccatcttttagctgaccatcgctcaacgttgcaccttCCAGCGCTGGCGTGATCATTTCAGCGCTGGAAGGTGCAACGTTGAGCGa                                                                                              >  1:314337‑M1/103‑151 (MQ=255)                                                tCGCGACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGa                                                                                              <  1:193455/151‑1 (MQ=255)                                                 cgcgACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGAt                                                                                             <  1:508887/151‑1 (MQ=255)                                                           gAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTaa                                                                                   >  1:16430/1‑151 (MQ=255)                                                                       aaaGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCa                                                                       >  1:129463/1‑151 (MQ=255)                                                                          gTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAg                                                                    >  2:162101/1‑151 (MQ=255)                                                                           tGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAgg                                                                   <  1:447105/151‑1 (MQ=255)                                                                           tGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAgg                                                                   <  1:298354/151‑1 (MQ=255)                                                                                  tGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTaaa                                                            <  1:467602/151‑1 (MQ=255)                                                                                    caacttcatcgattttaatacctttgccgccatcttttagctgaccatcgctcaacgttgcaccttCCAGCGCTGGCGTGATCATTTCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAAt                                                          >  2:300067‑M1/67‑151 (MQ=255)                                                                                           gcGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGaa                                                   >  1:15992/1‑151 (MQ=255)                                                                                                        ccACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGaa                                      >  2:31639/1‑151 (MQ=255)                                                                                                               aaGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAg                                >  2:71224/1‑150 (MQ=255)                                                                                                                    tcttttagcttaccatcgctcaacgttgcaccttCCAGCGCTGGCGTGATCATTTCAGCGCTGGAAGGTGCAACGTTGAGCGATTGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTc                          >  1:516671‑M1/35‑151 (MQ=255)                                                                                                                     aaGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTc                          <  2:448769/150‑1 (MQ=255)                                                                                                                     aaGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCa                         >  2:46667/1‑151 (MQ=255)                                                                                                                        tGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAgg                       <  2:378521/150‑1 (MQ=255)                                                                                                                               cGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTg               >  1:533912/1‑151 (MQ=255)                                                                                                                                  tACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCaa            >  1:493111/1‑151 (MQ=255)                                                                                                                                        cACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGAc      <  1:180839/151‑1 (MQ=255)                                                                                                                                          acgttgcaccttCCAGCGCTGAAATGCTCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGa    <  1:437327‑M1/139‑1 (MQ=255)                                                                                                                                            tctTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGAtg  <  2:6636/151‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                             CCTCAACGGCAAACCGCTGAATAGCTTTGCTGAGTTGCGCTCTCGTATCGCGACCACCGAGCCGGGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATG  >  NC_000913/3381690‑3381979

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