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

GCGGAACTGCGTCGTCCGTTAGGTATCGGCATGGTCGGCGGTCTGATTGTCAGCCAGGTGCTGACGCTGTTTACCACGCCGGTGATTTATTTGCTGTTCGACCGCCTGGCATTGTGGACCAAAAGCCGCTTTGCC  >  NC_000913/2158230‑2158364
                                                                   |                                                                   
gcgGAACTGCGTCGTCCGTTAGGTATCGGCATGGTCGGCGGTCTGATTGTCAGCCAGGTGCTGACGCt                                                                     <  1:814427/68‑1 (MQ=255)
          GTCGTCCGTTAGGTATCGGCATGGTCGGCGGTCTGATTGTCAGCCAGGTGCTGACGCTGTTTACCACg                                                           <  1:848180‑M1/68‑1 (MQ=255)
                   tAGGTATCGGCAGGGTCGGCGGTCTGATTGTCAGCCAGGTGCTGACGCTGTTTACCTCGCCGGTGatt                                                  <  1:1004212/68‑1 (MQ=255)
                   tAGGTAGCGGCATGGTCGTCGGTCTGATTGTCAGCCAGGTGCTGACGCTGTTTACCACGCCGGTGatt                                                  <  1:1121589/68‑1 (MQ=255)
                     ggTATCGGCATGGTCGGCGGTCTGATTGTCAGCCAGGTGCTGACGCTGTTTACAACGACGGTGattta                                                >  1:668887/1‑68 (MQ=255)
                     gctgacaatcagaccgccgaccatgccgatacctaacggacgacgcTGTTTACCACGCCGGTGattta                                                >  1:201844‑M2/47‑68 (MQ=255)
                     gctgacaatcagaccgccgaccatgccgatacctaacggacgacgcTGTTTACCACGCCGGTGatt                                                  <  1:155501‑M2/20‑1 (MQ=255)
                     gctgacaatcagaccgccgaccatgccgatacctaacggacgacgcTGTTTACCACGCCGGTGatt                                                  >  2:155501‑M2/47‑66 (MQ=255)
                           ggCATGGGCGGCGGTCTGATTGTCAGCCAGGTGCTGACGCTGTTTACCACGCCGGTGATTTATTTGCt                                          <  2:533469/68‑1 (MQ=255)
                             cATGGTCGGCGGTCTGATTGTCAGCCTGGTGCTGACGCTGTTTACCACGCCGGTGATTTATTTGCTGt                                        >  1:1016316/1‑68 (MQ=255)
                                 accgccgaccatgccgatacctaacggacgacgcTGTTTACCACGCCGGTGATTTATTTGCTGTTCGa                                    >  2:446706‑M2/35‑68 (MQ=255)
                                   cggcggTCTGATTGTCAGCCAGGTGCTGACGCTGTTTACCACGCCGGTGATTTTTTTGCTGTTCGAcc                                  >  1:1113165/1‑68 (MQ=255)
                                          cTGATTGTCAGCCAGGTGCTGACGCTGTTTACCACGCCGGTGATTTATTTGCTGTTCGTCCGCCTGGc                           <  1:591905/68‑1 (MQ=255)
                                               tGTCAGCCAGGTGCTGACGCTGTTTACCACGCCGGTGATTTATTTGCTGTTCGACCGCCTGGCATtgt                      >  2:610554/1‑68 (MQ=255)
                                               tGTCAGCCAGGTGCTGACGCTGTTTACCACGCCGGTGATTTATTTGCTGTTCGACCGCCTGGCATtgt                      >  2:857181/1‑68 (MQ=255)
                                                           gCTGACGCTGTTTACCACGCCGGTGATTTATTTGCTGTTCGACCGCCTGGCATTGTGGACCAAAAGcc          >  1:512417/1‑68 (MQ=255)
                                                               acgcTGTTTACCACGCCGGTGATTTATTTGCTGTTCGACCGCCTGGCATTGTGGACCAAAAGCCGCtt      >  2:305945‑M2/5‑68 (MQ=255)
                                                                cgcTGTTTACCACGCCGGTGATTTATTTGCTGTTCGACCGCCTGGc                           <  1:284316‑M2/43‑1 (MQ=255)
                                                                cgcTGTTTACCACGCCGGTGATTTATTTGCTGTTCGACCGCCTGGc                           >  2:284316‑M2/4‑46 (MQ=255)
                                                                cgcTGTTTACCACGCCGGTGATTTAATTGCTGTTCGACCGCCTGGCATTGTGGACCAAAAGCCGCttt     >  2:161734‑M2/4‑68 (MQ=255)
                                                                   tGTATACCACGCCGGTGATTTAAATGCAGTTCGACAGCATGGCATAGAGGACCAAAAGCCGCTATGcc  >  2:278500/1‑68 (MQ=255)
                                                                   |                                                                   
GCGGAACTGCGTCGTCCGTTAGGTATCGGCATGGTCGGCGGTCTGATTGTCAGCCAGGTGCTGACGCTGTTTACCACGCCGGTGATTTATTTGCTGTTCGACCGCCTGGCATTGTGGACCAAAAGCCGCTTTGCC  >  NC_000913/2158230‑2158364

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.