New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	31970 =	25 (0.690)	3 (0.090)	3/250	3.9	12.1%	coding (1154/3222 nt)	carB	carbamoyl phosphate synthetase subunit beta
	?	NC_000913	= 170093	21 (0.640)					coding (316/798 nt)	fhuC	iron(III) hydroxamate ABC transporter ATP binding subunit
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CCGCCGTCGGAAGGGGAGATTCTTCTTGATGCCCAACCGCTGGAAAGCTGGAGCAGCAAAGCGTTTGCCCGCAAAGTGGCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGA  >  NC_000913/169961‑170232                                                                                                                                     |                                                                                                                                            ccgccgTCGGAAGGGGAGATTCTTCTTGATGCCCAACCGCTGGAAAGCTGGAGCAGCAAAGCGTTTGCCCGCAAAGTGGCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGAt                                                                                                                                      <  2:537292/140‑1 (MQ=255) ccgccgTCGGAAGGGGAGATTCTTCTTGATGCCCAACCGCTGGAAAGCTGGAGCAGCAAAGCGTTTGCCCGCAAAGTGGCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGAt                                                                                                                                      <  2:524500/140‑1 (MQ=255)                 agaTTCTTCTTGATGCCCAACCGCTGGAAAGCTGGAGCAGCAAAGCGTTTGCCCGCAAAGTGGCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTgg                                                                                                                      >  1:466827/1‑140 (MQ=255)                          ttGATGCCAAACCGCTGGAAAGCTGGAGCAGCAAAGCGTTTGCCCGCAAAGTGGCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGgcgcg                                                                                                             <  2:331994/140‑1 (MQ=255)                            gATGCCCAACCGCTGGAAAGCTGGAGCAGCAAAGCGTTTGCCCGCAAAGTGGCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCt                                                                                                           >  2:481812/1‑140 (MQ=255)                               gCCCAACCGCTGGAAAGCTGGAGCAGCAAAGCGTTTGCCCGCAAAGTGGCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGgcg                                                                                                               >  2:586823/1‑133 (MQ=255)                               gCCCAACCGCTGGAAAGCTGGAGCAGCAAAGCGTTTGCCCGCAAAGTGGCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGgcg                                                                                                               <  1:586823/133‑1 (MQ=255)                                cccAACCGCTGGAAAGCTGGAGCAGCAAAGCGTTTGCCCGCAAAGTGGCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCa                                                                                                                    <  2:313940/127‑1 (MQ=255)                                cccAACCGCTGGAAAGCTGGAGCAGCAAAGCGTTTGCCCGCAAAGTGGCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCa                                                                                                                    >  1:313940/1‑127 (MQ=255)                                                     gcagcaAAGCGTTTGCCCGCAAAGTGGCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATcgc                                                                                  >  1:500445/1‑140 (MQ=255)                                                              cGTTTGCCCGCAAAGTGGCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTc                                                                         <  2:368345/140‑1 (MQ=255)                                                                      cGCAAAGTGGCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGc                                                                 <  1:496340/140‑1 (MQ=255)                                                                                     ttGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTCGCTGGTTgg                                                  >  1:479294/1‑140 (MQ=255)                                                                                       gCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGCGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTCGCTGGTTggat                                                >  1:178065/1‑138 (MQ=255)                                                                                         cGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTCGCTGGTTGGCTTa                                              >  2:405658/1‑140 (MQ=255)                                                                                                      tccGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGgcgc                                 >  2:437009/1‑140 (MQ=255)                                                                                                              aaGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTg                         >  2:521927/1‑140 (MQ=255)                                                                                                                     gACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATa                  >  1:106919/1‑140 (MQ=255)                                                                                                                     gACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATa                  <  2:466827/140‑1 (MQ=255)                                                                                                                           gcgtgAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGtctct            <  1:332282/140‑1 (MQ=255)                                                                                                                                    ggtggCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGgcggcg   <  1:405658/140‑1 (MQ=255)                                                                                                                                     gtggCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGa  <  2:500445/140‑1 (MQ=255)                                                                                                                                     gtggCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGa  >  1:140628/1‑140 (MQ=255)                                                                                                                                     |                                                                                                                                            CCGCCGTCGGAAGGGGAGATTCTTCTTGATGCCCAACCGCTGGAAAGCTGGAGCAGCAAAGCGTTTGCCCGCAAAGTGGCTTATTTGCCGCAGCAGCTTCCTCCGGCAGAAGGGATGACCGTGCGTGAACTGGTGGCGATTGGTCGTTACCCGTGGCATGGCGCGCTGGGGCGCTTTGGGGCGGCAGATCGCGAAAAAGTCGAGGAAGCTATCTCGCTGGTTGGCTTAAAACCGCTGGCGCATCGGCTGGTCGATAGTCTCTCTGGCGGCGA  >  NC_000913/169961‑170232

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