New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 19727510 (0.000)32 (0.990) 19/164 0.4 100% intergenic (‑297/‑483) flhD/uspC flagellar transcriptional regulator FlhD/universal stress protein UspC
?NZ_CP009273 = 2283592 NA (NA)intergenic (‑68/‑387) BW25113_RS11435/narP IS5‑like element IS5 family transposase/nitrate/nitrite response regulator protein NarP

ACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAG  >  NZ_CP009273/1972664‑1972839
                                                                                       |                                                                                        
aCTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAgg                                                                                        >  2:830264‑M1/1‑88 (MQ=255)
aCTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAgg                                                                                        >  2:731434‑M1/1‑88 (MQ=255)
    gagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATaaa                                                                                    <  2:576706‑M2/11‑1 (MQ=255)
    gagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATaaa                                                                                    <  2:60975‑M2/11‑1 (MQ=255)
    gagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATaaa                                                                                    <  1:776037‑M2/11‑1 (MQ=255)
    gTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaagg                                                                                    <  1:830264‑M1/90‑7 (MQ=255)
        aaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATAAATGta                                                                                <  2:707160‑M2/15‑1 (MQ=255)
           cacGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaac                                                                             <  1:731434‑M1/90‑14 (MQ=255)
             cGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaa                                                                           <  2:931917‑M1/90‑16 (MQ=255)
             cGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaa                                                                           >  2:719168‑M1/1‑75 (MQ=255)
                gctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATAAATGTATAAGTCAt                                                                        <  2:662608‑M2/23‑1 (MQ=255)
                  aTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccc                                                                      >  2:630547‑M1/1‑70 (MQ=255)
                  aTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccc                                                                      <  2:318445‑M1/90‑21 (MQ=255)
                  aTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccc                                                                      <  1:630547‑M1/90‑21 (MQ=255)
                      acatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATAAATGTATAAGTCATACtttt                                                                  <  2:303402‑M2/29‑1 (MQ=255)
                       aCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgata                                                                 >  2:383065‑M1/1‑65 (MQ=255)
                         ggAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatg                                                               <  2:821071‑M1/90‑28 (MQ=255)
                          gaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATAAATGTATAAGTCATACttttgttt                                                              >  2:972059‑M2/58‑90 (MQ=255)
                           aGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgag                                                             >  2:847030‑M1/1‑61 (MQ=255)
                            agaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATAAATGTATAAGTCATACttttgttttg                                                            <  2:723092‑M2/35‑1 (MQ=255)
                                atgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATAAATGTATAAGTCATACTTTTGTTTTGGgtg                                                        <  1:511726‑M2/39‑1 (MQ=255)
                                 gCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatg                                                       >  2:436391‑M1/1‑55 (MQ=255)
                                 gCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatg                                                       <  1:642517‑M1/90‑36 (MQ=255)
                                    gCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgttt                                                    >  2:942672‑M1/1‑52 (MQ=255)
                                    gCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgttt                                                    >  2:629086‑M1/1‑52 (MQ=255)
                                                   gCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccat                                     <  1:115638‑M1/90‑54 (MQ=255)
                                                        agaatttccccgcttattcgcaccttcCTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGtttt                                <  1:928811‑M2/63‑1 (MQ=255)
                                                        aTGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaac                                <  1:256319‑M1/90‑59 (MQ=255)
                                                        aTGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaac                                <  1:918838‑M1/90‑59 (MQ=255)
                                                          gCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacag                              <  1:383065‑M1/90‑61 (MQ=255)
                                                          gCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacag                              <  1:352620‑M1/90‑61 (MQ=255)
                                                             gATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggt                           <  2:105616‑M1/90‑64 (MQ=255)
                                                             gATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggt                           <  2:673463‑M1/90‑64 (MQ=255)
                                                             gATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggt                           >  2:717108‑M1/1‑27 (MQ=255)
                                                                  cacaAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatc                      >  1:605147‑M1/1‑22 (MQ=255)
                                                                   gcttattcgcaccttcCTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgcta                     >  2:617540‑M2/17‑90 (MQ=255)
                                                                                    ttAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggtt                          <  2:582298‑M1/68‑65 (MQ=255)
                                                                                    ttAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggtt                          >  1:582298‑M1/1‑4 (MQ=255)
                                                                                     tAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttacct   <  2:123490‑M1/90‑88 (MQ=255)
                                                                                     tAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttacct   <  1:827963‑M1/90‑88 (MQ=255)
                                                                                      aaggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttacctt  >  1:725180‑M1/1‑2 (MQ=255)
                                                                                      aaggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaaca                               <  2:643903‑M1/61‑60 (MQ=255)
                                                                                      aaggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaaca                               >  1:643903‑M1/1‑2 (MQ=255)
                                                                                       |                                                                                        
ACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAG  >  NZ_CP009273/1972664‑1972839

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