New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 248296 =16 (0.520)3 (0.110) 3/156 3.0 13.8% coding (912/1056 nt) dinB DNA polymerase IV
?NZ_CP009273 248360 = 23 (0.820)coding (976/1056 nt) dinB DNA polymerase IV
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CAAACCACCCAGGAGCACGTCTGGCCGCGGCTGAATAAAGCTGATCTAATCGCCACCGCGCGTAAAACCTGGGATGAACGCCGCGGCGGGCGCGGTGTGCGTCTGGTGGGGCTGCATGTGACGTTGCTTGACCCGCAAATGGAAAGACAACTGGTGCTGGGATTATGATGTATAC  >  NZ_CP009273/248273‑248447
                                                                                       |                                                                                       
cAAACCACCCAGGAGCACGTCTGGCCGCGGCTGAATAAAGCTGATCTAATCGCCACCGCGCGTAAAACCTGGGATGAACGCCGCggcgggc                                                                                      >  2:743397/1‑91 (MQ=255)
 aaaCCACCCAGGAGCACGTCTGGCCGCGGCTGAATAAAGCTGATCTAATCGCCACCGCGCGTAAAACCTGGGATGAACGCCGCggcgggcg                                                                                     <  1:743397/91‑1 (MQ=255)
  aaCCACCCAGGAGCACGTCTGGCCGCGGCTGAATAAAGCTGATCTAATCGCCACCGCGCGTAAAACCTGGGATGAACGCCGCGGCGGgcgc                                                                                    <  2:470152/91‑1 (MQ=255)
  aaCCACCCAGGAGCACGTCTGGCCGCGGCTGAATAAAGCTGATCTAATCGCCACCGCGCGTAAAACCTGGGATGAACGCCGCGGCGGgcgc                                                                                    <  2:376092/91‑1 (MQ=255)
                 cGTCTGGCCGCGGCTGAATAAAGCTGATCTAATCGCCACCGCGCGTAAAACCTGGGATGAACGCCGCAGCGGGCGCGGTGTGCGTCtggtg                                                                     >  1:496226/1‑91 (MQ=255)
                     tGGCCGCGGCTGAATAAAGCTGATCTAATCGCCACCGCGCGTAAAACCTGGGATGAACGCCGCGGCGGGCGCGGTGTGCGTCTGGTGGGGc                                                                 <  2:836481/91‑1 (MQ=255)
                           cgGCTGAATAAAGCTGATCTAATCGCCACCGCGCGTAAAACCTGGGATGAACGCCGCGGCGGGCGCGGTGTGCGTCTGGTGGGGCTGCAtg                                                           <  2:752763/91‑1 (MQ=255)
                           cgGCTGAATAAAGCTGATCTAATCGCCACCGCGCGTAAAACCTGGGATGAACGCCGCGGCGGGCGCGGTGTGCGTCTGGTGGGGCTGCAtg                                                           <  1:239316/91‑1 (MQ=255)
                           cgGCTGAATAAAGCTGATCTAATCGCCACCGCGCGTAAAACCTGGGATGAACGCCGCGGCGGGCGCGGTGTGCGTCTGGTGGGGCTGCAtg                                                           <  1:386923/91‑1 (MQ=255)
                                          gATCTAATCGCCACCGCGCGTAAAACCTGGGATGAACGCCGCGGCGGGCGCGGTGTGCGTCTGGTGGGGCTGCATGTGACGTTGCTTGAcc                                            >  2:424218/1‑91 (MQ=255)
                                                    ccACCGCGCGTAAAACCTGGGATGAACGCCGCGGCGGGCGCGGTGTGCGTCTGGTGGGGCTGCATGTGACGTTGCTTGACCCGCAAATGGa                                  >  2:100995/1‑91 (MQ=255)
                                                    ccACCGCGCGTAAAACCTGGGATGAACGCCGCGGCGGGCGCGGTGTGCGTCTGGTGGGGCTGCATGGGACGTTGCTTGACCCGCAAATGGa                                  >  2:104967/1‑91 (MQ=255)
                                                    ccACCGCGCGTAAAACCTGGGATGAACGCCGCGGCGGGCGCGGGGTGCGTCTGGTGGGGCTGCATGTGACGTTGCTTGACCCGCAAATGGa                                  >  2:156417/1‑91 (MQ=255)
                                                      aCCGCGCGTAAAACCTGGGATGAACGCCGCGGCGGGCGCGGTGTGCGTCTGGTGGGGCTGCATGTGACGTTGCTTGACCCGCAAATGGaaa                                >  1:125705/1‑91 (MQ=255)
                                                         gcgcgTAAAACCTGGGATGAACGCCGCGGCGGGCGCGGAGGGCGTCTGGTGGGGCTGCATGTGACGGTGCTTGACCCGCAAATGGAAAGAc                             >  1:756210/1‑91 (MQ=255)
                                                                   ccTGGGATGAACGCCGCGGCGGGCGCGGTGTGCGTCTGGTGGGGCTGCATGTGACGTTGCTTGACCCGCAAATGGAAAGACAACTGGTGCt                   <  1:153611/91‑1 (MQ=255)
                                                                    cTGGGATGAACGCCGCGGCGGGCGCGGTGTGCGTCTGGTGGGGCTGCATGTGACGTTGCTTGACCCGCAAATGGAAAGACAACTGGTGCTg                  <  1:501093/91‑1 (MQ=255)
                                                                        gATGAACGCCGCGGCGGGCGCGGTGTGCGTCTGGTGGGGCTGCATGTGACGTTGCTTGACCCGCAAATGGAAAGACAACTGGTGCTGGGAt              <  2:587688/91‑1 (MQ=255)
                                                                        gATGAACGCCGCGGCGGGCGCGGTGTGCGTCTGGTGGGGCTGCATGTGACGTTGCTTGACCCGCAAATGGAAAGACAACTGGTGCTGGGAt              <  2:327251/91‑1 (MQ=255)
                                                                        gATGAACGCCGCGGCGGGCGCGGTGTGCGTCTGGTGGGGCTGCATGTGACGTTGCTTGACCCGCAAATGGAAAGACAACTGGTGCTGGGAt              <  1:254909/91‑1 (MQ=255)
                                                                           gAACGCCGCGGCGGGCGCGGTGTGCGTCTGGTGGGGc                                                                 <  1:814078/37‑1 (MQ=255)
                                                                           gAACGCCGCGGCGGGCGCGGTGTGCGTCTGGTGGGGc                                                                 >  2:814078/1‑37 (MQ=255)
                                                                           gAACGCCGCGGCGGGCGCGGTGTGCGTCTGGTGGGGCTGCATGTGACGTTGCTTGACCCGCAAATGGAAAGACAACTGGTGCTGGGATTat           <  1:424218/91‑1 (MQ=255)
                                                                               gccgcGGCGGGCGCGGTGTGCGTCTGGTGGGGCTGCATGTGACGTTGCTTGACCCGCAAATGGAAAGACAACTGGTGCTGGGATTatgatg       <  2:503840/91‑1 (MQ=255)
                                                                                  gcggcgGGCGCGGTGTGCGTCTGGTGGGGCTGCATGTGACGTTGCTTGACCCGCAAATGGAAAGACAACTGGTGCTGGGATTATGATGtat    <  1:762921/91‑1 (MQ=255)
                                                                                  gcggcgGGCGCGGTGTGCGTCTGGTGGGGCTGCATGTGACGTTGCTTGACCCGCAAATGGAAAGACAACTGGTGCTGGGATTATGATGtat    <  2:229223/91‑1 (MQ=255)
                                                                                  gcggcgGGCGCGGTGTGCGTCTGGTGGGGCTGCATGTGACGTTGCTTGACCCGCAAATGGAAAGACAACTGGTGCTGGGATTATGATGtat    <  2:100088/91‑1 (MQ=255)
                                                                                    ggcgggcgCGGTGTGCGTCTGGTGGGGCTGCATGTGACGTTGCTTGACCCGCAAATGGAAAGACAACTGGTGCTGGGATTATGATGTATAc  <  1:274780/91‑1 (MQ=255)
                                                                                    ggcgggcgCGGTGTGCGTCTGGTGGGGCTGCATGTGACGTTGCTTGACCCGCAAATGGAAAGACAACTGGTGCTGGGATTATGATGTATAc  <  1:322945/91‑1 (MQ=255)
                                                                                    ggcgggcgCGGTGTGCGTCTGGTGGGGCTGCATGTGACGTTGCTTGACCCGCAAATGGAAAGACAACTGGTGCTGGGATTATGATGTATAc  <  1:400963/91‑1 (MQ=255)
                                                                                    ggcgggcgCGGTGTGCGTCTGGTGGGGCTGCATGTGACGTTGCTTGACCCGCAAATGGAAAGACAACTGGTGCTGGGATTATGATGTATAc  <  2:100852/91‑1 (MQ=255)
                                                                                    ggcgggcgCGGTGTGCGTCTGGTGGGGCTGCATGTGACGTTGCTTGACCCGCAAATGGAAAGACAACTGGTGCTGGGATTATGATGTATAc  <  1:125850/91‑1 (MQ=255)
                                                                                    ggcgggcgCGGTGTGCGTCTGGTGGGGCTGCATGTGACGTTGCTTGACCCGCAAATGGAAAGACAACTGGTGCTGGGATTATGATGTATAc  <  1:122166/91‑1 (MQ=255)
                                                                                       |                                                                                       
CAAACCACCCAGGAGCACGTCTGGCCGCGGCTGAATAAAGCTGATCTAATCGCCACCGCGCGTAAAACCTGGGATGAACGCCGCGGCGGGCGCGGTGTGCGTCTGGTGGGGCTGCATGTGACGTTGCTTGACCCGCAAATGGAAAGACAACTGGTGCTGGGATTATGATGTATAC  >  NZ_CP009273/248273‑248447

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