New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 2283591NA (NA)13 (0.810) 10/166 0.4 100% intergenic (‑67/‑388) BW25113_RS11435/narP IS5‑like element IS5 family transposase/nitrate/nitrite response regulator protein NarP
?NZ_CP009273 = 3512441 0 (0.000)intergenic (‑18/+383) aroK/hofQ shikimate kinase AroK/DNA uptake porin HofQ

CGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGATAAACCCGCTTCGCTCAGAGAGCGCCGCAGGACATGAACTGAAACCTCATAAG  >  NZ_CP009273/3512359‑3512509
                                                                                  |                                                                    
cgcCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGgaaggtg                                                               >  1:409763‑M2/1‑83 (MQ=255)
         gTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGgaaggtgcgaacaagt                                                      >  1:115551‑M2/1‑74 (MQ=255)
           gCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGgaaggtgcgaacaagtcc                                                    <  1:15157‑M2/90‑19 (MQ=255)
            cTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGgaaggtgcgaacaagtccc                                                   >  2:284676‑M2/1‑71 (MQ=255)
                    gaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAGACTAt                                                                 >  2:63764‑M2/59‑68 (MQ=255)
                    gaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAGACTAt                                                                 <  1:63764‑M2/10‑1 (MQ=255)
                    gaagaaatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAGACTATTCGTTAATGATAAACCCGCTTc                                           >  1:196607‑M2/59‑90 (MQ=255)
                         aatgaaatgactgagtcagccgagaagaatttccccgcttattcgcaccttccCTAAGACTATTCGTTAATGATAAACCCGCTTCGCTCa                                      >  1:195718‑M2/54‑90 (MQ=255)
                                  cAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGgaaggtgcgaacaagtccctgatatgagatcatgtttgtca                             <  1:68269‑M2/90‑42 (MQ=255)
                                  cAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGgaaggtgcgaacaagtccctgatatgagatcatgtttgtca                             <  1:401419‑M2/90‑42 (MQ=255)
                                     gagtcagccgagaagaatttccccgcttattcgcaccttccCTAAGACTATTCGTTAATGATAAACCCGCTTCGCTCAGAGAGCGCCGCa                          <  1:421282‑M2/49‑1 (MQ=255)
                                        gAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGgaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctgga                       <  1:91189‑M2/90‑48 (MQ=255)
                                        gAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGgaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctgga                       <  2:190443‑M2/90‑48 (MQ=255)
                                        gAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGgaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctgga                       >  2:328159‑M2/1‑43 (MQ=255)
                                           gccgagaagaatttccccgcttattcgcaccttccCTAAGACTATTCGTTAATGATAAACCCGCTTCGCTCAGAGAGCGCCGCAGGAAAt                    >  2:399876‑M2/36‑90 (MQ=255)
                                             atatTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGgaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccat                  >  2:127244‑M2/1‑38 (MQ=255)
                                             atatTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGgaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccat                  >  2:135234‑M2/1‑38 (MQ=255)
                                              tatTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGgaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccgta                 <  1:226114‑M2/90‑54 (MQ=255)
                                                     atttccccgcttattcgcaccttccCTAAGACTATTCGTTAATGATAAACCCGCTTCGCTCAGAGAGCGCCGCAGGACATGAACTGAAAc          <  1:272319‑M2/65‑1 (MQ=255)
                                                      ttCTCTGCCATTTTTTCGGTACTACTAAGgaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacaggg         <  2:236802‑M2/90‑62 (MQ=255)
                                                           ccgcttattcgcaccttccCTAAGACTATTCGTTAATGATAAACCCGCTTCGCTCAGAGAGCGCCGCAGGACATGAACTGAAACCTCATa    <  1:123141‑M2/71‑1 (MQ=255)
                                                             gcttattcgcaccttccCTAAGACTATTCGTTAATGATAAACCCGCTTCGCTCAGAGAGCGCCGCAGGACATGAACTGAAACCTCATAAg  <  2:171224‑M2/73‑1 (MQ=255)
                                                                                  |                                                                    
CGCCCAATAGTGCTTTTTCCGGCACCCATAGGCCCAACCAGAAAGATATTGCGTTTCTCTGCCATTTTTTCGGTACTACTAAGACTATTCGTTAATGATAAACCCGCTTCGCTCAGAGAGCGCCGCAGGACATGAACTGAAACCTCATAAG  >  NZ_CP009273/3512359‑3512509

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