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.

TCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGGAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAG  >  NZ_CP009273/3574631‑3574803
                                                                                        |                                                                                    
gttatcgacaaacacggttagttcaatgcgatctacccgcaaccagttgtcgcacatttcaatcatctctcgcatcagggcgctggcgAc                                                                                     <  2:543935‑M1/2‑1 (MQ=255)
   atcgacaaacacggttagttcaatgcgatctacccgcaaccagttgtcgcacatttcaatcatctctcgcatcagggcgctggcgAcgcc                                                                                  >  2:357797‑M1/86‑90 (MQ=255)
     gATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGGAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCAcgccgc                                                                                <  2:567166/90‑1 (MQ=255)
                    aaGCAACTCGTCGCCTGTATTGATGGAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGAt                                                                 >  2:212446/1‑90 (MQ=255)
                      tcaatgcgatctacccgcaaccagttgtcgcacatttcaatcatctctcgcatcagggcgctggcgACGCCGCAGTCATGTTGCCGAttt                                                               <  1:329817‑M1/24‑1 (MQ=255)
                        aaCTCGTCGCCTGTATTGATGGAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATATTGCCGATTTTg                                                             <  2:605967/90‑1 (MQ=255)
                         aCTCGTCGCCTGTATTGATGGAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTgg                                                            >  1:138807/1‑90 (MQ=255)
                                        tGATGGAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGActc                                             <  2:138807/90‑1 (MQ=255)
                                           tGGAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCg                                          <  1:604222/90‑1 (MQ=255)
                                              agaCGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTgtgt                                                   <  2:537607/78‑1 (MQ=255)
                                              agaCGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTgtgt                                                   >  1:537607/1‑78 (MQ=255)
                                              agaCGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTg                                       >  2:488764/1‑90 (MQ=255)
                                                   tcgtGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGgaaga                                  >  1:403717/1‑90 (MQ=255)
                                                   tcgtGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGgaaga                                  >  1:345619/1‑90 (MQ=255)
                                                    cgtGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGgaagaa                                 >  2:411015/1‑90 (MQ=255)
                                                       gggCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACcg                              >  2:348120/1‑90 (MQ=255)
                                                            aTCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGTACCgcggcg                         >  1:501816/1‑90 (MQ=255)
                                                            aTCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCgcggcg                         >  1:353438/1‑90 (MQ=255)
                                                            aTCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCgcggcg                         >  2:504679/1‑90 (MQ=255)
                                                            aTCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCgcggcg                         >  1:364343/1‑90 (MQ=255)
                                                             tctcACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGt                        >  2:182587/1‑90 (MQ=255)
                                                              ctcACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTc                       <  1:201714/90‑1 (MQ=255)
                                                                  ccATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCa                   >  2:398782/1‑90 (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)
                                                                                 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)
                                                                                        |                                                                                    
TCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGGAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAG  >  NZ_CP009273/3574631‑3574803

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