New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 1046001NA (NA)18 (0.710) 14/168 0.5 54.5% intergenic (+15/+418) BW25113_RS05135/cspH IS1‑like element IS1D family transposase/cold shock‑like protein CspH
?NZ_CP009273 = 1999990 15 (0.590)coding (354/1488 nt) amyA alpha‑amylase

GGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTTGTGAAATCTCCCTTTCCTCATACCGCTTG  >  NZ_CP009273/1045929‑1046082
                                                                        |                                                                                 
ggTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTg                                                                  <  1:199039/90‑1 (MQ=255)
ggTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTg                                                                  <  1:281840/90‑1 (MQ=255)
 gTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTga                                                                 >  1:184754/1‑90 (MQ=255)
         gCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGa                                                         <  1:33669/90‑1 (MQ=255)
              aCAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTc                                                    <  2:296192/90‑1 (MQ=255)
                   gTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACtt                                               >  2:74027/1‑90 (MQ=255)
                     cATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTg                                             <  1:74027/90‑1 (MQ=255)
                      aTCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGa                                            >  1:96610/1‑90 (MQ=255)
                      aTCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGa                                            <  2:534947/90‑1 (MQ=255)
                               tATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACAt                                   <  1:688995/90‑1 (MQ=255)
                                    gAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTtg                              >  1:528830/1‑90 (MQ=255)
                                     aaCATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTtgt                             <  1:454701/90‑1 (MQ=255)
                                                cTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTTGTGAAATCTCCCt                  <  2:473750/90‑1 (MQ=255)
                                                cTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTTGTGAAATCTCCCt                  <  2:528830/90‑1 (MQ=255)
                                                 tATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTTGTGAAATCTCCCtt                 >  2:254730/1‑90 (MQ=255)
                                                         aGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTTGTGAAATCTCCCTTTCCTCATa         <  2:150287/90‑1 (MQ=255)
                                                                gTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTTGTGAAATCTCCCTTTCCTCATACCGCTTg  <  2:544927/90‑1 (MQ=255)
                                                                gTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTTGTGAAATCTCCCTTTCCTCATACCGCTTg  <  2:711335/90‑1 (MQ=255)
                                                                    ttACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGtat                                        >  2:193522/1‑48 (MQ=255)
                                                                    ttACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGtat                                        <  1:193522/48‑1 (MQ=255)
                                                                        |                                                                                 
GGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAATACGTACGCTTGAGATCAGGATTGTCTACTTTGAGTATAACATCCTTGTGAAATCTCCCTTTCCTCATACCGCTTG  >  NZ_CP009273/1045929‑1046082

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