New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 4320675 =32 (0.610)4 (0.080) 4/236 NT 12.0% coding (134/627 nt) PP_3791 site‑specific recombinase
?NC_002947 4320700 = 28 (0.560)coding (109/627 nt) PP_3791 site‑specific recombinase

CAAACTGTACGGGACGTTGGGTTTTCTGCTGCATCACTATCGCCCGCGATGACACGTGCTCCCCATGAGCAACGTCGCGTACGCGAAGCTTGGTCAGATCGCAGGCTCGTAGCTTGAGCTTGCTATCGATGGCCAGATCGAAAAGGGCTAAATCACGAGTTTCCTCCGCGATTTGGAGCCTTACCCGAATGGCCCAGATATCTTTCAGTCGGAGCGGAGCTTTCTGCCCGACCAATTTTCCTTTGTTCAAAGGCTGATGGCTGCTGGTCGTAATGGCGTTCATGGCTTGTCCTCCATGT  >  NC_002947/4320526‑4320824
                                                                                                                                                     |                                                                                                                                                     
ggttgatcatccttgtcctacctcacatggaggacaagccatgaacgccattacgaccagcagccatcagcctttgaacaaaggaaaattggtcgggcagaaagctccgctccgactgaaagatatctgggccattcgggtaaggctcca                                                                                                                                                       <  1:171180‑M1/1‑1 (MQ=255)
 gttgatcatccttgtcctacctcacatggaggacaagccatgaacgccattacgaccagcagccatcagcctttgaacaaaggaaaattggtcgggcagaaagctccgctccgactgaaagatatctgggccattcgggtaaggctccaa                                                                                                                                                      <  1:43915‑M1/2‑1 (MQ=255)
    cTGTACGGGACGTTGGGTTTTCTGCTGCATCACTATCGCCCGCGATGACACGTGCTCCCCATGAGCAACGTCGCGTACGCGAAGCTTGGTCAGATCGCAGGCTCGTAGCTTGAGCTTGCTATCGATGGCCAGATCGAAAAGGGCTAAATc                                                                                                                                                   <  2:657094/150‑1 (MQ=255)
     tGTACGGGACGTTGGGTTTTCTGCTGCATCACTATCGCCCGCGATGACACGTGCTCCCCATGAGCAACGTCGCGTACGCGAAGCTTGGTCAGATCGCAGGCTCGTAGCTTGAGCTTGCTATCGATGGCCAGATCGGAAAGGGCTAAATCa                                                                                                                                                  >  2:70008/1‑150 (MQ=255)
        aCGGGACGTTGGGTTTTCTGCTGCATCACTATCGCCCGCGATGACACGTGCTCCCCATGAGCAACGTCGCGTACGCGAAGCTTGGTCAGATCGCAGGCTCGTAGCTTGAGCTTGCTATCGATGGCCAGATCGAAAAGGGCTAAATCACGa                                                                                                                                               <  2:1138461/150‑1 (MQ=255)
                ttGGGTTTTCTGCTGCATCACTATCGCCCGCGATGACACGTGCTCCCCATGAGCAACGTCGCGTACGCGAAGCTTGGTCAGATCGCAGGCTCGTAGCTTGAGCTTGCTATCGATGGCCAGATCGAAAAGGGCTAAATCACGAGTTtcctc                                                                                                                                       >  2:171180/1‑150 (MQ=255)
                ttGGGTTTTCTGCTGCATCACTATCGCCCGCGATGACACGTGCTCCCCATGAGCAACGTCGCGTACGCGAAGCTTGGTCAGATCGCAGGCTCGTAGCTTGAGCTTGCTATCGATGGCCAGATCGAAAAGGGCTAAATCACGAGTTtcctc                                                                                                                                       >  1:1324072/1‑150 (MQ=255)
                 ctacctcacatggaggacaagccatgaacgccattacgaccagcagccatcagcctttgaacaaaggaaaattggtcgggcagaaagctccgctccgactgaaagatatctgggccattcgggtaaggctccAAATCACGAGTTtcctcc                                                                                                                                      <  1:209875‑M1/18‑1 (MQ=255)
                    gTTTTCTGCTGCATCACTATCGCCCGCGATGACACGTGCTCCCCATGAGCAACGTCGCGTACGCGAAGCTTCGTCAGAGCGCAGGCTCGTAGCACGACTTTGTCAACCCAGTCCAATGCTGAAAGGGCTAAATCACGAGTTTCCTCcgcg                                                                                                                                   >  1:1027869/1‑150 (MQ=255)
                         ctgctgCATCACTATCGCCCGCGATGACACGTGCTCCCCATGAGCAACGTCGCGTACGCGAAGCTTGGTCAGATCGCAGGCTCGTAGCTTGAGCTTGCTATCGATGGCCAGATCGAAAAGGGCTAAATCACGAGTTTCCTCCGCGATTTg                                                                                                                              <  2:1324072/150‑1 (MQ=255)
                           gctgcATCACTATCGCCCGCGATGACACGTGCTCCCCATGAGCAACGTCGCGTACGCGAAGCTTGGTCAGATCGCAGGCTCGTAGCTTGAGCTTGCTATCGATGGCCAGATCGAAAAGGGCTAAATCACGAGTTTCCTCCGCGATTTAGa                                                                                                                            <  2:1348702/150‑1 (MQ=255)
                                aTCACTATCGCCCGCGATGACACGTGCTCCCCATGAGCAACGTCGCGTACGCGAAGCTTGGTCAGATCGCAGGCTCGTAGCTTGAGCTTGCTATCGATGGCCAGATCGAAAAGGGCTAAATCACGAGTTTCCTCCGCGATTTGGAGCCtt                                                                                                                       >  1:1340349/1‑150 (MQ=255)
                                    agccatgaacgccattacgaccagcagccatcagcctttgaacaaaggaaaattggtcgggcagaaagctccgctccgactgaaagatatctggtccattcgggtaaggctccAAATCACGAGTTTCCTCCGCGATTTGGAGCCTTAccc                                                                                                                   >  2:209875‑M1/114‑150 (MQ=255)
                                     tATCGCCCGCGATGACACGTGCTCCCCATGAGCAACGTCGCGTACGCGAAGCTTGGTCAGATCGCAGGCTCGTAGCTTGAGCTTGCTATCGATGGCCAGATCGAAAAGGGCTAAATCACGAGTTTCCTCCGCGATTTGGAGCCTTACCCg                                                                                                                  >  1:820188/1‑150 (MQ=255)
                                            cgcgATGACACGTGCTCCCCATGAGCAACGTCGCGTACGCGCAGCTTGGTCACATCGCATGCTCGTAGCTTGAGCTTGTTCTCCACAGCCCGGTGACGAACTGCCAAGGCTCGAGTTTCCTCCGCGATTTGGAGCCTTACCCGAATGGcc                                                                                                           <  1:1096725/150‑1 (MQ=255)
                                                aTGACACGTGCTCCCCATGAGCAACGTCGCGTACGCGAAGCTTGGTCCGCCCGCAGGTTCGTAGCTTGAGCCTGCTATCGATGGCCAGATCGAAAAGGGCTAAATCACGAGTTTCCTCCGCGATTTGGAGCCTTACCCGAATGGCCCAGa                                                                                                       <  2:338095/150‑1 (MQ=255)
                                                 tGACACGTGCTCCCCATGAGCAACGTCGCGTACGCGAAGCTTGGTCAGATCGCAGGCTCGTAGCTTGAGCTTGCTATCGATGGCCAGATCGAAAAGGGCTAAATCACGAGTTTCCTCCGCGATTTGGAGCCTTACCCGAATGGCCCAGat                                                                                                      >  1:973935/1‑150 (MQ=255)
                                                    cacGTGCTCCCCATGAGCAACGTCGCGTACGCGAAGCTTGGTCAGATCGCAGGCTCGTAGCTTGAGCTTGCTATCGATGGCCAGATCGAAAAGGGCTAAATCACGAGTTTCCTCCGCGATTTGGAGCCTTACCCGAATGGCCCAGATATc                                                                                                   >  1:451096/1‑150 (MQ=255)
                                                     acGTGCTCCCCATGAGCAACGTCGCGTACGCGAAGCTTGGTCAGATCGCAGGTTCGTAGCTTGAGCTTGCTATCGATGGCCAGATCGAAAAGGGCTAAATCACGAGTTTCCTCCGCGATTTGGAGCCTTACCCGAATGGCCCAGATATCt                                                                                                  <  1:1149790/150‑1 (MQ=255)
                                                         gCTCCCCATGAGCAACGTCGCGTACGCGAAGCTTGGTCAGATCGCAGGCTCGTAGCTTGAGCTTGCTATCGATGGCCAGATCGAAAAGGGCTAAATCACGAGTTTCCTCCGCGATTTGGAGCCTTACCCGAATGGCCCAGATATCTTTCa                                                                                              <  2:294457/150‑1 (MQ=255)
                                                           tCCCCATGAGCAACGTCGCGTACGCGAAGCTTGGTCAGATCGCAGGCTCGTAGCTTGAGCTTGCTATCGATGGCCAGATCGAAAAGGGCTAAATCACGAGTTTCCTCCGCGATTTGGAGCCTTACCCGAATGGCCCAGATATCTTTCAGt                                                                                            >  1:44900/1‑150 (MQ=255)
                                                              ccATGAGCAACGTCGCGTACGCGAAGCTTGGTCAGATCGCAGGCTCGTAGCTTGAGCTTGCTATCGATGGCCAGATCGAAAAGGGCTAAATCACGAGTTTCCTCCGCGATTTGGAGCCTTACCCGAATGGCCCAGATATCTTTCAGTcgg                                                                                         >  2:1286763/1‑150 (MQ=255)
                                                                           cgcgTACGCGAAGCTTGGTCAGATCGCAGGCTCGTAGCTTGAGCTTGCTATCGATGGCCAGATCGAAAAGGGCTAAATCACGAGTTTCCTCCGCGATTTGGAGCCTTACCCGAATGGCCCAGATATCTTTCAGTCGGAGCGGAGCTTTCt                                                                            <  1:226331/150‑1 (MQ=255)
                                                                             acaaaggaaaattggtcgggcagaaagctccgctccgactgaaagatatctgggccattcgggtaaggctccAAATCACGAGTTTCCTCCGCGATTTg                                                                                                                              <  2:1211914‑M1/26‑1 (MQ=255)
                                                                             acaaaggaaaattggtcgggcagaaagctccgctccgactgaaagatatctgggccattcgggtaaggctccAAATCACGAGTTTCCTCCGCGATTTg                                                                                                                              >  1:1211914‑M1/73‑98 (MQ=255)
                                                                                  gcgAAGCTTGGTCAGATCGCAGGCTCGTAGCATGAGCTTGCTATCGATGGCCAGATCGAAAAGGGCTAAATCACGAGTTTCCTCCGCGATTTGGAGCCTTACCCGAATGGCCCAGATATCTTTCAGTCGGAGCGGAGCTTTCTGCCCGAc                                                                     >  2:873822/1‑150 (MQ=255)
                                                                                               aGATCGCAGGCTCGTAGCTTGAGCTTGCTATCGATGGCCAGATCGAAAAGGGCTAAATCACGAGTTTCCTCCGCGATTTGGAGCCTTACCCGAATGGCCCAGATATCTTTCAGTCGGAGCGGAGCTTTCTGCCCGACCAATTTTCCTTTg                                                        <  2:820188/150‑1 (MQ=255)
                                                                                                 aTCGCAGGCTCGTAGCTTGAGCTTGCTATCGATGGCCAGATCGAAAAGGGCTAAATCACGAGTTTCCTCCGCGATTTGGAGCCTTACCCGAATGGCCCAGATATCTTTCAGTCGGAGCGGAGCTTTCTGCCCGACCAATTTTCCTTTGtg                                                      <  2:973935/150‑2 (MQ=255)
                                                                                                    gCAGGCTCGTAGCTTGAGCTTGCTATCGATGGCCAGATCGAAAAGGGCTAAATCACGAGTTTCCTCCGCGATTTGGAGCCTTACCCGAATGGCCCAGATATCTTTCAGTCGGAGCGGAGCTTTCTGCCCGACCAATTTTCCTTTGTTCaa                                                   >  1:104603/1‑150 (MQ=255)
                                                                                                        gCTCGTAGCTTGAGCTTTCTATCGATGGCCAGATCGAAAAGGGCTAAATCACGAGTTTCCTCCGCGATTTGGAGCCTTACCCGAATGGCCCAGATATCTTTCAGTCGGAGCGGAGCTTTCTGCCCGACCAATTTTCCTTTGTTCAAAGGc                                               >  2:780863/1‑150 (MQ=255)
                                                                                                         cTCGTAGCTTGAGCTTGCTATCGATGGCCAGATCGAAAAGGGCTAAATCACGAGTTTCCTCCGCGATTTGGAGCCTTACCCGAATGGCCCAGATATCTTTCAGTCGGAGCGGAGCTTTCTGCCCGACCAATTTTCCTTTGTTCAAAGGCt                                              <  1:1286763/150‑1 (MQ=255)
                                                                                                          tCGTAGCTTGAGCTTGCTATCGATGGCCAGATCGAAAAGGGCTAAATCACGAGTTTCCTCCGCGATTTGGAGCCTTACCCGAATGGCCCAGATATCTTTCAGTCGGAGCGGAGCTTTCTGCCCGACCAATTTTCCTTTGTTCAAAGGCTg                                             <  2:44900/150‑1 (MQ=255)
                                                                                                                 ttgagctcgCTATCGATGGCCAGATCGAAAAGGGATAAATCACGAGTTTCCTCCGCGATTTGGAGCCTTACCCGAATGGCCCAGATATCTTTCAGTCGGAGCGGAGCTTTCTGCCCCACCAATTTTCCTTTGTTCAAAGGCTGATGgctg                                      <  2:296457/150‑1 (MQ=255)
                                                                                                                  tgagcttgCTATCGATGGCCAGATCGAAAAGGGCTAAATCACGAGTTTCCTCCGCGATTTGGAGCCTTACCCGAATGGCCCAGATATCTTTCAGTCGGAGCGGAGCTTTCTGCCCGACCAATTTTCCTTTGTTCAAAGGCTGATGgctgc                                     <  2:451096/150‑1 (MQ=255)
                                                                                                                           tATCGATGGCCAGATCGAAAAGGGCTAAATCACGAGTTTCCTCCGCGATTTGGAGCCTTACCCGAATGGCCCAGATATCTTTCAGTCGGAGCGGAGCTTTCTGCCCGACCAATTTTCCt                                                           <  2:230095/119‑1 (MQ=255)
                                                                                                                           tATCGATGGCCAGATCGAAAAGGGCTAAATCACGAGTTTCCTCCGCGATTTGGAGCCTTACCCGAATGGCCCAGATATCTTTCAGTCGGAGCGGAGCTTTCTGCCCGACCAATTTTCCt                                                           >  1:230095/1‑119 (MQ=255)
                                                                                                                                aTGGCCAGATCGAAAAGGGCTAAATCACGAGTTTCCTCCGCGATTTGGAGCCTTACCCGAATGGCTCAGATATCTTTCAGTCGGAGCGGAGCTTTCTGCCCCACCAATTTTCCTTTGTTCAAAGGCTGATGGCTGCTGGTCGTAATGGCg                       >  2:296124/1‑150 (MQ=255)
                                                                                                                                    ccAGATCGAAAAGGGCTAAATCACGAGTTTCCTCCGCGATTTGGAGCCTTACCCGAATGGCCCAGATATCTTTCAGTCGGAGCGGAGCTTTCTGccc                                                                        <  1:997111/97‑1 (MQ=255)
                                                                                                                                    ccAGATCGAAAAGGGCTAAATCACGAGTTTCCTCCGCGATTTGGAGCCTTACCCGAATGGCCCAGATATCTTTCAGTCGGAGCGGAGCTTTCTGccc                                                                        >  2:997111/1‑97 (MQ=255)
                                                                                                                                                gggCTAAATCACGAGTTTCCTCCGCGATTTGGAGCCTTACCCGAATGGCCCAGATATCTTTCAGTCGGAGCGGAGCTTTCTGCCCGACCAATTTTCCTTTGTTCAAAGGCTGATGGCTGCTGGTCGTAATGGCGTTCATGGCTTGtcctc       >  2:43915/1‑150 (MQ=255)
                                                                                                                                                   cTAAATCACGAGTTTCCTCCGCGATTTGGAGCCTTACCCGAATGGCCCAGATATCTTTCAGTCGGAGCGGAGCTTTCTGCCCGACCAATTTTCCTTTGTTCCAAGGCTGATGGCTGCTGGTCGTAATGGCGCTCATGGCTTGTCCTCCAt    <  1:409756/150‑1 (MQ=255)
                                                                                                                                                     AAATCATGAGTTTCCTCCTCGATTTTGAGCCTTACCCGAATGTCCCAGATATCTTTCAGTCGGGGCGAAGCTTTCTGCCCGACCAATTTTCCTTTGTTCAAAGGCTTATGGCTGCTGGTCCTAATGCCTTTCGTGGCTTGTCCTCcatat  >  2:83691‑M1/1‑148 (MQ=255)
                                                                                                                                                     AAATCACGAGTTTCCTCCGCGATTTGGAGCCTTACCCGAATGGCCCAGATATCTTTCAGTCGGAGCGGAGCTTTCTGCCCGACCAATTTTCCTTTGTTCAAAGGCTGATGGCTGCTGGTCGTAATGGCGTTCATGGCTTGTCCTCCAtgt  >  1:70008‑M1/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAAACTGTACGGGACGTTGGGTTTTCTGCTGCATCACTATCGCCCGCGATGACACGTGCTCCCCATGAGCAACGTCGCGTACGCGAAGCTTGGTCAGATCGCAGGCTCGTAGCTTGAGCTTGCTATCGATGGCCAGATCGAAAAGGGCTAAATCACGAGTTTCCTCCGCGATTTGGAGCCTTACCCGAATGGCCCAGATATCTTTCAGTCGGAGCGGAGCTTTCTGCCCGACCAATTTTCCTTTGTTCAAAGGCTGATGGCTGCTGGTCGTAATGGCGTTCATGGCTTGTCCTCCATGT  >  NC_002947/4320526‑4320824

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