New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 3920778 =24 (0.880)3 (0.130) 3/144 2.0 12.3% coding (264/993 nt) asnA aspartate‑‑ammonia ligase
?NZ_CP009273 3920817 = 22 (0.940)coding (303/993 nt) asnA aspartate‑‑ammonia ligase
Rejected: Frequency below/above cutoff threshold.

AAGCTCTGCCTGATGCCCAGTTCGAAGTGGTTCATTCACTGGCGAAGTGGAAACGTCAGACCTTAGGGCAACACGACTTCAGCGCGGGCGAAGGGCTGTACACGCACATGAAAGCCCTTCGCCCCGATGAAGACCGTCTTTCTCCGTTGCACTCGGTCTATGTTGACCAGTG  >  NZ_CP009273/3920693‑3920864
                                                                                     |                                                                                      
aaGCTCTGCCTGATGCCCAGTTCGAAGTGGTTCATTCACTGGCGAAGTGGAAACGTCAGACCTTAGGGCAACACGACTTCAGCGCGGGCg                                                                                    >  2:431471/1‑90 (MQ=255)
 tgacgctcaccgtcgcccattacgcgttcccagtcccactggtcaacatagaccgagtgcaacggagaaagacggtcttcatcgGGGCGa                                                                                   >  1:200396‑M1/85‑90 (MQ=255)
 aGCTCTGCCTGATGCCCAGTTCGAAGTGGTTCATTCACTGGCGAAGTGGAAACGTCAGACCTTAGGGCAACACGACTTCAGCGCGGGCGa                                                                                   <  1:209018/90‑1 (MQ=255)
 aGCTCTGCCTGATGCCCAGTTCGAAGTGGTTCATTCACTGGCGAAGTGGAAACGTCAGACCTTAGGGCAACACGACTTCAGCGCGGGCGa                                                                                   <  1:389469/90‑1 (MQ=255)
     ctGCCTGATGCCCAGTTCGAAGTGGTTCATTCACTGGCGAAGTGGAAACGTCAGACCTTAGGGCAACACGACTTCAGCGCGGGCGAAggg                                                                               <  1:155992/90‑1 (MQ=255)
      tGCCTGATGCCCAGTTCGAAGTGGTTCATTCACTGGCGAAGTGGAAACGTCAGACCTTAGGGCAACACGACTTCAGCGCGGGCGAAGGGc                                                                              <  1:117269/90‑1 (MQ=255)
      ctcaccgtcgcccattacgcgttcccagtcccactggtcaacatagaccgagtgcaacggagaaagacggtcttcatcgGGGCGAAGGGc                                                                              <  2:623139‑M1/11‑1 (MQ=255)
       caccgtcgcccattacgcgttcccagtcccactggtcaacatagaccgagtgcaacgggagaaagacggtcttcatcgGGGCGaa                                                                                  >  1:636766‑M1/79‑85 (MQ=255)
       caccgtcgcccattacgcgttcccagtcccactggtcaacatagaccgagtgcaacgggagaaagacggtcttcatcgGGGCGaa                                                                                  <  2:636766‑M1/7‑1 (MQ=255)
          tGATGCCCAGTTCGAAGTGGTTCATTCACTGGCGAAGTGGAAACGTCAGACCTTAGGGCAACACGACTTCAGCGCGGGCGAAGGGCTGTa                                                                          <  1:123941/90‑1 (MQ=255)
                ccAGTTCGAAGTGGTTCATTCACTGGCGAAGTGGAAACGTCAGACCTTAGGGCAACACGACTTCAGCGCGGGCGAAGGGCTGTACACGca                                                                    <  1:217322/90‑1 (MQ=255)
                ccAGTTCGAAGTGGTTCATTCACTGGCGAAGTGGAAACGTCAGACCTTAGGGCAACACGACTTCAGCGCGGGCGAAGGGCTGTACACGca                                                                    <  2:754980/90‑1 (MQ=255)
                            tcccagtcccactggtcaacatagaccgagtgcaacggagaaagacggtcttcatcgGGGCGAAGGGCTGTACACGCACATGAAAGCCCt                                                        <  1:345498‑M1/33‑1 (MQ=255)
                                     cactggtcaacatagaccgagtgcaacggagaaagacggtcttcatcgGGGCGAAGGGCTGTACACGCACATGAAAGCCCTtcgcccgcg                                               <  2:165421‑M1/42‑4 (MQ=255)
                                       tGGCGAAGTGGAAACGTCAGACCTTAGGGCAACACGACTTCAGCGCGGGCGAAGGGCTGTACACGCACATGAAAGCCCTTCGccc                                                  <  2:73303/85‑1 (MQ=255)
                                       tGGCGAAGTGGAAACGTCAGACCTTAGGGCAACACGACTTCAGCGCGGGCGAAGGGCTGTACACGCACATGAAAGCCCTTCGccc                                                  >  1:73303/1‑85 (MQ=255)
                                       tGGCGAAGTGGAAACGTCAGACCTTAGGGCAACACGACTTCAGCGCGGGCGAAGGGCTGTACACGCACATGAAAGCCCTTCGCCCCGATg                                             <  1:40005/90‑1 (MQ=255)
                                          cGAAGTGGAAACGTCAGACCTTAGGGCAACACGACTTCAGCGCGGGCGAAGGGCTGTACACGCACATGAAAGCCCTTCGCCCCGATGAAg                                          <  1:561941/90‑1 (MQ=255)
                                          cGAAGTGGAAACGTCAGACCTTAGGGCAACACGACTTCAGCGCGGGCGAAGGGCTGTACACGCACATGAAAGCCCTTCGCCCCGATGAAg                                          >  1:409492/1‑90 (MQ=255)
                                                ggAAACGTCAGACCTTAGGGCAACACGACTTCAGCGCGGGCGAAGGGCTGTACACGCACATGAAAGCCCTTCGCCCCGATGAAGACCGTc                                    >  2:661524/1‑90 (MQ=255)
                                                      gTCAGACCTTAGGGCAACACGACTTCAGCGCGGGCGAAGGGCTGTACACGCACATGAAAGCCCTTCGCCCCGATGAAGACCGTCTTtctc                              <  1:661524/90‑1 (MQ=255)
                                                       tCAGACCTTAGGGCAACACGACTTCAGCGCGGGCGAAGGGCTGTACACGCACATGAAAGCCCTTCGCCCCGATGAAGACCGTCTTTCTcc                             >  2:90199/1‑90 (MQ=255)
                                                               tAGGGCAACACGACTTCAGCGCGGGCGAAGGGCTGTACACGCACATGAAAGCCCTTCGCCCCGATGAAGACCGTCTTTCTCCGTTGCACt                     >  1:51789/1‑90 (MQ=255)
                                                               tAGGGCAACACGACTTCAGCGCGGGCGAAGGGCTGTACACGCACATGAAAGCCCTTCGCCCCGATGAAGACCGTCTTTCTCCGTTGCACt                     <  1:512642/90‑1 (MQ=255)
                                                               tAGGGCAACACGACTTCAGCGCGGGCGAAGGGCTGTACACGCACATGAAAGCCCTTCGCCCCGATGAAGACCGTCTTTCTCCGTTGCACt                     <  1:153593/90‑1 (MQ=255)
                                                                   gaaagacggtcttcatcgGGGCGAAGGGCTGTACACGCACATGAAAGCCCTTCGCCCCGATGAAGACCGTCTTTCTCCGTTGCACTCGGt                 <  1:284454‑M1/72‑1 (MQ=255)
                                                                    cAACACGACTTCAGCGCGGGCGAAGGGCTGTACACGCACATGAAAGCCCTTCGCCCCGATGAAGACCGTCTTTCTCCGTTGCACTCGGTc                <  2:324102/90‑1 (MQ=255)
                                                                    cAACACGACTTCAGCGCGGGCGAAGGGCTGTACACGCACATGAAAGCCCTTCGCCCCGATGAAGACCGTCTTTCTCCGTTGCACTCGGTc                <  2:359623/90‑1 (MQ=255)
                                                                      acacGACTTCAGCGCGGGCGAAGGGCTGTACACGCACATGAAAGCCCTTCGCCCCGATGAAGACCGTCTTTCTCCGTTGCACTCGGTCTa              >  2:127878/1‑90 (MQ=255)
                                                                            cTTCAGCGCGGGCGAAGGGCTGTACACGCACATGAAAGCCCTTCGCCCCGATGAAGACCGTCTTTCTCCGTTGCACTCGGTCTATGTTGa        <  2:107093/90‑1 (MQ=255)
                                                                                  cgcgGGCGAAGGGCTGTACACGCACATGAAAGCCCTTCGCCCCGATGAAGACCGTCTTTCTCCGTTGCACTCGGTCTATGTTGACCAGTg  <  1:431471/90‑1 (MQ=255)
                                                                                     |                                                                                      
AAGCTCTGCCTGATGCCCAGTTCGAAGTGGTTCATTCACTGGCGAAGTGGAAACGTCAGACCTTAGGGCAACACGACTTCAGCGCGGGCGAAGGGCTGTACACGCACATGAAAGCCCTTCGCCCCGATGAAGACCGTCTTTCTCCGTTGCACTCGGTCTATGTTGACCAGTG  >  NZ_CP009273/3920693‑3920864

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.