New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 = 559956638 (1.290)6 (0.210) 6/234 NT 14.1% coding (3967/5418 nt) PP_4924 subtilase family serine protease
?NC_002947 = 5599598 37 (1.320)coding (3935/5418 nt) PP_4924 subtilase family serine protease

CATGAACCCATATTTGAACGCCCCTGAACGCGAATCGCCCAGATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGCCGGCAGGGTCAAAGACAAAGTCGCAACCTGCCCCGCCGCAACCTGTTGAGCCACTGGTTCACTGCTGGCTGTCAGCAG  >  NC_002947/5599462‑5599731
                                                                                                                                        |                                                                                                                                     
cATGAACCCATATTTGAACGCCCCTGAACGCGAATCGCCCAGATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCAGGTCGCCAAGGCGGCCGACCAt                                                                                                                          >  2:409529/1‑150 (MQ=255)
      cccATATTTGAACGCCCCTGAACGCGAATCGCCCAGATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTg                                                                                                                    <  2:379515/150‑1 (MQ=255)
       ccATATTTGAACGCCCCTGAACGCGAATCGCCCAGATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTg                                                                                                                    <  2:544621/149‑1 (MQ=255)
       ccATATTTGAACGCCCCTGAACGCGAATCGCCCAGATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTg                                                                                                                    >  1:544621/1‑149 (MQ=255)
            tttGAACGCCCCTGAACGCGAATCGCCCAGATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCa                                                                                                                           >  1:257471/1‑137 (MQ=255)
            tttGAACGCCCCTGAACGCGAATCGCCCAGATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCa                                                                                                                           <  2:257471/137‑1 (MQ=255)
            tttGAACGCCCCTGAACGCGAATCGCCCAGATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTg                                                                                                                    >  2:546141/1‑144 (MQ=255)
            tttGAACGCCCCTGAACGCGAATCGCCCAGATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTg                                                                                                                    <  1:546141/144‑1 (MQ=255)
               gAACGCCCCTGTGCGCTAATCGCCCAGATCGTCATCACCAGAGCCCGCCGGCTTGGCTGTGCGATCGAGCAGTACGCTGGGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTGGACATCAAc                                                                                                           <  2:187796/150‑1 (MQ=255)
               gAACGCCCCTGAACGCGAATCGCCCAGATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTGGACATCAAc                                                                                                           >  1:101075/1‑150 (MQ=255)
                    ccccTGAACGCGAATCGCCCAGATCGTAGTCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTGGACATCAACCTCca                                                                                                      >  1:619792/1‑150 (MQ=255)
                    ccccTGAACGCAACTAGCCGATAGCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTGGACATCAACCTCca                                                                                                      >  2:793668/1‑150 (MQ=255)
                     cccTGAACGCGAATCGCCCAGATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGAcc                                                                                                                            <  2:385767/127‑1 (MQ=255)
                     cccTGAACGCGAATCGCCCAGATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGAcc                                                                                                                            >  1:385767/1‑127 (MQ=255)
                               gAATCGCCCAGATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTGGACATCAAc                                                                                                           <  1:240402/134‑1 (MQ=255)
                               gAATCGCCCAGATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTGGACATCAAc                                                                                                           >  2:240402/1‑134 (MQ=255)
                                     cccAGATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCaaa                                                                                     <  2:156582/150‑1 (MQ=255)
                                       cAGATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCACCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATc                                                                                   <  2:450517/150‑1 (MQ=255)
                                           tCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGcc                                                                               <  1:669731/150‑1 (MQ=255)
                                                             ggCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGcc                                                                               >  1:712336/1‑132 (MQ=255)
                                                             ggCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGcc                                                                               <  2:712336/132‑1 (MQ=255)
                                                                  cTTGCCCGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAg                                                                                                                                      >  1:776557‑M2/1‑71 (MQ=255)
                                                                  cTTGCCCGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAg                                                                                                                                      <  2:776557‑M2/72‑2 (MQ=255)
                                                                   ttGCCTGTGCGATAGAGCAGTACTCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGCCGGCAGGGTCAAAGACAAAGTCGCa                                                       >  2:669879/1‑150 (MQ=255)
                                                                          tgCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATc                                                                                                                         <  1:47209/77‑1 (MQ=255)
                                                                          tgCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATc                                                                                                                         >  2:47209/1‑77 (MQ=255)
                                                                           gCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGc                                                                                                 >  2:44094/1‑100 (MQ=255)
                                                                           gCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGc                                                                                                 <  1:44094/100‑1 (MQ=255)
                                                                               tCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTAATGTCGCCAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGCCGGCAGGGTCAAAGACAAAGTCGCAACCTGCCccgcc                                           <  2:619792/150‑1 (MQ=255)
                                                                                    cAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCACCAAATCAGCCGGCAGGGTCAAAGACAAAGTCGCAACCTGCCCCGCCGCAAc                                      <  1:409529/150‑1 (MQ=255)
                                                                                     aGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGCCGGCAGGGTCAAAGACAAAGTCGCAACCTGCCCCGCCGCAAcc                                     >  1:94658/1‑150 (MQ=255)
                                                                                       tACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGCCGGCAGGGTCAAAGACAAAGTCGCAACCTGCCCCGCCGCAACCTg                                   >  1:362011/1‑150 (MQ=255)
                                                                                               gcatGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGCCGGCAGGGTCAAAGACAAAGTCGCAACCTGCCCCGCCGCAACCTGTTGAGCCa                           <  1:465569/147‑1 (MQ=255)
                                                                                                cgTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGCCGGCAGGGTCAAAGACAAAGTCGCAACCTGCCCCGCCGCAACCTGTTGAGCCAc                          <  1:793668/150‑1 (MQ=255)
                                                                                                      gccAGGCGACACCAGCGTGAGCGTCATGTCGCCAAggcggcacgcgcagcgtactgctcgatcgcacaggcaaggcgccgggctctggtgatgacgatctgggcgattcgcgttcaggggcgttcaaatatgggttcatgtcgacccatc                    <  2:702763‑M2/150‑116 (MQ=255)
                                                                                                         aGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGCCGGCAGGGTCAAAGACAAAGTCGCAACCTGCCCCGCCGCAACCTGTTGAGCCACTGGTTCAct                 <  1:190294/150‑1 (MQ=255)
                                                                                                           gCGACACCAGCGTGAGCGTCATGTCGCCTgcggcacgcgcagcgtactgctcgatcgcacaggcaaggcgccgggctctggtgatgacgatctgggcgattcgc                                                             <  1:546308‑M2/105‑76 (MQ=255)
                                                                                                           gCGACACCAGCGTGAGCGTCATGTCGCCTgcggcacgcgcagcgtactgctcgatcgcacaggcaaggcgccgggctctggtgatgacgatctgggcgattcgc                                                             >  2:546308‑M2/1‑30 (MQ=255)
                                                                                                             gACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGCCGGCAGGGTCAAAGACAAAGTCGCAACCTGCCCcgccgc                                         >  2:410580/1‑122 (MQ=255)
                                                                                                             gACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGCCGGCAGGGTCAAAGACAAAGTCGCAACCTGCCCcgccgc                                         <  1:410580/122‑1 (MQ=255)
                                                                                                               cacCAGCGTGAGCGTCATGTCGCCAAggcggcacgcgcagcgtactgctcgatcgcacaggcaaggcgccgggctctggtgatgacgatctg                                                                     <  1:133520‑M2/92‑67 (MQ=255)
                                                                                                               cacCAGCGTGAGCGTCATGTCGCCAAggcggcacgcgcagcgtactgctcgatcgcacaggcaaggcgccgggctctggtgatgacgatctg                                                                     >  2:133520‑M2/1‑26 (MQ=255)
                                                                                                                  cAGCGTGAGCGTCATGTCGCCTgcggcacgcgcagcgtactgctcgatcgcacaggcaaggcgccgggctctggtgatgacgatctgggcgattcgcgttcaggggcgttcaaatatgggttcatgtcgacccatcatcgcgcagagcg         <  2:786791‑M2/150‑128 (MQ=255)
                                                                                                                       ggAGCGTCATGTCGCCAAGGCGCCCGACCAACGAGGAGAAATCAACAACAACAAGCTCCACCAGAAAAAAAGCCGACAAGCTCAAAGACAAACTCGCAACCTAACCCCCCGCAACATGTCGCGCCACCCGTTCACTTCTGGCTATaagca   >  2:109668/2‑150 (MQ=255)
                                                                                                                        gAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGCCGGCAGGGTCAAAGACAAAGCCGCAACCTGCCCCGCCGCAACCTGTTGAGCCACTGGTTCACTGCTGGCTGTcagcag  >  2:466065/1‑150 (MQ=255)
                                                                                                                          gCGTCATGTCGCCAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGCCGGCAGGGTCAAAGACAAAGTCGCAACCTGCCCCGCCGCAACCTGTTGAGCCACTGGTTCActgc               >  1:796356/1‑135 (MQ=255)
                                                                                                                          gCGTCATGTCGCCAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGCCGGCAGGGTCAAAGACAAAGTCGCAACCTGCCCCGCCGCAACCTGTTGAGCCACTGGTTCActgc               <  2:796356/135‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                     
CATGAACCCATATTTGAACGCCCCTGAACGCGAATCGCCCAGATCGTCATCACCAGAGCCCGGCGCCTTGCCTGTGCGATCGAGCAGTACGCTGCGCGTGCCGCCAGGCGACACCAGCGTGAGCGTCATGTCGCCAAGGCGGCCGACCATCGAGTGGACATCAACCTCCACATGCTCCACCAGCAAATCAGCCGGCAGGGTCAAAGACAAAGTCGCAACCTGCCCCGCCGCAACCTGTTGAGCCACTGGTTCACTGCTGGCTGTCAGCAG  >  NC_002947/5599462‑5599731

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