New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 4577950	39 (0.720)	3 (0.060)	3/252	NT	6.7%	coding (2161/2211 nt)	glgB	1,4‑alpha‑glucan branching enzyme
	?	NC_002947	= 4577984	46 (0.900)					coding (2127/2211 nt)	glgB	1,4‑alpha‑glucan branching enzyme

GCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACC  >  NC_002947/4577839‑4578131                                                                                                                                                  |                                                                                                                                                    gccgGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCActg                                                                                                                                                  >  7:280025‑M2/1‑146 (MQ=255)                    tAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGc                                                                                                                               <  1:147728/149‑1 (MQ=255)                       gAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGa                                                                                                                            <  5:230500/149‑1 (MQ=255)                           gggTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGACAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGTCCACCGCCCCGAGATTGCCGACAtt                                                                                                                        <  4:5008/149‑1 (MQ=255)                            ggTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGa                                                                                                                            >  8:40565/1‑144 (MQ=255)                            ggTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGa                                                                                                                            <  7:40565/144‑1 (MQ=255)                                tcCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCCCCGCCCCGACATTGCCGACATTCGAcc                                                                                                                   >  8:281157/1‑149 (MQ=255)                                tcCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTCGAcc                                                                                                                   >  7:263440/1‑149 (MQ=255)                                  ccGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCActg                                                                                                                                                  <  4:202653‑M2/116‑4 (MQ=255)                                  ccGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCActg                                                                                                                                                  >  3:202653‑M2/1‑113 (MQ=255)                                  ccGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTCGACCCg                                                                                                                 >  2:120923/1‑149 (MQ=255)                                  ccGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGACTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTCGACCCg                                                                                                                 >  2:161317/1‑149 (MQ=255)                                      gAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGcc                                                                                                                              <  6:12143/132‑1 (MQ=255)                                      gAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGcc                                                                                                                              >  5:12143/1‑132 (MQ=255)                                              ctgTACAGCCAACCGACAAGGCGGGCTTCATGCTCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTCGACCCGGCATACAGCTCg                                                                                                     <  7:250883/147‑1 (MQ=255)                                                  aCAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTc                                                                                                                       <  2:62016/127‑1 (MQ=255)                                                      ccAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTCGACCCGGCATACAGCTCGGCATCACt                                                                                             <  4:136816/149‑1 (MQ=255)                                                      ccAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTCGACCCGGCATACAGCTCGGCATCACt                                                                                             <  4:172523/149‑1 (MQ=255)                                                       cAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTCGACCCGGCATACAGCTCGGCATCACTg                                                                                            <  2:8981/149‑1 (MQ=255)                                                         aCCGATCAGGCGGGCTTCATGATCAACACCCCCACCGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTCGACCCGGCATACAGCTCGGCATCACTGtt                                                                                          <  8:137045/149‑1 (MQ=255)                                                            gATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTCGACCCGGCATACAGCTCGGCATCACTGtt                                                                                          >  5:157596/1‑146 (MQ=255)                                                              tCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTCGACCCGGCATACAGCTCGGCATCACTGTTCAGaa                                                                                     <  8:263440/149‑1 (MQ=255)                                                                  gCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTc                                                                                 <  5:347605/149‑1 (MQ=255)                                                                        ttCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTCGACCCGGCATAc                                                                                                           <  6:78833/117‑1 (MQ=255)                                                                        ttCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGcc                                                                            >  5:100703/1‑148 (MQ=255)                                                                           aTGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACg                                                                        >  5:275526/1‑149 (MQ=255)                                                                           aTGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACg                                                                        >  8:105436/1‑149 (MQ=255)                                                                             gATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCt                                                                      >  3:138847/1‑149 (MQ=255)                                                                              aTCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTc                                                                     <  4:138847/149‑1 (MQ=255)                                                                                      ccccAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGc                                                                             <  8:203189/133‑1 (MQ=255)                                                                                      ccccAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCAGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGc                                                                             >  7:203189/1‑133 (MQ=255)                                                                                          aaCGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCAc                                                         <  6:100703/149‑1 (MQ=255)                                                                                                       tCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTaa                                            >  8:45661/1‑149 (MQ=255)                                                                                                        cAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAAc                                           >  6:143434/1‑149 (MQ=255)                                                                                                          gggCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCActggccctgaacctgccgccgttgggggtgttgatcatgaagcccgcctgatcggttggct                                                                                        <  6:35932‑M2/102‑62 (MQ=255)                                                                                                          gggCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCActggccctgaacctgccgccgttgggggtgttgatcatgaagcccgcctgatcggttggct                                                                                        >  5:35932‑M2/1‑41 (MQ=255)                                                                                                           ggCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCActggccctgaacctgccgccgttgggggtgttgatcatgaagcccgcctgatcggttggctgtactggccctttcgcgggagacccgttcctacgcaggcaccgcaccc                                        <  2:329377‑M2/149‑110 (MQ=255)                                                                                                             ccAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCa                                      >  4:279837/1‑149 (MQ=255)                                                                                                                          cTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCgggg                         >  2:40171/1‑149 (MQ=255)                                                                                                                              ccGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTg                       >  7:184533/1‑147 (MQ=255)                                                                                                                              ccGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTg                       >  7:161595/1‑147 (MQ=255)                                                                                                                               cGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGaa                     <  1:120923/148‑1 (MQ=255)                                                                                                                                 tggctggcCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGtt                  <  1:161317/149‑1 (MQ=255)                                                                                                                                   gctggcCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTg                 <  6:157596/148‑1 (MQ=255)                                                                                                                                    ctggcCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGc               <  7:45661/149‑1 (MQ=255)                                                                                                                                    ctggcCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGc               <  1:40171/149‑1 (MQ=255)                                                                                                                                       gcCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGcc              <  7:105436/147‑1 (MQ=255)                                                                                                                                          aCTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGccac            <  7:168506/146‑1 (MQ=255)                                                                                                                                          aCTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGccac            >  8:168506/1‑146 (MQ=255)                                                                                                                                           cTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGccaccacc        <  5:143434/149‑1 (MQ=255)                                                                                                                                           cTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGccaccacc        <  6:275526/149‑1 (MQ=255)                                                                                                                                            ttCGTCACTGGCCACCGCCCCGAGATTGCCGACATTCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGccaccacca       >  2:60386/1‑149 (MQ=255)                                                                                                                                                  actggccctgaacctgccgccgttgggggtgttgatcatgaagcccgcctgatcggttggctgtactggccctttcgcgggagacccgttcctacgcaggcaccgcacccggcggtgacccgtgtaggagcgggtttacccgcgaagg  >  1:329377‑M2/1‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                    GCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACC  >  NC_002947/4577839‑4578131

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.