New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	961784 =	57 (1.230)	3 (0.070)	3/484	4.0	6.4%	coding (584/684 nt)	cmk	cytidylate kinase
	?	NC_000913	= 4134337	31 (0.680)					coding (503/2181 nt)	katG	catalase/hydroperoxidase HPI
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CAACCGTGGTGGCCAGCCGACTGGGGCAGTTACGCCGGTCTGTTTATTCGTATGGCCTGGCACGGCGCGGGGACTTACCGTTCAATCGATGGACGCGGTGGCGCGGGTCGTGGTCAGCAACGTTTTGCACCGCTGAACTCCTGGCCGGATAACGTAAGCCTCGATAAAGCGCGTCGCCTGTTGTGGCCAATCAAACAGAAATATGGTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCG  >  NC_000913/4134090‑4134582                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      cAACCGTGGTGGCCAGCCGACTGGGGCAGTTACGCCGGTCTGTTTATTCGTATGGCCTGGCACGGCGCGGGGACTTACCGTTCAATCGATGGACGCGGTGGCGCGGGTCGTGGTCAGCAACGTTTTGCACCGCTGAACTCCTGGCCGGATAACGTAAGCCTCGATAAAGCGCGTCGCCTGTTGTGGCCAATCAAACAGCAATATGGTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACgg                                                                                                                                                                                                                                                    >  2:197460/1‑250 (MQ=255)           ggCCAGCCGACTGGGGCAGTTACGCCGGTCTGTTTATTCGTATGGCCTGGCACGGCGCGGGGACTTACCGTTCAATCGATGGACGCGGTGGCGCGGGTCGTGGTCAGCAACGTTTTGCACCGCTGAACTCCTGGCCGGATAACGTAAGCCTCGATAAAGCGCGTCGCCTGTTGTGGCCAATCAAACAGAAATATGGTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTGACGTGGCGCTAGaa                                                                                                                                                                                                                                          <  1:354705/251‑1 (MQ=255)                                 cGCCGGTCTGTTTATTCGTATGGCCTGCCACGGCGCGGGGACTTACCGCTCACTCGCTGGACGCGGTGGCGCGGGTCGTGGTCAGCAACGTTTTGCCCCGCTGAACTCCTGGCCGGATAACGTAAGCCTCGATAAAGCGCGTCGCCTGTTGTGGCCAATCAAACAGAAATATGGTCAGAAAATCTCCTGGGCCGCCCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAACCTCCGGCTTCCGTACCTTCg                                                                                                                                                                                                                    >  2:273438/1‑251 (MQ=255)                                     ggTCTGTTTATTCGTATGGCCTGGCACGGCGCGGGGACTTACCGTTCAATCGATGGACGCGGTGGCGCGGGTCGTGGTCAGCAACGTTTTGCACCGCTGAACTCCTGGCCGGATAACGTAAGCCTCGATAAAGCGCGTCGCCTGTTGTGGCCAATCAAACAGAAATATGGTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGttt                                                                                                                                                                                                                <  1:63090/251‑1 (MQ=255)                                             tATTCGTATGGCCTGGCACGGCGCGGGGACTTACCGTTCAATCGATGGACGCGGTGGCGCGGGTCGTGGTCAGCAACGTTTTGCACCGCTGAACTCCTGGCCGGATAACGTAAGCCTCGATAAAGCGCGTCGCCTGTTGTGGCCAATCAAACAGAAATATGGTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCg                                                                                                                                                                                                        >  2:357637/1‑251 (MQ=255)                                                                 gcgcGGGGACTTACCGTTCAATCGATGGACGCGGTGGCGCGGGTCGTGGTCAGCAACGTTTTGCACCGCTGAACTCCTGGCCGGATAACGTAAGCCTCGATAAAGCGCGTCGCCTGTTGTGGCCAATCAAACAGAAATATGGTCAGAAAATCTCCTGGGCCGACCTGTTTAGCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGaa                                                                                                                                                                                    <  2:26271/251‑1 (MQ=255)                                                                     ggggACTTACCGTTCAATCGATGGACGCGGTGGCGCGGGTCGTGGTCAGCAACGTTTTGCACCGCTGAACTCCTGGCCGGATAACGTAAGCCTCGATAAAGCGCGTCGCCTGTTGTGGCCAATCAAACAGAAATATGGTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGCACCgg                                                                                                                                                                                >  2:172143/1‑251 (MQ=255)                                                                         aCTTACCTTTCAATCGATGGACGCGGTGGCGCGGGTCGTGGACAGCAACGTTTTGCACCGCTGAACTCCTGGACGGATAACGTAAGCCTCGATAAAGCGCGTCGCCTGTTGTGGCCAATCAAACAGAAATATGGTCAGAAAATCTCCTGGGCGGACCTGTTTATCCTCGCGGGTAACGTGGCTCTATAAAACGCCGGCTTCCTTACCTGCGGTTGTGGTGCCGGTCGTGATGACGTCTGGTAACCGGATCt                                                                                                                                                                            <  2:48615/251‑1 (MQ=255)                                                                         aCTTACCGTTCAATCGATGGACGCGGTGGCGCGGGTCGTGGTCAGCAACGTTTTGCACCGCTGAACTCCTGGCCGGATAACGTAAGCCTCGATAAAGCGCGTCGCCTGTTGTGGCCAATCAAACAGAAATATGGTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCt                                                                                                                                                                            <  1:91150/251‑1 (MQ=255)                                                                         aCTTACCGTTCAATCGATGGACGCGGTGGCGCGGGTCGTGGTCAGCAACGTTTTGCACCGCTGAACTCCTGGCCGGATAACGTAAGCCTCGATAAAGCGCGTCGCCTGTTGTGGCCAATCAAACAGAAATATGGTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCt                                                                                                                                                                            <  1:314209/251‑1 (MQ=255)                                                                            tACCGTTCAATCGATGGACGCGGTGGCGCGGGTCGTGGTCAGCAACGTTTTGCACCGCTGAACTCCTGGCCGGATAACGTAAGCCTCGATAAAGCGCGTCGCCTGTTGTGGCCAATCAAACAGAAATATGGTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTg                                                                                                                                                                           <  1:415337/249‑1 (MQ=255)                                                                                gTTCAATCGATGGACGCGGTGGCGCGGGTCGTGGTCAGCAACGTTTTGCACCGCTGAACTCCTGGCCGGATAACGTAAGCCTCGATAAAGCGCGTCGCCTGTTGTGGCCAATCAAACAGAAATATGGTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGtt                                                                                                                                                                     >  1:355807/1‑251 (MQ=255)                                                                                                                ggTCAGCAACGTTTTGCACCGCTGACCTCCTGGCCGGATAACGTAAGCCTCGATAAAGCGCGTCGCCTGTTGTGGCCAATCAAACAGAAATATGGTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGTTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCa                                                                                                                                     <  1:436371/251‑1 (MQ=255)                                                                                                                    aGCAACGTTTTGCAACGCTGAACTCCTGGCCGGATAACGTAAGCCTCGATAAAGCGCGTCGCATGTTGTGGCCAATCAAACAGCAATATGGTCAGAAAATCTCCTGGGCCGCCCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCCTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCgc                                                                                                                                 >  1:180909/1‑250 (MQ=255)                                                                                                                          gTTTTGCCCCGCTGAACTACTGTCCGGAAAACCTAAGCCTCGATAAAGCGCGCCACCTGTTGTGGCCCATAAAACCGAACTATGGTCAGAACATCTCCTGGGCCGACCTGTTTATCCTCGCGGGAAACGTGGCGCTAGACACCTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGACGTGAAGACGTCTGGGAACCGGATCTGGATGTGAAATGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCg                                                                                                                           >  2:299357/1‑251 (MQ=255)                                                                                                                          gTTTTGCACCGCTGAACTCCAGGCCGGAGAACGTAAGCCTCGATAAAGCGCGTCGCCTGTTGTGGACAATCAAACAGAAATATGGTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGTTGCCGGTCGTGAAGACTTCTGTGAACCGGATCTGGATGTGAACTGGGTTGATGAAAAAGCCTGGCTGACTCACCGTCATCCg                                                                                                                           <  2:426695/251‑1 (MQ=255)                                                                                                                                   cGCTGAACTCCTGGCCGGATAACGTAAGCCTCGATAAAGCGCGTCGCCTGTTGTGGCCAATCAAACAGAAATATGGTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTg                                                                                                                  <  1:197460/251‑1 (MQ=255)                                                                                                                                         aCTCCTGGCCGGATAACGTAAGCCTCGATAAAGCGCGTAGCCTGATGTGCCCAATCAAACAGAAATATGGTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGaa                                                                                                             >  1:425364/1‑250 (MQ=255)                                                                                                                                                gCCGGATAACGTAAGCCTCGATAAAGCGCGTCGCCTGTTGTGGCCAATCAAACAGAAATATGGTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGCGCCGGTCGTGAAGACGACTGGGAACCGGATCTGGATGGTAACTGGGGCGATGAAAAAGCCTGTATGACTCACCGTCATCCGGAAGCGCGGGCGAAAGCACCGc                                                                                                     >  2:1279/1‑251 (MQ=255)                                                                                                                                                  cGGATAACGTAAGCCTCGATAAAGCGCGTCGCCTGTTGTGGCCAATCAAACAGAAATATGGTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCGCCGCTg                                                                                                   <  1:357637/251‑1 (MQ=255)                                                                                                                                                                  cGATAAAGCGCGTCGCCTGTTGTGGCCAATCACACAGAAATATGGTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGCGATgg                                                                                   >  2:45397/1‑251 (MQ=255)                                                                                                                                                                                    gttgtGGCCAATCAAACAGAAATATGGTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAAcc                                                                 <  1:233967/251‑1 (MQ=255)                                                                                                                                                                                         ggCCAATCAAACAGAAATATGGTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGa                                                             >  1:259717/1‑250 (MQ=255)                                                                                                                                                                                                  aaCAGAAATATGGTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGAt                                                   >  1:389085/1‑251 (MQ=255)                                                                                                                                                                                                    cAGAAATATGGTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAACAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGCTGGGGCTGATTTACGTTAACCCGGAAGGCCCGGaccc                                                 >  2:365349/1‑248 (MQ=255)                                                                                                                                                                                                             ggTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAgcggcg                                         <  1:78713/250‑1 (MQ=255)                                                                                                                                                                                                               tCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGAGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAAGACCGGATCACAGCGGCGaa                                       >  2:328246/1‑250 (MQ=255)                                                                                                                                                                                                                    aaaTCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCtt                                 <  2:355807/251‑1 (MQ=255)                                                                                                                                                                                                                             gggCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGgc                         >  1:119560/1‑250 (MQ=255)                                                                                                                                                                                                                                        tttATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCAGCTCTCGAACCTCCGCCTTTGGTGCCGGTCGAGAAAACGTCTGGGAACCGGATCTGGATGCTACCTGGGGTGATGCACAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCCCCGCTGGGTGCAACCGAGATGGGTCTGATTTCCGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACAGCTTTCTGCGGCAGCAGCTATCcg             >  1:437111/1‑251 (MQ=255)                                                                                                                                                                                                                                                   gggTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCg  <  1:1279/251‑1 (MQ=255)                                                                                                                                                                                                                                                   gggTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCg  >  1:25613/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      CAACCGTGGTGGCCAGCCGACTGGGGCAGTTACGCCGGTCTGTTTATTCGTATGGCCTGGCACGGCGCGGGGACTTACCGTTCAATCGATGGACGCGGTGGCGCGGGTCGTGGTCAGCAACGTTTTGCACCGCTGAACTCCTGGCCGGATAACGTAAGCCTCGATAAAGCGCGTCGCCTGTTGTGGCCAATCAAACAGAAATATGGTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCG  >  NC_000913/4134090‑4134582

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