New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 = 457795039 (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

CCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCAC  >  NC_002947/4577802‑4578096
                                                                                                                                                    |                                                                                                                                                  
ccccTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAG                                                                                                                                                    >  5:201731‑M1/1‑149 (MQ=255)
catcttcGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAG                                                                                                                                                    <  2:247192‑M1/146‑1 (MQ=255)
  ccTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGt                                                                                                                                                   >  1:329377‑M1/1‑147 (MQ=255)
   cTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGa                                                                                                                                                 >  7:137045/1‑149 (MQ=255)
       gcgGGTAAACCGGCTCCTACACGGGTCACCGCGGGGTGCGGTGCCTGCGTAGGAACGGGTCTCTCGCGAAAGGTCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAAt                                                                                                                                             <  5:188292/149‑1 (MQ=255)
        ggggTAAAGCCGCTCCTTCACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGTGCCAGTACAGCCATCCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATg                                                                                                                                            <  3:15210/148‑1 (MQ=255)
               aCCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGccc                                                                                                                                     >  8:280025/1‑149 (MQ=255)
                       cTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGtggct                                                                                                                               <  6:201731/147‑1 (MQ=255)
                         acacGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGtggctggc                                                                                                                            >  6:230500/1‑148 (MQ=255)
                             gggTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACtt                                                                                                                       >  3:136816/1‑149 (MQ=255)
                              gctCAGCGGCGGGTGCGTTGCCTGCGTAGTAACGGGTCTCCCGCGAAAGTGCCAGTACAGCCAACCGATGAGGCGGGCTTCATGATCAAGACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCGACTTc                                                                                                                      <  8:216929/147‑1 (MQ=255)
                                tCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGt                                                                                                                    <  2:152390/149‑1 (MQ=255)
                                     gccgGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCActg                                                                                                               >  7:280025‑M2/1‑146 (MQ=255)
                                        gggTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGtgacgaagtggccagccacgggcagccattatcgctggcc                                                                                                            <  2:329377‑M1/149‑41 (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                                                                                           >  5:12143/1‑132 (MQ=255)
                                                                          gAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGcc                                                                                           <  6:12143/132‑1 (MQ=255)
                                                                                  ctgTACAGCCAACCGACAAGGCGGGCTTCATGCTCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTCGACCCGGCATACAGCTCg                                                                  <  7:250883/147‑1 (MQ=255)
                                                                                      aCAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTc                                                                                    <  2:62016/127‑1 (MQ=255)
                                                                                        aGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGtgacgaagtggccagccacgggcagccattatcgctggccc                                                                                                           >  5:35932‑M1/1‑61 (MQ=255)
                                                                                        aGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGtgacgaagtggccagccacgggcagccattatcgctggccc                                                                                                           <  6:35932‑M1/102‑42 (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                                     >  8:105436/1‑149 (MQ=255)
                                                                                                               aTGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACg                                     >  5:275526/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)
                                                                                                                                                  cAGtgacgaagtggccagccacgggcagccattatcgctggccctgaacctgccgccgttgggggtgttgatcatgaagcccgcctgatcggttggctgtactggccctttcgcgggagacccgttcctacgcaggcaccgcacccggc  >  7:280025‑M1/1‑3 (MQ=255)
                                                                                                                                                  cAGtgacgaagtggccagccacgggcagccattatcgctggccctgaacctgccgccgttgggggtgttgatcatgaagcccgcctgatcggttggctgtactggccctttcgcgg                                   <  4:202653‑M1/116‑114 (MQ=255)
                                                                                                                                                  cAGtgacgaagtggccagccacgggcagccattatcgctggccctgaacctgccgccgttgggggtgttgatcatgaagcccgcctgatcggttggctgtactggccctttcgcgg                                   >  3:202653‑M1/1‑3 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CCCCTTCGCGGGTAAACCCGCTCCTACACGGGTCACCGCCGGGTGCGGTGCCTGCGTAGGAACGGGTCTCCCGCGAAAGGGCCAGTACAGCCAACCGATCAGGCGGGCTTCATGATCAACACCCCCAACGGCGGCAGGTTCAGGGCCAGCGATAATGGCTGCCCGTGGCTGGCCACTTCGTCACTGGCCACCGCCCCGAGATTGCCGACATTCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCAC  >  NC_002947/4577802‑4578096

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