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

GACGGGACCTGCCTGATACCCGCCGCCCAGCGATTTCATCAACTGAATGCTTTGGATCACCCGGCGGCTGTCCAGCATCAGTAATGACATCTCTT  >  NC_000913/4299579‑4299673
                                                |                                              
gacgGGACCTGCCTGATACCCGCCGCCCAGCGATTTCATCAACTGAATg                                                >  1:475804/1‑49 (MQ=255)
gacgGGACCTGCCTGATACCCGCCGCCCAGCGATTTCATCAACTGAATg                                                >  1:271425/1‑49 (MQ=255)
gacgGGACCTGCCTGATACCCGCCGCCCAGCGATTTCATCAACTGAATg                                                >  1:4689407/1‑49 (MQ=255)
gacgGGACCTGCCTGATACCCGCCGCCCAGCGATTTCATCAACTGAATg                                                >  1:45925/1‑49 (MQ=255)
  cGGGACCTGCCTGATACCCGCCGCCCAGCGATTTCATCAACTGAATGgat                                             >  1:2798897‑M1/1‑47 (MQ=255)
   gggACCTGCCTGATACCCGCCGCCCAGCGATTTCATCAACTGAATGgatg                                            >  1:4442954‑M1/1‑46 (MQ=255)
   gggACCTGCCTGATACCCGCCGCCCAGCGATTTCATCAACTGAATGgatg                                            >  1:3588283‑M1/1‑46 (MQ=255)
   gggACCTGCCTGATACCCGCCGCCCAGCGATTTCATCAACTGAATGgatg                                            >  1:3486690‑M1/1‑46 (MQ=255)
    ggACCTGCCTGATACCCGCCGCCCAGCGATTTCATCAACTGAATGgatgc                                           <  1:4002206‑M1/50‑6 (MQ=255)
     gACCTGCCTGATACCCGCCGCCCAGCGATTTCATCAACTGAATGgatgca                                          <  1:2920127‑M1/50‑7 (MQ=255)
       cctgcctgATACCCGCCGCCCAGCGATTTCATCAACTGAATGgatgcagc                                        >  1:5550074‑M1/1‑42 (MQ=255)
       cctgcctgATACCCGCCGCCCAGCGATTTCATCAACTGAATGgatgcagc                                        >  1:2351006‑M1/1‑42 (MQ=255)
        ctgcctgATACCCGCCGCCCAGCGATTTCATCAACTGAATGgatgcagcc                                       <  1:2486424‑M1/50‑10 (MQ=255)
        ctgcctgATACCCGCCGCCCAGCGATTTCATCAACTGAATGgatgcagcc                                       >  1:5297951‑M1/1‑41 (MQ=255)
         tgcctgATACCCGCCGCCCAGCGATTTCATCAACTGAATGgatgcagcct                                      >  1:5316181‑M1/1‑40 (MQ=255)
         tgcctgATACCCGCCGCCCAGCGATTTCATCAACTGAATGgatgcagcct                                      >  1:4827331‑M1/1‑40 (MQ=255)
          gcctgATACCCGCCGCCCAGCGATTTCATCAACTGAATGgatgcagcctc                                     >  1:5487643‑M1/1‑39 (MQ=255)
           cctgATACCCGCCGCCCAGCGATTTCATCAACTGAATGgatgcagcctcc                                    <  1:4115245‑M1/50‑13 (MQ=255)
            ctgATACCCGCCGCCCAGCGATTTCATCAACTGAATGgatgcagcctccc                                   >  1:23497‑M1/1‑37 (MQ=255)
            ctgATACCCGCCGCCCAGCGATTTCATCAACTGAATGgatgcagcctccc                                   >  1:675760‑M1/1‑37 (MQ=255)
            ctgAGACCCGCCGCCCAGCGCTTTCATCAACTGAATGgatgcagcctccc                                   >  1:5351458‑M1/1‑37 (MQ=255)
             tgATACCCGCCGCCCAGCGATTTCATCAACTGAATGgatgcagcctccc                                   >  1:3737297‑M1/1‑36 (MQ=255)
             tgATACCCGCCGCCCAGCGATTTCATCAACTGAATGgatgcagcctccc                                   >  1:4772004‑M1/1‑36 (MQ=255)
               aTACCCGCCGCCCAGCGATTTCATCAACTGAATGgatgcagcctccctaa                                <  1:4303351‑M1/50‑17 (MQ=255)
               aTACCCGCCGCCCAGCGATTTCATCAACTGAATGgatgcagcctcccta                                 <  1:4048301‑M1/49‑16 (MQ=255)
               aTACCCGCCGCCCAGCGATTTCATCAACTGAATGgatgcagcatccct                                  <  1:1875822‑M1/48‑15 (MQ=255)
                  cccGCCGCCCAGCGATTTCATCAACTGAATGgatgcagcctccctaaaa                              <  1:4248155‑M1/49‑19 (MQ=255)
                  cccGCCGCCCAGCGATTTCATCAACTGAATGgatgcagcctccctaaaa                              >  1:4018147‑M1/1‑31 (MQ=255)
                   ccgccgCCCAGCGATTTCATCAACTGAATGgatgcagcctccctaaaaaa                            <  1:3517158‑M1/50‑21 (MQ=255)
                   ccgccgCCCAGCGATTTCATCAACTGAATGgatgcagcctccctaaaaaa                            <  1:4997297‑M1/50‑21 (MQ=255)
                    cgccgcCCAGCGATTTCATCAACTGAATGgatgcagcctccctaaaaaag                           >  1:1399255‑M1/1‑29 (MQ=255)
                      ccgccCAGCGATTTCATCAACTGAATGgatgcagcctccctaaaaaagtg                         >  1:3651218‑M1/1‑27 (MQ=255)
                      ccgccCAGCGATTTCATCAACTGAATGgatgcagcctccctaaaaaagtg                         <  1:3765180‑M1/50‑24 (MQ=255)
                       cgccCAGCGATTTCATCAACTGAATGgatgcagcctccctaaaaaagtgg                        <  1:2391621‑M1/50‑25 (MQ=255)
                        gccCAGCGATTTCATCAACTGAATGgatgcagcctccctaaaaaagtggc                       >  1:5169494‑M1/1‑25 (MQ=255)
                         cccAGCGATTTCATCAACTGAATGgatgcagcctccctaaaaaagtggct                      >  1:2979664‑M1/1‑24 (MQ=255)
                          ccAGCGATTTCATCAACTGAATGgatgcagcctccctaaaaaagtggctc                     <  1:2974622‑M1/50‑28 (MQ=255)
                          ccAGCGATTTCATCAACTGAATGgatgcagcctccctaaaaaagtggctc                     >  1:1446197‑M1/1‑23 (MQ=255)
                           cAGCGATTTCATCAACTGAATGgatgcagcctccctaaaaaagtggctcc                    <  1:2203400‑M1/50‑29 (MQ=255)
                           cAGCGATTTCATCAACTGAATGgatgcagcctccctaaaaaagtggctc                     >  1:5250321‑M1/1‑22 (MQ=255)
                            aGCGATTTCATCAACTGAATGgatgcagcctccctaaaaaagtggctcca                   <  1:530185‑M1/50‑30 (MQ=255)
                             gCGATTTCATCAACTGAATGgatgcagcctccctaaaaaagtggctccat                  <  1:685679‑M1/50‑31 (MQ=255)
                             gCGATTTCATCAACTGAATGgatgcagcctccctaaaaaagtggctccat                  >  1:137431‑M1/1‑20 (MQ=255)
                              cGATTTCATCAACTGAATGgatgcagcctccctaaaaaagtggctccatc                 <  1:2621233‑M1/50‑32 (MQ=255)
                              cGATTTCATCAACTGAATGgatgcagcctccctaaaaaagtggctccatc                 >  1:3352945‑M1/1‑19 (MQ=255)
                               gATTTCATCAACTGAATGgatgcagcctccctaaaaaagtggctccatcg                >  1:1596908‑M1/1‑18 (MQ=255)
                                aTTTCATCAACTGAATGgatgcagcctccctaaaaaagtggctccatcgt               <  1:4509949‑M1/50‑34 (MQ=255)
                                aTTTCATCAACTGAATGgatgcagcctccctaaaaaagtggctccatcgt               <  1:3519823‑M1/50‑34 (MQ=255)
                                aTTTCATCAACTGAATGgatgcagcctccctaaaaaagtggctccatcgt               >  1:332698‑M1/1‑17 (MQ=255)
                                aTTTCATCAACTGAATGgatgcagcctccctaaaaaagtggctccatcgt               >  1:1144733‑M1/1‑17 (MQ=255)
                                 tttCATCAACTGAATGgatgcagcctccctaaaaaagtggctccatcgtg              <  1:1394064‑M1/50‑35 (MQ=255)
                                  ttCATCAACTGAATGgatgcagcctccctaaaaaagtggctccatcgtgg             >  1:2346198‑M1/1‑15 (MQ=255)
                                   tcatcaACTGAATGgatgcagcctccctaaaaaagtggctccatcgtggc            >  1:4574439‑M1/1‑14 (MQ=255)
                                   tcatcaACTGAATGgatgcagcctccctaaaaaagtggctccatcgtggc            >  1:931540‑M1/1‑14 (MQ=255)
                                     atcaACTGAATGgatgcagcctccctaaaaaagtggctccatcgtggcgg          <  1:4535743‑M1/50‑39 (MQ=255)
                                     atcaACTGAATGgatgcagcctccctaaaaaagtggctccatcgtggcg           >  1:4530812‑M1/1‑12 (MQ=255)
                                     atcaACTGAATGgatgcagcctccctaaaaaagtggctccatcgtggcg           <  1:2602228‑M1/49‑38 (MQ=255)
                                      tcaACTGAATGgatgcagcctccctaaaaaagtggctccatcgtggcgga         >  1:4075550‑M1/1‑11 (MQ=255)
                                       caACTGAATGgatgcagcctccctaaaaaagtggctccatcgtggcggat        >  1:4791837‑M1/1‑10 (MQ=255)
                                        aaCTGAATGgatgcagcctccctaaaaaagtggctccatcgtggcggatt       <  1:3353759‑M1/50‑42 (MQ=255)
                                        aaCTGAATGgatgcagcctccctaaaaaagtggctccatcgtggcggatt       <  1:538637‑M1/50‑42 (MQ=255)
                                        aaCTGAATGgatgcagcctccctaaaaaagtggctccatcgtggcggat        <  1:2388732‑M1/49‑41 (MQ=255)
                                        aaCTGAATGgatgcagcctccctaaaaaagtggctccatcgtggcggat        <  1:475431‑M1/49‑41 (MQ=255)
                                        aaCTGAATGgatgcagcctccctaaaaaagtggctccatcgtggcggat        >  1:3546861‑M1/1‑9 (MQ=255)
                                        aaCTGAATGgatgcagcctccctaaaaaagtggctccatcgtggcggat        >  1:369535‑M1/1‑9 (MQ=255)
                                         aCTGAATGgatgcagcctccctaaaaaagtggctccatcgtggcggatta      >  1:4479791‑M1/1‑8 (MQ=255)
                                          cTGAATGgatgcagcctccctaaaaaagtggctccatcgtggcggattat     <  1:5359338‑M1/50‑44 (MQ=255)
                                          cTGAATGgatgcagcctccctaaaaaagtggctccatcgtggcggattat     <  1:2030711‑M1/50‑44 (MQ=255)
                                          cTGAATGgatgcagcctccctaaaaaagtggctccatcgtggcggattat     <  1:5010690‑M1/50‑44 (MQ=255)
                                             aaTGgatgcagcctccctaaaaaagtggctccatcgtggcggattattgt  <  1:174210‑M1/50‑47 (MQ=255)
                                             aaTGgatgcagcctccctaaaaaagtggctccatcgtggcggattattgt  >  1:5124589‑M1/1‑4 (MQ=255)
                                             aaTGgatgcagcctccctaaaaaagtggctccatcgtggcggattatt    <  1:2138543‑M1/48‑45 (MQ=255)
                                                |                                              
GACGGGACCTGCCTGATACCCGCCGCCCAGCGATTTCATCAACTGAATGCTTTGGATCACCCGGCGGCTGTCCAGCATCAGTAATGACATCTCTT  >  NC_000913/4299579‑4299673

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑