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 71 [0] [0] 71 [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)69 (0.870) 51/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

CTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGAC  >  NC_000913/4325118‑4325214
                                                |                                                
acgggacctgcctgatacccgccgcccagcgatttcatcaactgaatgGa                                                 <  1:110873‑M2/2‑1 (MQ=255)
 cgggacctgcctgatacccgccgcccagcgatttcatcaactgaatgGAt                                                >  1:2840302‑M2/48‑50 (MQ=255)
  gggacctgcctgatacccgccgcccagcgatttcatcaactgaatgGATg                                               <  1:3626063‑M2/4‑1 (MQ=255)
   ggacctgcctgatacccgccgcccagcgatttcatcaactgaatgGATg                                               >  1:1552997‑M2/46‑49 (MQ=255)
    gacctgcctgatacccgccgcccagcgatttcatcaactgaatgGATGCa                                             <  1:985051‑M2/6‑1 (MQ=255)
      cctgcctgatacccgccgcccagcgatttcatcaactgaatgGATGCAGc                                           <  1:3172709‑M2/8‑1 (MQ=255)
      cctgcctgatacccgccgcccagcgatttcatcaactgaatgGATGCAGc                                           >  1:219385‑M2/43‑50 (MQ=255)
       ctgcctgatacccgccgcccagcgatttcatcaactgaatgGATGCAGcc                                          >  1:4850748‑M2/42‑50 (MQ=255)
        tgcctgatacccgccgcccagcgatttcatcaactgaatgGATGCAGCCt                                         >  1:3765968‑M2/41‑50 (MQ=255)
          cctgatacccgccgcccagcgatttcatcaactgaatgGATGCAGCCt                                         <  1:2214812‑M2/10‑1 (MQ=255)
          cctgatacccgccgcccagcgatttcatcaactgaatgGATGCAGCCTcc                                       >  1:1589557‑M2/39‑50 (MQ=255)
           ctgatacccgccgcccagcgatttcatcaactgaatgGATGCAGCCTccc                                      <  1:5961415‑M2/13‑1 (MQ=255)
           ctgatacccgccgcccagcgatttcatcaactgaatgGATGCAGCCTccc                                      >  1:3884415‑M2/38‑50 (MQ=255)
              atacccgccgcccagcgatttcatcaactgaatgGATGCAGCCTCCCTa                                    <  1:2883628‑M2/15‑1 (MQ=255)
              atacccgccgcccagcaatttcatcaactgaatgGATGCAGCCTCCCTa                                    >  1:292690‑M2/35‑49 (MQ=255)
               tacccgccgcccagcgatttcatcaactgaatgGATGCAGCCTCCCTaa                                   >  1:5982931‑M2/34‑49 (MQ=255)
                 cccgccgcccagcgatttcatcaactgaatgGATGCAGCCTCCCTaaaa                                 <  1:1700779‑M2/18‑1 (MQ=255)
                  ccgccgcccagcgatttcatcaactgaatgGATGCAGCCTCCCTaaaaaa                               <  1:4330476‑M2/20‑1 (MQ=255)
                   cgccgcccagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAg                              <  1:5845769‑M2/21‑1 (MQ=255)
                   cgccgcccagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAg                              >  1:298990‑M2/30‑50 (MQ=255)
                    gccgcccagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGt                             <  1:5764958‑M2/22‑1 (MQ=255)
                     ccgcccagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTg                            <  1:5666988‑M2/23‑1 (MQ=255)
                     ccgcccagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTg                            >  1:2628284‑M2/28‑50 (MQ=255)
                      cgcccagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTgg                           >  1:2491128‑M2/27‑50 (MQ=255)
                       gcccagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGc                          >  1:3192987‑M2/26‑50 (MQ=255)
                       gcccagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGc                          <  1:1507407‑M2/25‑1 (MQ=255)
                        cccagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCt                         <  1:4816765‑M2/26‑1 (MQ=255)
                         ccagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTc                        >  1:790052‑M2/24‑50 (MQ=255)
                         ccagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTc                        >  1:651835‑M2/24‑50 (MQ=255)
                         ccagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTc                        <  1:4511954‑M2/27‑1 (MQ=255)
                          cagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTcc                       <  1:5991002‑M2/28‑1 (MQ=255)
                           agcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCa                      <  1:4305608‑M2/29‑1 (MQ=255)
                           agcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCa                      <  1:1361719‑M2/29‑1 (MQ=255)
                           agcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCa                      <  1:1958839‑M2/29‑1 (MQ=255)
                            gcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCAt                     <  1:5005129‑M2/30‑1 (MQ=255)
                            gcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCAt                     >  1:1757997‑M2/21‑50 (MQ=255)
                             cgatttcatcaactgaatgGATGCGGCCTCCCTAAAAAAGTGGCTCCATc                    >  1:3945924‑M2/20‑50 (MQ=255)
                             cgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATc                    >  1:3693591‑M2/20‑50 (MQ=255)
                             cgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATc                    <  1:3975774‑M2/31‑1 (MQ=255)
                             cgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATc                    >  1:4616194‑M2/20‑50 (MQ=255)
                             cgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATc                    <  1:4603993‑M2/31‑1 (MQ=255)
                              gatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCg                   >  1:4192134‑M2/19‑50 (MQ=255)
                              gatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCg                   >  1:362064‑M2/19‑50 (MQ=255)
                               atttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGt                  >  1:2695829‑M2/18‑50 (MQ=255)
                               atttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGt                  >  1:5317006‑M2/18‑50 (MQ=255)
                               atttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGt                  >  1:5015672‑M2/18‑50 (MQ=255)
                                tttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTg                 <  1:1423064‑M2/34‑1 (MQ=255)
                                tttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTg                 >  1:4335527‑M2/17‑50 (MQ=255)
                                tttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTg                 >  1:1539451‑M2/17‑50 (MQ=255)
                                  tcatcaactggatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTgg                <  1:3357204‑M2/35‑1 (MQ=255)
                                  tcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTgg                <  1:1155951‑M2/35‑1 (MQ=255)
                                  tcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTgg                <  1:3046987‑M2/35‑1 (MQ=255)
                                   catcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCg              >  1:6019060‑M2/14‑50 (MQ=255)
                                   catcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCg              >  1:288912‑M2/14‑50 (MQ=255)
                                   catcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCg              >  1:4341674‑M2/14‑50 (MQ=255)
                                    atcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCg              <  1:5102327‑M2/37‑1 (MQ=255)
                                    atcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCg              <  1:550698‑M2/37‑1 (MQ=255)
                                    atcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCg              >  1:3308871‑M2/13‑49 (MQ=255)
                                      caaccgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGat           >  1:1518735‑M2/11‑50 (MQ=255)
                                       aactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGatt          <  1:1764849‑M2/41‑1 (MQ=255)
                                       aactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGatt          <  1:1236944‑M2/41‑1 (MQ=255)
                                       aactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGatt          <  1:4273583‑M2/41‑1 (MQ=255)
                                       aactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGat           <  1:1554528‑M2/40‑1 (MQ=255)
                                         ctgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGattat        <  1:5907864‑M2/43‑1 (MQ=255)
                                           gaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATtg      >  1:1109441‑M2/6‑50 (MQ=255)
                                            aatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATtgt     <  1:5404919‑M2/46‑1 (MQ=255)
                                             atgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATtgt     >  1:2026277‑M2/4‑49 (MQ=255)
                                               gGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGAc  >  1:401545‑M2/2‑50 (MQ=255)
                                               gGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGAc  >  1:2959705‑M2/2‑50 (MQ=255)
                                                gatgCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGAc  <  1:344108/49‑1 (MQ=255)
                                                gatgCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGAc  <  1:4851086/49‑1 (MQ=255)
                                                |                                                
CTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGAC  >  NC_000913/4325118‑4325214

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: