New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1541977 =	36 (0.640)	3 (0.050)	3/272	5.6	6.2%	coding (614/3741 nt)	narZ	nitrate reductase Z subunit alpha
	?	NC_000913	= 4406770	55 (0.990)					coding (117/2442 nt)	rnr	RNase R
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GATGACAACGGAGGAACCGAGATGTCACAAGATCCTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCGACCTGGTGA  >  NC_000913/4406633‑4406909                                                                                                                                          |                                                                                                                                            gaTGACAACGGAGGAACCGAGATGTCACAAGATCCTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCt                                                                                                                                           <  2:359892/140‑1 (MQ=255)      cAACGGAGGAACCGAGATGTCACAAGATCCTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCgg                                                                                                                                      <  1:647148/140‑1 (MQ=255)           gaggaACCGAGATGTCACAAGATCCTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGaa                                                                                                                                 <  2:651005/140‑1 (MQ=255)            aggaACCGAGATGTCACAAGATCCTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCgg                                                                                                                                      <  2:500534/134‑1 (MQ=255)            aggaACCGAGATGTCACAAGATCCTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCgg                                                                                                                                      >  1:500534/1‑134 (MQ=255)            aggaACCGAGATGTCACAAGATCCTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGAAc                                                                                                                                >  1:891047/1‑140 (MQ=255)            aggaACCGAGATGTCACAAGATCCTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGAAc                                                                                                                                >  1:650458/1‑140 (MQ=255)            aggaACCGAGATGTCACAAGATCCTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGAAc                                                                                                                                >  1:355761/1‑140 (MQ=255)              gaACCGAGATGTCACAAGATCCTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGAACTg                                                                                                                              <  2:795080/140‑1 (MQ=255)                           acaAGATCCTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCg                                                                                                                 >  1:449746/1‑140 (MQ=255)                           acaAGATCCTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCg                                                                                                                 >  1:821827/1‑140 (MQ=255)                                aTCCTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAagag                                                                                                            <  2:367745/140‑1 (MQ=255)                                   cTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGagcag                                                                                                         >  2:737359/1‑140 (MQ=255)                                   cTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGagcag                                                                                                         >  1:460304/1‑140 (MQ=255)                                       ccAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTg                                                                                                     <  2:688025/140‑1 (MQ=255)                                           gAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAgg                                                                                                 >  1:94310/1‑140 (MQ=255)                                            aaCGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGc                                                                                                <  2:891047/140‑1 (MQ=255)                                                cgAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGc                                                                                            >  2:853169/1‑140 (MQ=255)                                                    gCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTcg                                                                                        >  2:408140/1‑140 (MQ=255)                                                        aaaaaTACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTcgccgc                                                                                    <  2:460304/140‑1 (MQ=255)                                                             tACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCt                                                                                              >  2:184119/1‑125 (MQ=255)                                                             tACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCt                                                                                              <  1:184119/125‑1 (MQ=255)                                                                gcgAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTgcgcgc                                                                            >  2:110246/1‑140 (MQ=255)                                                                    atcccatcccTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATg                                                                        <  1:383334/140‑1 (MQ=255)                                                                        catcccTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAgc                                                                    >  2:875157/1‑140 (MQ=255)                                                                         atcccTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAgcg                                                                   <  1:737359/140‑1 (MQ=255)                                                                         atcccTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCACCGCCTGCGCGCGATGGAgcg                                                                   <  2:897817/140‑1 (MQ=255)                                                                          tcccTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAgcgc                                                                  >  2:641056/1‑140 (MQ=255)                                                                          tcccTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCAGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAgcgc                                                                  >  2:738759/1‑140 (MQ=255)                                                                            ccTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGa                                                                <  1:408140/140‑1 (MQ=255)                                                                              tAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGa                                                                          <  2:440749/128‑1 (MQ=255)                                                                              tAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGa                                                                          >  1:440749/1‑128 (MQ=255)                                                                                   gggAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCaa                                                         >  1:119839/1‑140 (MQ=255)                                                                                     gAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCaa                                                         <  1:151014/138‑1 (MQ=255)                                                                                     gAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCaa                                                         >  2:151014/1‑138 (MQ=255)                                                                                     gAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACt                                                       <  1:853169/140‑1 (MQ=255)                                                                                           aTCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCtt                                                 >  2:517296/1‑140 (MQ=255)                                                                                           aTCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCtt                                                 <  1:600401/140‑1 (MQ=255)                                                                                                 gAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACtcg                                           >  1:27141/1‑140 (MQ=255)                                                                                                         aaCCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTcgtcgt                                       <  1:890946/136‑1 (MQ=255)                                                                                                         aaCCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTcgtcgt                                       >  2:890946/1‑136 (MQ=255)                                                                                                               aCGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTcgtcgt                                       >  2:183728/1‑130 (MQ=255)                                                                                                               aCGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTcgtcgt                                       <  1:183728/130‑1 (MQ=255)                                                                                                                cGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATgc                            >  1:856997/1‑140 (MQ=255)                                                                                                                        aCCGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCgg                    <  2:829755/140‑1 (MQ=255)                                                                                                                         ccGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGa                   <  2:94310/140‑1 (MQ=255)                                                                                                                         ccGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGa                   >  2:215089/1‑140 (MQ=255)                                                                                                                               aGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCt             <  1:641056/140‑1 (MQ=255)                                                                                                                               aGCAGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCt             <  1:738759/140‑1 (MQ=255)                                                                                                                                 ccGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCg           >  1:680315/1‑140 (MQ=255)                                                                                                                                   gTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCGAc         <  1:215089/140‑1 (MQ=255)                                                                                                                                      atgaGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCGACCTg      <  2:27141/140‑1 (MQ=255)                                                                                                                                         aGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCGACCTGGTg   <  2:856997/140‑1 (MQ=255)                                                                                                                                         aGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCGACCTGGTg   <  2:821827/140‑1 (MQ=255)                                                                                                                                          gCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCGACCTGGTGa  >  2:335453/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                            GATGACAACGGAGGAACCGAGATGTCACAAGATCCTTTCCAGGAACGCGAAGCTGAAAAATACGCGAATCCCATCCCTAGTCGGGAATTTATCCTCGAACATTTAACCAAACGTGAAAAACCGGCCAGCCGTGATGAGCTGGCGGTAGAACTGCACATTGAAGGCGAAGAGCAGCTTGAAGGCCTGCGTCGCCGCCTGCGCGCGATGGAGCGCGATGGTCAACTGGTCTTCACTCGTCGTCAGTGCTATGCGCTGCCGGAACGCCTCGACCTGGTGA  >  NC_000913/4406633‑4406909

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