New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 19727511 (0.080)22 (1.710) 16/170 0.1 95.7% 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

TCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGG  >  NZ_CP009273/1972669‑1972840
                                                                                  |                                                                                         
tCACTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggt                                                                                    >  1:32427‑M1/1‑83 (MQ=255)
   cTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcg                                                                                 >  1:299076‑M1/1‑80 (MQ=255)
   cTACACGCACATACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcg                                                                                 >  2:265972‑M1/1‑80 (MQ=255)
          cacaTACAACGGGGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagt                                                                          <  2:299076‑M1/90‑18 (MQ=255)
                        gggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagat                                                            >  2:158266‑M1/1‑59 (MQ=255)
                          gggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatca                                                          >  1:90045‑M1/1‑57 (MQ=255)
                                         ataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctgga                                           <  2:78598‑M1/90‑49 (MQ=255)
                                         ataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctgga                                           <  1:150926‑M1/90‑49 (MQ=255)
                                              gCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccat                                      <  2:102572‑M1/90‑54 (MQ=255)
                                                gTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatag                                    >  1:85670‑M1/1‑35 (MQ=255)
                                                   aTGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaac                                 <  2:73484‑M1/90‑59 (MQ=255)
                                                     gCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacag                               <  2:108595‑M1/90‑61 (MQ=255)
                                                        gATCCCACAAAACACTCAATTCCTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggt                            <  2:129233‑M1/90‑64 (MQ=255)
                                                         aTCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggtt                           <  2:331338‑M1/90‑65 (MQ=255)
                                                                  aaCACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgag                  <  2:363856‑M1/90‑74 (MQ=255)
                                                                     acTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtca               <  2:85670‑M1/90‑77 (MQ=255)
                                                                              cTTAAggaaggtgcgaacaagtccctgatatgagatc                                                           <  2:61476‑M1/37‑33 (MQ=255)
                                                                              cTTAAggaaggtgcgaacaagtccctgatatgagatc                                                           >  1:61476‑M1/1‑5 (MQ=255)
                                                                               ttAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAAc     <  2:93369/90‑1 (MQ=255)
                                                                                tAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatct                                              >  1:71933‑M1/1‑3 (MQ=255)
                                                                                tAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatct                                              <  2:71933‑M1/48‑46 (MQ=255)
                                                                                 aaggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttacctt   >  2:57943‑M1/1‑2 (MQ=255)
                                                                                  aggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttc  >  2:309356‑M1/1‑1 (MQ=255)
                                                                                  |                                                                                         
TCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGG  >  NZ_CP009273/1972669‑1972840

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.