New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 71350	402 (1.060)	136 (0.370)	59/188	10.9	41.1%	intergenic (+85/‑1)	araC/yabI	DNA‑binding transcriptional dual regulator AraC/DedA family protein YabI
	?	NC_000913	= 2857042	0 (0.000)					intergenic (‑236/‑51)	ygbA/mutS	protein YgbA/DNA mismatch repair protein MutS
Rejected: Coverage evenness skew score above cutoff.

ACAGCATCTTTCCCGGAACCAGCATCAAGAACTCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAA  >  NC_000913/2856942‑2857050                                                                                                     |         aCAGCATCTTTCCCGGAACCAGCATCAAGAACTCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGaa          <  1:5205358/101‑1 (MQ=255) aCAGCATCTTTCCCGGAACCAGCATCAAGAACTCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGaa          >  1:5205770/1‑101 (MQ=255)  cAGCATCTTTCCCGGAACCAGCATCAAGAACTCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGaaa         <  1:5205266/101‑1 (MQ=255)  cAGCATCTTTCCCGGAACCAGCATCAAGAACTCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGaaa         <  1:5205303/101‑1 (MQ=255)  cAGCATCTTTCCCGGAACCAGCATCAAGAACTCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGaaa         >  2:5205771/1‑101 (MQ=255)  cAGCATCTTTCCCGGAACCAGCATCAAGAACTCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGaaa         <  1:5205428/101‑1 (MQ=255)   aGCATCTTTCCCGGAACCAGCATCAAGAACTCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGaaaa        >  1:5205772/1‑101 (MQ=255)   aGCATCTTTCCCGGAACCAGCATCAAGAACTCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGaaaa        >  1:5205773/1‑101 (MQ=255)   aGCATCTTTCCCGGAACCAGCATCAAGAACTCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGaaaa        >  1:5205774/1‑101 (MQ=255)   aGCATCTTTCCCGGAACCAGCATCAAGAACTCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGaaaa        >  1:5205775/1‑101 (MQ=255)    gCATCTTTCCCGGAACCAGCATCAAGAACTCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAAt       <  1:5205441/101‑1 (MQ=255)    gCATCTTTCCCGGAACCAGCATCAAGAACTCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAAt       <  1:5205552/101‑1 (MQ=255)     cATCTTTCCCGGAACCAGCATCAAGAACTCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGaaa         >  1:5205778/1‑98 (MQ=255)     cATCTTTCCCGGAACCAGCATCAAGAACTCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAAtg      >  2:5205776/1‑100 (MQ=255)     cATCTTTCCCGGAACCAGCATCAAGAACTCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAAtg      >  1:5205777/1‑100 (MQ=255)      aTCTTTCCCGGAACCAGCGTCAAGAACTCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAAtgt     <  2:5204950/101‑3 (MQ=255)      aTCTTTCCCGGAACCAGCATCAAGAACTCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAAtgt     >  1:5205779/1‑99 (MQ=255)      aTCTTTCCCGGAACCAGCATCAAGAACTCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAAtgt     <  2:5205564/101‑3 (MQ=255)       tCTTTCCCGGAACCAGCATCAAGAACTCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAAtgtg    >  1:5205780/1‑98 (MQ=255)        cTTTCCCGGAACCAGCATCAAGAACTCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAAtgtgc   >  2:5205782/1‑97 (MQ=255)        cTTTCCCGGAACCAGCATCAAGAACTCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAAtgtgc   <  1:5205218/101‑5 (MQ=255)        cTTTCCCGGAACCAGCATCAAGAACTCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAAtgtgc   >  1:5205781/1‑97 (MQ=255)        cTTTCCCGGAACCAGCATCAAGAACTCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAAtgtgc   <  1:5205348/101‑5 (MQ=255)        cTTTCCCGGAACCAGCATCAAGAACTCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAAtgtgc   <  1:5205317/101‑5 (MQ=255)         tttCCCGGAACCAGCATCAAGAACTCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGaaa         <  2:5205273/94‑1 (MQ=255)         tttCCCGGAACCAGCATCAAGAACTCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAAtgtgcc  >  1:5205785/1‑96 (MQ=255)         tttCCCGGAACCAGCATCAAGAACTCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAAtgtgcc  <  2:5204911/101‑6 (MQ=255)         tttCCCGGAACCAGCATCAAGAACTCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAAtgtgcc  >  1:5205783/1‑96 (MQ=255)         tttCCCGGAACCAGCATCAAGAACTCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAAtgtgcc  <  1:5204881/101‑6 (MQ=255)         tttCCCGGAACCAGCATCAAGAACTCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAAtgtgcc  >  2:5205784/1‑96 (MQ=255)                                cTCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGaaa         <  2:5205155/71‑1 (MQ=255)                                  cGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGaaaa        >  1:5205824/1‑70 (MQ=255)                                                                                                     |         ACAGCATCTTTCCCGGAACCAGCATCAAGAACTCGCCGTTCGCTTCTTCCCCTGAAATGATTAACTCCGGTATCATGTGCGCCTTATGTGATTACAACGAAAATAAAAA  >  NC_000913/2856942‑2857050

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