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.

TGGCGAAGTGGAAACGTCAGACCTTAGGGCAACACGACTTCAGCGCGGGCGAAGGGCTGTACACGCACATGAAAGCCCTTCGCCCCGATGAAGACCGTCTTTCTCCGTTGCACTCGGTCTATGTTGACCAGTGGGACTGGGAACGCGTAATGGGCGACGGTGAGCGTCAATTCTC  >  NZ_CP009273/3920732‑3920906
                                                                                     |                                                                                          
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)
             accgagtgcaacggagaaagacggtcttcatcggggcgaagggctttcatgtgcgtgtacagcccttcgcccCGATGAAGACCGTCTTtc                                                                           <  1:284454‑M2/18‑1 (MQ=255)
               gTCAGACCTTAGGGCAACACGACTTCAGCGCGGGCGAAGGGCTGTACACGCACATGAAAGCCCTTCGCCCCGATGAAGACCGTCTTtctc                                                                        <  1:661524/90‑1 (MQ=255)
                tCAGACCTTAGGGCAACACGACTTCAGCGCGGGCGAAGGGCTGTACACGCACATGAAAGCCCTTCGCCCCGATGAAGACCGTCTTTCTcc                                                                       >  2:90199/1‑90 (MQ=255)
                        tAGGGCAACACGACTTCAGCGCGGGCGAAGGGCTGTACACGCACATGAAAGCCCTTCGCCCCGATGAAGACCGTCTTTCTCCGTTGCACt                                                               <  1:153593/90‑1 (MQ=255)
                        tAGGGCAACACGACTTCAGCGCGGGCGAAGGGCTGTACACGCACATGAAAGCCCTTCGCCCCGATGAAGACCGTCTTTCTCCGTTGCACt                                                               >  1:51789/1‑90 (MQ=255)
                        tAGGGCAACACGACTTCAGCGCGGGCGAAGGGCTGTACACGCACATGAAAGCCCTTCGCCCCGATGAAGACCGTCTTTCTCCGTTGCACt                                                               <  1:512642/90‑1 (MQ=255)
                            gaaagacggtcttcatcgGGGCGAAGGGCTGTACACGCACATGAAAGCCCTTCGCCCCGATGAAGACCGTCTTTCTCCGTTGCACTCGGt                                                           <  1:284454‑M1/72‑1 (MQ=255)
                             cAACACGACTTCAGCGCGGGCGAAGGGCTGTACACGCACATGAAAGCCCTTCGCCCCGATGAAGACCGTCTTTCTCCGTTGCACTCGGTc                                                          <  2:359623/90‑1 (MQ=255)
                             cAACACGACTTCAGCGCGGGCGAAGGGCTGTACACGCACATGAAAGCCCTTCGCCCCGATGAAGACCGTCTTTCTCCGTTGCACTCGGTc                                                          <  2:324102/90‑1 (MQ=255)
                               acacGACTTCAGCGCGGGCGAAGGGCTGTACACGCACATGAAAGCCCTTCGCCCCGATGAAGACCGTCTTTCTCCGTTGCACTCGGTCTa                                                        >  2:127878/1‑90 (MQ=255)
                                     cTTCAGCGCGGGCGAAGGGCTGTACACGCACATGAAAGCCCTTCGCCCCGATGAAGACCGTCTTTCTCCGTTGCACTCGGTCTATGTTGa                                                  <  2:107093/90‑1 (MQ=255)
                                           cgcgggcgaagggctttcatgtgcgtgtacagcccttcgcccCGATGAAGACCGTCTTTCTCCGTTGCACTCGGTCTATGTTGACCAGTg                                            <  2:165421‑M2/48‑1 (MQ=255)
                                           cgcgGGCGAAGGGCTGTACACGCACATGAAAGCCCTTCGCCCCGATGAAGACCGTCTTTCTCCGTTGCACTCGGTCTATGTTGACCAGTg                                            <  1:431471/90‑1 (MQ=255)
                                                GCGAAGGGCTGTACACGCACATGAAAGCCCTTCGCCCCGATGAAGACCGTCTTTCTCCGTTGCACTCGGTc                                                          <  1:584807‑M1/71‑1 (MQ=255)
                                                GCGAAGGGCTGTACACGCACATGAAAGCCCTTCGCCCCGATGAAGACCGTCTTTCTCCGTTGCACTCGGTc                                                          <  1:295658‑M1/71‑1 (MQ=255)
                                                GCGAAGGGCTGTACACGCACATGAAAGCCCTTCGCCCCGATGAAGACCGTCTTTCTCCGTTGCACTCGGTc                                                          >  2:295658‑M1/1‑71 (MQ=255)
                                                GCGAAGGGCTGTACACGCACATGAAAGCCCTTCGCCCCGATGAAGACCGTCTTTCTCCGTTGCACTCGGTc                                                          >  2:584807‑M1/1‑71 (MQ=255)
                                                    agggctttcatgtgcgtgtacagcccttcgcccCGATGAAGACCGTCTTTCTCCGTTGCACTCGGTCTATGTTGACCAGTGGGACTGGGa                                   <  1:345498‑M2/57‑1 (MQ=255)
                                                    AGGGCTGTACACGCACATGAAAGCCCTTCGCCCCGATGAAGACCGTCTTTCTCCGTTGCACTCGGTCTATGTTGACCAGTGGGACTGGGa                                   >  1:310903‑M1/1‑90 (MQ=255)
                                                    AGGGCTGTACACGCACATGAAAGCCCTTCGCCCCGATGAAGACCGTCTTTCTCCGTTGCACTCGGTCTATGTTGACCAGTGGGACTGGGa                                   >  2:255561‑M1/1‑90 (MQ=255)
                                                         TGTACACGCACATGAAAGCCCTTCGCCCCGATGAAGACCGTCTTTCTCCGTTGCACTCGGTCTATGTTGACCAGTGGGACTGGGAAcgcg                              <  2:756382‑M1/90‑1 (MQ=255)
                                                         TGTACACGCACATGAAAGCCCTTCGCCCCGATGAAGACCGTCTTTCTCCGTTGCACTCGGTCTATGTTGACCAGTGGGACTGGGAAcgcg                              <  2:719421‑M1/90‑1 (MQ=255)
                                                         TGTACACGCACATGAAAGCCCTTCGCCCCGATGAAGACCGTCTTTCTCCGTTGCACTCGGTCTATGTTGACCAGTGGGACTGGGAAcgcg                              <  2:53360‑M1/90‑1 (MQ=255)
                                                          gTACACGCACATGAAAGCCCTTCGCCCCGATGAAGACCGTCTTTCTCCGTTGCACTCGGTCTATGTTGACCAGTGGGACTGGGAACGCGt                             >  2:664995/1‑90 (MQ=255)
                                                            acacGCACATGAAAGCCCTTCGCCCCGATGAAGACCGTCTTTCTCCGTTGCACTCGGTCTATGTTGACCAGTGGGACTGGGAACGCGTaa                           <  2:41042/90‑1 (MQ=255)
                                                               cGCACATGAAAGCCCTTCGCCCCGATGAAGACCGTCTTTCTCCGTTGCACTCGGTCTATGTTGACCAGTGGGACTGGGAACGCGTAATgg                        >  1:685460/1‑90 (MQ=255)
                                                                gCACATGAAAGCCCTTCGCCCCGATGAAGACCGTCTTTCTCCGTTGCACTCGTCTATGTTGACCAGTGGGACTGGGAACGCGTAATGGGc                      >  2:409766/1‑90 (MQ=255)
                                                                 cacaTGAAAGCCCTTCGCCCCGATGAAGACCGTCTTTCTCCGTTGCACTCGGTCTATGTTGACCAGTGGGACTGGGAACGCGTAATGGGc                      >  2:536095/1‑90 (MQ=255)
                                                                  acaTGAAAGCCCTTCGCCCCGATGAAGACCGTCTTTCTCCGTTGCACTCGGTCTATGTTGACCAGTGGGACTGGGAACGCGTAATGGGCg                     <  1:756382/90‑1 (MQ=255)
                                                                    aTGAAAGCCCTTCGCCCCGATGAAGACCGTCTTTCTCCGTTGCACTCGGTCTATGTTGACCAGTGGGACTGGGAACGCGTAATGGGCGAc                   >  2:622220/1‑90 (MQ=255)
                                                                     tGAAAGCCCTTCGCCCCGATGAAGACCGTCTTTCTCCGTTGCACTCGGTCTATGTTGACCAGTGGGACTGGGAACGCGTAATGGGCGACg                  <  1:127878/90‑1 (MQ=255)
                                                                          gcccttcgcccCGATGAAGACCGTCTTTCTCCGTTGCACTCGGTCTATGTTGACCAGTGGGACTGGGAACGCGTAATGGGCGACGGTGAg             <  2:623139‑M2/79‑1 (MQ=255)
                                                                              ttcgcccCGATGAAGACCGTCTTTCTCCCGTTGCACTCGGTCTATGTTGACCAGTGGGACTGGGAACGCGTAATGGGCGACGGTg               <  2:636766‑M2/78‑1 (MQ=255)
                                                                              ttcgcccCGATGAAGACCGTCTTTCTCCCGTTGCACTCGGTCTATGTTGACCAGTGGGACTGGGAACGCGTAATGGGCGACGGTg               >  1:636766‑M2/8‑85 (MQ=255)
                                                                               tcgcccCGATGAAGACCGTCTTTCTCCGTTGCACTCGGTCTATGTTGACCAGTGGGACTGGGAACGCGTAATGGGCGACGGTGAGCGTCa        >  1:200396‑M2/7‑90 (MQ=255)
                                                                                     cGATGAAGACCGTCTTTCTCCGTTGCACTCGGTCTATGTTGACCAGTGGGACTGGGAACGCGTAATGGGCGACGGTGAGCGTCAATtctc  <  1:31386/90‑1 (MQ=255)
                                                                                     |                                                                                          
TGGCGAAGTGGAAACGTCAGACCTTAGGGCAACACGACTTCAGCGCGGGCGAAGGGCTGTACACGCACATGAAAGCCCTTCGCCCCGATGAAGACCGTCTTTCTCCGTTGCACTCGGTCTATGTTGACCAGTGGGACTGGGAACGCGTAATGGGCGACGGTGAGCGTCAATTCTC  >  NZ_CP009273/3920732‑3920906

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