New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3955668 =	35 (0.610)	9 (0.170)	6/306	4.5	26.0%	coding (338/1545 nt)	ilvA	l‑threonine dehydratase, biosynthetic; also known as threonine deaminase
	?	NC_000913	3955897 =	18 (0.330)					coding (567/1545 nt)	ilvA	l‑threonine dehydratase, biosynthetic; also known as threonine deaminase
Rejected: Coverage evenness skew score above cutoff.

AAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTG‑CTTCTGCGGGTAACCACGCGCAGGGCGTCGCGTTTTCTTCTGCGCGGTTAGGCGTGAAGGCCCTGATCGTTATGCCAACCGCCACCGCCGACATCAAAGTCGACGCGGTGCGCGGCTTCGGCGGCGAAGTGCTGCTCCACGGCGCGAACTTTGATGAAGCGAAAGCCAAAGCGATCGAACTGTCACAGCAGCAGGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGA  >  NC_000913/3955443‑3955872                                                                                                                                                                                                                                   |                                                                                                                                                                                                             aaaaaCTGTCGCCGCGTTTGGATAACGTCATGCGGGTGAAGCGCGAAGATCGCCAGCCAGGGCACAGGTTTAAGCGGCGCGGCGCATCCGCCATGAGGGCGGGCCGGACGAAAGAACAGAAAGCGCACGGCGTGATCACTG‑CTTCTGCGGGTAACCCCGCGCGGGGCGTCGCGTTTTCTTCTGCGCGGTTGGGCGTGAAGGCCCTGATCGTTATGCCAACCGCCACCGCCGACATCaa                                                                                                                                                                                                  <  2:414262/238‑1 (MQ=255)                          cGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTG‑CTTCTGCGGGTAACCACGCGCAGGGCGTCGCGTTTTCTTCTGCGCGGTTAGGCGTGAAGGCCCTGATCGTTATGCCAACCGCCACCGCCGACATCAAAGTCGACGCGGTGCGCGGCTTcggcg                                                                                                                                                                        >  1:252313/1‑239 (MQ=255)                          cGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTG‑CTTCTGCGGGTAACCACGCGCAGGGCGTCGCGTTTTCTTCTGCGCGGTTAGGCGTGAAGGCCCTGATCGTTATGCCAACCGCCACCGCCGACATCAAAGTCGACGCGGTGCGCGGCTTcggcg                                                                                                                                                                        >  1:7196/1‑239 (MQ=255)                             cttgttGGTGACGCGCAAAGACCGCCGGCCGGTGCACAGTTTGAGGCTGCGCGGCGATAACGCCGTGTGGGGGGCCCGGACGGAGAAACAGAAAG‑GCCCGGGGTTGTCACCGCCTTTGGGGGGTAACCACGCGCGGGGCGTGGGGTTTTCTTCGGCGGGTTTGGGGGTGAAGCCCCTGGTCGTTATGCCACCCGCCACCGCCGACATCAAAGTCGACGCGGTGCGCGGCTTCGgcggcg                                                                                                                                                                     <  2:519607/234‑1 (MQ=255)                                  tGGAGAGGCGGGAGGACCGCCAGCCAGGCCACAGATTTAGGCGCCGGGGCGCAGACGCCATGATGGCGGGCTGGACGGAAGAACAAAAGGCGCCCGGCGTGATCCCTG‑CTTCTGCGGGTAACCACGCGCGGGGCGTCGCTTTTTCTTCTGCGCGGTTGGGCGTGAAGGCCCTGATGTTTATGCCAACCGCCACCGCCGACATCAAAGTCGACGCGGTGCGCGGCt                                                                                                                                                                              <  2:236658/225‑1 (MQ=255)                                            cgAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTG‑CTTCTGCGGGTAACCACGCGCAGGGCGTCGCGTTTTCTTCTGCGCGGTTAGGCGTGAAGGCCCTGATCGTTATGCCAACCGCCACCGCCGACATCAAAGTCGACGCGGTGCGCGGCTTCGGCGGCGa                                                                                                                                                                    >  1:236658/1‑225 (MQ=255)                                                             tGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTG‑CTTCTGCGGGTAACCACGCGCAGGGCGTCGCGTTTTCTTCTGCGCGGTTAGGCGTGAAGGCCCTGATCGTTATGCCAACCGCCACCGCCGACATCAAAGTCGACGCGGTGCGCGGCTTCGGCGGCGAAGTGCTGCTCCACGGCGCGAACTTTGATGaa                                                                                                                                     >  1:530047/1‑239 (MQ=255)                                                                       tAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTG‑CTTCTGCGGGTAACCACGCGCAGGGCGTCGCGTTTTCTTCTGCGCGGTTAGGCGTGAAGGCCCTGATCGTTATGCCAACCGCCACCGCCGACATCAAAGTCGACGCGGTGCGCGGCTTCGGCGGCGAAGTGCTGCTCCACGGCGCGAACTTTGATGAAGCGAAAGCCa                                                                                                                           >  1:398539/1‑239 (MQ=255)                                                                              cgcgGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTG‑CTTCTGCGGGTAACCACGCGCAGGGCGTCGCGTTTTCTTCTGCGCGGTTAGGCGTGAAGGCCCTGATCGTTATGCCAACCGCCACCGCCGACATCAAAGTCGACGCGGTGCGCGGCTTCGGCGGCGAAGTGCTGCTCCACGGCGCGAACTTTGATGAAGCGAAAGCCAAAGCGa                                                                                                                     >  1:360920/1‑238 (MQ=255)                                                                                      aTACGCCATGATGGCGGGCCGGACGGAAAAACAGAAAGCGCACGGCGTGATCACTG‑CTTCTGCGGGTAACCACGCGCAGGGCGTCGCTTTTTCTTCTGCGCGGTTAGGCGTGAAGGCCCTGATCGTTATGCCAACCGCCACCGCCGACATCAAAGTCGACGCGGTGCGCGGCTTCGGCGGCGAAGTGCTGCTCCACGGCGCGAACTTTGATGAAGCGAAAGCCAAAGCGATCGAACTGt                                                                                                            <  1:462984/239‑1 (MQ=255)                                                                                              tgatgGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTG‑CTTCTGCGGGTAACCACGCGCAGGGCGTCGCGTTTTCTTCTGCGCGGTTAGGCGTGAAGGCCCTGATCGTTATGCCAACCGCCACCGCCGACATCAAAGTCGACGCGGTGCGCGGCTTCGgcggcg                                                                                                                                                                     <  1:286772/174‑1 (MQ=255)                                                                                              tgatgGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTG‑CTTCTGCGGGTAACCACGCGCAGGGCGTCGCGTTTTCTTCTGCGCGGTTAGGCGTGAAGGCCCTGATCGTTATGCCAACCGCCACCGCCGACATCAAAGTCGACGCGGTGCGCGGCTTCGgcggcg                                                                                                                                                                     >  2:286772/1‑174 (MQ=255)                                                                                                                            cgcACGGCGTGATCACTG‑CTTCTGCGGGTAACCACGCGCAGGGCGTCGCGTTTTCTTCTGCGCGGTTAGGCGTGAAGGCCCTGATCGTTATGCCAACCGCCACCGCCGACATCAAAGTCGACGCGGTGCGCGGCTTCGGCGGCGAAGTGCTGCTCCACGGCGCGAACTTTGATGAAGCGAAAGCCAAAGCGATCGAACTGTCAcag                                                                                                       >  2:577338/1‑206 (MQ=255)                                                                                                                            cgcACGGCGTGATCACTG‑CTTCTGCGGGTAACCACGCGCAGGGCGTCGCGTTTTCTTCTGCGCGGTTAGGCGTGAAGGCCCTGATCGTTATGCCAACCGCAACCGCCGACATCAAAGTCGACGCGGTGCGCGGCTTCGGCGGCGAAGTGCTGCTCCACGGCGCGAACTTTGATGAAGCGAAAGCCAAAGCGATCGAACTGTCAcag                                                                                                       <  1:577338/206‑1 (MQ=255)                                                                                                                             gcACGGCGTGATCACTG‑CTTCTGCGGGTAACCACGCGCAGGGCGTCGCGTTTTCTTCTGCGCGGTTAGGCGTGAAGGCCCTGATCGTTATGCCAACCGCCACCGCCGACATCAAAGTCGACGCGGTGCGCGGCTTCGGCGGCGAAGTGCTGCTCCACGGCGCGAACt                                                                                                                                             <  1:494780/167‑1 (MQ=255)                                                                                                                             gcACGGCGTGATCACTG‑CTTCTGCGGGTAACCACGCGCAGGGCGTCGCGTTTTCTTCTGCGCGGTTAGGCGTGAAGGCCCTGATCGTTATGCCAACCGCCACCGCCGACATCAAAGTCGACGCGGTGCGCGGCTTCGGCGGCGAAGTGCTGCTCCACGGCGCGAACt                                                                                                                                             >  2:494780/1‑167 (MQ=255)                                                                                                                                   cgggATCGCTG‑CTTCGGCGGGAAACCACGCGCAGGGGGTCGCGTTTTCTTCTGCGCGGTTAGGCGTGAGGGCCCTGACCTTTATGCCAACCGCCACCGCCGACATCAAAGTCGACGCGGTGCGCGGc                                                                                                                                                                               <  2:87265/124‑1 (MQ=255)                                                                                                                                     tGATCACTG‑CTTCTGCGGGTAACCACGCGCAGGGCGTCGCGTTTTCTTCTGCGCGGTTAGGCGTGAAGGCCCTGATCGTTATGCCAACCGCCACCGCCGACATCAAAGTCGACGCGGTGCGCGGCTTCGGCGGCGAAGTGCTGCTCCACGGCGCGAACTTTGATGAAGCGAAAGCCAAAGCGATCGAACTGTCACAGCAGCACGGGTTCACCTGGGTGCCGCCGTTCGACAACCCGATg                                                             >  2:463281/1‑239 (MQ=255)                                                                                                                                             g‑CTTCTGCGGGTAACCACGCGCAGGGCGTCGCGTTTTCTTCTGCGCGGTTAGGCGTGAAGGCCCTGATCGTTATGCCAACCGCCACCGCCGACATCAAAGTCGACGCGGTGCGCGGCTTCGgcggcg                                                                                                                                                                     >  1:87265/1‑127 (MQ=255)                                                                                                                                             g‑CTTCTGCGGGTAACCACGCGCAGGGCGTCGCGTTTTCTTCTGCGCGGTTAGGCGTGAAGGCCCTGATCGTTATGCCAACCGCCACCGCCGACATCAAAGTCGACGCGGTGCGCGGCTTCGGCGGCGAAGTGCTGCTCCACGGCGCGAACTTTGATGAAGCGAAAGCCAAAGCGATCGAACTGTCACAGCAGCAGGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTg                                                      >  1:456109/1‑238 (MQ=255)                                                                                                                                                  cTGCGGGTAACCACGCGCAGGGCGTCGCGTTTTCTTCTGCGCGGTTAGGCGTGAAGGCCCTGATCGTTATGACAACCGCCACCGCCGACATCAAAGTCGACGCGGTGCGCGGCTTCGGCGGCGAAGTGCTGCTCCACGGCGCGAACTTTGATGAAGCGAAAGCCAAAGCGATCGAACTGTCACAGCAGCAGGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCgg                                                  >  1:867243/1‑238 (MQ=255)                                                                                                                                                         tAACCACGCGCAGGGCGTCGCGTTTTCTTCTGCGCGGTTAGGCGTGAAGGCCCTGATCGTTATGCCAACCGCCACCGCCGACATCAAAGTCGACGCGGTGCGCGGCTTCGGCGGCGAAGTGCTGCTCCACGGCGCGAACTTTGATGAAGCGAAAGCCAAAGCGATCGAACTGTCACAGCAGCAGGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGc                                           >  1:640091/1‑238 (MQ=255)                                                                                                                                                             cACGCGCGGGGGGTCGCGTTTTCTTCGGCGCGGTTGGGCGAGAAGGCCGTGATCGTTATGCCAACCGCCACCGCCGACATCAAAGTCGACGCGGTGCGCGGCTTCGGCGGCGAAGTGCTGCTCCACGGCGCGAACTTTGATGAAGCGAAAGCCAAAGCGATCGAACTGTCACAGCAGCAGGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCAcgct                                      <  2:456109/239‑1 (MQ=255)                                                                                                                                                                  gcgggggttggcgTTTTCTGCTGCGGGGTTGGGCGTGAAGGCCCTGATCGTTAGGCCACCCGCCACCGCCGACATCAAAGTCGACGCGGTGCGCGGCTTCGGCGGCGAAGTGCTGCTCCACGGCGCGAACTTTGATGAAGCGAAAGCCAAAGCGATCGAACTGTCACAGCAGCAGGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCAcgctggcg                                  <  2:360920/228‑1 (MQ=255)                                                                                                                                                                  gcAGGGCGTCGCGTTTTCTTCTGCGCGGTTAGGCGTGAAGGCCCTGATCGTTATGCCAACCGCCACCGCCGACATCAAAGTCGACGCGGTGCGCGGCTTCGGCGGCGAAGTGCTGCTCCACGGCGCGAACTTTGATGAAGCGAAAGCCAAAGCGATCGAACTGTCACAGCAGCAGGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCAc                                         >  1:402551/1‑231 (MQ=255)                                                                                                                                                                                   gtgcggcaggtttAGGCGAGAAGCCCGTGCTGGGTACGCCACCCGCCACCGCCGACAACGAAGGCGGCGGGGTGCGCGGCTTCGGCGGCGAAGTGCTGCTCCACGGCGCGAACTTTGATGAAGCGAAAGCCAAAGCGATCGAACTGTCACAGCAGCAGGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCcagcag                <  2:398539/228‑1 (MQ=255)                                                                                                                                                                                                        aGGCCCTGATCTTTAGGCCAACCGCCACCGCCGACATCAAAGTCGACGCGGTGCGCGGCTTCGGCGGCGAAGTGCTGCTCCACGGCGCGAACTTTGATGAAGCGAAAGCCAAAGCg                                                                                                                      <  2:729646/116‑1 (MQ=255)                                                                                                                                                                                                        aGGCCCTGATCGTTATGCCAACCGCCACCGCCGACATCAAAGTCGACGCGGTGCGCGGCTTCGGCGGCGAAGTGCTGCTCCACGGCGCGAACTTTGATGAAGCGAAAGCCAAAGCg                                                                                                                      >  1:729646/1‑116 (MQ=255)                                                                                                                                                                                                                     tATGCCAACCGTCACCGCCGACATCAAAGTCGACGCGGTGCGCGGCTTCGGCGGCGAAGTGCTGCTCCACGACGCGAACTTTGATGAAGCGAAAGCCAAAGCGATCGAACTGTCACAGCAGCAGGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGa  >  1:825032/1‑219 (MQ=255)                                                                                                                                                                                                                     tATGCCAACCGCCACCGCCGACATCAAAGTCGACGCGGTGCGCGGCTTCGGCGGCGAAGTGCTGCTCCACGGCGCGAACTTTGATGAAGCGAAAGCCAAAGCGATCGAACTGTCACAGCAGCAGGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGa  >  1:531966/1‑219 (MQ=255)                                                                                                                                                                                                                     tATGCCAACCGCCACCGCCGACATCAAAGTCGACGCGGTGCGCGGCTTCGGCGGCGAAGTGCTGCTCCACGGCGCGAACTTTGATGAAGCGAAAGCCAAAGCGATCGAACTGTCACAGCAGCAGGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGa  <  2:531966/219‑1 (MQ=255)                                                                                                                                                                                                                            accgccaccgccGACATCAAAGTCGACGCGGTGCGCGGCTTCGGCGGCGAAGTGCTGCTCCACGGCGCGAACTTTGATGAAGCGAAAGCCAAAGCGATCGAACTGTCAcagcagca                                                                                                  >  2:819412/1‑116 (MQ=255)                                                                                                                                                                                                                            accgccaccgccGACATCAAAGTCGACGCGGTGCGCGGCTTCGGCGGCGAAGTGCTGCTCCACGGCGCGAACTTTGATGAAGCGAAAGCCAAAGCGATCGAACTGTCACAGCAGCAGGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACg            <  2:469345/202‑1 (MQ=255)                                                                                                                                                                                                                            accgccaccgccGACATCAAAGTCGACGCGGTGCGCGGCTTCGGCGGCGAAGTGCTGCTCCACGGCGCGAACTTTGATGAAGCGAAAGCCAAAGCGATCGAACTGTCACAGCAGCAGGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACg            >  1:469345/1‑202 (MQ=255)                                                                                                                                                                                                                            accgccaccgccGACATCAAAGTCGACGCGGTGCGCGGCTTCGGCGGCGAAGTGCTGCTCCACGGCGCGAACTTTGATGAAGCGAAAGCCAAAGCGATCGAACTATCAcagcagca                                                                                                  <  1:819412/116‑1 (MQ=255)                                                                                                                                                                                                                                   |                                                                                                                                                                                                             AAAAACTGTCGTCGCGTCTTGATAACGTCATTCTGGTGAAGCGCGAAGATCGCCAGCCAGTGCACAGCTTTAAGCTGCGCGGCGCATACGCCATGATGGCGGGCCTGACGGAAGAACAGAAAGCGCACGGCGTGATCACTG‑CTTCTGCGGGTAACCACGCGCAGGGCGTCGCGTTTTCTTCTGCGCGGTTAGGCGTGAAGGCCCTGATCGTTATGCCAACCGCCACCGCCGACATCAAAGTCGACGCGGTGCGCGGCTTCGGCGGCGAAGTGCTGCTCCACGGCGCGAACTTTGATGAAGCGAAAGCCAAAGCGATCGAACTGTCACAGCAGCAGGGGTTCACCTGGGTGCCGCCGTTCGACCATCCGATGGTGATTGCCGGGCAAGGCACGCTGGCGCTGGAACTGCTCCAGCAGGACGCCCATCTCGA  >  NC_000913/3955443‑3955872

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