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 149 [0] [0] 146 [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)137 (0.920) 72/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

Only 100 of 149 total aligned reads displayed.

CGACGGGACCTGCCTGATACCCGCCGCCCAGCGATTTCATCAACTGAATGCTTTGGATCACCCGGCGGCTGTCCAGCATCAGTAATGACATCTCTTCGG  >  NC_000913/4299578‑4299676
                                                 |                                                 
cgacgGGACCTGCCTGATACCCGCCGCCCAGCGATTTCATCAACTGAATg                                                   >  2:1373746/1‑50 (MQ=255)
cgacgGGACCTGCCTGATACCCGCCGCCCAGCGATTTCATCAACTGAATg                                                   <  2:709785/50‑1 (MQ=255)
cgacgGGACCTGCCTGATACCCGCCGCCCAGCGATTTCATCAACTGAATg                                                   >  1:133839/1‑50 (MQ=255)
cgacgGGACCTGCCTGATACCCGCCGCCCAGCGATTTCATCAACTGAATg                                                   >  1:2585833/1‑50 (MQ=255)
cgacgGGACCTGCCTGATACCCGCCGCCCAGCGATTTCATCAACTGAATg                                                   >  2:4450117/1‑50 (MQ=255)
cgacgGGACCTGCCTGATACCCGCCGCCCAGCGATTTCATCAACTGAATg                                                   >  1:4366130/1‑50 (MQ=255)
 gacgGGACCTGCCTGATACCCGCCGCCCAGCGATTTCATCAACTGAATg                                                   >  2:5643048/1‑49 (MQ=255)
 gacgGGACCTGCCTGATACCCGCCGCCCAGCGATTTCATCAACTGAATg                                                   >  2:3577832/1‑49 (MQ=255)
 gacgGGACCTGCCTGATACCCGCCGCCCAGCGATTTCATCAACTGAATg                                                   >  2:4180069/1‑49 (MQ=255)
  aCGGGACCTGCCTGATACCCGCCGCCCAGCGATTTCATCAACTGAATGga                                                 >  2:3819011‑M1/1‑48 (MQ=255)
   cGGGACCTGCCTGATACCCGCCGCCCAGCGATTTCATCAACTGAATGgat                                                <  1:757531‑M1/50‑4 (MQ=255)
    gggACCTGCCTGATACCCGCCGCCCAGCGATTTCATCAACTGAATGgatg                                               >  1:3671464‑M1/1‑46 (MQ=255)
    gggACCTGCCTGATACCCGCCGCCCAGCGATTTCATCAACTGAATGgatg                                               >  1:3382495‑M1/1‑46 (MQ=255)
    gggACCTGCCTGATACCCGCCGCCCAGCGATTTCATCAACTGAATGgatg                                               <  2:4024859‑M1/50‑5 (MQ=255)
    gggACCTGCCTGATACCCGCCGCCCAGCGATTTCATCAACTGAATGgatg                                               >  1:4934119‑M1/1‑46 (MQ=255)
    gggACCTGCCTGATACCCGCCGCCCAGCGATTTCATCAACTGAATGgatg                                               >  2:1888266‑M1/1‑46 (MQ=255)
     ggACCTGCCTGATACCCGCCGCCCAGCGATTTCATCAACTGAATGgatgc                                              <  2:1958110‑M1/50‑6 (MQ=255)
     ggACCTGCCTGATACCCGCCGCCCAGCGATTTCATCAACTGAATGgatgc                                              <  1:4159378‑M1/50‑6 (MQ=255)
      gACCTGCCTGATACCCGCCGCCCAGCGCTTTCATCCACTGAATGgttg                                               >  1:1882272‑M1/1‑44 (MQ=255)
        cctgcctgATACCCGCCGCCCAGCGATTTCATGAACTGAATGgatgcagc                                           <  1:5060060‑M1/50‑9 (MQ=255)
         ctgcctgATACCCGCCGCCCAGCGATTTCATCAACTGAATGgatgcagcc                                          >  2:178116‑M1/1‑41 (MQ=255)
         ctgcctgATACCCGCCGCCCAGCGATTTCATCAACTGAATGgatgcagcc                                          <  1:3910975‑M1/50‑10 (MQ=255)
          tgcctgATACCCGCCGCCCAGCGATTTCATCAACTGAATGgatgcagcct                                         >  2:469253‑M1/1‑40 (MQ=255)
           gcctgATACCCGCCGCCCAGCGATTTCATCAACTGAATGgatgcagcctc                                        >  1:3967395‑M1/1‑39 (MQ=255)
            cctgATACCCGCCGCCCAGCGATTTCATCAACTGAATGgatgcagcctcc                                       <  2:3295498‑M1/50‑13 (MQ=255)
            cctgATACCCGCCGCCCAGCGATTTCATCAACTGAATGgatgcagcctcc                                       <  1:2261932‑M1/50‑13 (MQ=255)
            cctgATACCCGCCGCCCAGCGATTTCATCAACTGAATGgatgcagcctcc                                       <  1:1183876‑M1/50‑13 (MQ=255)
            cctgATACCCGCCGCCCAGCGATTTCATCAACTGAATGgatgcagcct                                         <  2:2391823‑M1/48‑11 (MQ=255)
             ctgATACCCGCCGCCCAGCGATTTCATCAACTGAATGgatgcagcctccc                                      >  1:4768638‑M1/1‑37 (MQ=255)
              tgATACCCGCCGCCCAGCGATTTCATCAACTGAATGgatgcagcctccc                                      >  1:2173402‑M1/1‑36 (MQ=255)
               gATACCCGCCGCCCAGCGATTTCATCAACTGAATGgatgcagcctccct                                     >  2:2684509‑M1/1‑35 (MQ=255)
               gATACCCGCCGCCCAGCGATTTCATCAACTGAATGgatgcagcctccct                                     >  1:2542396‑M1/1‑35 (MQ=255)
                aTACCCGCCGCCCAGCGATTTCATCAACTGAATGgatgcagcctccctaa                                   <  1:317575‑M1/50‑17 (MQ=255)
                aTACCCGCCGCCCAGCGATTTCATCAACTGAATGgatgcagcctcccta                                    <  2:1123084‑M1/49‑16 (MQ=255)
                aTACCCGCCGCCCAGCGATTTCATCAACTGAATGgatgcagcctcccta                                    <  1:796862‑M1/49‑16 (MQ=255)
                aTACCCGCCGCCCAGCGATTTCATCAACTGAATGgatgcagcctccct                                     <  2:3554366‑M1/48‑15 (MQ=255)
                   cccGCCGCCCAGCGATTTCATCAACTGAATGgatgcagcctccctaaaaa                                <  2:2506439‑M1/50‑20 (MQ=255)
                   cccGCCGCCCAGCGATTTCATCAACTGAATGgatgcagcctccctaaaa                                 <  2:3216371‑M1/49‑19 (MQ=255)
                   cccGCCGCCCAGCGATTTCATCAACTGAATGgatgcagcctccctaaaa                                 <  1:1114592‑M1/49‑19 (MQ=255)
                    ccgccgCCCAGCGATTTCATCAACTGAATGgatgcagcctccctaaaaaa                               <  1:4546987‑M1/50‑21 (MQ=255)
                     cgccgcCCAGCGATTTCATCAACTGAATGgatgcagcctccctaaaaaag                              >  1:4413347‑M1/1‑29 (MQ=255)
                     cgccgcCCAGCGATTTCATCAACTGAATGgatgcagcctccctaaaaaag                              >  2:5128163‑M1/1‑29 (MQ=255)
                     cgccgcCCAGCGATTTCATCAACTGAATGgatgcagcctccctaaaaaag                              >  2:5158995‑M1/1‑29 (MQ=255)
                      gccgccCAGCGATTTCATCAACTGAATGgatgcagcctccctaaaaaagt                             <  1:5551623‑M1/50‑23 (MQ=255)
                      gccgccCAGCGATTTCATCAACTGAATGgatgcagcctccctaaaaaagt                             >  1:3813372‑M1/1‑28 (MQ=255)
                       ccgccCAGCGATTTCATCAACTGAATGgatgcagcctccctaaaaaagtg                            >  2:3323618‑M1/1‑27 (MQ=255)
                        cgccCAGCGATTTCATCAACTGAATGgatgcagcctccctaaaaaagtgg                           >  2:4736709‑M1/1‑26 (MQ=255)
                         gccCAGCGATTTCATCAACTGAATGgatgcagcctccctaaaaaagtggc                          <  1:2627176‑M1/50‑26 (MQ=255)
                         gccCAGCGATTTCATCAACTGAATGgatgcagcctccctaaaaaagtggc                          <  2:1153524‑M1/50‑26 (MQ=255)
                          cccAGCGATTTCATCAACTGAATGgatgcagcctccctaaaaaagtggct                         <  2:210349‑M1/50‑27 (MQ=255)
                          cccAGCGATTTCATCAACTGAATGgatgcagcctccctaaaaaagtggct                         >  2:5223316‑M1/1‑24 (MQ=255)
                            cAGCGATTTCATCAACTGAATGgatgcagcctccctaaaaaagtggctc                        >  1:1293682‑M1/1‑22 (MQ=255)
                            cAGCGATTTCATCAACTGAATGgatgcagcctccctaaaaaagtggctc                        >  2:3432931‑M1/1‑22 (MQ=255)
                             aGCGATTTCATCAACTGAATGgatgcagcctccctaaaaaagtggctcca                      <  2:157254‑M1/50‑30 (MQ=255)
                              gCGATTTCATCAACTGAATGgatgcagcctccctaaaaaagtggctccat                     <  1:932899‑M1/50‑31 (MQ=255)
                              gCGATTTAATCAACTGAATGgatgcagcctccctaacaaagtggctccat                     >  2:5505045‑M1/1‑20 (MQ=255)
                               cGATTTCATCAACTGAATGgatgcagcctccctaaaaaagtggctccatc                    >  1:4776645‑M1/1‑19 (MQ=255)
                               cGATTTCATCAACTGAATGgatgcagcctccctaaaaaagtggctccatc                    <  1:3635577‑M1/50‑32 (MQ=255)
                               cGATTTCATCAACTGAATGgatgcagcctccctaaaaaagtggctccatc                    <  1:1217118‑M1/50‑32 (MQ=255)
                               cGATTTCATCAACTGAATGgatgcagcctccctaaaaaagtggctccatc                    <  1:2040522‑M1/50‑32 (MQ=255)
                                gATTTCATCAACTGAATGgatgcagcctccctaaaaaagtggctccatcg                   >  1:466915‑M1/1‑18 (MQ=255)
                                gATTTCATCAACTGAATGgatgcagcctccctaaaaaagtggctccatcg                   >  2:5510742‑M1/1‑18 (MQ=255)
                                 aTTTCATCAACTGAATGgatgcagcctccctaaaaaagtggctccatcgt                  >  1:5555508‑M1/1‑17 (MQ=255)
                                 aTTTCATCAACTGAATGgatgcagcctccctaaaaaagtggctccatcgt                  <  1:3142099‑M1/50‑34 (MQ=255)
                                 aTTTCATCAACTGAATGgatgcagcctccctaaaaaagtggctccatc                    >  1:519025‑M1/1‑17 (MQ=255)
                                  tttCATCAACTGAATGgatgcagcctccctaaaaaagtggctccatcgtg                 <  2:2131123‑M1/50‑35 (MQ=255)
                                  tttCATCAACTGAATGgatgcagcctccctaaaaaagtggctccatcgtg                 >  1:4114869‑M1/1‑16 (MQ=255)
                                  tttCATCAACTGAATGgatgcagcctccctaaaaaagtggctccatcgtg                 >  1:1286791‑M1/1‑16 (MQ=255)
                                  tttCATCAACTGAATGgatgcagcctccctaaaaaagtggctccatcgtg                 >  1:3237632‑M1/1‑16 (MQ=255)
                                    tcatcaACTGAATGgatgcagcctccctaaaaaagtggctccatcgtggc               >  2:493942‑M1/1‑14 (MQ=255)
                                      atcaACTGAATGgatgcagcctccctaaaaaagtggctccatcgtggcg              <  1:2786177‑M1/49‑38 (MQ=255)
                                       tcaACTGAATGgatgcagcctccctaaaaaagtggctccatcgtggcgga            <  2:2743326‑M1/50‑40 (MQ=255)
                                       tcaACTGAATGgatgcagcctccctaaaaaagtggctccatcgtggcgga            >  1:2788011‑M1/1‑11 (MQ=255)
                                        caACTGAATGgatgcagcctccctaaaaaagtggctccatcgtggcggat           >  1:2312511‑M1/1‑10 (MQ=255)
                                        caACTGAATGgatgcagcctccctaaaaaagtggctccatcgtggcggat           >  2:3739221‑M1/1‑10 (MQ=255)
                                         aaCTGAATGgatgcagcctccctaaaaaagtggctccatcgtggcggatt          <  1:3830182‑M1/50‑42 (MQ=255)
                                         aaCTGAATGgatgcagcctccctaaaaaagtggctccatcgtggcggatt          <  1:156428‑M1/50‑42 (MQ=255)
                                         aaCTGAATGgatgcagcctccctaaaaaagtggctccatcgtggcggatt          <  1:3600872‑M1/50‑42 (MQ=255)
                                         aaCTGAATGgatgcagcctccctaaaaaagtggctccatcgtggcggat           >  1:3399338‑M1/1‑9 (MQ=255)
                                           cTGAATGgatgcagcctccctaaaaaagtggctccatcgtggcggattat        <  1:2525173‑M1/50‑44 (MQ=255)
                                           cTGAATGgatgcagcctccctaaaaaagtggctccatcgtggcggattat        >  2:3607402‑M1/1‑7 (MQ=255)
                                           cTGAATGgatgcagcctccctaaaaaagtggctccatcgtggcggattat        <  1:466639‑M1/50‑44 (MQ=255)
                                             gAATGgatgcagcctccctaaaaaagtggctccatcgtggcggattattg      >  1:3666729‑M1/1‑5 (MQ=255)
                                              aaTGgatgcagcctccctaaaaaagtggctccatcgtggcggattattgt     >  2:5617483‑M1/1‑4 (MQ=255)
                                              aaTGgatgcagcctccctaaaaaagtggctccatcgtggcggattattgt     <  1:5408506‑M1/50‑47 (MQ=255)
                                              aaTGgatgcagcctccctaaaaaagtggctccatcgtggcggattattg      <  2:953653‑M1/49‑46 (MQ=255)
                                              aaTGgatgcagcctccctaaaaaagtggctccatcgtggcggattattg      <  1:2250251‑M1/49‑46 (MQ=255)
                                              aaTGgatgcagcctccctaaaaaagtggctccatcgtggcggattatt       <  2:5019453‑M1/48‑45 (MQ=255)
                                              aaTGgatgcagcctccctaaaaaagtggctccatcgtggcggattatt       <  1:4922286‑M1/48‑45 (MQ=255)
                                              aaTGgatgcagccgccctaaaaaggtgggtccttcgtggcggattattgt     <  2:5199111‑M1/50‑47 (MQ=255)
                                               aTGgatgcagcctccctaaaaaagtggctccatcgtggcggattattgtg    <  1:1590332‑M1/50‑48 (MQ=255)
                                               aTGgatgcagcctccctaaaaaagtggctccatcgtggcggattattgtg    <  2:3211160‑M1/50‑48 (MQ=255)
                                               aTGgatgcagcctccctaaaaaagtggctccatcgtggcggattattgtg    <  1:926576‑M1/50‑48 (MQ=255)
                                                tGgatgcagcctccctaaaaaagtggctccatcgtggcggattattgtga   <  2:4539474‑M1/50‑49 (MQ=255)
                                                tGgatgcagcctccctaaaaaagtggctccatcgtggcggattattgtga   <  2:3406899‑M1/50‑49 (MQ=255)
                                                tGgatgcagcctccctaaaaaagtggctccatcgtggcggattattgtga   <  1:3103667‑M1/50‑49 (MQ=255)
                                                tGgatgcagcctccctaaaaaagtggctccatcgtggcggattattgtga   <  2:1967763‑M1/50‑49 (MQ=255)
                                                 ggatgcagcctccctaaaaaagtggctccatcgtggcggattattgtgac  >  1:1379419‑M1/1‑1 (MQ=255)
                                                 ggatgcagcctccctaaaaaagtggctccatcgtggcggattattgtgac  >  2:3845375‑M1/1‑1 (MQ=255)
                                                 ggatgcagcctccctaaaaaagtggctccatcgtggcggattattgtgac  >  1:4290424‑M1/1‑1 (MQ=255)
                                                 |                                                 
CGACGGGACCTGCCTGATACCCGCCGCCCAGCGATTTCATCAACTGAATGCTTTGGATCACCCGGCGGCTGTCCAGCATCAGTAATGACATCTCTTCGG  >  NC_000913/4299578‑4299676

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