New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 1955597 =25 (0.840)3 (0.100) 3/168 2.4 7.7% coding (1279/1734 nt) argS arginine‑‑tRNA ligase
?NZ_CP009273 1955639 = 47 (1.570)coding (1321/1734 nt) argS arginine‑‑tRNA ligase
Rejected: Frequency below/above cutoff threshold.

CTGGAACGTGCACGCCGTCTGGTGGCAGAAAAGAACCCGGATATGCCAGCCGACGAGCTGGAAAAACTGGCTAACGCGGTTGGTATTGGTGCGGTGAAATATGCGGATCTCTCCAAAAACCGCACCACGGACTACATCTTCGACTGGGACAACATGCTGGCGTTTGAGGGTAATAC  >  NZ_CP009273/1955510‑1955685
                                                                                       |                                                                                        
cTGGAACGTGCACGCCGTCTGGTGGCAGAAAAGAACCCGGATATGCCAGCCGACGAGCTGGAAAAACTAGCTAACGCGGTTGGTATTGGt                                                                                        >  2:138673/1‑90 (MQ=255)
       gTGCACGCCGTCTGGTGGCAGAAAAGAACCCGGATATGCCAGCCGACGAGCTGGAAAAACTGGCTAACGCGGTTGGTATTGGTGCGGTGa                                                                                 >  1:414292/1‑90 (MQ=255)
       gTGCACGCCGTCTGGTGGCAGAAAAGAACCCGGATATGCCAGCCGACGAGCTGGAAAAACTGGCTAACGCGGTTGGTATTGGTGCGGTGa                                                                                 >  2:71337/1‑90 (MQ=255)
         gCACGCCGTCTGGTGGCAGAAAAGAACCCGGATATGCCAGCCGACGAGCTGGAAAAACTGGCTAACGCGGTTGGTATTGGTGCGGTGaaa                                                                               <  2:668403/90‑1 (MQ=255)
         gCACGCCGTCTGGTGGCAGAAAAGAACCCGGATATGCCAGCCGACGAGCTGGAAAAACTGGCTAACGCGGTTGGTATTGGTGCGGTGaaa                                                                               <  2:707111/90‑1 (MQ=255)
         gCACGCCGTCTGGTGGCAGAAAAGAACCCGGATATGCCAGCCGACGAGCTGGAAAAACTGGCTAACGCGGTTGGTATTGGTGCGGTGaaa                                                                               <  2:828417/90‑1 (MQ=255)
                 tCTGGTGGCAGAAAAGAACCCGGATATGCCAGCCGACGAGCTGGAAAAACTGGCTAACGCGGTTGGTATTGGTGCGGTGAAATATGCGGa                                                                       >  2:11218/1‑90 (MQ=255)
                 tCTGGTGGCAGAAAAGAACCCGGATATGCCAGCCGACGAGCTGGAAAAACTGGCTAACGCGGTTGGTATTGGTGCGGTGAAATATGCGGa                                                                       >  2:456294/1‑90 (MQ=255)
                  cTGGTGGCAGAAAAGAACCCGGATATGCCAGCCGACGAGCTGGAAAAACTGGCTAACGCGGTTGGTATTGGTGCGGTGAAATATGCGGAt                                                                      >  1:762772/1‑90 (MQ=255)
                           gtatggcgcggtattaccctcaaacgccagcatgttgtcccagtcgaagatgtagtccgtGGTGCGGTGAAATATGCGGATCTCTCCaaa                                                             <  2:421293‑M1/30‑1 (MQ=255)
                            aaaaGAACCCGGATATGCCAGCCGACGAGCTGGAAAAACTGGCTAACGCGGTTGGTATTGGTGCGGTGAAATATGCGGATCTCTCCaaaa                                                            >  2:7669/1‑90 (MQ=255)
                              aaGAACCCGGATATGCCAGCCGACGAGCTGGAAAAACTGGCTAACGCGGTTGGTATTGGTGCGGTGAAATATGCGGATCTCTCCAAAAAc                                                          <  2:395593/90‑1 (MQ=255)
                                  ggcgcggtattaccctcaaacgccagcatgttgtcccagtcgaagatgtagtcGGTGCGGTGAAATATGCGGATCTCTCCAAAAACCGca                                                      <  1:553150‑M1/37‑1 (MQ=255)
                                     cGGATATGCCAGCCGACGAGCTGGAAAAACTGGCTAACGCGGTTGGTATTGGTGCGGTGAAATATGCGGATCTCTCCAAAAACCGcacca                                                   >  2:777714/1‑90 (MQ=255)
                                      ggATATGCCAGCCGACGAGCTGGAAAAACTGGCTAACGCGGTTGGTATTGGTGCGGTGAAATATGCGGATCTCTCCAAAAACCGcaccac                                                  <  1:789231/90‑1 (MQ=255)
                                           tGCCAGCCGACGAGCTGGAAAAACTGGCTAACGCGGTTGGTATTGGTGCGGTGAAATATGCGGATCTCTCCAAAAACCGCACCACGGACt                                             >  1:506524/1‑90 (MQ=255)
                                             ccAGCCGACGAGCTGGAAAAACTGGCTAACGCGGTTGGTATTGGTGCGGTGAAATATGCGGATCTCTCCAAAAACCGCACCACGGACTAc                                           >  1:816091/1‑90 (MQ=255)
                                              cAGCCGACGAGCTGGAAAAACTGGCTAACGCGGTTGGTATTGGTGCGGTGAAATATGCGGATCTCTCCAAAAACCGCACCACGGACTACa                                          >  2:720894/1‑90 (MQ=255)
                                              cAGCCGACGAGCTGGAAAAACTGGCTAACGCGGTTGGTATTGGTGCGGTGAAATATGCGGATCTCTCCAAAAACCGCACCACGGACTACa                                          <  2:833983/90‑1 (MQ=255)
                                                 ccGACGAGCTGGAAAAACTGGCTAACGCGGTTGGTATTGGTGCGGTGAAATATGCGGATCTCTCCAAAAACCGCACCACGGACTACATCt                                       >  1:526778/1‑90 (MQ=255)
                                                 ccGACGAGCTGGAAAAACTGGCTAACGCGGTTGGTATTGGTGCGGTGAAATATGCGGATCTCTCCAAAAACCGCACCACGGACTACATCt                                       >  1:120035/1‑90 (MQ=255)
                                                         cTGGAAAAACTGGCTAACGCGGTTGGTATTGGTGCGGTGAAATATGCGGATCTCTCCAAAAACCGCACCACGGACTACATCTTCGACTgg                               <  1:156450/90‑1 (MQ=255)
                                                              tgttgtcccagtcgaagatgtagtcGGTGCGGTGAAATATGCGGATCTCTCCAAAAACCGCACCACGGACTACATCTTCGACTGGGacaa                          >  2:553150‑M1/26‑90 (MQ=255)
                                                                      cTAACGCGGTTGGTATTGGTGCGGTGAAATATGCGGATCTCTCCAAAAACCGCACCACGGACTACATCTTCGACTGGGACAACATGCTgg                  >  2:193937/1‑90 (MQ=255)
                                                                            cgGTTGGTATTGGTGCGGTGAAATATGCGGATCTCTCCAAAAACCGCACCACGGACTACATCTTCGACTGGGACAACATGCTGGCGTTTg            <  2:609528/90‑1 (MQ=255)
                                                                            cgGTTGGTATTGGTGCGGTGAAATATGCGGATCTCTCCAAAAACCGCACCACGGACTACATCTTCGACTGGGACAACATGCTGGCGTTTg            <  1:777714/90‑1 (MQ=255)
                                                                                      tGGTGCGGTGAAATATGCGGATCTCTCCAAAAACCGCACCACGGACTACATCTTCGACTGGGACAACATGCTGGCGTTTGAGGGTAATAc  >  1:538553/1‑90 (MQ=255)
                                                                                      tGGTGCGGTGAAATATGCGGATCTCTCCAAAAACCGCACCACGGACTACATCTTCGACTGGGACAACATGCTGGCGTTTGAGGGTAATAc  <  2:816091/90‑1 (MQ=255)
                                                                                       ggTGCGGTGAAATATGCGGATCTCTCCAAAAACCGCACCACGGACTACATCTTCGACTGGGACAACATGCTGGCGTTTg            <  1:723722/79‑1 (MQ=255)
                                                                                       ggTGCGGTGAAATATGCGGATCTCTCCAAAAACCGCACCACGGACTACATCTTCGACTGGGACAACATGCTGGCGTTTg            >  2:723722/1‑79 (MQ=255)
                                                                                       |                                                                                        
CTGGAACGTGCACGCCGTCTGGTGGCAGAAAAGAACCCGGATATGCCAGCCGACGAGCTGGAAAAACTGGCTAACGCGGTTGGTATTGGTGCGGTGAAATATGCGGATCTCTCCAAAAACCGCACCACGGACTACATCTTCGACTGGGACAACATGCTGGCGTTTGAGGGTAATAC  >  NZ_CP009273/1955510‑1955685

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