New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 379962	46 (0.880)	4 (0.080)	4/256	NT	8.1%	coding (259/1107 nt)	gbcB	glycine‑betaine dioxygenase subunit
	?	NC_002947	= 380031	47 (0.940)					coding (328/1107 nt)	gbcB	glycine‑betaine dioxygenase subunit

AGATCGAAGGCAAGCCCGTGATGCGGTCCTACACCATCTCCAGTTCGCCGTCAGTGCCCTACAGCTTCTCGATCACCGTCAAGCGCGTGCCGGGTGGCCTGGTGTCGAACTTCCTGCATGACACCATGCACGAAGGCGCCGAGCTGCCGGTGCACGGCCCGGTGGGGCTGTTCAACGCCATCGACTTCCCTGCTGGCAAGGTGCTGTACCTGTCGGGCGGCGTCGGTATTACCCCAGTGATGTCGATGGCGCGCTGGTTCTACGACACCAATGCCAACGTCGACATGGTGTTCG  >  NC_002947/379885‑380178                                                                                                                                                   |                                                                                                                                                    agaTCGAAGGCAAGCCCGTGATGCGGTCCTACACCATCTCCAGTTCGCCGTCAGTGCCCTACAGCTTCTCGATCACCGTCAAGCGCGTGCCGGGTGGCCTGGTGTCGAACTTCCTGCATGACACCATGCACGAAGGCGCCGAGCTGCcg                                                                                                                                                   >  4:282262‑M2/1‑147 (MQ=255)    tCGAAGGCAAGCCCGTGATGCGGTCCTACACCATCTCCAGTTCGCCGTCAGTGCCCTACAGCTTCTCGATCACCGTCAAGCGCGTGCCGGGTGGCCTGGTGTCGAACTTCCTGCATGACACCATGCACGAAGGCGCCGAGCTGCcggtg                                                                                                                                                <  2:238366‑M2/149‑6 (MQ=255)       aaGGCAAGCCCGTGATGCGGTCCTACACCATCTCCAGTTCGCCGTCAGTGCCCTACAGCTTCTCGATCACCGTCAAGCGCGTGCCGGGTGGCCTGGTGTCGAACTTCCTGCATGACACCATGCACGAAGGCGCCGAGCTGCcg                                                                                                                                                   >  6:248003‑M2/1‑141 (MQ=255)       aaGGCAAGCCCGTGATGCGGTCCTACACCATCTCCAGTTCGCCGTCAGTGCCCTACAGCTTCTCGATCACCGTCAAGCGCGTGCCGGGTGGCCTGGTGTCGAACTTCCTGCATGACACCATGCACGAAGGCGCCGAGCTGCcg                                                                                                                                                   <  5:248003‑M2/143‑3 (MQ=255)                 cGTGATGCGGTCCTACACCATCTCCAGTTCGCCGTCAGTGCCCTACAGCTTCTCGATCACCGTCAAGCGCGTGCCGGGTGGCCTGGTGTCGAACTTCCTGCATGACACCATGCACGAAGGCGCCGAGCTGCCGGTGCACGGCCCGGTgg                                                                                                                                   <  1:238366/149‑1 (MQ=255)                             cTACACCATCTCCAGTTCGCCGTCAGTGCCCTACAGCTTCTCGATCACCGTCAAGCGCGTGCCGGGTGGCCTGGTGTCGAACTTCCTGCATGACACCATGCACGAAGGCGCCGAGCTGCCGGTGCACGGCCCGGTGGGGCTGTTCAAc                                                                                                                        <  8:126599/148‑1 (MQ=255)                                cacCATCTCCAGTTCGCCGTCAGTGCCCTACAGCTTCTCGATCACCGTCAAGCGCGTGCCGGGTGGCCTGGTGTCGAACTTCCTGCATGACACCATGCACGAAGGCGCCGAGCTGCcggt                                                                                                                                                 <  7:249114‑M2/120‑5 (MQ=255)                                cacCATCTCCAGTTCGCCGTCAGTGCCCTACAGCTTCTCGATCACCGTCAAGCGCGTGCCGGGTGGCCTGGTGTCGAACTTCCTGCATGACACCATGCACGAAGGCGCCGAGCTGCcggt                                                                                                                                                 >  8:249114‑M2/1‑116 (MQ=255)                                   cATCTCCAGTTCGCCGTCAGTGCCCTACAGCTTCTCGATCACCGTCAAGCGCGTGCCGGGTGGCCTGGTGTCGAACTTCCTGCATGACACCATGCACGAAGGCGCCGAGCTGCCGGTGCACGGCCCGGTGGGGCTGTTCAACGCCATc                                                                                                                  >  4:30840/1‑148 (MQ=255)                                    aTCTCCAGTTCGCCGTCAGTGCCCTACAGCTTCTCGATCACCGTCAAGCGCGTGCCGGGTGGCCTGGTGTCGAACTTCCTGCATGACACCATGCACGAAGGCGCCGAGCTGCCGGTGCACGGCCCGGTGGGGCTGTTCAACGCCATCGa                                                                                                                >  6:249651/1‑149 (MQ=255)                                     tctcCAGTTCGCCGTCTGTGCCCTACAGCTTCTCGATCACCGTCAAGCGCGTGCCGGGTGGCCTGGTGTCGAACTTCCTGCATGACACCATGCACGAAGGCGCCGAGCTGCCGGTGCACGGCCCGGTGGGGCTGTTCAACGCCATCGAc                                                                                                               >  8:170649/1‑149 (MQ=255)                                          gttctgCCGTCAGTGCCCTACAGCTTCTCGATCACCGTCAAGCGCGTGCCGGGTGGCCTGGTGTCGAACTTCCTGCATGACACCATGCACGAAGGCGCCGAGCTGCCGGTGCACGGCCCGGTGGGGCTGTTCAACGCCATCGACTTccc                                                                                                          <  7:18582/144‑1 (MQ=255)                                           gTTCGCCGTCAGTGCCCTACAGCTTCTCGATCACCGTCAAGCGCGTGCCGGGTGGCCTGGTGTCGAACTTCCTGCATGACACCATGCACGAAGGCGCCGAGCTGCCGGTGCACGGCCCGGTGGGGCTGTTCAACGCCATCGACTTCCct                                                                                                         >  2:315327/1‑149 (MQ=255)                                              cGCCGTCAGTGCCCTACAGCTTCTCGATCACCGTCAAGCGCGTGCCGGGTGGCCTGGTGTCGAACTTCCTGCATGACACCATGCACGAAGGCGCCGAGCTGCCGGTGCACGGCCCGGTGGGGCTGTTCAACGCCATCGACTTCCctgct                                                                                                      >  5:95084/1‑149 (MQ=255)                                               gCCGTCAGTGCCCTACAGCTTCTCGATCACCGTCAAGCGCGTGCCGGGTGGCCTGGTGTCGAACTTCCTGCATGACACCATGCACGAAGGCGCCGAGCTGCCGGTGCACGGCCCGGTGGGGCTGTTCAACGCCATCGACTTCCctgctg                                                                                                     <  3:135036/149‑1 (MQ=255)                                                ccGTCAGTGCCCTACAGCTTCTCGATCACCGTCAAGCGCGTGCCGGGTGGCCTGGTGTCGAACTTCCTGCATGACCCCATGCACGAAGGCGCCGAGCTGCCGGTGCACGGCCCGGTGGGGCTGTTCAACGCCATCGACTTccc                                                                                                          >  8:18582/1‑143 (MQ=255)                                                 cGTCAGTGCCCTACAGCTTCTCGATCACCGTCAAGCGCGTGCCGGGTGGCCTGGTGTCGAACTTCCTGCATGACACCATGCACGAAGGCGCCGAGCTGCCGGTGCACGGCCCGGTGGGGCTGTTCAACGCCATCGACTTCCCTGCTgg                                                                                                    >  5:281398/1‑148 (MQ=255)                                                   tCAGTGCCCTACAGCTTCTCGATCACCGTCAAGCGCGTGCCGGGTGGCCTGGTGTCGAACTTCCTGCATGACACCATGCACGAAGGCGCCGAGCTGCCGGTGCACGGCCCGGTGGGGCTGTTCAACGCCATCGACTTCCCTGCTAGCaa                                                                                                 <  5:249651/149‑1 (MQ=255)                                                      gTGCCCTACAGCTTCTCGATCACCGTCAAGCGCGTGCCGGGTGGCCTGGTGTCGAACTTCCTGCATGACACCATGCACGAAGGCGCCGAGCTGCCGGTGCACGGCCCGGTGGGGCTGTTCAACGCCATCGACTTCCCTGCTGGCAAGGt                                                                                              <  8:67649/149‑1 (MQ=255)                                                       tGCCCTACAGCTTCTCGATCACCGTCAAGCGCGTGCCGGGTGGCCTGGTGTCGAACTTCCTGCATGACACCATGCACGAAGGCGCCGAGCTGCcggtgatcgagaagctgtagggcactgacggcgaactggagatggtgtaggaccg                                                                                              <  4:141936‑M2/148‑56 (MQ=255)                                                           cTACAGCTTCTCGATCACCGTCAAGCGCGTGCCGGGTGGCCTGGTGTCGAACTTCCTGCATGACACCATGCACGAAGGCGCCGAGCTGCCGGTGCACGGCCCGGTGGGGCTGTTCAACGCCATCGACTTCCCTGCTGGCAAGGTGCTg                                                                                          >  1:246792/1‑148 (MQ=255)                                                                  ttCTCGATCACCGTCAAGCGCGTGCCGGGTGGCCTGGTGTCGAACTTCCTGCATGACACCATGCACGAAGGCGCCGAGCTGCCGGTGCACGGCCCGGTGGGGCTGTTc                                                                                                                           >  3:304041/1‑108 (MQ=255)                                                                  ttCTCGATCACCGTCAAGCGCGTGCCGGGTGGCCTGGTGTCGAACTTCCTGCATGACACCATGCACGAAGGCGCCGAGCTGCCGGTGCACGGCCCGGTGGGGCTGTTc                                                                                                                           <  4:304041/108‑1 (MQ=255)                                                                  ttCTCGATCACCGTCAAGCGCGTGCCGGGTGGCCTGGTGTCGAACTTCCTGCATGACACCATGCACGAAGGCGCCGAGCTGCCGGTGCACGGCCCGGTGGGGCTGTTCAACGCCATCGACTTCCCTGCTGGCAAGGTGCTGTACCTGTc                                                                                  >  3:247741/1‑149 (MQ=255)                                                                    ctcGATCACCGTCAAGCGCGTGCCGGGTGGCCTGGTGTCGAACTTCCTGCATGACACCATGCACGAAGGCGCCGAGCTGCCGGTGCACGGCCCGGTGGGGCTGTTCAACGCCATCGACTTCCCTGCTGGCAAGGTGCTGTACCTGTc                                                                                  <  1:315327/147‑1 (MQ=255)                                                                       gATCACCGTCAAGCGCGTGCCGGGTGGCCTGGTGTCGAACTTCCTGCATGACACCATGCACGAAGGCGCCGAGCTGCCGGTGCACGGCCCGGTGGGGCTGTTCAACGCCATCGACTTCCCTGCTGGCAAGGTGCTGTACCTGTCGggc                                                                              <  6:95084/148‑1 (MQ=255)                                                                         tCACCGTCAAGCGCGTGCCGGGTGGCCTGGTGTCGAACTTCCTGCATGACACCATGCACGAAGGCGCCGAGCTGCcggtgatcgagaagctgtagggcactgacggcgaactggagatggtgtaggaccgcatcacgggcttgccttc                                                                            <  6:318063‑M2/148‑74 (MQ=255)                                                                               tCAAGCGCGTGCCGGGTGGCCTGGTGTCGAACTTCCTGCATGACACCATGCACGAAGGCGCCGAGCTGCCGGTGCACGGCCCGGTGGGGCTGTTCAACGCCATCGACTTCCCTGCTGGCAAGGTGCTGTACCTGTCGGGCGGCGTCGGt                                                                     <  5:109626/149‑1 (MQ=255)                                                                                 aaGCGCGTGCCGGGTGGCCTGGTGTCGAACTTCCTGCATGACACCATGCACGAAGGCGCCGAGCTGCCGGTGCACGGCCCGGTGGGGCTGTTCAACGCCATCGACTTCCCTGCTGGCAAGGTGCTGTACCTGTCGGGCGGCGTCGGTa                                                                    <  3:30840/148‑1 (MQ=255)                                                                                 aaGCGCGTGCCGGGTGGCCTGGTGTCGAACTTCCTGCATGACACCATGCACGAAGGCGCCGAGCTGCCGGTGCACGGCCCGGTGGGGCTGTTCAACGCCATCGACTTCCCTGCTGGCAAGGTGCTGTACCTGTCGGGCGGCGTCGGTAt                                                                   >  1:245323/1‑149 (MQ=255)                                                                                 aaGCGCGTGCCGGGTGGCCTGGTGTCGAACTTCCTGCATGACACCATGCACGAAGGCGCCGAGCTGCCGGTGCACGGCCCGGTGGGGCTGTTCAACGCCATCGACTTCCCTGCTGGCAAGGTGCTGTACCTGTCGGGCGGCGTCGGTAt                                                                   >  5:21067/1‑149 (MQ=255)                                                                                   gcgcGTGCCGGGTGGCCTGGTGTCGAACTTCCTGCATGACACCATGCACGAAGGCGCCGAGCTGCCGGTGCACGGCCCGGTGGGGCTGTTCAACGCCATCGACTTCCCTGCTGGCAAGGTGCTGTACCTGTCGGGCGGCGTCGGTATTa                                                                 <  2:246792/149‑1 (MQ=255)                                                                                     gcgTGCCGGGTGGCCTGGTGTCGAACTTCCTGCATGACACCATGCACGAAGGCGCCGAGCTGCCGGTGCACGGCCCGGTGGGGCTGTTCAACGCCATCGACTTCCCTGCTGGCAAGGTGCTGTACCTGTCGGGCGGCGTCGGTATTAcc                                                               >  8:45868/1‑149 (MQ=255)                                                                                           cGGGTGGCCTGGTGTCGAACTTCCTGCATGACACCATGCACGAAGGCGCCGAGCTGCCGGTGCACGGCCCGGTGGGGCTGTTCAACGCCATCGACTTCCCTGCTGGCAAGGTGCTGTACCTGTCGGGCGGCGTCGGTATTACCCCAGTg                                                         <  4:155932/149‑1 (MQ=255)                                                                                             ggTTGCCTGGTGTCGAACTTCCTGCATGACACCATGCACGAAGGCGCCGAGCTGCCGGTGCACGGCCCGGTGGGGCTGTTCAACGCCATCGACTTCCCTGCTGGCAAGGTGCTGTACCTGTCGGGCGGCGTCGGTATTACCCCAGTGAt                                                       <  7:45868/149‑1 (MQ=255)                                                                                                      gtgtCGAACTTCCTGCATGACACCATGCACGAAGGCGCCGAGCTGCCGGTGCACGGCCCGGTGGGGCTGTTCAACGCCATCGACTTCCCTGCTGGCAAGGTGCTGTACCTGTCGGGCGGCGTCGGTATTACCCCAGTGATGTCGATg                                                <  6:21067/147‑1 (MQ=255)                                                                                                         tCGAACTTCCTGCATGACACCATGCACGAAGGCGCCGAGCTGCCGGTGCACGGCCCGGTGGGGCTGTTCAACGCCATCGACTTCCCTGCTGGCAAGGTGCTGTACCTGTCGGGCGGCGTCGGTATTACCCCAGTGATGTCGATGgcgcg                                           >  4:48977/1‑149 (MQ=255)                                                                                                           gAACTTCCTGCATGACACCATGCACGAAGGCGCCGAGCTGCCGGTGCACGGCCCGGTGGGGCTGTTCAACGCCATCGACTTCCCTGCTGGCAAGGTg                                                                                             <  2:206595/97‑1 (MQ=255)                                                                                                           gAACTTCCTGCATGACACCATGCACGAAGGCGCCGAGCTGCCGGTGCACGGCCCGGTGGGGCTGTTCAACGCCATCGACTTCCCTGCTGGCAAGGTg                                                                                             >  1:206595/1‑97 (MQ=255)                                                                                                           gAACTTCCTGCATGACACCATGCACGAAGGCGCCGAGCTGCCGGTGCACGGCCCGGTGGGGCTGTTCAACGCCATCGACTTCCCTGCTGGCAAGGTGCTGTACCTGTCGGGCGGCGTCGGTATTACCCCAGTGATGTCGATGgcgcgc                                          <  7:170649/148‑1 (MQ=255)                                                                                                              cTTCCTGCATGACACCATGCACGAAGGCGCCGAGCTGCcggtgatcgagaagctgtagggcactgacggcgaactggagatggtgtaggaccgcatcacgggcttgccttcgatctccagctccagggtaacgaactgccctggcttga                                      >  3:141936‑M2/1‑38 (MQ=255)                                                                                                                  cTGCATGACACCATGCACGAAGGCGCCGAGCTGCCGGTGCACGGCCCGGTGGGGCTGTTCAACGCCATCGACTTCCCTGCTGGCAAGGTGCTGTACCTGTCGGGCGGCGTCGGTATTACCCCAGTGATGTCGATGGCGCGCTGGTTCt                                   >  2:122742/1‑148 (MQ=255)                                                                                                                   tGCATGACACCATGCACGAAGGCGCCGAGCTGCCGGTGCACGGCCCGGTGGGGCTGTTCAACGCCATCGACTTCCCTGCTGGCaa                                                                                                 <  1:201492/85‑1 (MQ=255)                                                                                                                   tGCATGACACCATGCACGAAGGCGCCGAGCTGCCGGTGCACGGCCCGGTGGGGCTGTTCAACGCCATCGACTTCCCTGCTGGCaa                                                                                                 >  2:201492/1‑85 (MQ=255)                                                                                                                   tGCATGACACCATGCACGAAGGCGCCGAGCTGCCGGTGCACGGCCCGGTGGGGCTGTTCAACGCCATCGACTTCCCTGCTGGCAAGGTGCTGTACCTGTCGGGCGGCGTCGGTATTACCCCAGTGATGTCGATGGCGCGCTGGTTCTAc                                 >  1:59778/1‑149 (MQ=255)                                                                                                                   tGCATGACACCATGCACGAAGGCGCCGAGCTGCCGGTGCACGGCCCGGTGGGGCTGTTCAACGCCATCGACTTCCCTGCTGGCAAGGTGCTGTACCTGTCGGGCGGCGTCGGTATTACCCCAGTGATGTCGATGGCGCGCTGGTTCTAc                                 >  3:34471/1‑149 (MQ=255)                                                                                                                     cATGACACCATGCACGAAGGCGCCGAGCTGCcggtgatcgagaagctgtagggcactgacggcgaactggagatggtgtaggaccgcatcacgggcttgccttcgatctccagctccagggtaacgaactgccctggcttgaagaaga                                >  5:318063‑M2/1‑31 (MQ=255)                                                                                                                     cATGACACCATGCACGAAGGCGCCGAGCTGCCGGTGCACGGCCCGGTGGGGCTGTTCAACGCCATCGACTTCCCTGCTGGCAAGGTGCTGTACCTGTCGGGCGGCGTCGGTATTACCCCAGTGATGTCGATGGCGCGCTGGTTCTAc                                 <  4:34471/147‑1 (MQ=255)                                                                                                                          cacCATGCACGAAGGCGCCGAGCTGCCGGTGCACGGCCCGGTGGGGCTGTTCAACGCCATCGACTTCCCTGCTGGCAAGGTGCTGTACCTGTCGGGCGGCGTCGGTATTACCCCAGTGATGTCGATGGCGCGCTGGTTCTa                                  <  1:253894/141‑1 (MQ=255)                                                                                                                          cacCATGCACGAAGGCGCCGAGCTGCCGGTGCACGGCCCGGTGGGGCTGTTCAACGCCATCGACTTCCCTGCTGGCAAGGTGCTGTACCTGTCGGGCGGCGTCGGTATTACCCCAGTGATGTCGATGGCGCGCTGGTTCTa                                  >  2:253894/1‑141 (MQ=255)                                                                                                                          cacCATGCACGAAGGCGCCGAGCTGCCGGTGCACGGCCCGGTGGGGCTGTTCAACGCAATCGACTTCCCTGCTGGCAAGGTGCTGTACCTGTCGGGCGGCGTCGGTATTACCCCAGTGATGTCGATGGCGCGCTGGTTCTACGACAcc                           >  4:151161/1‑148 (MQ=255)                                                                                                                           acCATGTACGAAGGCGCCGAGCTGCCGGTGCACGGCCCGGTGGGGCTGTTCAACGCCATCGACTTCCCTGCTGGCAAGGTGCTGTACCTGTCGGGCGGCGTCGGTATTACCCCAGTGATGTCGATGGCGCGCTgg                                       <  5:234070/135‑1 (MQ=255)                                                                                                                           acCATGAACGAAGGCGCCGAGCTGCCGGTGCACGGCCCGGTGGGGCTGTTCAACGCCATCGACTTCCCTGCTGGCAAGGTGCTGTACCTGTCGGGCGGCGTCGGTATTACCCCAGTGATGTCGATGGCGCGCTgg                                       >  6:234070/1‑135 (MQ=255)                                                                                                                                 cACGAAGGCGCCGAGCTGCCGGTGCACGGCCCGGTGGGGCTGTTCAACGCCATCGACTTCCCTGCTGGCAAGGTGCTGTACCTGTCGGGCGGCGTCGGTATTACCCCAGTGATGTCGATGGCGCGCTGGTTCTACGACACCAATGCCaa                   <  3:48977/149‑1 (MQ=255)                                                                                                                                         cgcCGAGCTGCCGGTGCACGGCCCGGTGGGGCTGTTCAACGCCATCGACTTCCCTGCTGGCAAGGTGCTGTACCTGTCGGGCGGCGTCGGTATTACCCCAGTGATGTCGATGGCGCGCTGGTTCTAGGACACCAATGCCAACGTCGACa           <  2:245323/149‑1 (MQ=255)                                                                                                                                           ccGAGCTGCGGGTGCACGGCCCGGTGGGGCTGTTCAACGCCATCGACTTCCCTGCTGGCAAGGTGCTGTACCTGTCGGGCGGCGTCGGTATTACCCCAGTGATGTCGATGGCGCGCTGGTTCTACGACACCAATGCCAACGTCGACATg         <  2:112978/149‑1 (MQ=255)                                                                                                                                             gAGCTGCCGGTGCACGGCCCGGTGGGGCTGTTCAACGCCATCGACTTCCCTGCTGGCAAGGTGCTGTACCTGTCGGGCGGCGTCGGTATTACCCCAGTGATGTCGATGGCGCGCTGGTTCTACGACACCAATGCCAACGTCGACATGgt       <  5:317282/149‑1 (MQ=255)                                                                                                                                                  gCCGGTGCACGGCCCGGTGGGGCTGTTCAACGCCATCGACTTCCCTGCTGGCAAGGTGCTGTACCTGTCGGGCGGCGTCGGTATTACCCCAGTGATGTCGATGGCGCGCTGGTTCTACGACACCAATGCCAACGTCGACATGGTGTTCg  >  8:25182/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                    AGATCGAAGGCAAGCCCGTGATGCGGTCCTACACCATCTCCAGTTCGCCGTCAGTGCCCTACAGCTTCTCGATCACCGTCAAGCGCGTGCCGGGTGGCCTGGTGTCGAACTTCCTGCATGACACCATGCACGAAGGCGCCGAGCTGCCGGTGCACGGCCCGGTGGGGCTGTTCAACGCCATCGACTTCCCTGCTGGCAAGGTGCTGTACCTGTCGGGCGGCGTCGGTATTACCCCAGTGATGTCGATGGCGCGCTGGTTCTACGACACCAATGCCAACGTCGACATGGTGTTCG  >  NC_002947/379885‑380178

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