New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 19212522 (0.810)4 (0.160) 4/154 1.8 16.4% coding (733/759 nt) ispU (2E,6E)‑farnesyl‑diphosphate‑specific ditrans,polycis‑undecaprenyl‑diphosphate synthase
?NZ_CP009273 = 192230 21 (0.850)coding (67/858 nt) cdsA phosphatidate cytidylyltransferase
Rejected: Frequency below/above cutoff threshold.

CAGATGTTCTCTGGCCCGATTTCGATGAACAAGACTTTGAAGGGGCGTTAAATGCCTTTGCTAATCGAGAGCGTCGTTTCGGCGGCACCGAGCCCGGTGATGAAACAGCCTGATGGGGGTCGCTTTTGCTGAAGTATCGCCTGATATCTGCTTTTGTGTTAATACCCGTCGTC  >  NZ_CP009273/192039‑192211
                                                                                      |                                                                                      
cAGATGTTCTCTGGCCCGATTTCGATGAACAAGACTTTGAAGGGGCGTTAAATGCCTTTGCTAATCGAGAGCGTCGTTTCGGCGGCACCg                                                                                     >  1:686474/1‑90 (MQ=255)
cAGATGTTCTCTGGCCCGATTTCGATGAACAAGACTTTGAAGGGGCGTTAAATGCCTTTGCTAATCGAGAGCGTCGTTTCGGCGGCACCg                                                                                     >  1:133263/1‑90 (MQ=255)
      ttCTCTGGCCCGATTTCGATGAACAAGACTTTGAAGGGGCGTTAAATGCCTTTGCTAATCGAGAGCGTCGTTTCGGCGGCACCGAGCCCg                                                                               <  1:391635/90‑1 (MQ=255)
                cGATTTCGATGAACAAGACTTTGAAGGGGCGTTAAATGCCTTTGCTAATCGAGAGCGTCGTTTCGGCGGCACCGAGCCCGGTGATGAAAc                                                                     <  2:764261/90‑1 (MQ=255)
                cGATTTCGATGAACAAGACTTTGAAGGGGCGTTAAATGCCTTTGCTAATCGAGAGCGTCGTTTCGGCGGCACCGAGCCCGGTGATGAAAc                                                                     <  2:121263/90‑1 (MQ=255)
                 gATTTCGATGAACAAGACTTTGAAGGGGCGTTAAATGCCTTTGCTAATCGAGAGCGTCGTTTCGGCGGCAacagaaacaacgccgcgatg                                                                    >  1:406161‑M1/1‑70 (MQ=255)
                 gATTTCGATGAACAAGACTTTGAAGGGGCGTTAAATGCCTTTGCTAATCGAGAGCGTCGTTTCGGCGGCACCGAGCCCGGTGATGAAACa                                                                    <  2:714512/90‑1 (MQ=255)
                               aGACTTTGAAGGGGCGTTAAATGCCTTTGCTAATCGAGAGCGTCGTTTCGGCGGCACCGAGCCCGGTGATGAAACAGCCTGATGGGGGTc                                                      >  2:169399/1‑90 (MQ=255)
                                   tttGAAGGGGCGTTAAATGCCTTTGCTAATCGAGAGCGTCGTTTCGGCGGCACCGAGCCCGGTGATGAAACAGCCTGATGGGGGTCGCtt                                                  <  2:221500/90‑1 (MQ=255)
                                    ttGAAGGGGCGTTAAATGCCTTTGCTAATCGAGAGCGTCGTTTCGGCGGCACCGAGCCCGGTGATGAAACAGCCTGATGGGGGTCGCttt                                                 >  1:325944/1‑90 (MQ=255)
                                         ggggCGTTAAATGCCTTTGCTAATCGAGAGCGTCGTTTCGGCGGCACCGAGCCCGGTGATGAAACAGCCTGATGGGGGTCGCTTTTGCTg                                            <  2:29084/90‑1 (MQ=255)
                                         ggggCGTTAAATGCCTTTGCTAATCGAGAGCGTCGTTTCGGCGGCACCGAGCCCGGTGATGAAACAGCCTGATGGGGGTCGCTTTTGCTg                                            <  1:392518/90‑1 (MQ=255)
                                         ggggCGTTAAATGCCTTTGCTAATCGAGAGCGTCGTTTCGGCGGCACCGAGCCCGGTGATGAAACAGCCTGATGGGGGTCGCTTTTGCTg                                            <  2:133263/90‑1 (MQ=255)
                                            gCGTTAAATGCCTTTGCTAATCGAGAGCGTCGTTTCGGCGGCACCGAGCCCGGTGATGAAACAGCCTGATGGGGGTCGCTTTTGCTGAAg                                         >  2:13840/1‑90 (MQ=255)
                                                tAAATGCCTTTGCTAATCGAGAGCGTCGTTTCGGCGGCACCGAGCCCGGTGATGAAACAGCCTGATGGGGGTCGCTTTTGCTGAAGTATc                                     <  2:22064/90‑1 (MQ=255)
                                                tAAATGCCTTTGCTAATCGAGAGCGTCGTTTCGGCGGCACCGAGCCCGGTGATGAAACAGCCTGATGGGGGTCGCTTTTGCTGAAGTATc                                     <  2:24355/90‑1 (MQ=255)
                                                   aTGCCTTTGCTAATCGAGAGCGTCGTTTCGGCGGCACCGAGCCCGGTGATGAAACAGCCTGATGGGGGTCGCTTTTGCTGAAGTATCGcc                                  <  2:612478/90‑1 (MQ=255)
                                                         ttGCTAATCGAGAGCGTCGTTTCGGCGGCAacagaaacaacgccgcgatgacgacgggtatttacacaaaagcagatatcaggcgatact                            >  1:298293‑M1/1‑30 (MQ=255)
                                                          tGCTAATCGAGAGCGTCGTTTCGGCGGCAacagaaacaacgccgcgatgacgacgggtattaacacaaaagcagatatcaggcgatactt                           >  2:214240‑M1/1‑29 (MQ=255)
                                                          tGCTAATCGAGAGCGTCGTTTCGGCGGCACCGAGCCCGGTGATGAAACAGCCTGATGGGGGTCGCTTTTGCTGAAGTATCGCCTGATATc                           <  1:570385/90‑1 (MQ=255)
                                                          tGCTAATCGAGAGCGTCGTTTCGGCGGCACCGAGCCCGGTGATGAAACAGCCTGATGGGGGTCGCTTTTGCTGAAGTATCGCCTGATATc                           <  1:243678/90‑1 (MQ=255)
                                                           gCTAATCGAGAGCGTCGTTTCGGCGGCACCGAGCCCGGTGATGAAACAGCCTGATGGGGGTCGCTTTTGCTGAAGTATCGCCTGATATCt                          <  2:686474/90‑1 (MQ=255)
                                                                 cGAGAGCGTCGTTTCGGCGGCACCGAGCCCGGTGATGAAACAGCCTGATGGGGGTCGCTTTTGCTGAAGTATCGCCTGATATCTGCtttt                    <  1:380643/90‑1 (MQ=255)
                                                                   agagCGTCGTTTCGGCGGCACCGAGCCCGGTGATGAAACAGCCTGATGGGGGTCGCTTTTGCTGAAGTATCGCCTGATATCTGCTTTtgt                  <  2:423576/90‑1 (MQ=255)
                                                                   agagCGTCGTTTCGGCGGCACCGAGCCCGGTGATGAAACAGCCTGATGGGGGTCGCTTTTGCTGAAGTATCGCCTGATATCTGCTTTtgt                  <  2:505281/90‑1 (MQ=255)
                                                                             ttCGGCGGCAacagaaacaacgccgcgatgacgacgggtattaacacaaaagcagatatcaggcgatacttcagcaaaagcgacccccat        >  1:545130‑M1/1‑10 (MQ=255)
                                                                                 gcggcAacagaaacaacgccgcgatgacgacgggtattaacacaaaagcagatatcaggcgatacttcagcaaaagcgacccccatcagg    <  2:13359‑M1/90‑85 (MQ=255)
                                                                                   ggcAacagaaacaacgccgcgatgacgacgggtattaacacaaaagcagatatcaggcgatacttcagcaaaagcgacccccatcaggct  <  2:315592‑M1/90‑87 (MQ=255)
                                                                                      |                                                                                      
CAGATGTTCTCTGGCCCGATTTCGATGAACAAGACTTTGAAGGGGCGTTAAATGCCTTTGCTAATCGAGAGCGTCGTTTCGGCGGCACCGAGCCCGGTGATGAAACAGCCTGATGGGGGTCGCTTTTGCTGAAGTATCGCCTGATATCTGCTTTTGTGTTAATACCCGTCGTC  >  NZ_CP009273/192039‑192211

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