New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1539752	38 (0.880)	3 (0.070)	3/580	5.6	8.7%	coding (2839/3741 nt)	narZ	nitrate reductase 2 (NRZ), alpha subunit
	?	NC_000913	= 1539851	26 (0.620)					coding (2740/3741 nt)	narZ	nitrate reductase 2 (NRZ), alpha subunit
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below cutoff threshold.

GGCTGGAGATAATTTTACGCGGCTGCGCCTGAATATCGCGAAAGCGAATCTTCTCGTCCTCTTTGTGCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTGGTGTAATTGAGTTTGCCGAGGAAATCGATTTCATCCTGGGTATTCCACGAAATCCCTTTACCGCCATTGCCAAGTTTGTCCATCAATGGCCCGAGCGAGGTAAAGCGTTCATACGTAGCAGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCGATACCTTTGTAGATTTCCCAGTCGCTGCGTGATTCCCACGCCGGAT  >  NC_000913/1539563‑1540150                                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                                            ggCTGGAGATAATTTTACGCGGCTGCGCCTGAATATCGCGAAAGCGAATCTTCTCGTCCTCTTTGTGCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTGGTGTAATTGAGTTTGCCGAGGAAATCAATTTCATCCTGGGTATTCCACGAAATCCCTTTACCGCCAttac                                                                                                                                                                                                                                                                                                  >  2:139336/1‑298 (MQ=255)                          cgcCTGAATATCGCGAAAGCGAATCTTCTCGTCCTCTTTGTGCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTCGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGTCGGGCCATCACGCTTGGTTTAATTGAGTTTGCCGAGGCAAACGATTTCATCCTCGGTATTCCACGAAATCCCTTTACCCCCATTGCCAAGTTTGTCCATCAATGGCCcgag                                                                                                                                                                                                                                                                         >  2:10461/1‑300 (MQ=255)                                     ctggAAAGCGTATCTTCTCGGCCTCGTTGGGCAGCGCCAGGTGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTGGTGTAATTGAGTTTGCCGAGGAAATCGATTTCATCCTGGGTATTCCACGAAATCCCTTTACCGCCATTGCCAAGTTTGTCCATCAATGGCCCGAGCGAGGTAAAGCg                                                                                                                                                                                                                                                             <  1:45317/298‑1 (MQ=255)                                             cGAATCTTCTCGTCCTCTTTGTGCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTGGTGTAATTGAGTTTGCCGAGGAAATCGATTTCATCCTGGGTATTCCACGAAATCCCTTTACCGCCATTGCCAAGTTTGTCCATCAATGGCCCGAGCGAGGTAAAGCGTTCATa                                                                                                                                                                                                                                                       <  1:211020/299‑1 (MQ=255)                                                   ttctCGTCCTCTTTGTGCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTGTTGTAATTGAGTTTGTCGAGGAAATCGATTTCATCCTGGGTATTCCACGAAATCCCTTTACCGCAATTGCCTAGTTTGTCCATTAAGGGCCCGAGCGAGGTAAAGCGGTCACACGTatacg                                                                                                                                                                                                                                               >  2:157433/1‑297 (MQ=255)                                                                   gCAGGGCCAAAGGGGTAGGTTCGGGCCGGGTTATCTCGCCGAGCGCCTGCCACGCGTTGACTGCAACATGACCGTTGGTTTCCGGGGCCAGTGCCCGAATCACTTCTGAAGCGTCAATGGCGGGGTAAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTGGTGTAATTGAGTTTGCCGAGGAAATCGATTTCATCCTGGGTATTCCACGAAATCCCTTTACCGCCATTGCCAAGTTTGTCCATCAATGGCCCGAGCGAGGTAAAGCGTTCATACGTAGCAGGGTAGTCGCGCTCCa                                                                                                                                                                                                                                <  2:381785/300‑1 (MQ=255)                                                                                            cccGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTGGTGTAATTGAGTTTGCCGAGGAAATCGATTTCATCCTGGGTATTCCACGAAATCCCTTTACCGCCATTGCCAAGTTTGTCCATCAATGGCCCGAGCGAGGTAAAGCGTTCATACGTAACCGGGTAGTCGCGCTCCCCCGCCACAATTTTCGGCCCGGttttt                                                                                                                                                                                                      >  2:91069/1‑300 (MQ=255)                                                                                                gTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTGGTGTAATTGAGTTTGCCGAGGAAATCGATTTCATCCTGGGTATTCCACGAAATCCCTTTACCGCCATTGCCAAGTTTGTCCATCAATGGCCCGAGCGAGGTAAAGCGTTCATACGTAGCAGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGCGGTTTTGCCt                                                                                                                                                                                                   <  1:70049/300‑1 (MQ=255)                                                                                                   aacgcGCCCAGCGCCGTCGCCGCTGTGACTGCAACATGACCGTGGGTTTCCGGTGCCATTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTGGTGTAATTGAGTTTGCCGAGGAAATCGATTTCATCCTGGGTATTCCACGAAATCCCTTTACCGCCATTGCCAAGTTTGTCCATCAATGGCCCGAGCGAGGTAAAGCGTTCATACGTAGCAGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGCGGTTTTGCCTGGaa                                                                                                                                                                                               <  2:108893/297‑1 (MQ=255)                                                                                                          ccAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTGGTGTAATTGAGTTTGCCGAGGAAATCGATTTCATCCTGGGTATTCCACGAAATCCCTTTACCGCCATTGCCAAGTTTGTCCATCAATGGCCCGAGCGAGGTAAAGCGTTCATACGTAGCAGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGCGGTTTTGCCTGGAATCAGATc                                                                                                                                                                                        <  1:157439/301‑1 (MQ=255)                                                                                                                  tGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTGGTGTAATTGAGTTTGCCGAGGAAATCGATTTCATCCTGGGTATTCCACGAAATCCCTTTACCGCCATTGCCAAGTTTGTCCATCAATGGCCCGAGCGAGGTAAAGCGTTCATACGTAGCAGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGc                                                                                                                                                                                <  1:139341/301‑1 (MQ=255)                                                                                                                    ccACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTGGTGTAATTGAGTTTGCCGAGGAAATCGATTTCAGCCTGGGTATTCCACGAAATCCCTTTACCGCCATTGCCAAGTTTGTCCATCAATGGCCCGAGCGAGGTAAAGCGTTCATACGTAGCAGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGc                                                                                                                                                                                >  1:147033/1‑299 (MQ=255)                                                                                                                      ccgcgtgaaCTGCAACAGTACCGTTGGTTTCCGGTGCNAGTGCCAGAATCACTTCTGAAGATTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTGGTGTAATTGAGTTTGCCGAGGAAATCGATTTCATCCTGGGTATTCCACGAAATCCCTTTACCGCCATTGCCAAGTTTGTCCATCAATGGCCCGAGCGAGGTAAAGCGTTCATACGTAGCAGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCttt                                                                                                                                                                            <  2:56336/293‑1 (MQ=255)                                                                                                                             gACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTGGTGTAATTGAGTTTGCCGAGGAAATCGATTTCATCCTGGGTATTCCACGAAATCCCTTTACCGCCATTGCCACGTTTGTCAATCAATGGCCCGAGCGAGGTAACGCGTTCATACGTAGCAGGGTTGCCGCGCTCCACCGCCACACACTTCGGCGCGGCTTTGCTGGGAATCGGATCCCATTCGCCTTTGCGCCAg                                                                                                                                                                     >  2:29323/1‑301 (MQ=255)                                                                                                                                    aCATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTGGTGTAATTGAGTTTGCCGAGGAAATCGATTTCATCCTGGGTATTCCACGAAATCCCTTTACCGCCATTGCCAAGTTTGTCCATCAATGGCCCGAGCGAGGTAAAGCGTTCATACGTAGCAGGGTAGTCGCGCTCCACCGCCCCAATATTCGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCa                                                                                                                                                              >  2:394568/1‑301 (MQ=255)                                                                                                                                           tgTTGGTTTCCTGTGCCCTGGCCAGAATGCCTGCTTGAGTGTCAATTGTGGTGTCAATCAGCGGACGCCCCTGGTCCTGGCCATCACGGTTTGTGTAATTGAGTTTGCCGAGGAAATCGATTTCATCCTGGGTATTCCACGAAATCCCTTTACCGCCATTGCCAAGTTTGTCCATCAATGGCCCGAGCGAGGTAAAGCGTTCATACGTAGCAGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTcaca                                                                                                                                                       <  2:216437/300‑1 (MQ=255)                                                                                                                                                gTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTGGTGTAATTGAGTTTGCCGAGGAAATCGATTTCATCCTGGGTATTCCACGAAATCCCTTTACCGCCATTGCCAAGTTTGTCCATCAATGGCCCGAGCGAGGTAAAGCGTTCATACGTAGCAGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACgg                                                                                                                                                    >  1:263160/1‑299 (MQ=255)                                                                                                                                                     cGGGGCCAGTGGCAGAATCGCTTATGAAGCGTCAATGGCGGTTGCAATCAGCGGACACCCCTGTGGCGGGCCATCACGCTTGGTGTAATTGAGTTTGCCGAGGAAATCGATTTCATCCTGGGTATTCCACGAAATCACTTTACCGCCATTGCCAAGTTTGTCCATCAATGGCCCGGGCGAGGTAAAGCGTTCATACGTAGAAGAATAGTCGCGCTCCACCGCCACAATATTCGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCAAGCACTTCACACGGCTGTgag                                                                                                                                             <  1:10461/301‑1 (MQ=255)                                                                                                                                                               gccagAATCACTGCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTGGTGTAATTGAGTTTGCCGAGGAAATCGATTTCATCCTGGGTATTCCACGAAATCCCTTTACCGCCATTGCCAAGTTTGTCCATCAATGGCCCGAGCGAGGTAAAGCGTTCATACGTAGCAGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGcc                                                                                                                                    <  1:91072/300‑1 (MQ=255)                                                                                                                                                                         cTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTGGTGTAATTGAGTTTGCCGAGGAAATCGATTTCATCCTGGGTATTCCACGAAATCCCTTTACCGCCATTGCCAAGTTTGTCCATCAATGGCCCGAGCGAGGTAAAGCGTTCATACGTAGCAGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAGAGTCATg                                                                                                                         >  1:360806/1‑301 (MQ=255)                                                                                                                                                                                  cGTCAATTGCGGTGTAAAGCAGCGGACGCACCTGTGCCGGGCCATCACGCTTGGTGTAAGTGGGTTTGCCGAGGAAATCGATTTCATCCTGGGTATTCCACGAAATCCCTTTACCGCCATTGCCAAGTTTGTCCATCAATGGCCCGAGCGAGGTAAAGCGTTCATACGTAGCAGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAGAGTCATGCAGTAgggg                                                                                                                <  1:209574/301‑1 (MQ=255)                                                                                                                                                                                           cGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTGGTGTAATTGAGTTTGCCGAGGAAATCGATTTCATCCTGGGTATTCCACGAAATCCCTTTACCGCCATTGCCAAGTTTGTCCATCAATGGCCCGAGCGAGGTAAAGCGTTCATACGTAGCAGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAGAGTCATGCAGTAGGGGTTGTAATAc                                                                                                       <  1:394586/301‑1 (MQ=255)                                                                                                                                                                                                                ccTGTGCCGGGCCATCACGCTTGGTGTAATTGAGTTTGCCGAGGAAATCGATTTCATCCTGGGTATTCCACGAAATCCCTTTACCGCCATTGCCAAGTTTGTCCATCAATGGCCCGAGCGAGGTAAAGCGTTCATACGTAGCAGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGCGCACTTCACACGGCTGTGTGAGCTCTTCCGGAGAGTCATGCAGTAGGGGTTGTAATTCAACGTCGGTTTCTCTGCCCAg                                                                                  >  1:278698/1‑301 (MQ=255)                                                                                                                                                                                                                  ttgtccgtggcATCACGCTTGGTTTTATTGAGTTTGCCGCGTAAATCGATTGCATACTTGCTATTCCACGAACACCCTTTACCGCCATTGCCAAGTTTGTCCAGCAATGACCCGAACCAGGTAAAGCTTTCATTCGTACCCGGGTCGTCGCGCTCCCGCGCCACAATATTCGGCACGTTTTTGCCTGGACTCAGATCGCATTCGCCTGTGCGGCAGTAGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAGAGTCATGCAGTAGGGGTTGTGATACTACATCGGTTTCTTTGCCAatat                                                                                <  2:147028/291‑1 (MQ=255)                                                                                                                                                                                                                                                      ccGAGGAAATCGATTTCATCCTGGGTATTCCACGAAATCCCTTTACCGCCATTGCCAAGTTTGTCCATCAATGGCCCGAGCGAGGTAAAGCGTTCATACGTAGCAGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCCCACGGCTGTGAGCGCTCTGCCGGAGAGTCATGCAGTCGGGGTTGTAATCCCACCTCGGTTTCTTTGCCAAGATTACCCACGCACCCTTCCTAAAATGCttttgacata                                             >  2:382224/1‑293 (MQ=255)                                                                                                                                                                                                                                                          ggAAATCGATTTCATCCTGGGTATTCCACGAAATCCCTTTACCGCCATTGCCAAGTTTTTCCATCAATGGCCCGAGCGAGGTAAAGCGTTCATACGTAGCAGGGTAGTCACGCTCCACCGCCACAATATTCGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCGATACCtt                                         <  1:29323/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                                ttACCGCCATTGCCAAGTTTGTCCATCAATGGCCCGAGCGAGGTAAAGCGTTCATACGTAGCAGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCGATACCTTTGTAGATTTCCCAGGCGCTGCTTGATTCCCACGCCGGaa  >  1:238284/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                                            GGCTGGAGATAATTTTACGCGGCTGCGCCTGAATATCGCGAAAGCGAATCTTCTCGTCCTCTTTGTGCAGCGCCAGATGGGTATGTTCGCGCCCGGTGATCTCGCCCAGCGCCTGCCACGCTTTGACTGCAACATGACCGTTGGTTTCCGGTGCCAGTGCCAGAATCACTTCTGAAGCGTCAATGGCGGTGTCAATCAGCGGACGCCCCTGTGCCGGGCCATCACGCTTGGTGTAATTGAGTTTGCCGAGGAAATCGATTTCATCCTGGGTATTCCACGAAATCCCTTTACCGCCATTGCCAAGTTTGTCCATCAATGGCCCGAGCGAGGTAAAGCGTTCATACGTAGCAGGGTAGTCGCGCTCCACCGCCACAATATTCGGCGCGGTTTTGCCTGGAATCAGATCGCATTCGCCTTTGCGCCAGTCGAGCACTTCACACGGCTGTGAGAGCTCTGCCGGAGAGTCATGCAGTAGGGGTTGTAATACCACGTCGGTTTCTTTGCCAAGATGACCCACGCACACTTGCGAAAATGCTTTGGCGATACCTTTGTAGATTTCCCAGTCGCTGCGTGATTCCCACGCCGGAT  >  NC_000913/1539563‑1540150

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