New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 258951928 (0.770)6 (0.170) 6/232 2.4 16.4% coding (1925/3114 nt) acrD aminoglycoside/multidrug efflux system
?NC_000913 = 2589539 34 (0.970)coding (1945/3114 nt) acrD aminoglycoside/multidrug efflux system
Rejected: Frequency below/above cutoff threshold.

ATGTCGGTGTTTGCCACCGTTGGTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATCGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGTTCTGCAGGTTTTGATATGGAGTTGCAGGACCACGCTGGAGCGGGTCACGATG  >  NC_000913/2589410‑2589671
                                                                                                                                 |                                                                                                                                    
atGTCGGTGTTTGCCACCGTTGGTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTTGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATCGAGCGTGCAACGAAGGCGTTTaa                                                                                                                   <  4:57036/149‑1 (MQ=255)
      gtgtTTGCCACCGTTGGTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATCGAGCGTGCAACGAAGGc                                                                                                                         <  6:137834/137‑1 (MQ=255)
      gtgtTTGCCACCGTTGGTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATCGAGCGTGCAACGAAGGc                                                                                                                         >  5:137834/1‑137 (MQ=255)
              cACCGTTGGTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATCGAGCGTGCAACGAAGGCGTTTAACCAAAt                                                                                                             <  4:55544/141‑1 (MQ=255)
              cACCGTTGGTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATCGAGCGTGCAACGAAGGCGTTTAACCAAAt                                                                                                             >  3:55544/1‑141 (MQ=255)
                      gTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATCGAGCGTGCAACGAAGGCGTTTAACCAAATTAaag                                                                                                        <  4:101512/138‑1 (MQ=255)
                      gTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATCGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTcggaag                                                                                             >  3:101512/1‑145 (MQ=255)
                              ccTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATCGAGCGTGCAACGAAGGCGtt                                                                                                                      <  3:13896/116‑1 (MQ=255)
                              ccTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATCGAGCGTGCAACGAAGGCGtt                                                                                                                      >  4:13896/1‑116 (MQ=255)
                                      tAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATCGAGCGTGCAACGAAGGCGTTTAACCAAATTaaa                                                                                                         >  8:121755/1‑121 (MQ=255)
                                      tAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATCGAGCGTGCAACGAAGGCGTTTAACCAAATTaaa                                                                                                         <  7:121755/121‑1 (MQ=255)
                                              aaaaCGTGGCGCGAATGTTTATCCACCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATCGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAAt                                                                      >  6:148886/1‑148 (MQ=255)
                                                                   tCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATCGAGCGTGCAACGAAGGCGTTTAACCaaa                                                                                                              <  8:147944/87‑1 (MQ=255)
                                                                   tCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATCGAGCGTGCAACGAAGGCGTTTAACCaaa                                                                                                              >  7:147944/1‑87 (MQ=255)
                                                                      gCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATCGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGTTCTGc                                                 >  3:164317/1‑145 (MQ=255)
                                                                      gCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATCGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGTTCTGc                                                 <  4:164317/145‑1 (MQ=255)
                                                                        cTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATCGAGCGTGCAACGAAGGCGTTTAACCAAATTAaag                                                                                                        >  2:78896/1‑88 (MQ=255)
                                                                        cTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATCGAGCGTGCAACGAAGGCGTTTAACCAAATTAaag                                                                                                        <  1:78896/88‑1 (MQ=255)
                                                                          gAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATCGAGCGTGCAACGACGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGTTCTGc                                                 <  8:122827/141‑1 (MQ=255)
                                                                          gAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATCGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGTTCTGc                                                 >  7:122827/1‑141 (MQ=255)
                                                                           aaaGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATCGAGCGTGCAACGAAGGCGTTTAACCAAATTAaagaa                                                                                                      <  2:6992/87‑1 (MQ=255)
                                                                           aaaGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATCGAGCGTGCAACGAAGGCGTTTAACCAAATTAaagaa                                                                                                      >  1:6992/1‑87 (MQ=255)
                                                                                   gAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATCGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGTTCTGCAGGTTTTGATATGGa                                  <  5:148886/147‑1 (MQ=255)
                                                                                        aaCGCGACAGTAAGACCGGCACCTCGTTTGCCATTATCGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGTTCTGCAGGTTTTGATATgg                                   >  7:158785/1‑141 (MQ=255)
                                                                                        aaCGCGACAGTAAGACCGGCACCTCGTTTGCCATTATCGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGGCTTGGTAGTTCTGCAGGTTTTGATATgg                                   <  8:158785/141‑1 (MQ=255)
                                                                                             gACAGTAAGACCGGCACCTCGTTTGCCATTATCGAGCgtgccggtcttactgtcgcgttcgctccagtctttcaggcggataaacattcgcgccacgttttgcccgttacccccagggcca                                                  <  8:154761‑M2/121‑85 (MQ=255)
                                                                                             gACAGTAAGACCGGCACCTCGTTTGCCATTATCGAGCgtgccggtcttactgtcgcgttcgctccagtctttcaggcggataaacattcgcgccacgttttgcccgttacccccagggcca                                                  >  7:154761‑M2/1‑37 (MQ=255)
                                                                                                aGTAAGACCGGCACCTCGTTTGCCATTATCGAGCgtgccggtcttactgtcgcgttcgctccagtctttcaggcggataaacattcgcgccacgttttgcccgttac                                                             <  3:5232‑M2/107‑74 (MQ=255)
                                                                                                aGTAAGACCGGCACCTCGTTTGCCATTATCGAGCgtgccggtcttactgtcgcgttcgctccagtctttcaggcggataaacattcgcgccacgttttgcccgttac                                                             >  4:5232‑M2/1‑34 (MQ=255)
                                                                                                          gCACCTCGTTTGCCATTATCGAGCgtgccggtcttactgtcgcgttcgctccagtctttcaggcggataaacattcgcgccacgttttgccc                                                                  >  5:113594‑M2/1‑24 (MQ=255)
                                                                                                          gCACCTCGTTTGCCATTATCGAGCgtgccggtcttactgtcgcgttcgctccagtctatcaggcggataaacattcgcgccacgttttgccc                                                                  <  6:113594‑M2/92‑69 (MQ=255)
                                                                                                             ccTCGTTTGCCATTATCGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGTTCTGCAGGttt                                           <  5:95863/112‑1 (MQ=255)
                                                                                                             ccTCGTTTGCCATTATCGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGTTCTGCAGGttt                                           >  6:95863/1‑112 (MQ=255)
                                                                                                              cTCGTTTGCCATTATCGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAAt                                                                      <  4:42005/84‑1 (MQ=255)
                                                                                                              cTCGTTTGCCATTATCGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAAt                                                                      >  3:42005/1‑84 (MQ=255)
                                                                                                                 gTTTGCCATTATCGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGTTCTGCAGGTTTTGATATGGAGTTGCAGGACCACGCTGGAGCGGGTCACGATg  >  2:124872/1‑149 (MQ=255)
                                                                                                                      ccATTATCGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGTTCTGCAGGTTTTGATATGGAGTTGc                             >  6:4690/1‑117 (MQ=255)
                                                                                                                      ccATTATCGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGTTCTGCAGGTTTTGATATGGAGTTGc                             <  5:4690/117‑1 (MQ=255)
                                                                                                                       cATTATCGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGTTCTGCAGGTTTTGATATGGAGTTGCAGGACCACGCTGGAGCgg           <  7:39188/134‑1 (MQ=255)
                                                                                                                       cATTATCGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGTTCTGCAGGTTTTGATATGGAGTTGCAGGACCACGCTGGAGCgg           >  8:39188/1‑134 (MQ=255)
                                                                                                                                 |                                                                                                                                    
ATGTCGGTGTTTGCCACCGTTGGTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATCGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGTTCTGCAGGTTTTGATATGGAGTTGCAGGACCACGCTGGAGCGGGTCACGATG  >  NC_000913/2589410‑2589671

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