New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 1972747 =0 (0.000)9 (0.450) 5/156 1.3 58.1% intergenic (‑293/‑487) flhD/uspC flagellar transcriptional regulator FlhD/universal stress protein UspC
?NZ_CP009273 2095231 = 13 (0.650)intergenic (‑1/+145) wbbL/BW25113_RS10605 N‑acetylglucosaminyl‑diphospho‑decaprenol L‑rhamnosyltransferase WbbL/IS5‑like element IS5 family transposase

AACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACG  >  NZ_CP009273/1972658‑1972817
                                                                                         |                                                                      
gcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcc                                                                        >  1:122912‑M1/90‑90 (MQ=255)
gcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcc                                                                        >  2:129267‑M1/90‑90 (MQ=255)
gcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcc                                                                        >  2:361172‑M1/90‑90 (MQ=255)
gcatttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcc                                                                        >  1:49005‑M1/90‑90 (MQ=255)
 catttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcct                                                                       >  1:250279‑M1/89‑90 (MQ=255)
  atttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcctt                                                                      <  1:166356‑M1/3‑1 (MQ=255)
   tttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTa                                                                     <  2:467102‑M1/4‑1 (MQ=255)
   tttacgggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTa                                                                     <  2:498786‑M1/4‑1 (MQ=255)
        gggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATa                                                                >  2:191055‑M1/82‑90 (MQ=255)
        gggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATa                                                                >  1:625484‑M1/82‑90 (MQ=255)
        gggagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATa                                                                >  1:23057‑M1/82‑90 (MQ=255)
         cGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaag                                                               >  2:593040‑M1/1‑85 (MQ=255)
          gagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATaaa                                                              >  2:526056‑M1/80‑90 (MQ=255)
          gagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATaaa                                                              >  1:503612‑M1/80‑90 (MQ=255)
          gagaaaaaatcggctcaaacatgaagaaatgaaatgactgagtaagccgagaagaatttccccgcttattcgcaccttcCTTAACATaaa                                                              >  1:148178‑M1/80‑90 (MQ=255)
             aaaaaatcggctcaaacatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATAAATGt                                                           <  2:524221‑M1/14‑1 (MQ=255)
              cTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcg                                                          >  2:262519‑M1/1‑80 (MQ=255)
                   cGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaa                                                     >  1:205978‑M1/1‑75 (MQ=255)
                     cacaTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagt                                                   >  1:28483‑M1/1‑73 (MQ=255)
                        aTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccc                                                <  1:543674‑M1/90‑21 (MQ=255)
                          aCAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctg                                              <  1:382262‑M1/90‑23 (MQ=255)
                            acatgaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttcCTTAACATAAATGTATAAGTCATACtttt                                            >  2:6210‑M1/62‑90 (MQ=255)
                                      ggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcat                                  <  2:425666‑M1/90‑35 (MQ=255)
                                                 tCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatct                       <  1:450216‑M1/90‑46 (MQ=255)
                                                    ataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctgga                    <  1:262519‑M1/90‑49 (MQ=255)
                                                    ataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcaaacaagtccctgatatgagatcatgtttgtcatctgga                    <  1:291696‑M1/90‑49 (MQ=255)
                                                      aatGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagc                  <  1:593040‑M1/90‑51 (MQ=255)
                                                      aatGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagc                  <  2:560902‑M1/90‑51 (MQ=255)
                                                         gCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccat               >  1:134017‑M1/1‑37 (MQ=255)
                                                         gCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccat               <  1:110454‑M1/90‑54 (MQ=255)
                                                                gCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacag        <  1:86005‑M1/90‑61 (MQ=255)
                                                                   gATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggt     <  2:538362‑M1/90‑64 (MQ=255)
                                                                   gATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggt     <  2:616570‑M1/90‑64 (MQ=255)
                                                                      cccgcttattcgcaccttcCTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacg  <  1:577028‑M1/71‑1 (MQ=255)
                                                                                         |                                                                      
AACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACG  >  NZ_CP009273/1972658‑1972817

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