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.

CTACTTCACCCATGAAAAAGACAACATCATGTCGGTGTTTGCCACCGTTGGTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATCGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGTTCTGCAGGTTTTGATATGGA  >  NC_000913/2589382‑2589639
                                                                                                                                         |                                                                                                                        
ctactTCACCCATGAAAAAGACAACATCATGTCGGTGTTTGCCACCGTTGGTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGcc                                                                                                                <  1:102411/148‑1 (MQ=255)
                            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)
                                                      tGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACgctcgataatggcaaacgaggtgccggtcttactgtc                                                                                     <  8:154761‑M1/121‑38 (MQ=255)
                                                      tGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACgctcgataatggcaaacgaggtgccggtcttactgtc                                                                                     >  7:154761‑M1/1‑84 (MQ=255)
                                                          ccTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATCGAGCGTGCAACGAAGGCGtt                                                                                      <  3:13896/116‑1 (MQ=255)
                                                          ccTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATCGAGCGTGCAACGAAGGCGtt                                                                                      >  4:13896/1‑116 (MQ=255)
                                                                 gTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACgctcgataatggcaaacgaggtgccggtcttact                                                                                        >  4:5232‑M1/1‑73 (MQ=255)
                                                                 gTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACgctcgataatggcaaacgaggtgccggtcttact                                                                                        <  3:5232‑M1/107‑35 (MQ=255)
                                                                  tAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATCGAGCGTGCAACGAAGGCGTTTAACCAAATTaaa                                                                         >  8:121755/1‑121 (MQ=255)
                                                                  tAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATCGAGCGTGCAACGAAGGCGTTTAACCAAATTaaa                                                                         <  7:121755/121‑1 (MQ=255)
                                                                      gggCAAAACGTGGCGCGAATGTTTATCCGCCTGATAGACTGGAGCGAACGCGACAGTAAGACCGGCACgctcgataatggcaaacgaggtgc                                                                                                  <  6:113594‑M1/92‑25 (MQ=255)
                                                                      gggCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACgctcgataatggcaaacgaggtgc                                                                                                  >  5:113594‑M1/1‑68 (MQ=255)
                                                                          aaaaCGTGGCGCGAATGTTTATCCACCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATCGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAAt                                      >  6:148886/1‑148 (MQ=255)
                                                                                               tCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATCGAGCGTGCAACGAAGGCGTTTAACCaaa                                                                              <  8:147944/87‑1 (MQ=255)
                                                                                               tCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATCGAGCGTGCAACGAAGGCGTTTAACCaaa                                                                              >  7:147944/1‑87 (MQ=255)
                                                                                                  gCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATCGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGTTCTGc                 <  4:164317/145‑1 (MQ=255)
                                                                                                  gCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATCGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGTTCTGc                 >  3:164317/1‑145 (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                                                                      >  1:6992/1‑87 (MQ=255)
                                                                                                       aaaGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATCGAGCGTGCAACGAAGGCGTTTAACCAAATTAaagaa                                                                      <  2:6992/87‑1 (MQ=255)
                                                                                                               gAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATCGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGTTCTGCAGGTTTTGATATGGa  <  5:148886/147‑1 (MQ=255)
                                                                                                                    aaCGCGACAGTAAGACCGGCACCTCGTTTGCCATTATCGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGTTCTGCAGGTTTTGATATgg   >  7:158785/1‑141 (MQ=255)
                                                                                                                    aaCGCGACAGTAAGACCGGCACCTCGTTTGCCATTATCGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGGCTTGGTAGTTCTGCAGGTTTTGATATgg   <  8:158785/141‑1 (MQ=255)
                                                                                                                         gACAGTAAGACCGGCACCTCGTTTGCCATTATCGAGCgtgccggtcttactgtcgcgttcgctccagtctttcaggcggataaacattcgcgccacgttttgcccgttacccccagggcca                  >  7:154761‑M2/1‑37 (MQ=255)
                                                                                                                         gACAGTAAGACCGGCACCTCGTTTGCCATTATCGAGCgtgccggtcttactgtcgcgttcgctccagtctttcaggcggataaacattcgcgccacgttttgcccgttacccccagggcca                  <  8:154761‑M2/121‑85 (MQ=255)
                                                                                                                            aGTAAGACCGGCACCTCGTTTGCCATTATCGAGCgtgccggtcttactgtcgcgttcgctccagtctttcaggcggataaacattcgcgccacgttttgcccgttac                             >  4:5232‑M2/1‑34 (MQ=255)
                                                                                                                            aGTAAGACCGGCACCTCGTTTGCCATTATCGAGCgtgccggtcttactgtcgcgttcgctccagtctttcaggcggataaacattcgcgccacgttttgcccgttac                             <  3:5232‑M2/107‑74 (MQ=255)
                                                                                                                                      gCACCTCGTTTGCCATTATCGAGCgtgccggtcttactgtcgcgttcgctccagtctttcaggcggataaacattcgcgccacgttttgccc                                  >  5:113594‑M2/1‑24 (MQ=255)
                                                                                                                                      gCACCTCGTTTGCCATTATCGAGCgtgccggtcttactgtcgcgttcgctccagtctatcaggcggataaacattcgcgccacgttttgccc                                  <  6:113594‑M2/92‑69 (MQ=255)
                                                                                                                                         ccTCGTTTGCCATTATCGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGTTCTGCAGGttt           >  6:95863/1‑112 (MQ=255)
                                                                                                                                         ccTCGTTTGCCATTATCGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGTTCTGCAGGttt           <  5:95863/112‑1 (MQ=255)
                                                                                                                                         |                                                                                                                        
CTACTTCACCCATGAAAAAGACAACATCATGTCGGTGTTTGCCACCGTTGGTTCTGGCCCTGGGGGTAACGGGCAAAACGTGGCGCGAATGTTTATCCGCCTGAAAGACTGGAGCGAACGCGACAGTAAGACCGGCACCTCGTTTGCCATTATCGAGCGTGCAACGAAGGCGTTTAACCAAATTAAAGAAGCTCGCGTTATCGCCAGCAGCCCGCCAGCAATTAGCGGTCTTGGTAGTTCTGCAGGTTTTGATATGGA  >  NC_000913/2589382‑2589639

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.