New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 884909 =22 (0.590)3 (0.090) 3/246 3.4 11.2% intergenic (+4/+37) mdfA/ybjH multidrug efflux system protein/uncharacterized protein
?NC_000913 884940 = 27 (0.780)intergenic (+35/+6) mdfA/ybjH multidrug efflux system protein/uncharacterized protein
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATGGCGAACATAAGGATGCGGGCGGTAAGCCGGAGCTGGCGCAGGAGCCACATAGACGGTTTGCGGCACAGTCGCGACG  >  NC_000913/884827‑885084
                                                                                                                 |                                                                                                                                                
tGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTc                                                                                                               >  2:38346/1‑149 (MQ=255)
    aaCGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGAt                                                                                                            >  7:128136/1‑148 (MQ=255)
    aaCGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGAt                                                                                                            <  5:58904/148‑1 (MQ=255)
     aCGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATTGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATaa                                                                                                          <  2:168900/149‑1 (MQ=255)
      cGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATtctc                                                                                                                                                <  8:66178/110‑1 (MQ=255)
      cGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATtctc                                                                                                                                                >  7:66178/1‑110 (MQ=255)
              tgtgGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCtt                                                                                                 >  6:127591/1‑149 (MQ=255)
                    tGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTg                                                                                                                     <  6:90395/123‑1 (MQ=255)
                    tGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTg                                                                                                                     >  5:90395/1‑123 (MQ=255)
                                ttATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGAt                                                                               >  3:78893/1‑149 (MQ=255)
                                  aTCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATgg                                                                             >  8:14429/1‑149 (MQ=255)
                                  aTCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATgg                                                                             >  6:51197/1‑149 (MQ=255)
                                       tttAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATGGCGaa                                                                         <  7:14429/148‑1 (MQ=255)
                                         tAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATGGCGAACa                                                                       >  6:57928/1‑148 (MQ=255)
                                             aGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTg                                                                                      >  8:33704/1‑129 (MQ=255)
                                             aGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTg                                                                                      <  7:33704/129‑1 (MQ=255)
                                                     aGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATGGCGAACATAAGGATgcgggc                                                          <  5:57928/149‑1 (MQ=255)
                                                      gATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATGGCGAACATAAGGATgcgggcg                                                         <  8:128136/149‑1 (MQ=255)
                                                                    aCGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATGGCGAACATAAGGATGCGGGCGGTAAGCCGGAGCTg                                           >  3:98661/1‑149 (MQ=255)
                                                                               aaaaaTGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATGGCGAACATAAGGATGCGGGCGGTAAGCCGGAGCTGGCGCATGAGcc                                <  5:51197/149‑1 (MQ=255)
                                                                                         gACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATGGCGAACATAAGGATGCGGGCGGTAAGCCGGAGCTGGCGCAGGAGCCACATAGACgg                      >  3:129703/1‑149 (MQ=255)
                                                                                               tttGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATGGCGAACATAAGGATGCGGGCGGTAAGCCGGCGCTGGCGCAGGAGCCACATAGCCGGTGTGCg                >  8:100161/1‑149 (MQ=255)
                                                                                                  gtgCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATGGCGAACATAAGGATGCGGGCGGTAAGCCGGAGCTGGCGCAGGAGCca                               <  7:113678/131‑1 (MQ=255)
                                                                                                  gtgCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATGGCGAACATAAGGATGCGGGCGGTAAGCCGGAGCTGGCGCAGGAGCca                               >  8:113678/1‑131 (MQ=255)
                                                                                                   tgCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATGGCGAACATAAGGATGCGGGCGGTAAGCCGGAGCTGGCGCAGGAGCCACATAGACGGTTTGCGGca             >  4:647/1‑148 (MQ=255)
                                                                                                     cGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATGGCGAACATAAGGATGCGGGCGGTAAGCCGGAGCTGGCGCAGGAGCCACATAGACGGTTTGCGGcaca           <  4:129703/148‑1 (MQ=255)
                                                                                                       aTCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATGGCGAACATAAGGATGCGGGCGGTAAGCCGGAGCTGGCGCAGGAGCCACATAGACGGTTTGCGGCACAGTc        >  2:179241/1‑149 (MQ=255)
                                                                                                          aGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATGGCGAACATAAGGATGCGGGCGGTAAGCCGGAg                                              >  3:103856/1‑108 (MQ=255)
                                                                                                          aGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATGGCGAACATAAGGATGCGGGCGGTAAGCCGGAg                                              <  4:103856/108‑1 (MQ=255)
                                                                                                          aGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATGGCGAACATAAGGATGCGGGCGGTAAGCCGGAGCTGGCGCAGGAGCCACATAGACGGTTTGCGGCACAGTcgc      <  5:127591/148‑1 (MQ=255)
                                                                                                          aGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATGGCGAACATAAGGATGCGGGCGGGAAGCCGGAGCTGGCGCAGGAGCCACATAGACGGTTTGCGGCACAGTcgcg     <  4:78893/149‑1 (MQ=255)
                                                                                                              aTTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATGGCGAACATAAGGATGCGGGCGGTAAGCCGGAGCTGGCGCAGGAGCCACATAGACGGTTTGCGGCACAGTCGCGACg  >  3:47452/1‑148 (MQ=255)
                                                                                                              aTTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATGGCGAACATAAGGATGCGGGCGGTAAGCCGGAGCTGGCGCAGGAGCCACATAGACGGTTTGCGGCACAGTCGCGACg  <  4:47452/148‑1 (MQ=255)
                                                                                                                 |                                                                                                                                                
TGTCAACGGAATTTTGTGGCTGTCGCTGATGGTTATCTTTTTAAAAGATAAACAGATGGGAAATTCTCACGAAGGGTAAAAAAATGCCTGACTGCTTTGTGCGATCAGGCATTCTCGAATTAATGGTGATGGTCGTCAATCTGGTGTTCGATAACCATCCCTTCACCTACGCTGGCAAGATGGCGAACATAAGGATGCGGGCGGTAAGCCGGAGCTGGCGCAGGAGCCACATAGACGGTTTGCGGCACAGTCGCGACG  >  NC_000913/884827‑885084

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