New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 3574719 =29 (1.330)4 (0.200) 4/156 1.4 12.6% coding (222/489 nt) yhhY N‑acetyltransferase
?NZ_CP009273 3574782 = 29 (1.450)coding (285/489 nt) yhhY N‑acetyltransferase
Rejected: Frequency below/above cutoff threshold.

CTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTAACCGTGTTTGTCGATAAC  >  NZ_CP009273/3574693‑3574869
                                                                                         |                                                                                       
ctcACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTc                                                                                         <  1:201714/90‑1 (MQ=255)
    ccATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCa                                                                                     >  2:398782/1‑90 (MQ=255)
     tctcgcatcagggcgctggcgacgccgcggttcttccaacgagagtcgacacagataccaaaatcggcaacatgactgcggcgtCGCCAg                                                                                    <  2:124527‑M2/6‑1 (MQ=255)
     cATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAg                                                                                    <  1:615915/90‑1 (MQ=255)
      ctcgcatcagggcgctggcgACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAgc                                                                                   <  1:44455‑M1/70‑1 (MQ=255)
       ttGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAgcg                                                                                  >  2:29277/1‑90 (MQ=255)
       tcgcatcagggcgctggcgACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAgcg                                                                                  <  2:156546‑M1/71‑1 (MQ=255)
         gACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGcc                                                                                >  2:259450/1‑90 (MQ=255)
         gACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGcc                                                                                >  1:282967/1‑90 (MQ=255)
              gCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGAt                                                                           >  2:219199/1‑90 (MQ=255)
                  cgctggcgacgccgcggttcttccaacgagagtcgacacagataccaaaatcggcaacatgactgcggcgtCGCCAGCGCCCTGATGCga                                                                       <  2:156546‑M2/19‑1 (MQ=255)
                   gctggcgacgccgcggttcttccaacgagagtcgacacagataccaaaatcggcaacatgactgcggcgtCGCCAGCGCCCTGATGCgag                                                                      <  1:44455‑M2/20‑1 (MQ=255)
                   aGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCgag                                                                      >  1:412329/1‑90 (MQ=255)
                   aGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCgag                                                                      <  1:102601/90‑1 (MQ=255)
                    ctggcgACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCgaga                                                                     <  2:124527‑M1/84‑1 (MQ=255)
                     cgcCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCgagag                                                                    <  1:437724/90‑1 (MQ=255)
                     cgcCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCCGATGCgagag                                                                    <  1:418204/90‑1 (MQ=255)
                             CCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGa                                                            >  1:348120‑M1/1‑90 (MQ=255)
                                CAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAAt                                                         >  2:437059‑M1/1‑90 (MQ=255)
                                               ttttGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGtt                                          <  2:385385/90‑1 (MQ=255)
                                                tttGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTg                                         <  2:412329/90‑1 (MQ=255)
                                                 ttGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGGGAGATGATTGAAATGTGCGACAACTGGTTGc                                        >  1:463989/1‑90 (MQ=255)
                                                      aTCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCgagag                                                                    >  1:239994/1‑57 (MQ=255)
                                                      aTCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCgagag                                                                    <  2:239994/57‑1 (MQ=255)
                                                        cTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGa                                 >  1:213465/1‑90 (MQ=255)
                                                            gtCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGc                             >  2:312884/1‑90 (MQ=255)
                                                            gtCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGc                             >  1:342009/1‑90 (MQ=255)
                                                              cGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCAt                           <  2:389289/90‑1 (MQ=255)
                                                                 aaatcggcaacatgactgcggcgtCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGa                        <  1:329817‑M2/66‑1 (MQ=255)
                                                                  tctcGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGaa                       >  2:214134/1‑90 (MQ=255)
                                                                   ctcGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAAc                      <  2:235678/90‑1 (MQ=255)
                                                                   ctcGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAAc                      <  2:345619/90‑1 (MQ=255)
                                                                   ctcGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAAc                      <  2:403717/90‑1 (MQ=255)
                                                                    tcGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACt                     >  2:353014/1‑90 (MQ=255)
                                                                    tcGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACt                     >  1:329098/1‑90 (MQ=255)
                                                                               aCCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTTAACTAACCGTGTTTg          <  1:312884/90‑1 (MQ=255)
                                                                               aCCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTAACCGTGTTTg          >  2:593106/1‑90 (MQ=255)
                                                                                 cgcgGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTAACCGTGTTTGTc        <  2:11697/90‑1 (MQ=255)
                                                                                    ggcgtCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTAACCGTGTTTGTCGAt     >  2:357797‑M2/6‑90 (MQ=255)
                                                                                       gtCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTAACCGTGTTTGTCGATAAc  <  2:543935‑M2/88‑1 (MQ=255)
                                                                                         |                                                                                       
CTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTAACCGTGTTTGTCGATAAC  >  NZ_CP009273/3574693‑3574869

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