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

CTGCTGCATCACTATCGCCCGCGATGACACGTGCTCCCCATGAGCAACGTCGCGTACGCGAAGCTTGGTCAGATCGCAGGCTCGTAGCTTGAGCTTGCTATCGATGGCCAGATCGAAAAGGGCTAAATCACGAGTTTCCTCCGCGATTTGGAGCCTTACCCGAATGGCCCAGATATCTTTCAGTCGGAGCGGAGCTTTCTGCCCGACCAATTTTCCTTTGTTCAAAGGCTGATGGCTGCTGGTCGTAATGGCGTTCATGGCTTGTCCTCCATGTGAGGTAGGACAAGGATGATCAACCC  >  NC_002947/4320551‑4320849
                                                                                                                                                     |                                                                                                                                                     
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                                                                                    >  1:230095/1‑119 (MQ=255)
                                                                                                  tATCGATGGCCAGATCGAAAAGGGCTAAATCACGAGTTTCCTCCGCGATTTGGAGCCTTACCCGAATGGCCCAGATATCTTTCAGTCGGAGCGGAGCTTTCTGCCCGACCAATTTTCCt                                                                                    <  2:230095/119‑1 (MQ=255)
                                                                                                       aTGGCCAGATCGAAAAGGGCTAAATCACGAGTTTCCTCCGCGATTTGGAGCCTTACCCGAATGGCTCAGATATCTTTCAGTCGGAGCGGAGCTTTCTGCCCCACCAATTTTCCTTTGTTCAAAGGCTGATGGCTGCTGGTCGTAATGGCg                                                >  2:296124/1‑150 (MQ=255)
                                                                                                           ccAGATCGAAAAGGGCTAAATCACGAGTTTCCTCCGCGATTTGGAGCCTTACCCGAATGGCCCAGATATCTTTCAGTCGGAGCGGAGCTTTCTGccc                                                                                                 >  2:997111/1‑97 (MQ=255)
                                                                                                           ccAGATCGAAAAGGGCTAAATCACGAGTTTCCTCCGCGATTTGGAGCCTTACCCGAATGGCCCAGATATCTTTCAGTCGGAGCGGAGCTTTCTGccc                                                                                                 <  1:997111/97‑1 (MQ=255)
                                                                                                                gggtaaggctccaaatcgcggaggaaactcgtgatttGGAGCCTTACCCGAATGGACCAGATATCTTTCAGTCGGAGCGGAGCTTTCTGCCCGACCAATTTTCCTTTGTTCAAAGGCTGATGGCTGCTGGTCGTAATGGCGTTCATGGCt                                       >  2:209875‑M2/38‑150 (MQ=255)
                                                                                                                       gggCTAAATCACGAGTTTCCTCCGCGATTTGGAGCCTTACCCGAATGGCCCAGATATCTTTCAGTCGGAGCGGAGCTTTCTGCCCGACCAATTTTCCTTTGTTCAAAGGCTGATGGCTGCTGGTCGTAATGGCGTTCATGGCTTGtcctc                                >  2:43915/1‑150 (MQ=255)
                                                                                                                          cTAAATCACGAGTTTCCTCCGCGATTTGGAGCCTTACCCGAATGGCCCAGATATCTTTCAGTCGGAGCGGAGCTTTCTGCCCGACCAATTTTCCTTTGTTCCAAGGCTGATGGCTGCTGGTCGTAATGGCGCTCATGGCTTGTCCTCCAt                             <  1:409756/150‑1 (MQ=255)
                                                                                                                           caaatcgcggaggaaactcgtgatttGGAGCCTTACCCGAATGGCCCAGATATCTTTCAGTCGGAGCGGAGCTTTCTGCCCGACCAATTTTCCTTTGt                                                                                >  1:1211914‑M2/27‑98 (MQ=255)
                                                                                                                           caaatcgcggaggaaactcgtgatttGGAGCCTTACCCGAATGGCCCAGATATCTTTCAGTCGGAGCGGAGCTTTCTGCCCGACCAATTTTCCTTTGt                                                                                <  2:1211914‑M2/72‑1 (MQ=255)
                                                                                                                            AAATCATGAGTTTCCTCCTCGATTTTGAGCCTTACCCGAATGTCCCAGATATCTTTCAGTCGGGGCGAAGCTTTCTGCCCGACCAATTTTCCTTTGTTCAAAGGCTTATGGCTGCTGGTCCTAATGCCTTTCGTGGCTTGTCCTCcatat                           >  2:83691‑M1/1‑148 (MQ=255)
                                                                                                                            AAATCACGAGTTTCCTCCGCGATTTGGAGCCTTACCCGAATGGCCCAGATATCTTTCAGTCGGAGCGGAGCTTTCTGCCCGACCAATTTTCCTTTGTTCAAAGGCTGATGGCTGCTGGTCGTAATGGCGTTCATGGCTTGTCCTCCAtgt                           >  1:70008‑M1/1‑150 (MQ=255)
                                                                                                                                   ggaggaaactcgtgatttGGAGCCTTACCCGAATGGCCCAGATATCTTTCAGTCGGAGCGGAGCTTTCTGCCCGACCAATTTTCCTTTGTTCAAAGGCTGATGGCTGCTGGTCGTAATGGCGTTCATGGCTTGTCCTCCATGTGAGGTAg                    <  1:209875‑M2/132‑1 (MQ=255)
                                                                                                                                   gAGTTTTCTCTGCTATTTGGAGCCTTACCCGAATGGCCCATATATCTTTCAGTCGGAGCGGAGCTTTCTGCCCGACCAATTTTCCTTTGTTCAAAGGCTGATGGCTGCTGGTCGTAATGGCGTTCATGGCTTGTCCTCCATGTGAGGTAg                    <  2:104603/150‑1 (MQ=255)
                                                                                                                                           tccGCGATTTGGAGCCTTACCCGAATGGCCCAGATATCTTTCAGTCGGAGCGGAGCTTTCTGCCCGACCAATTTTCCTTTGTTCAAAGGCTGATGGCTGCTGGTCGTAATGGCGTTCATGGCTTGTCCTCCATGTGAGGTAGGACAAGga            >  1:1332424/1‑150 (MQ=255)
                                                                                                                                                   ttGGAGCCTTACCCGAATGGCCCAGATATCTTTCAGTCGGAGCGGAGCTTTCTGCCCGACCAATTTTCCTTTGTTCAAAGGCTGATGGCTGCTGGTCGTAATGGCGTTCATGGCTTGTCCTCCATGTGAGGTAGGACAAGGATGATCAAc    <  1:43915‑M2/148‑1 (MQ=255)
                                                                                                                                                    tGGAGCCTTACCCGAATGGCCCAGATATCTTTCAGTCGGAGCGGAGCTTTCTGCCCGACCAATTTTCCTTTGTTCAAAGGCTGATGGCTGCTGGTCGTAATGGCGTTCATGGCTTGTCCTCCATGTGAGGTAGGACAAGGATGATCAAcc   <  1:171180‑M2/149‑1 (MQ=255)
                                                                                                                                                     ggAGCCTTACCCGAATGGCCCAGATATCTTTCAGTCGGAGCGGAGCTTTCTGCCCGACCAATTTTCCTTTGTTCAAAGGCTGATGGCTGCTGGTCGTAATGGCGTTCATGGCTTGTCCTCCATGTGAGGTAGGACAAGGATGATCAAccc  <  2:1332424/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTGCTGCATCACTATCGCCCGCGATGACACGTGCTCCCCATGAGCAACGTCGCGTACGCGAAGCTTGGTCAGATCGCAGGCTCGTAGCTTGAGCTTGCTATCGATGGCCAGATCGAAAAGGGCTAAATCACGAGTTTCCTCCGCGATTTGGAGCCTTACCCGAATGGCCCAGATATCTTTCAGTCGGAGCGGAGCTTTCTGCCCGACCAATTTTCCTTTGTTCAAAGGCTGATGGCTGCTGGTCGTAATGGCGTTCATGGCTTGTCCTCCATGTGAGGTAGGACAAGGATGATCAACCC  >  NC_002947/4320551‑4320849

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.