New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 19727511 (0.040)19 (0.850) 11/164 0.6 95.0% 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

ACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGG  >  NZ_CP009273/1972676‑1972840
                                                                           |                                                                                         
acGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaaca                                                                             <  2:3169‑M1/90‑15 (MQ=255)
 cGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaa                                                                            <  1:50371‑M1/90‑16 (MQ=255)
           aCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgata                                                                  >  1:373779‑M1/1‑65 (MQ=255)
               aGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgag                                                              >  1:126416‑M1/1‑61 (MQ=255)
                ggggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgaga                                                             >  2:521453‑M1/1‑60 (MQ=255)
                 gggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagat                                                            <  2:373779‑M1/90‑32 (MQ=255)
                          gATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgt                                                   <  1:521453‑M1/90‑41 (MQ=255)
                          gATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgt                                                   <  1:94771‑M1/90‑41 (MQ=255)
                                                         aaaaCACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatg                    <  2:162822‑M1/90‑72 (MQ=255)
                                                           aaCACTCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgag                  >  2:429993‑M1/1‑17 (MQ=255)
                                                                tCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatc             >  1:670380‑M1/1‑12 (MQ=255)
                                                                tCAATTACTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatc             >  2:373410‑M1/1‑12 (MQ=255)
                                                                   cttccTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAAt          <  1:135885/86‑1 (MQ=255)
                                                                      aCTTAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaactta       >  2:102882‑M1/1‑6 (MQ=255)
                                                                         tAAggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcatcttacct    >  1:679708‑M1/1‑3 (MQ=255)
                                                                          aaggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttacctt   >  1:633598‑M1/1‑2 (MQ=255)
                                                                           aggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttc  >  2:378004‑M1/1‑1 (MQ=255)
                                                                           aggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagccatagaacagggttcatcatgagtcatcaacttaccttc  >  1:515721‑M1/1‑1 (MQ=255)
                                                                           aggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagc                                         >  1:546371‑M1/1‑1 (MQ=255)
                                                                           aggaaggtgcgaacaagtccctgatatgagatcatgtttgtcatctggagc                                         <  2:546371‑M1/51‑51 (MQ=255)
                                                                           |                                                                                         
ACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGG  >  NZ_CP009273/1972676‑1972840

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