Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 4,299,628 Δ25,538 bp [mdtP][phnC] 27 genes

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 4299628 4325165 25538 72 [0] [0] 70 [mdtP]–[phnC] [mdtP],mdtO,mdtN,ytcA,yjcS,alsK,alsE,alsC,alsA,alsB,alsR,rpiB,yjdP,phnP,phnO,phnN,phnM,phnL,phnK,phnJ,phnI,phnH,phnG,phnF,phnE,phnD,[phnC]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 42996270 (0.000)68 (0.940) 52/98 0.0 100% coding (1404/1467 nt) mdtP putative multidrug efflux pump outer membrane channel
?NC_000913 4325166 = 0 (0.000)intergenic (‑1/+132) phnC/yjdN phosphonate ABC transporter ATP binding subunit/conserved protein YjdN

TGATTGAAGGTTTTGGCGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGAC  >  NC_000913/4325119‑4325214
                                               |                                                
cgggacctgcctgatacccgccgcccagcgatttcatcaactgaatgGAt                                                >  1:2798897‑M2/48‑50 (MQ=255)
 gggacctgcctgatacccgccgcccagcgatttcatcaactgaatgGATg                                               >  1:3588283‑M2/47‑50 (MQ=255)
 gggacctgcctgatacccgccgcccagcgatttcatcaactgaatgGATg                                               >  1:4442954‑M2/47‑50 (MQ=255)
 gggacctgcctgatacccgccgcccagcgatttcatcaactgaatgGATg                                               >  1:3486690‑M2/47‑50 (MQ=255)
  ggacctgcctgatacccgccgcccagcgatttcatcaactgaatgGATGc                                              <  1:4002206‑M2/5‑1 (MQ=255)
   gacctgcctgatacccgccgcccagcgatttcatcaactgaatgGATGCa                                             <  1:2920127‑M2/6‑1 (MQ=255)
     cctgcctgatacccgccgcccagcgatttcatcaactgaatgGATGCAGc                                           >  1:5550074‑M2/43‑50 (MQ=255)
     cctgcctgatacccgccgcccagcgatttcatcaactgaatgGATGCAGc                                           >  1:2351006‑M2/43‑50 (MQ=255)
      ctgcctgatacccgccgcccagcgatttcatcaactgaatgGATGCAGcc                                          <  1:2486424‑M2/9‑1 (MQ=255)
      ctgcctgatacccgccgcccagcgatttcatcaactgaatgGATGCAGcc                                          >  1:5297951‑M2/42‑50 (MQ=255)
       tgcctgatacccgccgcccagcgatttcatcaactgaatgGATGCAGCCt                                         >  1:4827331‑M2/41‑50 (MQ=255)
       tgcctgatacccgccgcccagcgatttcatcaactgaatgGATGCAGCCt                                         >  1:5316181‑M2/41‑50 (MQ=255)
        gcctgatacccgccgcccagcgatttcatcaactgaatgGATGCAGCCTc                                        >  1:5487643‑M2/40‑50 (MQ=255)
         cctgatacccgccgcccagcgatttcatcaactgaatgGATGCAGCCTcc                                       <  1:4115245‑M2/12‑1 (MQ=255)
          ctgatacccgccgcccagcgatttcatcaactgaatgGATGCAGCCTccc                                      >  1:675760‑M2/38‑50 (MQ=255)
          ctgatacccgccgcccagcgatttcatcaactgaatgGATGCAGCCTccc                                      >  1:23497‑M2/38‑50 (MQ=255)
          ctgagacccgccgcccagcgctttcatcaactgaatgGATGCAGCCTccc                                      >  1:5351458‑M2/38‑50 (MQ=255)
           tgatacccgccgcccagcgatttcatcaactgaatgGATGCAGCCTccc                                      >  1:3737297‑M2/37‑49 (MQ=255)
           tgatacccgccgcccagcgatttcatcaactgaatgGATGCAGCCTccc                                      >  1:4772004‑M2/37‑49 (MQ=255)
             atacccgccgcccagcgatttcatcaactgaatgGATGCAGCCTCCCTaa                                   <  1:4303351‑M2/16‑1 (MQ=255)
             atacccgccgcccagcgatttcatcaactgaatgGATGCAGCCTCCCTa                                    <  1:4048301‑M2/15‑1 (MQ=255)
             atacccgccgcccagcgatttcatcaactgaatgGATGCAGCATCCCt                                     <  1:1875822‑M2/14‑1 (MQ=255)
                cccgccgcccagcgatttcatcaactgaatgGATGCAGCCTCCCTaaaa                                 <  1:4248155‑M2/18‑1 (MQ=255)
                cccgccgcccagcgatttcatcaactgaatgGATGCAGCCTCCCTaaaa                                 >  1:4018147‑M2/32‑49 (MQ=255)
                 ccgccgcccagcgatttcatcaactgaatgGATGCAGCCTCCCTaaaaaa                               <  1:4997297‑M2/20‑1 (MQ=255)
                 ccgccgcccagcgatttcatcaactgaatgGATGCAGCCTCCCTaaaaaa                               <  1:3517158‑M2/20‑1 (MQ=255)
                  cgccgcccagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAg                              >  1:1399255‑M2/30‑50 (MQ=255)
                    ccgcccagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTg                            <  1:3765180‑M2/23‑1 (MQ=255)
                    ccgcccagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTg                            >  1:3651218‑M2/28‑50 (MQ=255)
                     cgcccagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTgg                           <  1:2391621‑M2/24‑1 (MQ=255)
                      gcccagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGc                          >  1:5169494‑M2/26‑50 (MQ=255)
                       cccagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCt                         >  1:2979664‑M2/25‑50 (MQ=255)
                        ccagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTc                        >  1:1446197‑M2/24‑50 (MQ=255)
                        ccagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTc                        <  1:2974622‑M2/27‑1 (MQ=255)
                         cagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTcc                       <  1:2203400‑M2/28‑1 (MQ=255)
                         cagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTc                        >  1:5250321‑M2/23‑49 (MQ=255)
                          agcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCa                      <  1:530185‑M2/29‑1 (MQ=255)
                           gcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCAt                     <  1:685679‑M2/30‑1 (MQ=255)
                           gcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCAt                     >  1:137431‑M2/21‑50 (MQ=255)
                            cgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATc                    <  1:2621233‑M2/31‑1 (MQ=255)
                            cgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATc                    >  1:3352945‑M2/20‑50 (MQ=255)
                             gatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCg                   >  1:1596908‑M2/19‑50 (MQ=255)
                              atttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGt                  >  1:1144733‑M2/18‑50 (MQ=255)
                              atttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGt                  <  1:3519823‑M2/33‑1 (MQ=255)
                              atttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGt                  <  1:4509949‑M2/33‑1 (MQ=255)
                              atttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGt                  >  1:332698‑M2/18‑50 (MQ=255)
                               tttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTg                 <  1:1394064‑M2/34‑1 (MQ=255)
                                ttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTgg                >  1:2346198‑M2/16‑50 (MQ=255)
                                 tcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGc               >  1:931540‑M2/15‑50 (MQ=255)
                                 tcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGc               >  1:4574439‑M2/15‑50 (MQ=255)
                                   atcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCgg             <  1:4535743‑M2/38‑1 (MQ=255)
                                   atcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCg              >  1:4530812‑M2/13‑49 (MQ=255)
                                   atcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCg              <  1:2602228‑M2/37‑1 (MQ=255)
                                    tcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGa            >  1:4075550‑M2/12‑50 (MQ=255)
                                     caactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGat           >  1:4791837‑M2/11‑50 (MQ=255)
                                      aactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGatt          <  1:3353759‑M2/41‑1 (MQ=255)
                                      aactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGatt          <  1:538637‑M2/41‑1 (MQ=255)
                                      aactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGat           <  1:2388732‑M2/40‑1 (MQ=255)
                                      aactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGat           >  1:3546861‑M2/10‑49 (MQ=255)
                                      aactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGat           >  1:369535‑M2/10‑49 (MQ=255)
                                      aactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGat           <  1:475431‑M2/40‑1 (MQ=255)
                                       actgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGatta         >  1:4479791‑M2/9‑50 (MQ=255)
                                        ctgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGattat        <  1:5010690‑M2/43‑1 (MQ=255)
                                        ctgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGattat        <  1:2030711‑M2/43‑1 (MQ=255)
                                        ctgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGattat        <  1:5359338‑M2/43‑1 (MQ=255)
                                           aatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGattatt       <  1:2138543‑M2/44‑1 (MQ=255)
                                           aatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATtgt     >  1:5124589‑M2/5‑50 (MQ=255)
                                           aatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATtgt     <  1:174210‑M2/46‑1 (MQ=255)
                                               gatgCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGAc  <  1:4615322/49‑1 (MQ=255)
                                               gatgCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGAc  <  1:4772558/49‑1 (MQ=255)
                                               |                                                
TGATTGAAGGTTTTGGCGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGAC  >  NC_000913/4325119‑4325214

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