New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 19971110 (0.000)40 (0.450) 30/98 1.2 100% coding (704/753 nt) tcyN cystine ABC transporter ATP binding subunit
?NC_000913 2023679 = 0 (0.000)intergenic (+1/‑289) fliR/rcsA flagellar biosynthesis protein FliR/DNA‑binding transcriptional activator RcsA

TTTTTAATTTGCTGGCTGATATTATTAGTGAATTGCCATTAATATAATTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGAAAAATAATAC  >  NC_000913/2023632‑2023727
                                               |                                                
ttgcagcagaaacttctcgaggaactggcgggtgcgaggctgctcggTTc                                                >  1:4060274‑M2/48‑50 (MQ=255)
  gcagcagaaacttctcgaggaactggcgggtgcgaggctgctcggTTcc                                               <  1:871268‑M2/4‑1 (MQ=255)
   cagcagaaacttctcgaggaactggcgggtgcgaggctgcccggTTCCg                                              >  1:2332531‑M2/45‑49 (MQ=255)
      cagaaacttctcgaggaactggcgggcgcgaggccgctcggTTCCGGAAc                                          >  1:6318070‑M2/42‑50 (MQ=255)
       agaaacttctcgaggaactggcgggtgcgaggctgctcggTTCCGTAACg                                         <  1:6346775‑M2/10‑1 (MQ=255)
       agaaacttctcgaggaactggcgggtgcgaggctgctcggTTCCGTAACg                                         >  1:2480136‑M2/41‑50 (MQ=255)
       agaaacttctcgaggaactggcgggtgcgaggctgctcggTTCCGTAACg                                         <  1:3417350‑M2/10‑1 (MQ=255)
        gaaacttctcgaggaactggcgggtgcgaggctgctcggTTCCGTAACGt                                        <  1:6493128‑M2/11‑1 (MQ=255)
          aacttctcgaggaactggcgggtgcgaggctgctcggTTCCGTAACGttt                                      <  1:2814480‑M2/13‑1 (MQ=255)
          aacttctcgaggaactggcgggggcgcggctgctcggTTCCCTAACGttt                                      >  1:3376462‑M2/38‑50 (MQ=255)
                  gaggaactggcgggtgcgaggctgctcggTTCCGTAACGTTTATCATg                                <  1:2228302‑M2/19‑1 (MQ=255)
                  gaggaactggcgggtgcgaggctgctcggTTCCGTAACGTTTATCATg                                <  1:5897722‑M2/19‑1 (MQ=255)
                  gaggaactggcgggtgcgaggctgctcggTTCCGTAACGTTTATCATGtt                              <  1:4321187‑M2/21‑1 (MQ=255)
                  gaggaactggcgggtgcgaggctgctcggTTCCGTAACGTTTATCATGtt                              <  1:3531707‑M2/21‑1 (MQ=255)
                    ggaactggcgggtgcgaggctgctcggTTCCGTAACGTTTATCATGTTAt                            <  1:3440703‑M2/23‑1 (MQ=255)
                      aactggcgggtgcgaggctgctcggTTCCGTAACGTTTATCATGTTATcc                          >  1:1832077‑M2/26‑50 (MQ=255)
                        ctggcgggtgcgaggctgctcggTTCCGTAACGTTTATCATGTTATCCTa                        <  1:321619‑M2/27‑1 (MQ=255)
                             gggtgcgaggctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGAt                   >  1:1873997‑M2/19‑50 (MQ=255)
                              ggtgcgaggctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGAtt                  <  1:622103‑M2/33‑1 (MQ=255)
                              ggtgcgaggctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGAtt                  <  1:5048318‑M2/33‑1 (MQ=255)
                               gtgcgaggctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGATTa                 <  1:5319186‑M2/34‑1 (MQ=255)
                                 gcgaggctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATc               >  1:1149002‑M2/15‑50 (MQ=255)
                                  cgaggctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATcc              <  1:265621‑M2/37‑1 (MQ=255)
                                  cgaggctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATc               <  1:6191177‑M2/36‑1 (MQ=255)
                                    aggctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGa            <  1:5301358‑M2/39‑1 (MQ=255)
                                     ggctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGa            >  1:3055859‑M2/11‑49 (MQ=255)
                                     ggctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGa            >  1:6389009‑M2/11‑49 (MQ=255)
                                        tgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGaaaaa        >  1:1875434‑M2/8‑50 (MQ=255)
                                         gctcggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGaaaaa        >  1:4994281‑M2/7‑49 (MQ=255)
                                         gctcggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGaaaaa        >  1:1292002‑M2/7‑49 (MQ=255)
                                         gctcggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGaaaaa        >  1:6628212‑M2/7‑49 (MQ=255)
                                          ctcggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGaaaaa        <  1:455803‑M2/43‑1 (MQ=255)
                                           tcggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGAAAaataa     >  1:4465853‑M2/5‑50 (MQ=255)
                                            cggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGAAAaataa     >  1:3597082‑M2/4‑49 (MQ=255)
                                             ggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGAAAAataata   >  1:5761685‑M2/3‑50 (MQ=255)
                                              gTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGAAAAATAATAc  >  1:5686599‑M2/2‑50 (MQ=255)
                                              gTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGAAAAATAATAc  <  1:1368661‑M2/49‑1 (MQ=255)
                                              gTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGAAAAATAATAc  >  1:791649‑M2/2‑50 (MQ=255)
                                              gTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGAAAAATAATAc  <  1:272940‑M2/49‑1 (MQ=255)
                                              gTTCCGTAACGTTTATCATCTTATCCTAAGGATTATCCGAAAAATAATAc  <  1:2293883‑M2/49‑1 (MQ=255)
                                               |                                                
TTTTTAATTTGCTGGCTGATATTATTAGTGAATTGCCATTAATATAATTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGAAAAATAATAC  >  NC_000913/2023632‑2023727

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