New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	4464118 =	43 (0.800)	4 (0.080)	4/250	NT	8.8%	coding (1837/1974 nt)	betT‑III	choline/carnitine/betaine transporter BetTC
	?	NC_002947	4464151 =	42 (0.820)					coding (1804/1974 nt)	betT‑III	choline/carnitine/betaine transporter BetTC

CCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCGTAGTCCTGGCCCCCTTCTCGCAGATGCACTTCCGCGCGGAAATACTTGCGGGAATCGGCACCGTCTTCGGTGTCCCGCATGACGAAGCTTGGCGTATTGAAGGCCCGTGGCCGCACCTCGTAGAGGAAACGTCCTTCACCGGCATGGGAGAGCTCCAGAGATACCCGCCCATCATCACGCTCATTCACGCTCAC  >  NC_002947/4463970‑4464265                                                                                                                                                     |                                                                                                                                                    gtcgcccaggaactgcgcaagcaaggctacgaggtgagcgtgaatgagcgtgatgatgggcgggtatctctggagctctcccatgccggtgaaggacgtttcctctacgaggtgcggccacgggccttcaatacgccaagcttcgtcat                                                                                                                                                     <  5:285974‑M1/1‑1 (MQ=255) ccGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCGTAGTCCTGGCCCCCTTCTCGCAGATGCACTTCCGCGCGGAAATACtt                                                                                                                                                     <  7:92758/149‑1 (MQ=255)  cGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCGTAGTCCTGGCCCCCTTCTCGCAGATGCACTTCCGCGCGGAAATACTTg                                                                                                                                                    <  1:74546/149‑1 (MQ=255)    aaTCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCGTAGTCCTGGCCCCCTTCTCGCAGATGCACTTCCGCGCGGAAATACTTGCg                                                                                                                                                  >  3:56891/1‑149 (MQ=255)     cccaggaactgcgcaagcaaggctacgaggtgagcgtgaatgagcgtgatgatgggcgggtatctctggagctctcccatgccggtgaaggacgtttcctctacgaggtgcggccacgggccttcaatacgccaagcttcgtcaTGCgg                                                                                                                                                 >  7:161546‑M1/145‑149 (MQ=255)     cccaggaactgcgcaagcaaggctacgaggtgagcgtgaatgagcgtgatgatgggcgggtatctctggagctctcccatgccggtgaaggacgtttcctctacgaggtgcggccacgggccttcaatacgccaagcttcgtcaTGCgg                                                                                                                                                 >  3:197014‑M1/145‑149 (MQ=255)     aTCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCGTAGTCCTGGCCCCCTTCTCGCAGATGCACTTCCGCGCGGAAATACTTGCgg                                                                                                                                                 <  4:56891/149‑1 (MQ=255)         aCTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCGTAGTCCTGGCCCCCTTCTCGCAGATGCACTTCCGCGCGGAAATACTTGCGGGAAt                                                                                                                                             <  8:120292/149‑1 (MQ=255)             cAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCCCATCGTCACGGCTCGCCCCCATGATGTCGTAGTCCTGGCCCCCTGCTCGCAGATGCACTTCCGCGCGGAAATACTTGCGGGAATCgg                                                                                                                                          >  6:204976/1‑148 (MQ=255)             cAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCGTAGTCCTGGCCCCCTTCTCGCAGATGCACTTCCGCGCGGAAATACTTGCGGGAATCgg                                                                                                                                          <  6:245855/148‑1 (MQ=255)                  caccacGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGTATGTCGCCAATCACATCGTCACGGATCCACCCCATGATGTCGTAGTCCTGGCCCCCTTCTCGCAGATGCACTTACGAGCGGAAATACTTGCGGGAATCGGCACCGt                                                                                                                                    >  6:167474/1‑149 (MQ=255)                   accacGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCGTAGTCCTGGCCCCCTTCTCGCAGATGCACTTCCGCGCGGAAATACTTGCGGGAATCGGCACCGTc                                                                                                                                   >  4:87740/1‑149 (MQ=255)                     cacGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCGTAGTCCTGGCCCCCTTCTCGCAGATGCACTTCCGCGCGGAAATACTTGCGGGAATCGGCACCGTCtt                                                                                                                                 <  2:263950/149‑1 (MQ=255)                       cGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCGTAGTCCTGGCCCCCTTCTCGCAGATGCACTTCCGCGCGGAAATACTTGCGGGAATCGGCACCGTCTTCg                                                                                                                               <  2:128713/149‑1 (MQ=255)                                   gCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCGTAGTCCTGGCCCCCTTCTCGCAGATGCACTTCCGCGCGGAAATACTTGCGGGAATCGGCACCGTCTTCGGTGTCCCGCAt                                                                                                                    >  8:192367/1‑148 (MQ=255)                                    cAGGTGGCTTTCGTACTGATCCAGGATGTCGACTATCAGATCGTCACGGCTCGACCCCATGATGTGGTAGTCGTGGCCCGCTTCTCTCAGTTGCTCTTCCGCGCGGAAATACTTGCGGGATTCGGCACCGTCTTCGGTGTCCCTCATGa                                                                                                                  <  3:125589/149‑1 (MQ=255)                                            gTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCGTAGTCCTGGCCCCCTTCTCGCAGATGCACTTCCGCGCGGAAATACTTGCGGGAATCGGCACCGTCTTCGGTGTCCCGCATGACGAAGCtt                                                                                                          <  5:184137/149‑1 (MQ=255)                                                 tACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCGTAGTCCTGGCCCCCTTCTCGCAGATGCACTTCCGCGCGGAAATACTTGCGGGAATCGGCACCGTCTTCGGTGTCCCGCATGACGAAGCTTGGCGt                                                                                                     >  7:243712/1‑149 (MQ=255)                                                 tACTGATCAAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCGTAGTCCTGGCCCCCTTCTCGCAGATGCACTTCCGCGCGGAAATACTTGCGGGAATCGGCACCGTCTTCGGTGTCCCGCATGACGAAGCTTGGCg                                                                                                      >  5:317983/1‑148 (MQ=255)                                                    tGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCGTAGTCCTGGCCCCCTTCTCGCAGATGCACTTCCGCGCGGAAATACTTGCGGGAATCGGCACCGTCTTCGGTGTCCGGCAGGACGACGCTTGGCGTAtt                                                                                                  >  6:104034/1‑149 (MQ=255)                                                    tGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCGTAGTCCTGGCCCCCTTCTCGCAGATGCACTTCCGCGCGGAAATACTTGCGGGAATCGGCACCGTCTTCGGTGTCCCGCATGACGAAGCTTGGCGTAtt                                                                                                  >  1:14191/1‑149 (MQ=255)                                                       tCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCGTAGTCCTGGCCCCCTTCTCGCAGATGCACTTCCGCGCGGAAATACTTGCGGGAATCGGCACCGTCTTCGGTGTCCCGCATGACGAAGCTTGGCGTATTGaa                                                                                               <  3:120926/149‑1 (MQ=255)                                                       tCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCGTAGTCCTGGCCCCCTTCTCGCAGATGCACTTCCGCGCGGAAATACTTGCGGGAATCGGCACCGTCTTCGGTGTCCCGCATGACGAAGCTTGGCGTATTGaa                                                                                               <  6:144475/149‑1 (MQ=255)                                                         tgggcgggtatctctggagctctcccatgccggtgaaggacgtttcctctacgaggtgcggccacgggccttcaatacgccaagcttcgtcaTGCGGGAATCGGCACCGTCTTCGGTGTCCCGCATGa                                                                                                                  <  3:237528‑M1/36‑1 (MQ=255)                                                         tgggcgggtatctctggagctctcccatgccggtgaaggacgtttcctctacgaggtgcggccacgggccttcaatacgccaagcttcgtcaTGCGGGAATCGGCACCGTCTTCGGTGTCCCGCATGa                                                                                                                  >  4:237528‑M1/93‑128 (MQ=255)                                                         cAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCGTAGTCCTGGCCCCCTTCTCGCAGATGCACTTCCGCGCGGAAATACTTGCGGGAATCGGCACCGTCTTCGGTGTCCCGCATGACGAAGCTTGGc                                                                                                       <  3:192016/139‑1 (MQ=255)                                                         cAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCGTAGTCCTGGCCCCCTTCTCGCAGATGCACTTCCGCGCGGAAATACTTGCGGGAATCGGCACCGTCTTCGGTGTCCCGCATGACGAAGCTTGGc                                                                                                       >  4:192016/1‑139 (MQ=255)                                                          aGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCGTAGTCCTGGCCCCCTTCTCGCACATGCACTTCCGCGCGGAAATACTTGCGGGAATCGGCACCGTCTTCGGAGTCCCGCATGACGAAGCTTGGCGTATTGAAGGc                                                                                            >  2:324635/1‑149 (MQ=255)                                                                            tCGTCACGGCTCCACCCCATGATGTCGTAGTCCTGGCCCCCTTCTCGCAGATGCACTTCCGCGCGGAAATACTTGCGGGAATCGGCACCGTCTTCGGTGTCCCGCATGACGAAGCTTGGCGTATTGAAGGCCCGTGGCCGCACCTCGTa                                                                          <  7:192367/149‑1 (MQ=255)                                                                             cctCACGGCTCCACCCCATGATGTCGTAGTCCTGGCCCCCTTCTCGCAGATGCACTTCCGCGCGGAAATACTTGCGGGAATCGGCACCGTCTTCGGTGTCCCGCATGACGAAGCTTGGCGTATTGAAGGCCCGTGGCCGCACCTCGTag                                                                         <  4:254592/147‑1 (MQ=255)                                                                               tcccatgccggtgaaggacgtttcctctacgaggtgcggccacgggccttcaatacgccaagcttcgtcaTGCGGGAATCGGCACCGTCTTCGGTGTCCCGCATg                                                                                                                   >  7:79311‑M1/71‑105 (MQ=255)                                                                               tcccatgccggtgaaggacgtttcctctacgaggtgcggccacgggccttcaatacgccaagcttcgtcaTGCGGGAATCGGCACCGTCTTCGGTGTCCCGCATg                                                                                                                   <  8:79311‑M1/35‑1 (MQ=255)                                                                                 aCGGCTCCACCCCATGATGTCGTAGTCCTGGCACCCTTCTCGCAGATGCACTTCCGCGCGGAAATACTTGCGGGAATCGTCACCGTCTTCGTTGTCCCGCATGACGAAGCTTGGCGTATTGAAGGCCCGTGGCCGCACCTCGTAGAgg                                                                      <  5:204976/148‑1 (MQ=255)                                                                                    gCTCCACCCCATGATGTCGTAGTCCTGGCCCCCTTCTCGCAGATGCACTTCCGCGCGGAAATACTTGCGGGAATCGGCACCGTCTTCGGTGTCCCGCATGACGAAGCTTGGCGTATTGAAGGCCCGTGGCCGCACCTCGTAGAGGAAAc                                                                  <  2:27919/149‑1 (MQ=255)                                                                                              atgatgTCGTAGTCCTGGCCCCCTTCTCGCAGATGCACTTCCGCGCGGAAATACTTGCGGGAATCGGCACCGTCTTCGGTGTCCCGCATGACGAAGCTTGGCGTATTGAAGGCCCGTGGCCGCACCTCGTAGAGGAAACGTCCTTCAcc                                                        >  1:56751/1‑149 (MQ=255)                                                                                                 atgTCGTAGTCCTGGCCCCCTTCTCGCAGATGCACTTCCGCCCGCAAATACTTGCCGGAATCGGCACCGTCATCGGTGTCCCGCAAGACGAAGATTCGAGTATTGAAGGCCCGTGCCCGCACCTCGTAGAGGAAACGTCCATCACCgg                                                      >  6:285974/1‑148 (MQ=255)                                                                                                     cGTAGTCCTGGCCCCCTTCTCGCAGATGCACTTCCGCGCGGAAATACTTGCGGGAATCGGCACCGTCTTCGGTGTCCCGCATGACGAAGCTTGGCGTATTGAAGGCCCGTGGCCGCACCTCGTAGAGGAAACGTCCTTCACCGGCATgg                                                 >  5:328620/1‑149 (MQ=255)                                                                                                      gTAGTCCTGGCCCCCTTCTCGCAGATGCACTTCCGCGCGGAAATACTTGCGGGAATCGGCACCGTCTTCGGTGTCCCGCATGACGAAGCTTGGCGTATTGAAGGCACGTGGCCGCACCTCGTAGAGGAAACGTCCTTCACCGGCATgg                                                 >  7:115516/1‑148 (MQ=255)                                                                                                       tAGTCCTGGCCCCCTTCTCGCAGATGCACTTCCGCGCGGAAATACTTGCGTGAATCGGCACCGTCTTCGGTGTCCCGCATGACGAAGCTTGGCGTATTGAAGGCCCGTGGCCGCACCTCGTAGAGGAAACGTCCTTCACCGGCATGGga                                               <  6:317983/149‑1 (MQ=255)                                                                                                       tAGTCCTGGCCCCCTTCTCGCAGATGCACTTCCGCGCGGAAATACTTGCGGGAATCGGCACCGTCTTCGGTGTCCCGCATGACGAAGCTTGGCGTATTGAAGGCCCGTGGCCGCACCTCGTAGAGGAAACGTCCTTCACCGGCATgg                                                 >  3:201732/1‑147 (MQ=255)                                                                                                       tAGTCCTGGCCCCCTTCTCGCAGATGCACTTCCGCGCGGAAATACTTGCGGGAATCGGCACCGTCTTCGGTGTCCCGCATGACGAAGCTTGGCGTATTGAAGGCCCGTGGCCGCACCTCGTAGAGGAAACGTCCTTCACCGGCATGGga                                               <  3:87740/149‑1 (MQ=255)                                                                                                         gTCCTGGCCCCCTTCTCGCAGATGCACTTCCGCGCGGAAATACTTGCGGGAATCGGCACCGTCTTCGGTGTCCCGCATGACGAAGCTTGGCGTATTGAAGGCCCGTGGCCGCACCTCGTAGAGGAAACGTCCTTCACCGGCATGGgaga                                             <  7:161151/149‑1 (MQ=255)                                                                                                         gTCCTGGCCCCCTTCTCGCAGATGCACTTCCGCGCGGAAATACTTGCGGGAATCGGCACCGTCTTCGGTGTCCCGCATGACGAAGCTTGGCGTATTGAAGGCCCGTGGCCGCACCTCGTAGAGGAAACGTCCTTCACCGGCATGGgag                                              <  2:14191/148‑1 (MQ=255)                                                                                                          tCCTGGCCCCCTTCTCGCAGATGCACTTCCGCGCGGAAATACTTGCGGGAATCGGCACCGTCTTCGGTGTCCCGCATGACGAAGCTTGGCGTATTGAAGGCCCGTGGCCGCACCTCGTAGAGGAAACGTCCTTCACCGGCATGGGagag                                            >  4:166103/1‑149 (MQ=255)                                                                                                           ccTGGCCCCCTTCTCGCACATGCACTTCAGCGCGGAAATACTTGCGGGAATCGGCACCGTCTTCGGTGTCCCGCATGCCGAAGCTTGGCGTATTGAAGGCCCGTGGCCGCACCTCGTAGAGGAAACGTCCGTCACCGGCATGGGagag                                            >  2:153024/1‑148 (MQ=255)                                                                                                                 ccccTTCTCGCAGATGCACTTCCGCGCGGAAATACTTGCGGGAATCGGCACCGTCTTCGGTGTCCCGCATGACGAAGCTTGGCGTATTGAAGGCCCGTGGCCGCACCTCGTAGAGGAAACGTCCTTCACCGGCATGGGAGAGCTCCaga                                     <  2:56751/149‑1 (MQ=255)                                                                                                                        tcGCAGATGCACTTCCGCGCGGAAATACTTGCGGGAATCGGCACCGTCTTCGGTGTCCCGCATGACGAAGCTTGGCGTATTGAAGGCCCGTGGCCGCACCTCGTAGAGGAAACGTCCTTCACCGGCATGGGAGAGCTCCAGAGATAccc                              >  2:51603/1‑149 (MQ=255)                                                                                                                          tcAGATGCACTTCCGCGCGGAAATACTTGCGGGAATCGGCACCGTCTTCGGTGTCCCGCATGACGAAGCTTGGCGTATTGAAGGCCCGTGGCCGCACCTCGTAGAGGAAACGTCCTTCACCGGCATGGGAGAGCTCCAGAGATAccc                              <  5:167474/146‑1 (MQ=255)                                                                                                                          gCAGATGCACTTCCGCGCGGAAATACTTGCGGGAATCGGCACCGTCTTCGGTGTCCCGCATGACGAAGCTTGGCGTATTGAAGGCCCGTTGCCGCACCTCGTAGAGGAAACGTCCTTCACCGGCATGGGAGAGCTCCAGAGATACCCGc                            >  7:249570/1‑149 (MQ=255)                                                                                                                          gCAGATGCACTTCCGCGCGGAAATACTTGCGGGAATCGGCACCGTCTTCGGTGTCCCGCATGACGAAGCTTGGCGTATTGAAGGCCCGTGGCCGCACCTCGTAGAGGAAACGTCCTTCACCGGCATGGGAGAGCTCCAGAGATACCCGc                            >  1:312143/1‑149 (MQ=255)                                                                                                                            aGATGCACTTCCGCGCGGAAATACTTGCGGGAATCGGCACCGTCTTCGGTGTCCCGCATGACGAAGCTTGGCGTATTGAAGGCCCGTGGCCGCACCTCGTAGAGGAAACGTCCTTCACCGGCATGGGAGAGCTCCAGAGATACCCGccc                          <  8:243712/149‑1 (MQ=255)                                                                                                                                gCACTTCCGCGCGGAAATACTTGCGGGAATCGGCACCGTCTTCGGTGTCCCGCATGACGAAGCTTGGCGTATTGAAGGCCCGTGGCCGCACCTCGTAGAGGAAACGTCCTTCACCGGCATgg                                                 >  7:183100/1‑122 (MQ=255)                                                                                                                                gCACTTCCGCGCGGAAATACTTGCGGGAATCGGCACCGTCTTCGGTGTCCCGCATGACGAAGCTTGGCGTATTGAAGGCCCGTGGCCGCACCTCGTAGAGGAAACGTCCTTCACCGGCATgg                                                 <  8:183100/122‑1 (MQ=255)                                                                                                                                             gAAATACTTGCGGGAATCGGCACCGTCTTCGGTGTCCCGCATGACGAAGCTTGGCGTATTGAAGGCCCGTGGCCGCACCTCGTAGAGGAAACGTCCTTCACCGGCATGGGAGAGCTCCAGAGATACCCGCCCATCATCACGCTCATTCa         >  6:212539/1‑149 (MQ=255)                                                                                                                                                  aCTTGCGGGAATCGGCACCGTCTTCGGTGTCCCGCATGACGAAGCTTGGCGTATTGAAGGCCCGTGGCCGCACCTCGTAGAGGAAACGTCCTTCACCGGCATGGGAGAGCTCCAGAGATACCCGcc                           >  6:68903/1‑126 (MQ=255)                                                                                                                                                  aCTTGCGGGAATCGGCACCGTCTTCGGTGTCCCGCATGACGAAGCTTGGCGTATTGAAGGCCCGTGGCCGCACCTCGTAGAGGAAACGTCCTTCACCGGCATGGGAGAGCTCCAGAGATACCCGcc                           <  5:68903/126‑1 (MQ=255)                                                                                                                                                    ttGCGGGAATCGGCACCGTCTTCGGTGTCCCGCATGACGAAGCTTGGCGTATTGAAGGCCCGTGGCCGCACCTCGTAGAGGAAACGTCCGTCACCGGCATGGGAGAGCTGCAGAGATACCCGCCCATCATCCCGCTCATTCACGCTCAc  >  2:279316/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                    CCGAATCCACTTCAGCCCACCACGTGCAAGTAATGCAGGTGGCGTTCGTACTGATCCAGGATGTCGCCAATCACATCGTCACGGCTCCACCCCATGATGTCGTAGTCCTGGCCCCCTTCTCGCAGATGCACTTCCGCGCGGAAATACTTGCGGGAATCGGCACCGTCTTCGGTGTCCCGCATGACGAAGCTTGGCGTATTGAAGGCCCGTGGCCGCACCTCGTAGAGGAAACGTCCTTCACCGGCATGGGAGAGCTCCAGAGATACCCGCCCATCATCACGCTCATTCACGCTCAC  >  NC_002947/4463970‑4464265

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