New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	4377070 =	24 (0.730)	3 (0.090)	3/588	4.1	12.4%	coding (196/318 nt)	sugE	multidrug efflux system protein
	?	NC_000913	4377134 =	19 (0.590)					coding (260/318 nt)	sugE	multidrug efflux system protein
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below cutoff threshold.

ATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCTTTAGTCGTTTGACGCCGAGTGTTATTACTGTGACGGCGATGATTGTCAGTATGGCGCTACTTGCCTGGGCGATGAAATCGTTACCAGTAGGGACGGCTTATGCCGTGTGGACGGGTATTGGCGCAGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCACTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGGAGAGTATCCACAGGTAGTCGCGGTCCGGGCCGCAAACCAGCGCATGGCGGTATTCCCGATCGAGTGCA  >  NC_000913/4376770‑4377366                                                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                                         aTTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCTTTAGTCGTTTGACGCCGAGTGTTATTACTGTGACGGCGATGATTGTCAGTATGGCGCTACTTGCCTGGGCGATGAAATCGTTACCAGTAGGGACGGCTTATGCCGTGTGGACGGGTATTg                                                                                                                                                                                                                                                                                                          <  1:107275/301‑1 (MQ=255)    ccccCGCCACCGGGCAAAAAAGCCGTTCACCCATCAGCTTGGAAGGGCGACCTGCACACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCTTTAGTCGTTTGACGCCGAGTGTTATTACTGTGACGGCGATGATTGTCAGTATGGCGCTACTTGCCTGGGCGATGAAATCGTTACCAGTAGGGACGGCTTATGCCGTGTGGACGGGTATTGgcg                                                                                                                                                                                                                                                                                                       <  2:266193/301‑1 (MQ=255)                       tatCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGTGCCGTTGGCCTGAAATATACCCACGGCTTTAGTCGTTTGACGCCGAGTGTTATTACTGTGACGGCGATGATTGTCAGTATGGCGCTACTTGCCTGGGCTATGAAATCGTTACCAGTAGGGACGGCTTATGCCGTGTGGATGGGTATTGGCGCAGTCGGCGCGGCCACaa                                                                                                                                                                                                                                                                                     >  1:261263/1‑300 (MQ=255)                           cGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGTCTGAAATATACCCACGGCTTTAGTCGTTTGACGCCGAGTGTTATTAATGTGACGGCGAGGATTGTCAGTATGGCGCTACTTGCCTGGGCGAGGAAAACGTTACATGGAGGGACGGCGTATGTCGTGGGGACGGGGATTGGCGCGGGCGGCGGGGcagttacgggg                                                                                                                                                                                                                                                                                >  2:166168/1‑291 (MQ=255)                               cAACCAGCCGCTTTGCAGGCGGACCTGCTACGGCGTCTTTTCGCCGGGGCCGGTCCCAATGTGTCGGTGCCTAAATTGTACTGGATTAACTTAGTTATTGCTGGTCTGCGGGAAGTGGTATGGGCCGTTGGCCTGAAGTATACCGACGGCTTTGGTCGTTTGACGCCGGTGGTTATTACTGTGACGGGGATGATTGTCAGTATGGCGCTACTTGCCTGGGCGATGAAATCGTTACCAGTAGGGACGGCTTATGCCGTGTGGACGGGTATTGGCGCAGTCGGCGCGGCCATAACCGGCATtg                                                                                                                                                                                                                                                                            <  2:108868/301‑1 (MQ=255)                                  ccATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCTTTAGTCGTTTGACGCCGAGTGTTATTACTGTGACGGCGATGATTGTCAGTATGGCGCTACTTGCCTGGGCGATGAAATCGTTACCAGTAGGGACGGCTTATGCCGTGTGGACGGGTATTGGCGCAGTCGGCGCGGCCATAACCGGCATTGtgc                                                                                                                                                                                                                                                                         <  1:272641/301‑1 (MQ=255)                                                        ctcacaccccTCTCTGCTCCGGGGACTGCCCCAATTCTCCGTATCCTTTTTTGTCCTGGTTTATCTTTGTTATTGCTGGTGTGCTGGAAGTTCGTTGGGCCGTTGGCATGAAATATACCCACGTCTTTAGTCGTTTTACTCCGAGGGTTATTGCTGTGACGGCGATGATTGTCAGTATGGCGCTACTTGCCTGGGCGATGAAATCGTTACCAGTAGGGACGGCTTATGCTGTGTGGACGGGTATTGGCGCAGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAAccc                                                                                                                                                                                                                                                   <  2:167690/293‑1 (MQ=255)                                                          cAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCTTTAGTCGTTTGACGCCGAGTGTTATTACTGTGACGGCGATGATTGTCAGTATGGCGCTACTTGCCTGGGCGATGAAATCGTTACCAGTAGGGACGGCTTATGCCGTGTGGACGGGTATTGGCGAAGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGGGCGTCCGCTataccgg                                                                                                                                                                                                                                                 >  1:154645/1‑295 (MQ=255)                                                                     tgtCACCGAGGGCGGCCCGATGTCTCGGGGGCCTGTTNGGTCTGGGATTTGCTTAGTTATGGCTGGTCTGCTGGAAGTGGTATGGGTCGTTGGCCTGAAATATACCCGCGGCTTTAGTCGTTTGACGCCGAGTGTTATTACTGTGACGGCGATGATTGTCAGTATGGCGCTACTTGCCTGGGCGATGAAATCGTTACCAGTAGGGACGGCTTATGCCGTGTGGACGGGTATTGGCGCAGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCg                                                                                                                                                                                                                                      <  2:46082/299‑1 (MQ=255)                                                                                                   gCCTTATATGTCCTGGTTTATCTTAGTTATTGTTGGTCTTCTGGAAGTGGTTTGGGCCGTTGGCCTGAAATATACCCACGGCTTTAGTCGTTTGACGCCTAGTGTTATTACTGTGACGGCGATGATTGTCAGTATGGCGCTACTTGCCTGGGCGATGAAATCGTTACCAGTAGGGACGGCTTATGCCGTGTGGACGGGTATTGGCGCAGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGattat                                                                                                                                                                                                         <  1:166178/300‑1 (MQ=255)                                                                                                           tGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCTTTAGTCGTTTGACGCCGAGTGTTATTACTGTGACGGCGATGATTGTCAGTATGGCGCTACTTGCCTGGGCGATGAAATCGTTACCAGTAGGGACGGCTTATGCCGTGTGGACGGGTATTGGCGCAGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGGGGCCGCTCACCCGGTGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTTGttggaa                                                                                                                                                                                                >  2:41523/1‑296 (MQ=255)                                                                                                                                   ttGGTCTGCTGAAAGTGGTATGGGCCGTTTGCCTGAAATATACCCACGGCTGTAGTCGTTTGACGCCTGGTGTTATTACTATGACGGCGATGATTGTCAGTATGGCGCTACTTGCCTGGGCGATGAAATCGTTACCAGTAGGGACGGCTTATGCCGTGTGGACGGGTATTGGCGCAGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCACTAACTACCAg                                                                                                                                                                         <  1:261975/299‑1 (MQ=255)                                                                                                                                   gtGGTCTGGTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGTTTTAGTCGTTTGACGCCGAGTGTTATTACTGTGACGGCGATGATTGTCAGTATGGCGCTACTTGCCTGGGCGATGAAATCGTTACCAGTAGGGACGGCTTATGCCGTGTGGACGGGTATTGGCGCAGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCACTAACTACCAg                                                                                                                                                                         <  1:41526/299‑1 (MQ=255)                                                                                                                                      gTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCTTTAGTCGTTTGACGCCGAGTGTTATTACTGTGACGGCGATGATTGTCAGTATGGCGCTACTTGCCTGGGCGATGAAATCGTTACCAGTAGGGACGGCTTATGCCGTGTGGACGGGTATTGGCGCGGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCACTAACTACCAGgctt                                                                                                                                                                     >  1:273470/1‑300 (MQ=255)                                                                                                                                                             gTTGGCCTGAAATATACCCACGGCTTTAGTCGTTTGACGCCGAGTGTTATTACTGTGACGGCGATGATTGTCAGTATGGCGCTACTTGCCTGGGCGATGAAATCGTTACCAGTAGGGACGGCTTATGCCGTGTGGACGGGTATTGGCGCAGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGGGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCAACCACTAACTACCAGGCTGGTGTACCCCAATAAATTTACTGa                                                                                                                                              >  2:33515/1‑301 (MQ=255)                                                                                                                                                             gTTGGCCCGAAATATACCCACGGCTTTAGTCGTTTGACGCCGAGTGTTATTACTGTGACGGCGATGATTGTCAGTATGGCGCTACTTGCCTGGGCGATGAAATCGTTACCAGTAGGGACGGCTTATGCCGTGTGGACGGGTATTGGCGCAGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCACTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGa                                                                                                                                              >  1:190256/1‑301 (MQ=255)                                                                                                                                                                                                  cGCCGAGTGTTATTACTGTGACGGCGATGATTGTCAGTATGGCGCTACTTGCCTGGGCGATGAAATCGTTACCAGTAGGGACGGCTTATGCCGTGTGGACGGGTATTGGCGCAGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCACTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCAtct                                                                                                          >  1:132814/1‑300 (MQ=255)                                                                                                                                                                                                                    tggggtcGATGATTGGCAGGATTGCTCTGCTTGCCTGGGCGATTAAATGTTTACCAGTGGGGACGGGTTATGCCGTGTGGACTGGTATGGGGGCATTCGGCGCGGCCATGACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCACTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAg                                                                                        <  2:132805/294‑1 (MQ=255)                                                                                                                                                                                                                                      gTATGGCGCTACTTGCCTGGGCGATGAAATCGTTACCAGTAGGGTCCGCTTATGCCGTGTGGACGGGTATTGGCGCAGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCTGCTGCGCCTGGCGAGTCTGGCGTTACTTGTATTGGGGATTATTGGTCTGATACTCAGCACTCACTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAACCCCTTCCCGTGTCGCGACTTCCAGCATCTCCTGTTTCACTTCGTCAGAAATTGTTGGCGTGCgggg                                                                     >  1:271421/1‑300 (MQ=255)                                                                                                                                                                                                                                              cgaCTGGCCTGGGCGAGGAAGTCGGTACCAGTGGGGACGGCGTATGGCGTGTGGACGGGTATTGGCGCAGTCGGCGAGGCCAGAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCACTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGGAGAGTATc                                                              <  2:13112/298‑1 (MQ=255)                                                                                                                                                                                                                                                                                               tgtgGACGGGTATTGGCGCAGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCACTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGGAGAGTATCCACAGGTAGTCGCGGTCAGGGCCGCAAACCAGCGCATGGCGGTATTCCCg            >  1:43715/1‑301 (MQ=255)                                                                                                                                                                                                                                                                                               tgagggcGGGTATTGGCGCAGTCGGCGCGGGCATAACCGGCATTGTGCGGGTCGGGGAGGCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCACGAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGGAGAGTATCCACAGGTAGTCGCGGTCCGGGCCGCAAACCAGCGCATGGCGGTATTccc             <  1:59509/294‑1 (MQ=255)                                                                                                                                                                                                                                                                                                  gaaCGGGTATGGGCGCCATGTGCGCGGCCATAAGCGGCATTAGGCTGCTCGTGGAGTGCGCTAACCCGATGCGCCTGGGGAGTCTGTCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCACTAACTACCAGGCTGCTGTACCCAGATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGGGGAGTATCCACAGGTAGTCGCGGTCCGGGCCGCAAACCAGCGCATGGCGGTATTCCCGAt          <  1:237764/298‑1 (MQ=255)                                                                                                                                                                                                                                                                                                         tATTGGCGCAGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCACTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGCCTTCCAGCATCCCCTTTTTCACTTCGTAAGCAATGGTTGGCGTGCGGACACGTATCCACCGTTCGTCGCGGTCGGGGCCGACAACCCGCCCATGTCGGTTTTCCCgtgcgctggca  >  2:181086/1‑291 (MQ=255)                                                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                                         ATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATGGGCCGTTGGCCTGAAATATACCCACGGCTTTAGTCGTTTGACGCCGAGTGTTATTACTGTGACGGCGATGATTGTCAGTATGGCGCTACTTGCCTGGGCGATGAAATCGTTACCAGTAGGGACGGCTTATGCCGTGTGGACGGGTATTGGCGCAGTCGGCGCGGCCATAACCGGCATTGTGCTGCTCGGTGAGTCCGCTAACCCGATGCGCCTGGCGAGTCTGGCGTTAATCGTATTGGGGATTATTGGTCTGAAACTCAGCACTCACTAACTACCAGGCTGCTGTACCCAAATAAATTTACTGACATCAAACCCTTCCCGGGTCGCGACTGCCAGCATCTCCTGTTTCACTTCGTCAGAAATGGTTGGCGTGCGGGAGAGTATCCACAGGTAGTCGCGGTCCGGGCCGCAAACCAGCGCATGGCGGTATTCCCGATCGAGTGCA  >  NC_000913/4376770‑4377366

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