New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 223597 =NA (NA)13 (0.090) 6/168 NT NA intergenic (+189/‑174) gmhB/rrsH D,D‑heptose 1,7‑bisphosphate phosphatase/16S ribosomal RNA of rrnH operon
?NC_000913 223624 = NA (NA)intergenic (+216/‑147) gmhB/rrsH D,D‑heptose 1,7‑bisphosphate phosphatase/16S ribosomal RNA of rrnH operon

CTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGCTCTTTAACAATTTATCAGA  >  NC_000913/223508‑223648
                                                                                         |                                                   
cTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAg                                           >  2:3958889/1‑100 (MQ=31)
 ttCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGc                                          <  1:1567987/100‑1 (MQ=32)
 ttCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGc                                          <  1:2322767/100‑1 (MQ=32)
 ttCACACAAACAGCCGGGTAAAAAAAAAAAAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGa                                            <  1:677860/98‑1 (MQ=11)
  tCACACAATCAGCCGGTTAAAAAAAAAAAAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGc                                          <  1:1985260/99‑1 (MQ=11)
   cacacaAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTg                                        <  2:3701119/100‑1 (MQ=33)
   cacacaAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTg                                        >  2:3513881/1‑100 (MQ=33)
   cacacaAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTg                                        <  1:3266057/100‑1 (MQ=33)
   cacacaAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTg                                        <  1:2338659/100‑1 (MQ=33)
     cacaAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGaa                                      >  1:1120800/1‑100 (MQ=38)
     cacaAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGaa                                      >  2:2612799/1‑100 (MQ=38)
     cacaAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGaa                                      >  2:2604040/1‑100 (MQ=38)
     cacaAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGaa                                      >  2:2603385/1‑100 (MQ=38)
     cacaAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGaa                                      >  1:761083/1‑100 (MQ=38)
       caAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAg                                    <  2:246690/100‑1 (MQ=40)
       caAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAg                                    >  1:502902/1‑100 (MQ=38)
         aaCAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGcc                                  >  2:3720260/1‑100 (MQ=255)
            aGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCcgcg                               <  1:1305249/100‑1 (MQ=255)
            aGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCcgcg                               <  1:1942349/100‑1 (MQ=255)
            aGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCcgcg                               <  1:3479688/100‑1 (MQ=255)
            aGCCGGCGCGATTGAAGATAAAAAACCTGAAAGTCAGGGGTGACGCTGAAAGAGAAAAGCTGAATATACGCCACCTCGCGACAGTAATCTTAAAGCcg                                 >  1:353416/1‑98 (MQ=2)
             gCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGt                              >  1:1475856/1‑100 (MQ=255)
             gCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGt                              >  2:216511/1‑100 (MQ=255)
               cggttcggttGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCg                            <  2:3769689/100‑1 (MQ=255)
               cggttcggttGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACATCGCGACAGTGAGCTGAAAGCCGCGTCg                            >  2:2337979/1‑100 (MQ=255)
                ggttcggttGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATATGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGc                           >  2:185602/1‑100 (MQ=255)
                ggttcggttGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGc                           <  2:4131976/100‑1 (MQ=255)
                 gttcggttGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATACACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCa                          <  1:3643548/100‑1 (MQ=255)
                  ttcggttGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCaa                         >  2:207088/1‑100 (MQ=255)
                   tcggttGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAAc                        <  1:2530162/100‑1 (MQ=255)
                   tcggttGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAAc                        <  2:2321959/100‑1 (MQ=255)
                   tcggttGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAAc                        <  2:2320921/100‑1 (MQ=255)
                     ggttGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGAGTCGCAACTg                      >  1:3416420/1‑100 (MQ=255)
                      gttGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGc                     <  1:2082861/100‑1 (MQ=40)
                      gttGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGc                     >  2:2963638/1‑100 (MQ=40)
                      gttGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGc                     >  1:198054/1‑100 (MQ=40)
                      gttGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGc                     <  2:619528/100‑1 (MQ=40)
                       ttGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGct                    >  2:2230558/1‑100 (MQ=38)
                       ttGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGct                    >  2:2196465/1‑100 (MQ=38)
                        tGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGctc                   <  1:3022487/100‑1 (MQ=35)
                         gAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGctct                  <  1:1860558/100‑1 (MQ=35)
                           agagAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGCTCttt                >  2:758747/1‑100 (MQ=32)
                           agagAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGCTCttt                <  1:2213993/100‑1 (MQ=33)
                           agagAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGCTCttt                >  1:1068388/1‑100 (MQ=32)
                           agagAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGCTCttt                >  1:2430504/1‑100 (MQ=32)
                            gagaAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGCTCttt                >  2:3837951/1‑99 (MQ=32)
                            gagaAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGCTCTTTa               <  2:2587412/100‑1 (MQ=32)
                            gagaAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGCTCTTTa               <  2:1831743/100‑1 (MQ=32)
                            gagaAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGCTCTTTa               <  1:191253/100‑1 (MQ=32)
                             agaAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGCTCTTTaa              >  1:1334319/1‑100 (MQ=31)
                              gaAAAATCCTGTAATGCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGCTCTTTAAc             <  1:1881203/100‑1 (MQ=2)
                              gaAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGCTCTTTAAc             <  2:2220350/100‑1 (MQ=11)
                              gaAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGCTCTTTAAc             <  1:3973607/100‑1 (MQ=11)
                                          aagtCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGCTCTTTAACAATTTATCAGa  <  1:3329750/96‑1 (MQ=1)
                                                                                         |                                                   
CTTCACACAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGAGCTGAAAGCCGCGTCGCAACTGCTCTTTAACAATTTATCAGA  >  NC_000913/223508‑223648

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