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

GTTTCGGGCAAACGCGGGTAACCGAGTTATGAATGCCGCCGCGCATGCCAGTTACTTCCGAACCAGGAATATTCATAATGCGTTCCTGGGCGTGATACATCATGGTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATG  >  NC_000913/1539011‑1539303                                                                                                                                                   |                                                                                                                                                   gtctCGGGAAAACGCGGGTAACCGAGTTATGAATGCCGCCGCGCATGCCAGTTACTTCCGATCCAGGAATATTCATAATGCGTTCCTGGGAGTGATACATCATGGTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTc                                                                                                                                                  <  4:420707/146‑1 (MQ=255) gctctctcgcggtggaccgattgtctggatcagcgaaacagatgcccgagaactgaccattgtcgataacgactgggtggaagtattcaacgccaacggcgcgctgacTGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTc                                                                                                                                                  >  1:6293‑M1/109‑149 (MQ=255)  tttCGGGCAAACGCGGGTAACCGAGTTATGAATGCCGCCGCGCATGCCAGTTACTTCCGAACCAGGAATATTCATAATGCGTTCCTGGGCGTGATACATCATGGTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCa                                                                                                                                                 >  4:194163/1‑149 (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                                                                                                                                  <  2:165942/149‑1 (MQ=255)                 ggTAACCGAGTTATGAATGCCGCCGCGCATGCCAGTTACTTCCGAACCAGGAATATTCATAATGCGTTCCTGGGCGTGATACATCATGGTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGt                                                                                                                                  <  4:337633/149‑1 (MQ=255)                  gTAACCGAGTTATGAATGCCGCCGCGCATGCCAGTTACTTCCGAACCAGGAATATTCATAATGCGTTCCTGGGCGTGATACATCATGGTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGtt                                                                                                                                 <  2:184543/149‑1 (MQ=255)                  gTAACCGAGTGATGAATGCCCCCGCGCATGCCTGTTACTTCCGAACCAGGAATATTCATAATGCGTTCCTAGGCGTGATACATCATGGTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCGGTCAGCGCGCCGTTGGCGtt                                                                                                                                 >  4:240142/1‑149 (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                                                                                                    >  3:598523/1‑149 (MQ=255)                                               gCCAGTTACTTCCGAACCAGGAATATTCATAATGCGTTCCTGGGCGTGATACATCATGGTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACaa                                                                                                    <  2:87819/149‑1 (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                                                                                                         <  3:85044/126‑1 (MQ=255)                                                                 aGGAATATTCATAATGCGTTCCTGGGCGTGATACATCATGGTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATc                                                                                                         >  4:85044/1‑126 (MQ=255)                                                                 aGGAATATTCATAATGCGTTCCTGGGCGTGATACATCATGGTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCAt                                                                                  >  1:108298/1‑149 (MQ=255)                                                                  ggAATATTCATAATGCGTTCCTGGGCGTGATACATCATGGTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATc                                                                                 <  2:174593/149‑1 (MQ=255)                                                                           aTAATGCGTTCCTGGGCGTGATACATCATGGTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCt                                                                        >  1:116717/1‑149 (MQ=255)                                                                             aaTGCGTTCCTGGGCGTGATACATCATGGTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAAt                                                                                                                             <  1:26787/94‑1 (MQ=255)                                                                             aaTGCGTTCCTGGGCGTGATACATCATGGTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAAt                                                                                                                             >  2:26787/1‑94 (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)                                                                                                          gactgcccgcgcggtggtcagccaacgtgtaccgccgggcaGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGc                                         >  1:6293‑M2/42‑149 (MQ=255)                                                                                                             tgcccgcgcggtggtcagccaacgtgtaccgccgggcaGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTc                                      >  2:321221‑M2/39‑149 (MQ=255)                                                                                                              GCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCa                                     <  2:643206‑M1/149‑1 (MQ=255)                                                                                                               cccgcgcggtggtcagccaacgtgtaccgccgggcaGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAg                                    >  4:366665‑M2/37‑149 (MQ=255)                                                                                                                        tacacttGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAg                           >  3:64343/6‑149 (MQ=255)                                                                                                                        acacGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAg                           <  2:108298/149‑1 (MQ=255)                                                                                                                          acGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGt                         <  4:660036/149‑1 (MQ=255)                                                                                                                            gTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGttt                       >  1:185244/1‑149 (MQ=255)                                                                                                                               ggCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGc                                   >  1:614300/1‑134 (MQ=255)                                                                                                                               ggCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGc                                   <  2:614300/134‑1 (MQ=255)                                                                                                                                    accaccGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTa               >  2:347687/1‑149 (MQ=255)                                                                                                                                    accaccGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTa               <  4:64343/149‑1 (MQ=255)                                                                                                                                       accGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGt            >  3:585678/1‑149 (MQ=255)                                                                                                                                       accGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGt            >  1:12353/1‑149 (MQ=255)                                                                                                                                               ggcaGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGaa    <  2:230819‑M2/145‑1 (MQ=255)                                                                                                                                                gcaGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAAt   <  4:598523‑M2/146‑1 (MQ=255)                                                                                                                                                 caGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATg  <  2:12353‑M2/147‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                   GTTTCGGGCAAACGCGGGTAACCGAGTTATGAATGCCGCCGCGCATGCCAGTTACTTCCGAACCAGGAATATTCATAATGCGTTCCTGGGCGTGATACATCATGGTCATGCCCGGCGGTACACGTTGGCTGACCACCGCGCGGGCAGTCAGCGCGCCGTTGGCGTTGAATACTTCCACCCAGTCGTTATCGACAATGGTCAGTTCTCGGGCATCTGTTTCGCTGATCCAGACAATCGGTCCACCGCGAGAGAGCGTCAGCATTAGCAGGTTTTCACTGTAGGTTGAGTGAATG  >  NC_000913/1539011‑1539303

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