New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 387401538 (0.940)3 (0.080) 3/150 2.8 8.1% coding (567/1074 nt) recF DNA replication/repair protein RecF
?NZ_CP009273 = 3874077 35 (0.990)coding (505/1074 nt) recF DNA replication/repair protein RecF
Rejected: Frequency below/above cutoff threshold.

AGTCAGAGAAAACTCAGGGAGAAATTGCTTACAGGTATCAGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGC  >  NZ_CP009273/3873931‑3874102
                                                                                    |                                                                                       
aGTCAGAGAAAACTCAGGGGGAAATTGCTTACAGGTATCAGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGAt                                                                                    >  2:306419/1‑90 (MQ=255)
aGTCAGAGAAAACTCAGGGAGAAATTGCTTACAGGTATCAGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGAt                                                                                    >  1:449475/1‑90 (MQ=255)
         aaaCTCAGGGAGAAATTGCTTACAGGTATCAGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGc                                                                           <  2:737055/90‑1 (MQ=255)
              cAGGGAGAAATTGCTTACAGGTATCAGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCg                                                                      <  2:990114/90‑1 (MQ=255)
               aGGGAGAAATTGCTTACAGGTATCAGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCgg                                                                     >  2:85861/1‑90 (MQ=255)
                    gaAATTGCTTACAGGTATCAGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCa                                                                <  1:322619/90‑1 (MQ=255)
                    gaAATTGCTTACAGGTATCAGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCa                                                                >  1:472775/1‑90 (MQ=255)
                      aaTTGCTTACAGGTATCAGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGacacgttacgaacagctacgcccgtgg                                                              >  2:123621‑M1/1‑63 (MQ=255)
                      aaTTGCTTACAGGTATCAGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGc                                                              <  2:472775/90‑1 (MQ=255)
                      aaTTGCTTACAGGTATCAGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCCCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGc                                                              <  1:333773/90‑1 (MQ=255)
                          gCTTACAGGTATCAGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCtt                                                          >  2:813359/1‑90 (MQ=255)
                                 ggTATCAGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCa                                                   <  2:370902/90‑1 (MQ=255)
                                    aTCAGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGacacgttacgaacagccacgcccgtgggataaagagctgat                                                >  1:684688‑M1/1‑49 (MQ=255)
                                    aTCAGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACgg                                                >  2:404030/1‑90 (MQ=255)
                                    aTCAGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACgg                                                >  1:472852/1‑90 (MQ=255)
                                    aTCAGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACgg                                                >  1:546050/1‑90 (MQ=255)
                                           atatCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGc                                                              >  2:958285/1‑69 (MQ=255)
                                           atatCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGc                                                              <  1:958285/69‑1 (MQ=255)
                                           atatCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGc                                         >  2:588784/1‑90 (MQ=255)
                                           atatCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGc                                         <  1:27551/90‑1 (MQ=255)
                                           atatCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGc                                         >  2:157216/1‑90 (MQ=255)
                                           atatCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGc                                         >  2:563457/1‑90 (MQ=255)
                                           atatCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGc                                         >  2:717313/1‑90 (MQ=255)
                                           atatCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGc                                         >  1:461233/1‑90 (MQ=255)
                                           atatCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGc                                         >  1:615947/1‑90 (MQ=255)
                                           atatCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCCGATTTGCTCCGCCAGCGGGATCCGCTCTTTATCCCACGGGCGTaac                                         >  1:271357/1‑88 (MQ=255)
                                                ggCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTc                                    >  2:334589/1‑90 (MQ=255)
                                                  ccGCGATACCGGCGCTATACTCCGCGCGCCAGGTGacacgttacgaacagctacgcccgtgggataaagagctgatcccgctgccggagc                                  <  1:123621‑M1/90‑56 (MQ=255)
                                                   cgcgATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTa                                 <  1:950128/90‑1 (MQ=255)
                                                              cgcTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTcacct                      >  1:90165‑M2/1‑85 (MQ=255)
                                                               gcTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTcacctg                     >  2:582591‑M2/1‑84 (MQ=255)
                                                               gcTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTcacctg                     <  1:743425‑M2/90‑7 (MQ=255)
                                                                 tataCTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTcacctggc                   <  2:208291‑M2/90‑9 (MQ=255)
                                                                  ataCTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTcacctggcg                  <  1:296758‑M2/90‑10 (MQ=255)
                                                                   taCTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTcacctggcgc                 >  2:4922‑M2/1‑80 (MQ=255)
                                                                       ccGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATc                                                                 >  1:234315/1‑38 (MQ=255)
                                                                       ccGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATc                                                                 <  2:234315/38‑1 (MQ=255)
                                                                       ccGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAgcg             >  2:647294/1‑90 (MQ=255)
                                                                         gcgcgcCAGGTGCTGATTTGCTCCGCCAGCGGGATc                                                                 <  1:73763/36‑1 (MQ=255)
                                                                         gcgcgcCAGGTGCTGATTTGCTCCGCCAGCGGGATc                                                                 >  2:73763/1‑36 (MQ=255)
                                                                         gcgcgcCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGcgcc           >  2:1147296/1‑90 (MQ=255)
                                                                         gcgcgcCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGcgcc           >  2:847203/1‑90 (MQ=255)
                                                                         gcgcgcCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGcgcc           >  2:831334/1‑90 (MQ=255)
                                                                         gcgcgcCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGcgcc           >  1:949197/1‑90 (MQ=255)
                                                                           gcgcCAGGTGacacgttacgaacagctacgcccgtgggataaagagctgatcccgctggcggagcaaatcagcacctggcgcgcggagta         >  2:4922‑M1/1‑10 (MQ=255)
                                                                            cgcCAGGTGacacgttacgaacagctacgcccgtgggataaagagctgatcccgctggcggagcaaatcagcacctggcgcgcggagtat        <  1:296758‑M1/90‑82 (MQ=255)
                                                                             gcCAGGTGacacgttacgaacagctacgcccgtgggataaagagctgatcccgctggcggagcaaatcagcacctggcgcgcggagtata       <  2:208291‑M1/90‑83 (MQ=255)
                                                                              ccAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCAtt      >  1:6009/1‑90 (MQ=255)
                                                                               cAGGTGacacgttacgaacagctacgcccgtgggataaagagctgatcccgctggcggagcaaatcagcacctggcgcgcggagtatagc     <  1:743425‑M1/90‑85 (MQ=255)
                                                                               cAGGTGacacgttacgaacagctacgcccgtgggataaagagctgatcccgctggcggagcaaatcagcacctggcgcgcggagtatagc     >  2:582591‑M1/1‑6 (MQ=255)
                                                                                aGGTGacacgttacgaacagctacgcccgtgggataaagagctgatcccgctggcggagcaaatcagcacctggcgcgcggagtatagcg    >  1:90165‑M1/1‑5 (MQ=255)
                                                                                 ggTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTgcg   <  2:1021101/90‑1 (MQ=255)
                                                                                  gTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTgcgc  <  1:336297/90‑1 (MQ=255)
                                                                                    |                                                                                       
AGTCAGAGAAAACTCAGGGAGAAATTGCTTACAGGTATCAGCCATATCGGCCGCGATACCGGCGCTATACTCCGCGCGCCAGGTGCTGATTTGCTCCGCCAGCGGGATCAGCTCTTTATCCCACGGGCGTAGCTGTTCGTAACGTGTCACCTGGCGCAGCGCCGCATTGCGC  >  NZ_CP009273/3873931‑3874102

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