New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 2158238 =13 (0.420)4 (0.140) 3/122 2.7 25.8% coding (2976/3123 nt) mdtB multidrug efflux system, subunit B
?NC_000913 2158297 = 11 (0.390)coding (3035/3123 nt) mdtB multidrug efflux system, subunit B

TGACCACTCTGGCGGCTCTGCTTGGCGCGCTGCCGCTGATGTTGAGTACCGGGGTCGGCGCGGAACTGCGTCGTCCGTTAGGTATCGGCATGGTCGGCGGTCTGATTGTCAGCCAGGTGCTGACGCT  >  NC_000913/2158171‑2158297
                                                                   |                                                           
tgaCCACTCTGGCGGCTCTGCTTGGCGCGCTGCCGCTGATGTTGAGTACCGGGGTCGGCGCGGAACTg                                                             >  2:672875/1‑68 (MQ=255)
  aaagcggcttttggtccacaatgccaggcggtcgaacagcaattaaatcaccggcgtggtaaacaGcg                                                           >  2:161734‑M1/66‑68 (MQ=255)
   aagcggcttttggtccacaatgccaggcggtcgaacagcaaataaatcaccggcgtggtaaacaGcgt                                                          >  2:305945‑M1/65‑68 (MQ=255)
          ggcggcTCTGCTTGGCGCGCTGCCGCTGATGTTGAGTACCGGGGTCGGCGCGGAACTGCGTCGTCCGt                                                   <  2:78134/68‑1 (MQ=255)
          ggcggcGCTGCTTGGCGCGCTGCCGCTGATGTTGAGTACCGGCGGCGGCGCGGAACTGCGTCGTCCGt                                                   <  1:1133615/68‑1 (MQ=255)
             ggcTCTGCTTGGCGCGCTGCCGCTGATGTTGAGTACCGGGGTCGGCGCGGAACTGCGTCGTCCGTTAg                                                >  2:369379/1‑68 (MQ=255)
               ctctGCTTGGCGCGCTGCCGCTGATGTTGAGTACCGGGGTCGGCGCGGAACTGCGTCGTCCGTTAGGt                                              <  1:805785/68‑1 (MQ=255)
                        gccaggcggtcgaacagcaaataaatcaccggcgtggtaaacaGcg                                                           <  1:284316‑M1/3‑1 (MQ=255)
                        gccaggcggtcgaacagcaaataaatcaccggcgtggtaaacaGcg                                                           >  2:284316‑M1/44‑46 (MQ=255)
                          gcgcTGCCGCTGATGTTGAGTACCGGGGTCGGCGCGGAACTGCGTCGTCCGTTAGGTATCGGCATGGt                                   >  2:1120388/1‑68 (MQ=255)
                                 tcgaacagcaaataaatcaccggcgtggtaaacaGCGTCGTCCGTTAGGTATCGGCATGGTCGGCGGt                            >  2:446706‑M1/35‑68 (MQ=255)
                                 cGCTGATGTTGAGTACCGGGGGCGGCGCGGAACTGCGTCGTCCGTTAGGTATCGGCATGGTCGGCGGt                            <  2:1140753/68‑1 (MQ=255)
                                     gATGTTGAGTTCCGGGGTCGGCGCGGAACTGCGTCGTCCGTTAGGTATCGGCATGGTCGGCGGTCTGa                        <  2:16717/68‑1 (MQ=255)
                                             taaatcaccggcgtggtaaacaGCGTCGTCCGTTAGGTATCGGCATGGTCGGCGGTCTGATTGTCAGc                >  1:201844‑M1/23‑68 (MQ=255)
                                               aatcaccggcgtggtaaacaGCGTCGTCCGTTAGGTATCGGCATGGTCGGCGGTCTGATTGTCAGc                >  2:155501‑M1/21‑66 (MQ=255)
                                               aatcaccggcgtggtaaacaGCGTCGTCCGTTAGGTATCGGCATGGTCGGCGGTCTGATTGTCAGc                <  1:155501‑M1/46‑1 (MQ=255)
                                               aCCGGGGTCGGCGCGGAACTGCGTCGTCCGTTAGGTATCGGCATGGTCGGCGGTCTGATTGTCAGCCa              >  2:714610/1‑68 (MQ=255)
                                                   gggTCGGCGCGGAACTGCGTCGTCCGTTAGGTATCGGCATGGTCGGCGGTCTGATTGTCAGCCAGGTg          <  2:340462/68‑1 (MQ=255)
                                                    ggTCGGCGCGGAACTGCGTCGTCCGTTAGGTATCGGCATGGTCGGCGGTCTGATTGTCAGCCAGGTGc         >  1:534411/1‑68 (MQ=255)
                                                     gTCGGCGCGGAACTGCGTCGTCCGTTAGGTATCGGCATGGTCGGCGGTCTGATTGTCAGCCAGGTGCt        <  1:613912/68‑1 (MQ=255)
                                                           gcgGAACTGCGTCGTCCGTTAGGTATCGGCATGGTCGGCGGTCTGATTGTCAGCCAGGTGCTGACGCt  <  1:814427/68‑1 (MQ=255)
                                                                   |                                                           
TGACCACTCTGGCGGCTCTGCTTGGCGCGCTGCCGCTGATGTTGAGTACCGGGGTCGGCGCGGAACTGCGTCGTCCGTTAGGTATCGGCATGGTCGGCGGTCTGATTGTCAGCCAGGTGCTGACGCT  >  NC_000913/2158171‑2158297

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