New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1539119 =	46 (0.600)	3 (0.040)	3/252	NT	6.3%	coding (3472/3741 nt)	narZ	nitrate reductase 2 (NRZ), alpha subunit
	?	NC_000913	1539157 =	46 (0.620)					coding (3434/3741 nt)	narZ	nitrate reductase 2 (NRZ), alpha subunit

CCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGTTTCGGGCAAACGCGGGTAACCGAGTTATGAATGCCGCCGCGCATGCCAGTTACTTCCGAACCAGGAATATTCATAATGCGTTCCTGGGCGTGATACATCATGGTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAG  >  NC_000913/1538972‑1539263                                                                                                                                                    |                                                                                                                                                 cattcactcaacctacagtgaaaacctgctaatgctgacgctctctcgcggtggaccgattgtctggatcagcgaaacagatgcccgagaactgaccattgtcgataacgactgggtggaagtattcaacgccaacggcgcgctgacTg                                                                                                                                                 <  2:12353‑M1/2‑1 (MQ=255)  attcactcaacctacagtgaaaacctgctaatgctgacgctctctcgcggtggaccgattgtctggatcagcgaaacagatgcccgagaactgaccattgtcgataacgactgggtggaagtattcaacgccaacggcgcgctgacTGc                                                                                                                                                <  4:598523‑M1/3‑1 (MQ=255)   ttcactcaacctacagtgaaaacctgctaatgctgacgctctctcgcggtggaccgattgtctggatcagcgaaacagatgcccgagaactgaccattgtcgataacgactgggtggaagtattcaacgccaacggcgcgctgacTGcc                                                                                                                                               <  2:230819‑M1/4‑1 (MQ=255)                cgTAACCGCCAATCATATGCGTTGGTTTCGGGCAAACGCGGGTAACCGAGTTATGAATGCCGCCGCGCATGCCAGTTACTTCCGAACCAGGAATATTCATAATGCGTTCCTGGGCGTGATACATCATGGTCATGCCCGGCGGTACACGt                                                                                                                                  <  3:104732/149‑1 (MQ=255)                 gTAACCGCCAATCATATGCGTTGGTTTCGGGCAAACGCGGGTAACCGAGTTATGAATGCCGCCGCGCATGCCAGTTACTTCCGAACCAGGAATATTCATAATGCGTTCCTGGGCGTGATACATCATGGTCATGCCCGGCGGTACACGtt                                                                                                                                 >  2:564652/1‑149 (MQ=255)                     ccGCCAATCATATGCGTTGGTTTCGGGCAAACGCGGGTAACCGAGTTATGAATGCCGCCGCGCATGCCAGTTACTTCCGAACCAGGAATATTCATAATGCGTTCCTGGGCGTGATACATCATGGTCATGCCCGGCGGTACACGTTGGCt                                                                                                                             <  4:99902/149‑1 (MQ=255)                       gCCAATCATATGCGTTGGTTTCGGGCAAACGCGGGTAACCGAGTTATGAATGCCGCCGCGCATGCCAGTTACTTCCGAACCAGGAATATTCATAATGCGTTCCTGGGCGTGATACATCATGGTCATGCCCGGCGGTACACGTTGGCTGa                                                                                                                           >  3:633138/1‑149 (MQ=255)                         cAATCATATGCGTTGGTTTCGGGCAAACGCGGGTAACCGAGTTATGAATGCCGCCGCGCATGCCAGTTACTTCCGAACCAGGAATATTCATAATGCGTTCCTGGGCGTGATACATCATGGTCATGCCCGGCGGTACACGTTGGCTGacc                                                                                                                         >  1:56498/1‑149 (MQ=255)                          aaTCATATGGGTTGGTTTCGGGCAAACGCGGGTAACCGAGTTATGAATTCCGCCGCGCATGCCAGTTACTTCCGAACCAGGAATATTCATAATGCGTTCCTGGGCGTGATACATCCTGGTCATGCCCGGCGGTACACGTTGGCTGacca                                                                                                                        >  2:429708/1‑149 (MQ=255)                                   ctgacgctctctcgcggtggaccgattgtctggatcagcgaaacagatgcccgagaactgaccattgtcgataacgactgggtggaagtattcaacgccaacggcgcgctgacTGCCCGGCGGTACACGTTGGCTGACCACCGCGCggg                                                                                                               >  4:366665‑M1/114‑149 (MQ=255)                                   cGTTGGTTTCGGGCAAACGCGGGTAACCGAGTTATGAATGCCGCCGCGCATGCCAGTTACTTCCGAACCAGGAATATTCATAATGCGTTCCTGGGCGTGATACATCATGGTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCggg                                                                                                               <  2:129592/149‑1 (MQ=255)                                     gacgctctctcgcggtggaccgattgtctggatcagcgaaacagatgcccgagaactgaccattgtcgataacgactgggtggaagtattcaacgccaacggcgcgctgacTGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCa                                                                                                             >  2:321221‑M1/112‑149 (MQ=255)                                        gtctCGGGAAAACGCGGGTAACCGAGTTATGAATGCCGCCGCGCATGCCAGTTACTTCCGATCCAGGAATATTCATAATGCGTTCCTGGGAGTGATACATCATGGTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTc                                                                                                          <  4:420707/146‑1 (MQ=255)                                        gctctctcgcggtggaccgattgtctggatcagcgaaacagatgcccgagaactgaccattgtcgataacgactgggtggaagtattcaacgccaacggcgcgctgacTGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTc                                                                                                          >  1:6293‑M1/109‑149 (MQ=255)                                         tttCGGGCAAACGCGGGTAACCGAGTTATGAATGCCGCCGCGCATGCCAGTTACTTCCGAACCAGGAATATTCATAATGCGTTCCTGGGCGTGATACATCATGGTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCa                                                                                                         >  4:194163/1‑149 (MQ=255)                                             gggCAAACGCGGGTAACCGAGTTATGAATGCCGCCGCGCATGCCAGTTACTTCCGAACCAGGAATATTCATAATGCGTTCCTGGGCGTGATACATCATGGTCATGCCCGGCGGTa                                                                                                                                       >  1:591773/1‑115 (MQ=255)                                             gggCAAACGCGGGTAACCGAGTTATGAATGCCGCCGCGCATGCCAGTTACTTCCGAACCAGGAATATTCATAATGCGTTCCTGGGCGTGATACATCATGGTCATGCCCGGCGGTa                                                                                                                                       <  2:591773/115‑1 (MQ=255)                                             gggCAAACGCGGGTAACCGAGTTATGAATGCCGCCGCGCATGCCAGTTACTTCCGAACCAGGAATATTCATAATGCGTTCCTGGGCGTGATACATCATGGTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAgcgc                                                                                                     >  1:222590/1‑149 (MQ=255)                                               gCAAACGCGGGTAACCGAGTTATGAATGCCGCCGCGCATGCCAGTTACTTCCGAACCAGGAATATTCATAATGCGTTCCTGGGCGTGATACATCATGGTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAgcgcgc                                                                                                   >  1:371947/1‑149 (MQ=255)                                               gCAAACGCGGGTAACCGAGTTATGAATGCCGCCGCGCATGCCAGTTACTTCCGAACCAGGAATATTCATAATGCGTTCCTGGGCGTGATACATCATGGTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAgcgcgc                                                                                                   >  3:270236/1‑149 (MQ=255)                                                 aaaCGCGGGTAACCGAGTTATGAATGCCGCCGCGCATGCCAGTTACTTCCGAACCAGGAATATTCATAATGCGTTCCTGGGCGTGATACATCATGGTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCg                                                                                                 <  1:323708/149‑1 (MQ=255)                                                       gggTAACCGAGTTATGAATGCCGCCGCGCATGCCAGTTACTTCCGAACCAGGAATATTCATAATGCGTTCCTGGGCGTGATACATCATGGTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCg                                                                                           >  3:42351/1‑149 (MQ=255)                                                        ggTAACCGAGTTATGAATGCCGCCGCGCATGCCAGTTACTTCCGAACCAGGAATATTCATAATGCGTTCCTGGGCGTGATACATCATGGTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGt                                                                                          <  4:337633/149‑1 (MQ=255)                                                        ggTAACCGAGTTATGAATGCCGCCGCGCATGCCAGTTACTTCCGAACCAGGAATATTCATAATGCGTTCCTGGGCGTGATACATCATGGTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGt                                                                                          <  2:165942/149‑1 (MQ=255)                                                         gTAACCGAGTTATGAATGCCGCCGCGCATGCCAGTTACTTCCGAACCAGGAATATTCATAATGCGTTCCTGGGCGTGATACATCATGGTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGtt                                                                                         <  2:184543/149‑1 (MQ=255)                                                         gTAACCGAGTGATGAATGCCCCCGCGCATGCCTGTTACTTCCGAACCAGGAATATTCATAATGCGTTCCTAGGCGTGATACATCATGGTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCGGTCAGCGCGCCGTTGGCGtt                                                                                         >  4:240142/1‑149 (MQ=255)                                                              cGAGTTATGAATGCCGCCGCGCATGCCAGTTACTTCCGAACCAGGAATATTCATAATGCGTTCCTGGGCGTGATACATCATGGTCATGc                                                                                                                                                >  3:314865/1‑89 (MQ=255)                                                              cGAGTTATGAATGCCGCCGCGCATGCCAGTTACTTCCGAACCAGGAATATTCATAATGCGTTCCTGGGCGTGATACATCATGGTCATGc                                                                                                                                                <  4:314865/89‑1 (MQ=255)                                                              cGAGTTATGAATGCCGCCGCGCATGCCAGTTACTTCCGAACCAGGAATATTCATAATGCGTTCCTCGGCGTGATACATCATGGTCATGCCCGGCGGTACACGTTGGCTTACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATa                                                                                    >  1:482342/1‑149 (MQ=255)                                                                     tGAATGCCGCCGCGCATGCCAGTTACTTCCGAACCAGGAATATTCATAATGCGTTCCTGGGCGTGATACATCATGGTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCAc                                                                             >  1:698586/1‑149 (MQ=255)                                                                     tGAATGCCGCCGCGCATGCCAGTTACTTCCGAACCAGGAATATTCATAATGCGTTCCTGGGCGTGATACATCATGGTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCAc                                                                             >  1:230819/1‑149 (MQ=255)                                                                                 cgcATGCCAGTTACTTCCGAACCAGGAATATTCATAATGCGTTCCTGGGCGTGATACATCATGGTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATc                                                                 <  1:506515/149‑1 (MQ=255)                                                                                      gCCAGTTACTTCCGAACCAGGAATATTCATAATGCGTTCCTGGGCGTGATACATCATGGTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACaa                                                            <  2:87819/149‑1 (MQ=255)                                                                                      gCCAGTTACTTCCGAACCAGGAATATTCATAATGCGTTCCTGGGCGTGATACATCATGGTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACaa                                                            >  3:598523/1‑149 (MQ=255)                                                                                       ccAGTTACTTCCGAACCAGGAATATTCATAATGCGTTCCTGGGCGTGATACATCATGGTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAAt                                                           <  3:24755/149‑1 (MQ=255)                                                                                        cAGTTACTTCCGAACCAGGAATATTCATAATGCGTTCCTGGGCGTGATACATCATGGTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATg                                                          >  4:456698/1‑149 (MQ=255)                                                                                               ttCCGAACCAGGAATATTCATAATGCGTTCCTGGGCGTGATACATCATGGTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGtt                                                   <  2:56498/149‑1 (MQ=255)                                                                                                   gAACCAGGAATATTCATAATGCGTTCCTGGGCGTGATACATCATGGTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCg                                               >  1:3860/1‑149 (MQ=255)                                                                                                    aaCCAGGAATATTCATAATGCGTTCCTGGGCGTGATACATCATGGTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCgg                                              <  2:222590/149‑1 (MQ=255)                                                                                                        aGGAATATTCATAATGCGTTCCTGGGCGTGATACATCATGGTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATc                                                                 >  4:85044/1‑126 (MQ=255)                                                                                                        aGGAATATTCATAATGCGTTCCTGGGCGTGATACATCATGGTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATc                                                                 <  3:85044/126‑1 (MQ=255)                                                                                                        aGGAATATTCATAATGCGTTCCTGGGCGTGATACATCATGGTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCAt                                          >  1:108298/1‑149 (MQ=255)                                                                                                         ggAATATTCATAATGCGTTCCTGGGCGTGATACATCATGGTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATc                                         <  2:174593/149‑1 (MQ=255)                                                                                                                  aTAATGCGTTCCTGGGCGTGATACATCATGGTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCt                                >  1:116717/1‑149 (MQ=255)                                                                                                                    aaTGCGTTCCTGGGCGTGATACATCATGGTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAAt                                                                                     >  2:26787/1‑94 (MQ=255)                                                                                                                    aaTGCGTTCCTGGGCGTGATACATCATGGTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAAt                                                                                     <  1:26787/94‑1 (MQ=255)                                                                                                                     aTGCGTTCCTGGGCGTGATACATCATGGTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTCCTCGGGCATCTGTTTCGCTGAt                             <  2:684880/149‑1 (MQ=255)                                                                                                                                ggCGTGATACATCATGGTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCg                  >  3:389557/1‑149 (MQ=255)                                                                                                                                 gCGTGATACATCATGGTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCgg                 <  4:633138/149‑1 (MQ=255)                                                                                                                                   gTGATACATCATGGTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTc               >  4:85987/1‑149 (MQ=255)                                                                                                                                         catcatGGTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACcgc         >  3:477919/1‑149 (MQ=255)                                                                                                                                          atcatGGTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACcgcg        >  3:195090/1‑149 (MQ=255)                                                                                                                                           tcatGGTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCga       >  4:111705/1‑149 (MQ=255)                                                                                                                                             atGGTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCgaga     <  2:3860/149‑1 (MQ=255)                                                                                                                                                gTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGagagag  >  1:72501/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                 CCCAGGCCAGCTGCGCGTAACCGCCAATCATATGCGTTGGTTTCGGGCAAACGCGGGTAACCGAGTTATGAATGCCGCCGCGCATGCCAGTTACTTCCGAACCAGGAATATTCATAATGCGTTCCTGGGCGTGATACATCATGGTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAG  >  NC_000913/1538972‑1539263

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