New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 5775531 =49 (0.930)4 (0.080) 4/236 NT 8.2% coding (1342/1473 nt) betB NAD‑dependent betaine aldehyde dehydrogenase
?NC_002947 5775547 = 42 (0.830)coding (1358/1473 nt) betB NAD‑dependent betaine aldehyde dehydrogenase

CGAAACCGAAGAAGAAGTGATCCGCCGTGCCAACGACACCGACTACGGCCTGGCCGCCGGCGTCTGCACCAACGACATCACCCGCGCCCACCGCATCATCCACAAGCTCGAAGCCGGTATCTGCTGGATCAACGCCTGGGGTGAATCGCCGGCCGAAATGCCGGTAGGTGGCTACAAGCAGTCGGGCGTCGGCCGTGAAAACGGCGTCAGCTCGCTGGCTCAATACACTCGCATCAAGTCGGTCCAGGTCGAGCTGGGCGGCTACAACTCGGTTTTCTGAACCCTGTCTAGCCCC  >  NC_002947/5775383‑5775677
                                                                                                                                                    |                                                                                                                                                  
cGAAACCGAAGAAGAAGTGATCCGCCGTGCCAACGACACCGACTACGGCCTGGCCGCCGGCGTCTGCACCAACGACATCACCCGCGCCCACCGCATCATCCACAAGCTCGAAGCCGGTATCTGCTGGATCAACGCCTGGGGTGAATCgcc                                                                                                                                                   <  2:896454/150‑1 (MQ=255)
    aCCGAAGAAGAAGTGATCCGCCGTGCCAACGACACCGACTACGGCCTGGCCGCCGGCGTCTGCACCAACGACATCACCCGCGCCCACCGCATCATCCACAAGCTCGAAGCCGGTATCTGCTGGATCAACGCCTGGGGTGAATCgccggcc                                                                                                                                               >  1:584666/1‑150 (MQ=255)
      cGAAGAAGAAGTGATCCGCCGTGCCAACGACACCGACTACGGCCTGGCCGCCGGCGTCTGCACCAACGACATCACCCGCGCCCACCGCATCATCCACAAGCTCGAAGCCGGTATCTGCTGGATCAACGCCTGGGGTGAATCGCCGGCCGa                                                                                                                                             >  1:1332891/1‑150 (MQ=255)
        aagaagaagTGATCCGCCGTGCCAACGACACCGACTACGGCCTGGCCGCCGGCGTCTGCACCAACGACATCACCCGCGCCCACCGCATCATCCACAAGCTCGAAGCCGGTATCTGCTGGATCAACGCCTGGGGTGAATCGCCGGCCGaaa                                                                                                                                           <  2:1324562/150‑1 (MQ=255)
        aagaagaagTGATCCGCCGTGCCAACGACACCGACTACGGCCTGGCCGCCGGCGTCTGCACCAACGACATCACCCGCGCCCACCGCATCATCCACAAGCTCGAAGCCGGTATCTGCTGGATCAACGCCTGGGGTGAATCGCCGGCCGaaa                                                                                                                                           <  2:1365361/150‑1 (MQ=255)
             gaagTGATCCGCCGTGCCAACGACACCGACTACGGCCTGGCCGCCGGCGTCTGCACCAACGACATCACCCGCGCCCACCGCATCATCCACAAGCTCGAAGCCGGTATCTGCTGGATCAACGCCTGGGGTGAATCGCCGGCCGAAATGCCg                                                                                                                                      >  1:1059004/1‑150 (MQ=255)
             gaagTGATCCGCCGTGCCAACGACACCGACTACGGCCTGGCCGCCGGCGTCTGCACCAACGACATCACCCGCGCCCACCGCATCATCCACAAGCTCGAAGCCGGTATCTGCTGGATCAACGCCTGGGGTGAATCGCCGGCCGAAATGCCg                                                                                                                                      <  1:880008/150‑1 (MQ=255)
                           tGCCAACGACACCGACTACGGCCTGGCCGCCGGCGTCTGCACCAACGACATCACCCGCGCCCACCGCATCATCCACAAGCTCGAAGCCGGTATCTGCTGGATCAACGCCTGGGGTGAATCGCCGGCCGAAATGCCGGTAGGTGGCTACaa                                                                                                                        <  2:1061885/150‑1 (MQ=255)
                           tGCCAACGACACCGACTACGGCCTGGCCGCCGGCGTCTGCACCAACGACATCACCCGCGCCCACCGCATCATCCACAAGCTCGAAGCCGGTATCTGCTGGATCAACGCCTGGGGTGAATCGCCGGCCGAAATGCCGGTAGGTGGCTACaa                                                                                                                        <  1:206226/150‑1 (MQ=255)
                           tGCCAACGACACCGACTACGGCCTGGCCGCCGGCGTCTGCACCAACGACATCACCCGCGCCCACCGCATCATCCACAAGCTCGAAGCCGGTATCTGCTGGATCAACGCCTGGGGTGAATCGCCGGCCGAAATGCCGGTAGGTGGCTACaa                                                                                                                        <  1:1402192/150‑1 (MQ=255)
                           tGCCAACGACACCGACTACGGCCTGGCCGCCGGCGTCTGCACCAACGACATCACCCGCGCCCACCGCATAATCCACATGCTCGAAGCCGGTATCTGCTGGATCAACGCCTGGGGTGAATCGCCGGCCGAAATGCCGGTAGGTGGCTACaa                                                                                                                        <  2:409207/150‑1 (MQ=255)
                              cAACGACACCGACTACGGCCTGGCCGCCGGCGTCTGCACCAACGACATCACCCGCGCCCACCGCATCATCCACAAGCTCGAAGCCGGTATCTGCTGGATCAACGCCTGGGGTGAATCGCCGGCCGAAATGCCGGTAg                                                                                                                                  <  1:350309/137‑1 (MQ=255)
                              cAACGACACCGACTACGGCCTGGCCGCCGGCGTCTGCACCAACGACATCACCCGCGCCCACCGCATCATCCACAAGCTCGAAGCCGGTATCTGCTGGATCAACGCCTGGGGTGAATCGCCGGCCGAAATGCCGGTAg                                                                                                                                  >  2:350309/1‑137 (MQ=255)
                               aaCGACACCGACTACGGCCTGGCCGCCGGCGTCTGCACCAACGACATCACCCGCGCCCACCGCATCATCCACAAGCTCGAAGCCGGTATCTGCTGGATCAACGCCTGGGGTGAATCGCCGGCCGAAATGCCGGTAGGTGGCTACAAGCAg                                                                                                                    >  2:721172/1‑150 (MQ=255)
                               aaCGACACCGACTACGGCCTGGCCGCCGGCGTCTGCACCAACGACATCACCCGCGCCCACCGCATCATCCACAAGCTCGAAGCCGGTATCTGCTGGATCAACGCCTGGGGTGAATCGCCGGCCGAAATGCCGGTAGGTGGCTACAAGCAg                                                                                                                    >  1:1354145/1‑150 (MQ=255)
                                aCGACACCGACTACGGCCTGGCCGCCGGCGTCTGCACCAACGACATCACCCGCGCCCACCGCATCATCCACAAGCTCGAAGCCGGTATCTGCTGGATCAACGCCTGGGGTGAATCGCCGGCCGAAATGCCGGTAGGTGGCTACAAGCAGt                                                                                                                   <  2:1354145/150‑1 (MQ=255)
                                         aCTACGGCCTGGCCGCCGGCGTCTGCACCAACGACATCACCCGCGCCCACCGCATCATCCACAAGCTCGAAGCCGGTATCTGCTGGATCAACGCCTGGGGTGAATCGCCGGCCGAAATGCCGGTAGGTGGCTATGAGCAGTCGGGCGTCg                                                                                                          <  2:1059545/150‑1 (MQ=255)
                                         aCTACGGCCTGGCCGCCGGCGTCTGCACCAACGACATCACCCGCGCCCACCGCATCATCCACAAGCTCGAAGCCAGTATCTGCTGGATCAACGCCTGGGGTGAATCGCCGGCCGAAATGCCGGTAGGTGGCTACAAGCAGTCGGGCGTCg                                                                                                          >  2:637228/1‑150 (MQ=255)
                                              ggCCTGGCCGCCGGCGTCTGCACCAACGACATCACCCGCGCCCACCGCATCATCCACAAGCTCGAAGCCGGTATCTGCTGGATCAACGCCTGGGGTGAATCGCCGGCCGAAATGCCGGTAGGTGGCTACAAGCAGTCGGGCGTCGGCCGt                                                                                                     <  1:721172/150‑1 (MQ=255)
                                                ccTGGCCGCCGGCGTCTGCACCAACGACATCACCCGCGCCCACCGCATCATCCACAAGCTCGAAGCCGGTATCTGCTTGATCAACGCCTGGGGTGAATCGCCGGCCGAAATGCCGGTAGGTGGCTACAAGCAGTCGGGCGTCGGCCGTGa                                                                                                   >  2:1276301/1‑150 (MQ=255)
                                                   ggCCGCCGGCGTCTGCACCAACGACATCACCCGCGCCCACCGCATCATCCACAAGCTCGAAGCCGGTATCTGCTGGATCAACGCCTGGGGTGAATCGCCGGCCGAAATGCCGGTAGGTGGCTACAAGCAGTCGGGCGTCGGCCGTGaaaa                                                                                                <  1:1124805/150‑1 (MQ=255)
                                                     ccgccgGCGTCTGCACCAACGACATCACCCGCGCCCACCGCATCATCCACAAGCTCGAAGCCGGTATCTGCTGGATCAACGCCTGGGGTGAATCGCCGGCCGAAATGCCGGTAGGTGGCTACAAGCAGTCGGGCGTCGGCCGTGAAAACg                                                                                              >  2:585526/1‑150 (MQ=255)
                                                     ccgccgGCGTCTGCACCAACGACATCACCCGCGCCCACCGCATCATCCACAAGCTCGAAGCCGGTATCTGCTGGATCAACGCCTGGGGTGAATCGCCGGCCGAAATGCCGGTAGGTGGCTACAAGCAGTCGGGCGTCGGCCGTGAAAACg                                                                                              >  2:901503/1‑150 (MQ=255)
                                                           gCGTCTGCACCAACGACATCACCCTCGCCTATCGCATCACCCACACGCTCGAAATCGGTACCCCGTGTATCACCCCCTGGGGTCAATCGCCGGCCGAAATGTCAGTAGCGGACCACAAGCAGTCGCCCGGCGGCCGTGAAAACAGCGtaa                                                                                        <  2:592537/150‑3 (MQ=255)
                                                                  cACCAACGACATCACCCGCGCCCACCGCATCATCCACAAGCTCGAAGCCGGTATCTGCTGGATCAACGCCTGGGGTGAATCGCCGGCCGAAATGCCGGTAGGTGGCTACAAGCAGTCGGGCGTCGGCCGTGAAAACGGCGTCAGCTCGCt                                                                                 <  2:816777/150‑1 (MQ=255)
                                                                    ccAACGACATCACCCGCGCCCACCGCATCATCCACAAGCTCGAAGCCGGTATCTGCTGGATCAACGCCTGGGGTGAATCGCCGGCCGAAATGCCGGTAGGTGGCTACAAGCAGTCGGGCGTCGGCCGTGAAAACGGCATCAGCTCGCTgg                                                                               >  2:129674/1‑150 (MQ=255)
                                                                           cATCACCCGCGCCCACCGCATCATCCACAAGCTCGAAGCCGGTATCTGCTGGATCAACGCCTGGGGTGAATCGCCGGCCGAAATGCCGGTAGGTCTGCACAATCAGTCGGGCGTCGGCCGTGAAAACGGCGTCAGCTCGCTGGCTCAATa                                                                        <  2:1059004/150‑1 (MQ=255)
                                                                              cACCCGCGCCCACCGCATCATCCACAAGCTCGAAGCCGGTATCTGCTGGATCAACGCCTGGGGTGAATCGCCGGCCGAAATGCCGGTAGGTGGCTACAAGCAGTCGGGCGTCGGCCGTGAAAACGGCGTCAGCTCGCTGGCTCAATacat                                                                     <  2:930137/150‑2 (MQ=255)
                                                                                cccGCGCCCACCGCATCATCCACAAGCTCGAAGCCGGTATCTGCTGGATCAACGCCTGGGGTGAATCGCCGGCCGAAATGCCGGTAGGTGGCTACAAGCAGTCGGGCGTCGGCCGTGAAAACGGCGTCAGCTCGCTGGCTCAATACACTc                                                                   <  1:585526/150‑1 (MQ=255)
                                                                                 ccGCGCCCACCGCGTCATCCACAAGCTCGAAGCCGGTATCTGCTGGATCAACGCCTGGGGTGAATCGCCGGCCGAAATGCCGGTAGGTGGCTACAAGCAGTCGGTCGTCGGCCGTGAAAAAGGCGTCCGCTCGCTGGCTCAATACACTcc                                                                  >  2:142308/1‑149 (MQ=255)
                                                                                 ccGCGCCCACCGCATCATCCACAAGCTCGAAGCCGGTATCTGCTGGATCAACGCCTGGGGTGAATCGCCGGCCGAAATGCCGGTAGGTGGCTACAAGCAGTCGGGCGTCGGCCGTGAAAACGGCGTCAGCTCGCTGGCTCAATACACTCg                                                                  >  2:835820/1‑150 (MQ=255)
                                                                                         aCCGCATACACCACAGGCTCGAAGCCGGTATCTGCTGGATCAACGCCTGGGGTGAATCGCCGGCCGAAATGCCGGTAGGTGGCTACAAGCAGTCGGGCGTCGGCCGTGAAAACGGCGTCAGCTCGCTGGCTCAATACACTCGCATCAAGt                                                          >  2:59540/1‑150 (MQ=255)
                                                                                           cGCATCATCCACAAGCTCGAAGCCGGTATCTGCTGGATTAACGCCTGGGGTGAATCGCCGGCCGAAATGCCGGTAGGTGGCTACAAGCAGTCGGGCGTCGGCCGTGAAAACGGCg                                                                                           >  1:629982/1‑115 (MQ=255)
                                                                                           cGCATCATCCACAAGCTCGAAGCCGGTATCTGCTGGATCAACGCCTGGGGTGAATCGCCGGCCGAAATGCCGGTAGGTGGCTACAAGCAGTCGGGCGTCGGCCGTGAAAACGGCg                                                                                           <  2:629982/115‑1 (MQ=255)
                                                                                            gCATCATCCACAAGCTCGAAGCCAGTATCTGCTGGATCAACGCCTGGGGTGAATCGCCGGCCGAAATGCCGGTAGGTGGCTACAAGCAGTCGGGCGTCGGCCGTGAAAACGGCGTCAGCTCGCTGGCTCAATACACTCGCATCAAGTCgg                                                       <  1:637228/150‑1 (MQ=255)
                                                                                                  tcCACAAGCTCGAAGCCGGTATCTGCTGGATCAACGCCTGTGGTGAATCGCCGGCCGAAATGCCGGTAGGTGGCTACAAGCAGTCGGGCGTCGGCCGTGAAAACGGCGTCAGCTCGCTGGCTCAATACACTCGCATCAAGTCGGTCCAgg                                                 >  1:1355073/1‑150 (MQ=255)
                                                                                                  tcCACAAGCTCGAAGCCGGTATCTGCTGGATCAACGCCTGGGGTGAATCGCCGGCCGAAATGCCGGTAGGTGGCTACAAGCAGTCGGGCGTCGGCCGTGAAAACgg                                                                                             <  2:1094562/106‑1 (MQ=255)
                                                                                                  tcCACAAGCTCGAAGCCGGTATCTGCTGGATCAACGCCTGGGGTGAATCGCCGGCCGAAATGCCGGTAGGTGGCTACAAGCAGTCGGGCGTCGGCCGTGAAAACgg                                                                                             >  1:1094562/1‑106 (MQ=255)
                                                                                                    cacaAGCTCGAAGCCGGTATCTGCTGGATCAACGCCTGGGGTGAATCGCCGGCCGAAATGCCGGTAGGTGGCTACAAGCAGTCGGGCGTCGGCCGTGAAAACGGCGTCAGCTCGCTGGCTCAATACACTCGCATCAAGTCGGTCCAGGTc                                               >  2:293830/1‑150 (MQ=255)
                                                                                                     acaAGCTCGAAGCCGGTATCTGCTGGATCAACGCCTGGGGTGAATCGCCGGCCGAAATGCCGGTAGGTGGCTACGAGCAGTCGGGCGTCGGCCGTGAAAAc                                                                                               >  2:1262751/1‑101 (MQ=255)
                                                                                                     acaAGCTCGAAGCCGGTATCTGCTGGATCAACGCCTGGGGTGAATCGCCGGCCGAAATGCCGGTAGGTGGCTACGAGCAGTCGGGCGTCGGCCGTGAAAAc                                                                                               <  1:1262751/101‑1 (MQ=255)
                                                                                                     acaAGCTCGAAGCCGGTATCTGCTGGATCAACGCCTGGGGTGAATCGCCGGCCGAAATGCCGGTAGGTGGCTACAAGCCGCAGGGCGTCGGCCATGAAAACGGCGTCAGCTCGCTGGCTCAATACACTCGCATCAAGTCGGTCCAGGTCg                                              >  1:247044/1‑150 (MQ=255)
                                                                                                           tCGAAGCCGGTATCTGCTGGATCAACGCCTGGGGTGAATCGCCGGCCGAAATGCCGGTAGGTGGCTACAAGCAGTCGGGCGTCGGCCGTGAAAACGGCGTCAGCTCGCTGGCTCAATACACTCGCATCAAGTCGGTCCAGGTCGAGCTgg                                        <  2:584666/150‑1 (MQ=255)
                                                                                                             gAAGCCGGTATCTGCTGGATCAACGCCTGGGGTGAATCGCCGGCCGAAATGCCGGTAGGTGGCTACAAGCAGTCGGGCGTCGGCCGTGAAAACGGCGTCAGCTCGCTGGCTCAATACACTCGCATCAAGTCGGTCCAGGTCGAGTTGggc                                      >  1:832231/1‑150 (MQ=255)
                                                                                                                       gccgacgcccgactgcttgtagccacctaCCGGCCGAAATGCCGGTAg                                                                                                                                  <  1:415325‑M1/19‑1 (MQ=255)
                                                                                                                       gccgacgcccgactgcttgtagccacctaCCGGCCGAAATGCCGGTAg                                                                                                                                  >  2:415325‑M1/30‑48 (MQ=255)
                                                                                                                         tgctgGATCAACGCCTGGGGTGAATCGCCGGCCGAAATGCCGGTAGGTGGCTACAAGCAGTCGGGCGTCGGCCGTGAAAACGGCGTCAGCTCGCTGGCTCAATACACTCGCATCAAGTCGGTCCAGGTCGAGCTGGGCGGCTACAACTCg                          <  1:142308/150‑1 (MQ=255)
                                                                                                                         tgctgGATCAACGCCTGGGGTGAATCGCCGGCCGAAATGCCGGTAGGTGGCTACAAGCAGTCGGGCGTCGGCCGTGAAAACGGCGTCAGCTCGCTGGCTCAATACACTCGCATCAAGTCGGTCCAGGTCGAGCTGGGCGGCTACAACTCg                          >  1:953635/1‑150 (MQ=255)
                                                                                                                                tCAACGCCTGGGGTGAATCGCCGGCCGAAATGCCGGTAGGTGGCTACAAGCAGTCGGGCGTCGGCCGTGAAAACGGCGTCAGCTCGCTGGCTCAATACACTCGCATCAAGTCGGTCCAGGTCGAGCTGGGCGGCTACAACTCGGTTTTCt                   <  1:835820/150‑1 (MQ=255)
                                                                                                                                 gactgcttgtagccacctaCCGGCATTTCGGCCGGTAGGTGGCTACAAGCAGTCGGGCGTCGGCCGTGAAAACGGCGTCAGCTCGCTGGCTCAATACACTCGCATCAAGTCGGTCCAGGTCGAGCTGGGCGGCTACAACTCGGTTTTc                    <  2:1341024‑M1/129‑1 (MQ=255)
                                                                                                                                 gactgcttgtagccacctaCCGGCATTTCGGCCGGTAGGTGGCTACAAGCAGTCGGGCGTCGGCCGTGAAAACGGCGTCAGCTCGCTGGCTCAATACACTCGCATCAAGTCGGTCCAGGTCGAGCTGGGCGGCTACAACTCGGTTTTc                    >  1:1341024‑M1/20‑148 (MQ=255)
                                                                                                                                       cTGGGGTGAATCGCCGGCCGAAATGCCGGTAGGTGGCTACAAGCAGTCGGGCGTCGGCCGTGAAAACGGCGTCAGCTCGCTGGCTCAATACACTCGCATCAAGTCGGTCCAGGTCGa                                             >  1:233784/1‑117 (MQ=255)
                                                                                                                                       cTGGGGTGAATCGCCGGCCGAAATGCCGGTAGGTGGCTACAAGCAGTCGGGCGTCGGCCGTGAAAACGGCGTCAGCTCGCTGGCTCAATACACTCGCATCAAGTCGGTCCAGGTCGa                                             <  2:233784/117‑1 (MQ=255)
                                                                                                                                            gTGAATCGCCGGCCGAAATGCCGGTAGGTGGCTACAAGCAGTCGGGCGTCGGCCATGAAAACGGCGTCAGCTCGCTGGCTCAATACACTCGCATCAAGTCGGTCCAGGTCGAGCTGGGCGGCTACAACTCGGTTTTCTGAACCCTGTCTa       <  2:247044/150‑1 (MQ=255)
                                                                                                                                                 tCGCCGGCCGAAATGCCGGTAGGTGGCTACAAGCAGTCGGGCGTCGGCCGTGAAAACGGCGTCAGCTCGCTGGCTCAATACACTCGCATCAAGTCGGTCCAGGTCGAGCTGGGCGGCTACAACTCGGTTTTCTGAACCCTGTCTAGcccc  <  1:129674/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CGAAACCGAAGAAGAAGTGATCCGCCGTGCCAACGACACCGACTACGGCCTGGCCGCCGGCGTCTGCACCAACGACATCACCCGCGCCCACCGCATCATCCACAAGCTCGAAGCCGGTATCTGCTGGATCAACGCCTGGGGTGAATCGCCGGCCGAAATGCCGGTAGGTGGCTACAAGCAGTCGGGCGTCGGCCGTGAAAACGGCGTCAGCTCGCTGGCTCAATACACTCGCATCAAGTCGGTCCAGGTCGAGCTGGGCGGCTACAACTCGGTTTTCTGAACCCTGTCTAGCCCC  >  NC_002947/5775383‑5775677

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