New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	2617476 =	27 (1.090)	4 (0.170)	4/258	2.7	13.4%	coding (1383/1464 nt)	bepA	OM protein maintenance and assembly metalloprotease and chaperone, periplasmic
	?	NC_000913	2617529 =	26 (1.090)					coding (1436/1464 nt)	bepA	OM protein maintenance and assembly metalloprotease and chaperone, periplasmic
Rejected: Frequency below/above cutoff threshold.

GATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAGCTGGCTGCGCGAGCAGAAGGTTATGCGCTCGCCGGACGACTCGATCAGGCCATTTCGCTGTTGAGTAGCGCCAGTTCGCAGGTGAAATTAGGCAGCCTGCAACAAGCGCGTTACGATGCGCGCATCGACCAGTTGCGCCAGCTGCAGGAACGCTTTAAGCCTTATACCAAGATGTAATTCGCATTAAAGGAAGAACAATGACCAAACAGGTAAAAATCTACCATAACCCACGCTGTTCAAAGAG  >  NC_000913/2617333‑2617624                                                                                                                                                |                                                                                                                                                     gATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAGCTGGCTGCGCGAGCAGAAGGTTATGCGCTCGCCGGACGACTCGATCAGGCCATTTCGCTGTTGAGTAGCGCCAGTTCGCAGGTGAAATTAGGCAGCCTGc                                                                                                                                                  >  6:119819/1‑148 (MQ=255)  aTTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAGCTGGCTGCGCGAGCAGAAGGTTATGCGCTCGCCGGACGACTCGATCAGGCCATTTCGCTGCTTAGTAGCGCCAGTTCGCAGGTGAAATTAGGCAGCCTGcaa                                                                                                                                                >  4:61978/1‑149 (MQ=255)            accggCGGAAGCCGCGCTAAATAACCGCGATCAGGAGCTGGCTGCGCGAGCAGAAGGTTATGCGCTCGCCGGACGCCTCGATCAGGCCCGTTCGCTGTGGAGTCGCGCCAGTTCGCAGGTGAAATTAGGCCGCCTGCAACAAGCGCGt                                                                                                                                       >  5:113857/4‑148 (MQ=255)               ggCGGAAGCCGCGCTAAATAACCGCGATCAGGAGCTGGCTGCGCGAGCAGAAGGTTATGCGCTCGCCGGACGACTCGATCAGGCCATTTCGCTGTTGAGTAGCGCCAGTTCGCAGGTGAAATTAGGCAGCCTGCAACAAGCGCGTTAc                                                                                                                                    >  5:69346/1‑148 (MQ=255)                            cTAAATAACCGCGATCAGGAGCTGGCTGCGCGAGCAGAAGGTTATGCGCTCGCCGGACGACTCGATCAGGCCATTTCGCTGTTGAGTAGCGCCAGTTCGCAGGTGAAATTAGGCAGCCTGCAACAAGCGCGTTACGATGCGCGCATCGa                                                                                                                      <  6:69346/149‑1 (MQ=255)                               aaTAACCGCGATCAGGAGCTGGCTGCGCGAGCAGAAGGTTATGCGCTCGCCGGACGACTCGATCAGGCCATTTCGCTGTTGAGTAGCGCCAGTTCGCAGGTGAAATTAGGCAGCCTGCAACAAGCGCGTTACGATGCGCGCATCGAc                                                                                                                     <  5:119819/147‑1 (MQ=255)                                         aTCAGGAGCTGGCTGCGCGAGCAGAAGGTTATGCGCTCGCCGGACGACTCGATCAGGCCATTTCGCTGTTGAGTAGCGCCAGTTCGCAGGTGAAATTAGGCAGCCTGCAACAAGCGCGTTACGATGCGCGCATCGACCAGTTGCGCCAg                                                                                                         >  4:123535/1‑149 (MQ=255)                                          tCAGGAGCTGGCTGCGCGAGCAGAAGGTTATGCGCTCGCCGGACGACTCGATCAGGCCATTTCGCTGGTGAGTAGCGCCAGTTCGCAGGTGAAATTAGGCCGCCTGCAACAAGCGCGTTACGATGCGCGCATCGACCAGTTGCGCCAGc                                                                                                        >  5:142393/1‑149 (MQ=255)                                            ccggctctttgaacagcgtgggttatggtagatttttacctgtttggtcattgttcttcctttaatgcgaattacatcttggtataaggcttaaagcgttCCTGCAACAAGCGCGTTACGATGCGCGCATCGACCAGTTGCGCCAGCTg                                                                                                      <  3:116685‑M1/49‑1 (MQ=255)                                                       aacagcgtgggttatggtagatttttacctgtttggtcattgttcttcctttaatgcgaattacatcttggtataaggcttaaagcgttCCTGCAACAAGCGCGTTACGATGCGCGCATCGACCAGTTGCGCCAGCTGCAGGAACGCtt                                                                                           >  4:116685‑M1/90‑149 (MQ=255)                                                        cgcgAGCAGAAGGTTATGCGCTCGCCGGACGACTCGATCAGGCCATTTCGCTGTTGAGTAGCGCCAGTTCGCAGGTGAAATTAGGCAGCCTGCAACAAGCGCGTTACGATGCGCGCATCGACCAGTTGCGCCAGCTGCAGGAACGCttt                                                                                          <  3:123535/149‑1 (MQ=255)                                                         gcgAGCAGAAGGTTATGCGCTCGCCGGACGACTCGATCAGGCCATTTCGCTGTTGAGTAGCGCCAGTTCGCAGGTGAAATTAGGCAGCCTGCAACAAgcgc                                                                                                                                         >  6:126966/1‑101 (MQ=255)                                                         gcgAGCAGAAGGTTATGCGCTCGCCGGACGACTCGATCAGGCCATTTCGCTGTTGAGTAGCGCCAGTTCGCAGGTGAAATTAGGCAGCCTGCAACAAgcgc                                                                                                                                         <  5:126966/101‑1 (MQ=255)                                                          cgAGCAGAAGGTTATGCGCTCGCCGGACGACTCGATCAGGCCATTTCGCTGTTGAGTAGCGCCAGTTCGCAGGTGAAATTAGGCAGCCTGCAACAAGCGCGTTACGATGCGCGCATCGACCAGTTGCGCCAGCTGCAGGAACGCTTTaa                                                                                        <  2:101697/149‑1 (MQ=255)                                                                   gggtaggcgcTCGCCGGACGACGCGATCAGGCCCATTCGCTGGTGAGTAGCGCCAGTTCGCAGGCGAAATTAGGCAGCCTGCAACAAGCGCGTTACGATGCGCGCATCGACCAGTTGCGCCAGCTGCAGGAACGCTTTAAGCCTTATAc                                                                               <  6:142393/143‑1 (MQ=255)                                                                     ttATGCGCTCGCCGGACGACTCGATCAGGCCATTTCGCTGTTGAGTAGCGCCAGTTCGCAGGTGAAATTAGGCAGCGTGCAACAAGCGCGTTACGATGCGCGCATCGACCAGTTGCGCCAGCTGCAGGAACGCTTTAAGCCTTATAcc                                                                              >  5:28257/1‑148 (MQ=255)                                                                          gatttttacctgtttggtcattgttcttcctttaatgcgaattacatcttggtataaggcttaaagcgttCCTGCAACAAGCGCGTTACGATGCGCGCATCGACCAGTTGCGCCAGCTGCAGGAAc                                                                                               >  2:123423‑M1/71‑126 (MQ=255)                                                                          gatttttacctgtttggtcattgttcttcctttaatgcgaattacatcttggtataaggcttaaagcgttCCTGCAACAAGCGCGTTACGATGCGCGCATCGACCAGTTGCGCCAGCTGCAGGAAc                                                                                               <  1:123423‑M1/56‑1 (MQ=255)                                                                                   gacgacTCGATCAGGCCATTTCGCTGTTGAGTAGCGCCAGTTCGCAGGTGAAATTAGGCAGCCTGCAACAAGCGCGTTACGATGCGCGCATCGACCAGTTGCGCCAGCTGCa                                                                                                    <  6:42558/112‑1 (MQ=255)                                                                                            aTCAGGCCATTTCGCTGTTGAGTAGCGCCAGTTCGCAGGTGAAATTAGGCAGCCTGCAACAAGCGCGTTACGATGCGCGCATCGACCAGTTGCGCCAGCTGCAGGAACGCTTTAAGCCTTATACCAAGATGTAATTCGCATTAAAgg                                                        >  5:63873/1‑147 (MQ=255)                                                                                                ggCCATTTCGCTGTTGAGTAGCGCCAGTTCGCAGGTGAAATTAGGCAGCCTGCAACAAGCGCGTTACGATGCGCGCATCGACCAGTTGCGCCAGCTGCAGGAACGCTTTAAGCCTTATACCAAGATGTAATTCGCATTAAAGgaagaa                                                   >  2:57823/1‑148 (MQ=255)                                                                                                             ctGAGTAGCTCCAGTTCTCAGTTGAAATTAGGCAGCCTGCAACAAGCTAGTTACTATGCGCGCATCGACCAGTTGCGCAAGCTGCAGGAACGCTTTAAGCCTTATACCAAGATGTAATTCGCATTAAAGGAAGAACAATGACCAAACa                                      <  3:61978/147‑1 (MQ=255)                                                                                                                      gcCAGTTCGCAGGTGAAATTAGGCAGCCTGCAACAAGCGCGTTACGATGCGCGCATCGACCAGTTGCGCCAGCTGCAGGAACGCTTTAAGCCTTATACCAAGATGTAATTCGCATTAAAGGAAGAACAATGACCAAACAGGTAAAAATc                            >  3:135647/1‑149 (MQ=255)                                                                                                                        cAGTTCGCAGGTGAAATTAGGCAGCGTGCAACAAGCGCGTTACGATGCGCGCATCGACCAGTTGCGCCAGCTGCAGGAACGCTTTAAGCCTTATACCAAGATGTAATTCGCATTAAAGGAAGAACAATGACCAAACAGGTAAAAATCTa                          <  6:28257/149‑1 (MQ=255)                                                                                                                         aGTTCGCAGGTGAAATTAGGCAGCCTGCAACAAGCGCGTTACGATGCGCGCATCGACCAGTTGCGCCAGCTGCAGGAACGCTTTAAGCCTTATACCAAGATGTAATTCGCATTAAAGGAAGAACAATGACCAAACAGGTAAAAATCTAc                         <  1:57823/149‑1 (MQ=255)                                                                                                                         aGTTCGCAGGTGAAATTAGGCAGCCTGCAACAAGCGCGTTACGATGCGCGCATCGACCAGTTGCGCCAGCTGCAGGAACGCTTTAAGCCTTATACCAAGATGTAATTCGCATTAAAGGAAGAACAATGACCAAACAGGTAAAAATCTAc                         <  6:63873/149‑1 (MQ=255)                                                                                                                          gTTCGCAGGTGAAATTAGGCAGCCTGCAACAAGCGCGTTACGATGCGCGCATCGACCAGTTGCGCCAGCTGCAGGAACGCTTTaa                                                                                        >  4:47211/1‑85 (MQ=255)                                                                                                                          gTTCGCAGGTGAAATTAGGCAGCCTGCAACAAGCGCGTTACGATGCGCGCATCGACCAGTTGCGCCAGCTGCAGGAACGCTTTaa                                                                                        <  3:47211/85‑1 (MQ=255)                                                                                                                                aGGTGAAATTAGGCAGCCTGCAACAAGCGCGTTACGATGCGCGCATCGACCAGTTGCGCCAGCTGCAGGAACGCTTTAAGCCTTATACCAAGATGTAATTCGCATTAAAGGAAGAACAATGACCAAACAGGTAAAAATCTACCATAAcc                  >  3:28931/1‑149 (MQ=255)                                                                                                                                 ggTGAAATTAGGCAGCCTGCAACAAGCGCGTTACGATGCGCGCATCGACCAGTTGCGCCAGCTGCAGGAACGCTTTAAGCCTTATACCAAGATGTAATTCGCATTAAAGGAAGAACAATGACCAAACAGGTAAAAATCTACCATAAccc                 <  4:28931/149‑1 (MQ=255)                                                                                                                                             cAGCCTGCAACAAGCGCGTTACGATGCGCGCATCGACCAGTTGCGCCAGCTGCAGGAACGCTTTAAGCCTTATACCAAGATGTAATTCGCATTAAAGGAAGAACAATGACCAAACAGGTAAAAATCTACCATAACCCACGCTGTTCa       >  3:48715/1‑147 (MQ=255)                                                                                                                                                CCTGCAACAAGCGCGTTACGATGCGCGCATCGACCAGTTGCGCCAGCTGCAGAAACGCTTTAAGCCTTATACCAAGATGTAATTCGCATTAAAGGAAGAACAATGACCAAACAGGTAAAAATCTACCATAACCCACGCTGTTCAAagag  <  2:67919‑M1/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                     GATTTGCTGGCACAGGCGGAAGCCGCGCTAAATAACCGCGATCAGGAGCTGGCTGCGCGAGCAGAAGGTTATGCGCTCGCCGGACGACTCGATCAGGCCATTTCGCTGTTGAGTAGCGCCAGTTCGCAGGTGAAATTAGGCAGCCTGCAACAAGCGCGTTACGATGCGCGCATCGACCAGTTGCGCCAGCTGCAGGAACGCTTTAAGCCTTATACCAAGATGTAATTCGCATTAAAGGAAGAACAATGACCAAACAGGTAAAAATCTACCATAACCCACGCTGTTCAAAGAG  >  NC_000913/2617333‑2617624

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