New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 303119335 (0.590)3 (0.060) 3/148 3.7 8.6% coding (273/1734 nt) recJ single‑stranded‑DNA‑specific exonuclease RecJ
?NZ_CP009273 = 3031201 32 (0.590)coding (265/1734 nt) recJ single‑stranded‑DNA‑specific exonuclease RecJ
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCC  >  NZ_CP009273/3031115‑3031289
                                                                                        |                                                                                        
tAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGt                                                                                       >  1:1119022/1‑90 (MQ=255)
            gTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGcc                                                                           >  1:1453079/1‑90 (MQ=255)
                    aaCGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGt                                                                   <  1:212518/90‑1 (MQ=255)
                                 cAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGAc                                                      <  1:1731241/90‑1 (MQ=255)
                                       gTCGATATTGGCGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGGGCTGGTCGCGCCGCCGGCGTCGAAATCACCGACCACaat                                                >  2:576537/1‑90 (MQ=255)
                                       gTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACaat                                                <  1:124367/90‑1 (MQ=255)
                                       gTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACaat                                                >  2:609075/1‑90 (MQ=255)
                                       gTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACaat                                                >  2:1781256/1‑90 (MQ=255)
                                       gTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACaat                                                >  2:1682252/1‑90 (MQ=255)
                                       gTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACaat                                                <  1:483139/90‑1 (MQ=255)
                                              ttGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCcgc                                         >  2:924728/1‑90 (MQ=255)
                                                gctgcAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGt                                       >  1:1165401/1‑90 (MQ=255)
                                                     aaCCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCg                                                        >  1:202394/1‑68 (MQ=255)
                                                     aaCCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCg                                                        <  2:202394/68‑1 (MQ=255)
                                                     aaCCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACaata                                               <  1:583251/77‑1 (MQ=255)
                                                     aaCCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACaata                                               >  2:583251/1‑77 (MQ=255)
                                                     aaCCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCtt                                  <  2:1217948/90‑1 (MQ=255)
                                                     aaCCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCtt                                  >  2:1319298/1‑90 (MQ=255)
                                                     aaCCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCtt                                  >  1:165446/1‑90 (MQ=255)
                                                        cAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTcgc                               >  1:815351/1‑90 (MQ=255)
                                                            cgagcgCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGaaa                           <  2:629687/90‑1 (MQ=255)
                                                              agcgCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAg                         <  1:982540/90‑1 (MQ=255)
                                                               gcgcATCGCCAGCACGGCTCTAAGCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTAGCGAAAAGc                        >  1:73103/1‑90 (MQ=37)
                                                               gcgcATCGCCAGCACGGCTCTAAGCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTAGCGAAAAGc                        <  2:73103/90‑1 (MQ=37)
                                                               gcgcATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGc                        >  1:460510/1‑90 (MQ=255)
                                                                cgcATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCg                       <  1:255154/90‑1 (MQ=255)
                                                                 gcATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGt                      <  1:1286930/90‑1 (MQ=255)
                                                                 gcATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGt                      <  2:1165401/90‑1 (MQ=255)
                                                                  cATCGCCAGCACGCTTAGAGCCGTGCTGGTTGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGtt                     >  2:1383604/1‑90 (MQ=255)
                                                                     cGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTa                  <  1:609075/90‑1 (MQ=255)
                                                                      gCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTaa                 <  1:1781256/90‑1 (MQ=255)
                                                                                   tAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAAc      <  2:1453079/90‑1 (MQ=255)
                                                                                       gCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGcc  >  1:946128/1‑90 (MQ=255)
                                                                                        |                                                                                        
TAAGCCGTAACCGTCTTCGAAACGGTTTGGTACCAGGTAGTCGATATTGCTGCAACCAAGCGAGCGCATCGCCAGCACGCTTAGAGCCGTGCTGGTCGCGCCGTCGGCGTCGAAATCACCGACCACAATAATCCGCGTTCCTTCGCGAAAAGCGTTGTAAAGGATCTCAACGGCC  >  NZ_CP009273/3031115‑3031289

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