New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 1238033 =36 (0.650)3 (0.070) 3/130 3.6 8.5% coding (504/915 nt) ldcA muramoyltetrapeptide carboxypeptidase
?NZ_CP009273 1238051 = 36 (0.830)coding (486/915 nt) ldcA muramoyltetrapeptide carboxypeptidase
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TTGGCATCCACGGTGTACCGATCAGTGAAATCAGCATCGCAAGATTGCCTCCCCACAACGTGCCTTCGGCCCGACATGTCGGGCCTTCCCCTTGCCATTCAATGGTGAAGGTTTCATTGCGTAACGCCAGCCAGAAATGGTGCTCGGTAA  >  NZ_CP009273/1237966‑1238115
                                                                                     |                                                                
ttGGCATCCACGGTGTACCGATCAGTGAAATCAGCATCGCAAGATTGCCTCCCCACAACGTGCCTTTGGCCCGACATGTCGGGCCTTccc                                                              <  1:710765/90‑1 (MQ=255)
 tGGCATCCACGGTGTACCGATCAGTGAAATCAGCATCGCAAGATTGCCTCCCCACAACGTGCCTTCGGCCCGACATGTCGGGCCTTcccc                                                             >  2:1546279/1‑90 (MQ=255)
    cATCCACGGTGTACCGATCAGTGAAATCAGCATCGCAAGATTGCCTCCCCACAACGTGCCTTCGGCCCGACATGTCGGGCCTTCCCCTTg                                                          >  2:948139/1‑90 (MQ=255)
     aTCCACGGTGTACCGATCAGTGAAATCAGCATCGCAAGATTGCCTCCCCACAACGTGCCTTCGGCCCGACATGTCGGGCCTTCCCCTTGc                                                         <  1:1262945/90‑1 (MQ=255)
     aTCCACGGTGTACCGATCAGTGAAATCAGCATCGCAAGATTGCCTCCCCACAACGTGCCTTCGGCCCGACATGTCGGGCCTTCCCCTTGc                                                         <  1:146102/90‑1 (MQ=255)
           ggTGTACCGATCAGTGAAATCAGCATCGCAAGATTGCCTCCCCACAACGTGCCTTCGGCCCGACATGTCGGGCCTTCCCCTTGCCATTCa                                                   <  1:1580466/90‑1 (MQ=255)
            gtgtACCGATCAGTGAAATCAGCATCGCAAGATTGCCTCCCCACAACGTGCCTTCGGCCCGACATGTCGGGCCTTCCCCTTGCCATTCaa                                                  >  2:436855/1‑90 (MQ=255)
                aCCGATCAGTGAAATCAGCATCGCAAGATTGCCTCCCCACAACGTGCCTTCGGCCCGACATGTCGGGCCTTCCCCTTGCCATTCAATggg                                              >  2:581120/1‑89 (MQ=255)
                   gATCAGTGAAATCAGCATCGCAAGATTGCCTCCCCACAACGTGCCTTCGGCCCGACATGTCGGGCCTTCCCCTTGCCATTCAATGGTGaa                                           >  1:1322422/1‑90 (MQ=255)
                   gATCAGTGAAATCAGCATCGCAAGATTGCCTCCCCACAACGTGCCTTCGGCCCGACATGTCGGGCCTTCCCCTTGCCATTCAATGGTGaa                                           >  2:1237722/1‑90 (MQ=255)
                   gATCAGTGAAATCAGCATCGCAAGATTGCCTCCCCACAACGTGCCTTCGGCCCGACATGTCGGGCCTTCCCCTTGCCATTCAATGGTGaa                                           <  1:581120/90‑1 (MQ=255)
                      cAGTGAAATCAGCATCGCAAGATTGCCTCCCCACAACGTGCCTTCGGCCCGACATGTCGGGCCTTCCCCTTGCCATTCAATGGTg                                             <  1:788557/85‑1 (MQ=255)
                      cAGTGAAATCAGCATCGCAAGATTGCCTCCCCACAACGTGCCTTCGGCCCGACATGTCGGGCCTTCCCCTTGCCATTCAATGGTg                                             >  2:788557/1‑85 (MQ=255)
                      cAGTGAAATCAGCATCGCAAGATTGCCTCCCCACAACGTGCCTTCGGCCCGACATGTCGGGCCTTCCCCTTGCCATTCAATGGTGAAGGt                                        >  1:76435/1‑90 (MQ=255)
                      cAGTGAAATCAGCATCGCAAGATTGCCTCCCCACAACGTGCCTTCGGCCCGACATGTCGGGCCTTCCCCTTGCCATTCAATGGTGAAGGt                                        >  1:1436645/1‑90 (MQ=255)
                      cAGTGAAATCAGCATCGCAAGATTGCCTCCCCACAACGTGCCTTCGGCCCGACATGTCGGGCCTTCCCCTTGCCATTCAATGGTGAAGGt                                        <  2:939043/90‑1 (MQ=255)
                      cAGTGAAATCAGCATCGCAAGATTGCCTCCCCACAACGTGCCTTCGGCCCGACATGTCGGGCCTTCCCCTTGCCATTCAATGGTGAAGGt                                        >  1:466362/1‑90 (MQ=255)
                      cAGTGAAATCAGCATCGCAAGATTGCCTCCCCACAACGTGCCTTCGGCCCGACATGTCGGGCCTTCCCCTTGCCATTCAATGGTGAAGGt                                        >  2:1488341/1‑90 (MQ=255)
                      cAGTGAAATCAGCATCGCAAGATTGCCTCCCCACAACGTGCCTTCGGCCCGACATGTCGGGCCTTCCCCTTGCCATTCAATGGTGAAGGt                                        >  2:1416548/1‑90 (MQ=255)
                        gTGAAATCAGCATCGCAAGATTGCCTCCCCACAACGTGCCTTCGGCCCGACATGTCGGGCCTTCCCCTTGCCATTCAATGGTGAAGGttt                                      >  2:303354/1‑90 (MQ=255)
                           aaaTCAGCATCGCAAGATTGCCTCCCCACAACGTGCCTTCGGCCCGACATGTCGGGCCTTCCCCTTGCCATTCAATGGTGAAGGTTTCAt                                   <  2:955768/90‑1 (MQ=255)
                           aaaTCAGCATCGCAAGATTGCCTCCCCACAACGTGCCTTCGGCCCGACATGTCGGGCCTTCCCCTTGCCATTCAATGGTGAAGGTTTCAt                                   <  2:76435/90‑1 (MQ=255)
                                    tCGCAAGATTGCCTCCCCACAACGTGCCTTCGGCCCGACATGTCGGGCCTTCCCCTTGCCATTc                                                    >  1:781254/1‑64 (MQ=255)
                                    tCGCAAGATTGCCTCCCCACAACGTGCCTTCGGCCCGACATGTCGGGCCTTCCCCTTGCCATTc                                                    <  2:781254/64‑1 (MQ=255)
                                     cGCAAGATTGCCTCCCCACAACGTGCCTTCGGCCCGACATGTCGGGCCTTCCCCTTGCCATTCAATGGTGAAGGTTTCATTGCGTAACgc                         <  2:63661/90‑1 (MQ=255)
                                     cGCAAGATTGCCTCCCCACAACGTGCCTTCGGCCCGACATGTCGGGCCTTCCCCTTGCCATTCAATGGTGAAGGTTTCATTGCGTAACgc                         <  1:948139/90‑1 (MQ=255)
                                              gCCTCCCCACAACGTGCCTACGGCCCGACATGTCGGGCCTTCCCCTTGCCATTCAATGGTGAAGGTTTCATTGCGTAACGCCAGCCAGaa                >  2:529700/1‑90 (MQ=255)
                                                   cccACAACGTGCCTTCGGCCCGACATGTCGGGCCTTCCCCTTGCCATTCAATGGTGAAGGTTTCATTGCGTAACGCCAGCCAGAAATGGt           >  1:44172/1‑90 (MQ=255)
                                                     cacaACGTGCCTTCGGCCCGACATGTCGGGCCTTCCCCTTGCCATTCAATGGTGAAGGTTTCATTGCGTAACGCCAGCCAGAAATGGTGc         <  2:12107/90‑1 (MQ=255)
                                                           gTGCCTTCGGCCCGACATGTCGGGCCTTCCCCTTGc                                                         >  2:593792/1‑36 (MQ=255)
                                                           gTGCCTTCGGCCCGACATGTCGGGCCTTCCCCTTGc                                                         <  1:593792/36‑1 (MQ=255)
                                                           gTGCCTTCGGCCCGACATGTCGGGCCTTCCCCTTGCCATTCAATGGTGAAGGTTTCATTGCGTAACGCCAGCCAGAAATGGTGCTCGGTa   >  2:1015768/1‑90 (MQ=255)
                                                            tGCCTTCGGCCCGACATGTCGGGCCTTCCCCTTGCCATTCAATGGTg                                             <  1:806585/47‑1 (MQ=255)
                                                            tGCCTTCGGCCCGACATGTCGGGCCTTCCCCTTGCCATTCAATGGTg                                             >  2:806585/1‑47 (MQ=255)
                                                            tGCCTTCGGCCCGACATGTCGGGCCTTCCCCTTGCCATTCAATGGTGAAGGTTTCATTGCGTAACGCCAGCCAGAAATGGTGCTCGGTaa  >  2:702921/1‑90 (MQ=255)
                                                            tGCCTTCGGCCCGACATGTCGGGCCTTCCCCTTGCCATTCAATGGTGAAGGTTTCATTGCGTAACGCCAGCCAGAAATGGTGCTCGGTaa  <  1:529700/90‑1 (MQ=255)
                                                                                     |                                                                
TTGGCATCCACGGTGTACCGATCAGTGAAATCAGCATCGCAAGATTGCCTCCCCACAACGTGCCTTCGGCCCGACATGTCGGGCCTTCCCCTTGCCATTCAATGGTGAAGGTTTCATTGCGTAACGCCAGCCAGAAATGGTGCTCGGTAA  >  NZ_CP009273/1237966‑1238115

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.