New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 = 559951544 (0.890)4 (0.090) 4/244 NT 8.6% coding (3967/5418 nt) PP_4924 subtilase family serine protease
?NC_002947 = 5599547 43 (0.910)coding (3935/5418 nt) PP_4924 subtilase family serine protease

ATGATGGGTCGACATGAACCCATATTTGAACGCCCCTGAACGCGAATCGCCCAGATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGC‑‑‑CAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGCCGGCAGGGTCAAAGACAAAGTCGCAACCTGCCCCGCCGCAACCTGTTGAGCCACTGGTTCACTGCTGGCTGTCAGCAGCCGTTCATTG  >  NC_002947/5599399‑5599690
                                                                                                                                                       |                                                                                                                                               
atgatgGGTCGACATGAACCCATATTTGAACGCCCCTGAACGCGAATCGCCCAGATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGC‑‑‑Caa                                                                                                                                                 <  7:13350/149‑1 (MQ=255)
 tgatgGGTCGACATGAACCCATATTTGAACGCCCCTGAACGCGAATCGCCCAGATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGC‑‑‑CAAg                                                                                                                                                <  1:239824‑M2/149‑2 (MQ=255)
                  cccATATTTGAACGCCCCTGAACGCGAATCGCCCAGATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGGGAGCGTCATGTCGC‑‑‑CAAGGCGGCCGACCATCGAGt                                                                                                                               <  5:688/149‑1 (MQ=255)
                     atatTTGAACGCCCCTGAACGCGAATCGCCCAGATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGC‑‑‑CAAGGCGGCCGACCATCGAGTGGa                                                                                                                            >  1:296901/1‑149 (MQ=255)
                      tatTTGAACGCCCCTGAACGCGAATCGCCCAGATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGC‑‑‑CAAGGCGGCCGACCATCGAGTGGAc                                                                                                                           <  5:70871/149‑1 (MQ=255)
                        tttGAACGCCCCTGAACGCGAATCGCCCAGATCGTCAACACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGC‑‑‑CAAGGCGGCCGACCATCGAg                                                                                                                                >  1:66489/1‑142 (MQ=255)
                        tttGAACGCCCCTGAACGCGAATCGCCCAGATCGTCAACACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGC‑‑‑CAAGGCGGCCGACCATCGAg                                                                                                                                <  2:66489/142‑1 (MQ=255)
                                 cccTGAACGCGAATCGCCCAGATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGC‑‑‑CAAGGCGGCCGACCATCGAGTGGACATCAACCTCca                                                                                                                <  2:2797/149‑1 (MQ=255)
                                 cccTGAACGCGAATCGCCCAGATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGC‑‑‑CAAGGCGGCCGACCATCGAGTGGACATCAACCTCca                                                                                                                >  6:269137/1‑149 (MQ=255)
                                           gAATCGCCCAGATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGC‑‑‑CAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTccac                                                                                                      >  5:195358/1‑149 (MQ=255)
                                           gAATCGCCCAGATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGC‑‑‑CAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTccac                                                                                                      <  5:13223/149‑1 (MQ=255)
                                           gAATCGCCCAGATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGC‑‑‑CAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTccac                                                                                                      >  2:94682/1‑149 (MQ=255)
                                              tCGCCCAGATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGC‑‑‑CAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAg                                                                                                   <  6:195358/149‑1 (MQ=255)
                                                gCCCAGATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGC‑‑‑CAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGc                                                                                                  >  6:83100/1‑148 (MQ=255)
                                                 cccAGATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGC‑‑‑CAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCa                                                                                                 >  7:87678/1‑148 (MQ=255)
                                                 cccAGATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTACCGCCAGGCGACACCAGCGTGAGCGTCATGTCGC‑‑‑CAAGGCGGCCGACCATc                                                                                                                                   >  7:152070/1‑114 (MQ=255)
                                                 cccAGATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTACCGCCAGGCGACACCAGCGTGAGCGTCATGTCGC‑‑‑CAAGGCGGCCGACCATc                                                                                                                                   <  8:152070/114‑1 (MQ=255)
                                                               aCCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGC‑‑‑CAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGCCGGCAgg                                                                                   <  2:296901/148‑1 (MQ=255)
                                                                ccAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGC‑‑‑CAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTc                                                                                                         <  7:114528/125‑1 (MQ=255)
                                                                ccAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGC‑‑‑CAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTc                                                                                                         >  8:114528/1‑125 (MQ=255)
                                                                   gagCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACTCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCCTCAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGCCGGCAGGGt                                                                                 <  1:94682/149‑1 (MQ=255)
                                                                   gagCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGC‑‑‑CAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATc                                                                                             <  3:151726/134‑1 (MQ=255)
                                                                   gagCCCGGCGCCTTGCCTGTGCGATCGAGCAGTAAGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGC‑‑‑CAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATc                                                                                             >  4:151726/1‑134 (MQ=255)
                                                                    agCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGC‑‑‑CAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGCCGGCAGGGTCaaa                                                                             <  2:239824/149‑1 (MQ=255)
                                                                        cGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGC‑‑‑CAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGCCGGCAGGGTCaaa                                                                             <  6:42402/145‑1 (MQ=255)
                                                                        cGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGC‑‑‑CAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGCCGGCAGGGTCaaa                                                                             >  5:42402/1‑145 (MQ=255)
                                                                         ggCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGC‑‑‑CAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGCCGGCAGGGTCAAAGACaa                                                                        <  5:269137/149‑1 (MQ=255)
                                                                         ggCGCCTTGCCTGTGCGATCAAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGC‑‑‑CAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGCCGGCAGGGTCAAAGACa                                                                         >  3:171697/1‑148 (MQ=255)
                                                                          gcgcCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGC‑‑‑CAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGCCGGCAGGGTCAAAGACaaa                                                                       <  5:83100/149‑1 (MQ=255)
                                                                                      tgCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGC‑‑‑CAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGCCGGCAGGGTCAAAGACAAAGTCGCAACCTGc                                                           >  4:120153/1‑149 (MQ=255)
                                                                                             gAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGC‑‑‑CAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGCCGGCAGGGTCAAAGACAAAGTCGCAACCTGCCccgcc                                                     <  3:120153/148‑1 (MQ=255)
                                                                                                  gTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGC‑‑‑CAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATCCTCCACCAGCAAATCAGCCGGCAGGGTCAAAGACAAAGTCGCAACCTGCCCCGCCGCAAcc                                               >  2:249549/1‑149 (MQ=255)
                                                                                                      gCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGC‑‑‑CAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGCCGGCAGGGTCAAAGACAAAGTCGCAACCTGCCCCGCCGCAACCTGtt                                           <  4:171697/149‑1 (MQ=255)
                                                                                                      gCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGC‑‑‑CAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGCCGGCAGGGTCAAAGACAAAGTCGCAACCTGCCCCGCCGCAACCTGtt                                           >  6:326534/1‑149 (MQ=255)
                                                                                                      gCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGC‑‑‑CAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGCCGGCAGGGTCAAAGACAAAGTCGCAACCTGCCCCGCCGCAACCTGtt                                           <  1:308379/149‑1 (MQ=255)
                                                                                                       cTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGC‑‑‑CAAggcggcacgcgcagcgtactgctcgatcgcacaggcaaggcgccgggctctggtgatgacgatctgggcgattcgcgttcaggggcgttcaaatatg                                                >  1:184450‑M2/1‑46 (MQ=255)
                                                                                                       cTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGC‑‑‑CAAggcggcacgcgcagcgtactgctcgatcgcacaggcaaggcgccgggctctggtgatgacgatctgggcgattcgcgttcaggggcgttcaaatatg                                                <  2:184450‑M2/143‑98 (MQ=255)
                                                                                                        tGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGC‑‑‑CAAggcggcacgcgcagcgtactgctcgatcgcacaggcaaggcgccgggctctggtgatgacgatctgggcgattcgcgttcaggggcgttcaaatatgggttcat                                         <  8:56006‑M2/149‑105 (MQ=255)
                                                                                                             gTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGC‑‑‑CAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGCCGGCAGGGTCAAAGACAAAGTCGCAACCTGCCCCGCCGCAACCTGTTGAGCCAc                                    <  8:87678/149‑1 (MQ=255)
                                                                                                                 cgccAGGCGACACCAGCGTGAGCGTCATGTCGC‑‑‑CAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGCCGGCAGGGTCAAAGACAAAGTCGCAACCTGCCCCGCCGCAACCTGTTGAGCCACTGGt                                >  6:57894/1‑149 (MQ=255)
                                                                                                                  gccAGGCGACACCAGCGTGAGCGTCATGTCGC‑‑‑CAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGCCGGCAGGGTCAAAGACAAAGTCGCAACCTGCCCCGCCGCAACCTGTTGAGCCACTGGtt                               >  8:89112/1‑149 (MQ=255)
                                                                                                                              cAGCGTGAGCGTCATGTCGC‑‑‑CAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGCCGGCAGGGTCAAAGACAAAGTCGCAACCTGCCCCGCCGCAACCTGTTGAGCCACTGGTTCACTgctggctg                   >  3:137222/1‑149 (MQ=255)
                                                                                                                               aGCGTGAGCGTCATGTCGC‑‑‑CAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGCCGGCAGGGTCAAAGACAAAGTCGCAACCTGCCCCTCCGCAACCTGTTGAGCCACTGGTTCACTGCTGGCTGt                  >  8:221536/1‑149 (MQ=255)
                                                                                                                                      gCGTCATGTCGC‑‑‑CAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGCCGGCAGGGTCAAAGACAAAGTCGCAACCTGCCCCGCCGCAACCTGTTGAGCCACTGGTTCACTGCTGACAGCcagcagt           <  7:89112/149‑2 (MQ=255)
                                                                                                                                        gTCATGTCGC‑‑‑CAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGCCGGCAGGGTCAAAGACAAAGTCGCAACCTGCCCCGCCGCAACCTGTTGAGCCACTGGTTCACTGCTGGCTGTCAGCAGCCg         >  8:184621/1‑149 (MQ=255)
                                                                                                                                         tCATGTCGC‑‑‑CAAggcggcacgcgcagcgtactgctcgatcgcacaggcaaggcgccgggctctggtgatgacgatctgggcgattcgcgttcaggggcgttcaaatatgggttcatgtcgacccatcatcgcgcagagcgctcggctgg        >  7:56006‑M2/1‑12 (MQ=255)
                                                                                                                                            tGTCGC‑‑‑CAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGCCGGCAGGGTCAAAGACAAAGTCGCAACCTGCCCCGCCGCAACCTGTTGAGCCACTGGTTCACTGCTGGCTGTCAGCAGCCGTTCa     >  2:81751/1‑149 (MQ=255)
                                                                                                                                             gTCGC‑‑‑CAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGCCGGCAGGGTCAAAGACAAAGTCGCAACCTGCCCCGCCGCAACCTGTTGAGCCACTGGTTCACTGCTGGCTGTCAGCAGCCGTTCAt    >  7:165655/1‑149 (MQ=255)
                                                                                                                                             gTCGC‑‑‑CAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGCCGGCAGGGTCAAAGACAAAGTCGCAACCTGCCCCGCCGCAACCTGTTGAGCCACTGGTTCACTGCTGGCTGTCAGCAGCCGTTCAt    >  8:173464/1‑149 (MQ=255)
                                                                                                                                                gC‑‑‑CAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGCCGGCAGGGTCAAAGACAAAGTCGCAACCTGCCCCGCCGCAACCTGTTGAGCCACTGGTTCACTGCTGGCTGTca                <  3:122505/134‑1 (MQ=255)
                                                                                                                                                gC‑‑‑CAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGCCGGCAGGGTCAAAGACAAAGTCGCAACCTGCCCCGCCGCAACCTGTTGAGCCACTGGTTCACTGCTGGCTGTca                >  4:122505/1‑134 (MQ=255)
                                                                                                                                                gC‑‑‑CAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGCCGGCAGGGTCAAAGACAAAGTCGCAACCTGCCCCGCCGCAACCTGTTGAGCCACTGGTTCACTGCTGGCTGTCAGAAGCCGTTCATTg  <  7:221536/148‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                               
ATGATGGGTCGACATGAACCCATATTTGAACGCCCCTGAACGCGAATCGCCCAGATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGC‑‑‑CAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGCCGGCAGGGTCAAAGACAAAGTCGCAACCTGCCCCGCCGCAACCTGTTGAGCCACTGGTTCACTGCTGGCTGTCAGCAGCCGTTCATTG  >  NC_002947/5599399‑5599690

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.