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 111 [0] [0] 109 [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)106 (0.980) 65/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 109 total aligned reads displayed.

CTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACG  >  NC_000913/4325118‑4325215
                                                |                                                 
acgggacctgcctgatacccgccgcccagcgatttcatcaactgaatgGa                                                  >  1:3056031‑M2/49‑50 (MQ=255)
acgggacctgcctgatacccgccgcccagcgatttcatcaactgaatgGa                                                  >  1:4326806‑M2/49‑50 (MQ=255)
acgggacctgcctgatacccgccgcccagcgatttcatcaactgaatgGa                                                  >  1:80706‑M2/49‑50 (MQ=255)
acgggacctgcctgatacccgccgcccagcgatttcatcaactgaatgGa                                                  <  1:6316673‑M2/2‑1 (MQ=255)
 cgggacctgcctgatacccgccgcccagcgatttcatcaactgaatgGAt                                                 >  1:1082503‑M2/48‑50 (MQ=255)
  gggacctgcctgatacccgccgcccagcgatttcatcaactgaatgGATg                                                <  1:6066324‑M2/4‑1 (MQ=255)
  gggacctgcctgatacccgccgcccagcgatttcatcaactgaatgGATg                                                <  1:1225788‑M2/4‑1 (MQ=255)
  gggacctgcctgatacccgccgcccagcgatttcatcaactgaatgGATg                                                <  1:7155847‑M2/4‑1 (MQ=255)
   ggacctgcctgatacccgccgcccagcgatttcatcaactgaatgGATg                                                >  1:5179763‑M2/46‑49 (MQ=255)
   ggacctgcctgatacccgccgcccagcgatttcatcaactgaatgGATGc                                               <  1:1051239‑M2/5‑1 (MQ=255)
   ggacctgcctgatacccgccgcccagcgatttcatcaactgaatgGATGc                                               <  1:3343370‑M2/5‑1 (MQ=255)
    gacctgcctgatacccgccgcccagcgatttcatcaactgaatgGATGCa                                              <  1:6868024‑M2/6‑1 (MQ=255)
    gacctgcctgatacccgccgcccagcgatttcatcaactgaatgGATGCa                                              <  1:2472194‑M2/6‑1 (MQ=255)
     acctgcctgatacccgccgcccagcgatttcatcaactgaatgGATGCAg                                             <  1:2207171‑M2/7‑1 (MQ=255)
      cctgcctgatacccgccgcccagcgatttcatcaactgaatgGATGCAGc                                            <  1:3943445‑M2/8‑1 (MQ=255)
       ctgcctgatacccgccgcccagcgatttcatcaactgaatgGATGCAGcc                                           <  1:7413750‑M2/9‑1 (MQ=255)
       ctgcctgatacccgccgcccagcgatttcatcaactgaatgGATGCAGcc                                           >  1:8204218‑M2/42‑50 (MQ=255)
       ctgcctgatacccgccgcccagcgatttcatcaactgaatgGATGCAGcc                                           <  1:4041211‑M2/9‑1 (MQ=255)
         gcctgatacccgccgcccagcgatttcatcaactgaatgGATGCAGCCTc                                         >  1:1821694‑M2/40‑50 (MQ=255)
         gcctgatacccgccgcccagcgatttcatcaactgaatgGATGCAGCCTc                                         >  1:5130560‑M2/40‑50 (MQ=255)
          cctgatacccgccgcccagcgatttcatcaactgaatgGATGCAGCCt                                          <  1:1484304‑M2/10‑1 (MQ=255)
          cctgatacccgccgcccagcgatttcatcaactgaatgGATGCAGCCTcc                                        <  1:5735795‑M2/12‑1 (MQ=255)
          cctgatacccgccgcccagcgatttcatcaactgaatgGATGCAGCCTc                                         <  1:7687539‑M2/11‑1 (MQ=255)
           ctgatacccgccgcccagcgatttcatcaactgaatgGATGCAGCCTccc                                       >  1:1083161‑M2/38‑50 (MQ=255)
           ctgatacccgccgcccagcgatttcatcaactgaatgGATGCAGCCTccc                                       >  1:318671‑M2/38‑50 (MQ=255)
            tgatacccgccgcccagcgatttcatcaactgaatgGATGCAGCCTccc                                       >  1:1566568‑M2/37‑49 (MQ=255)
            tgatacccgccgcccagcgatttcatcaaccgaatgGATGCAGCCTccc                                       >  1:3976992‑M2/37‑49 (MQ=255)
              atacccgccgcccagcgatttcatcaactgaatgGATGCAGCCTCCCTa                                     <  1:414728‑M2/15‑1 (MQ=255)
              atacccgccgcccagcgatttcatcaactgaatgGATGCAGCCTCCCTa                                     <  1:2037086‑M2/15‑1 (MQ=255)
                 cccgccgcccagcgatttcatcaactgaatgGATGCAGCCTCCCTaaaa                                  >  1:2433599‑M2/32‑49 (MQ=255)
                  ccgccgcccagcgatttcatcaactgaatgGATGCAGCCTCCCTaaaaaa                                <  1:3762216‑M2/20‑1 (MQ=255)
                   cgccgcccagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAg                               <  1:482246‑M2/21‑1 (MQ=255)
                   cgccgcccagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAg                               >  1:864762‑M2/30‑50 (MQ=255)
                    gccgcccagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGt                              >  1:466303‑M2/29‑50 (MQ=255)
                     ccgcccagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTg                             >  1:2359684‑M2/28‑50 (MQ=255)
                      cgcccagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTgg                            >  1:6989589‑M2/27‑50 (MQ=255)
                      cgcccagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTgg                            <  1:3202597‑M2/24‑1 (MQ=255)
                      cgcccagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTgg                            <  1:4532506‑M2/24‑1 (MQ=255)
                      cgcccagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTgg                            >  1:2537785‑M2/27‑50 (MQ=255)
                        cccagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCt                          <  1:3039127‑M2/26‑1 (MQ=255)
                        cccagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCt                          >  1:1022247‑M2/25‑50 (MQ=255)
                        cccagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCt                          >  1:7640347‑M2/25‑50 (MQ=255)
                        cccagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCt                          <  1:1931427‑M2/26‑1 (MQ=255)
                         ccagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTc                         >  1:3556451‑M2/24‑50 (MQ=255)
                         ccagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTc                         <  1:3002990‑M2/27‑1 (MQ=255)
                         ccagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTc                         >  1:41593‑M2/24‑50 (MQ=255)
                         ccagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTc                         >  1:4014292‑M2/24‑50 (MQ=255)
                         ccagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTc                         >  1:7170055‑M2/24‑50 (MQ=255)
                         ccagcgatttcatcaactgaatgGATGAAGCCTCCCTAAAAAAGGGGCTc                         >  1:5847645‑M2/24‑50 (MQ=255)
                          cagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTcc                        <  1:3369683‑M2/28‑1 (MQ=255)
                          cagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTcc                        <  1:1491872‑M2/28‑1 (MQ=255)
                          cagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTcc                        <  1:2510783‑M2/28‑1 (MQ=255)
                          cagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTcc                        <  1:3627629‑M2/28‑1 (MQ=255)
                           agcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCa                       >  1:6802032‑M2/22‑50 (MQ=255)
                           agcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCa                       >  1:1929770‑M2/22‑50 (MQ=255)
                           agcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCa                       <  1:1746859‑M2/29‑1 (MQ=255)
                            gcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCAt                      <  1:622443‑M2/30‑1 (MQ=255)
                            gcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCAt                      >  1:5980452‑M2/21‑50 (MQ=255)
                             cgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATc                     <  1:2869493‑M2/31‑1 (MQ=255)
                             cgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATc                     >  1:3578282‑M2/20‑50 (MQ=255)
                             cgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATc                     >  1:8299241‑M2/20‑50 (MQ=255)
                             cgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATc                     >  1:5692218‑M2/20‑50 (MQ=255)
                             cgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATc                     <  1:5657696‑M2/31‑1 (MQ=255)
                              gatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCg                    <  1:1575801‑M2/32‑1 (MQ=255)
                              gatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCg                    <  1:8073112‑M2/32‑1 (MQ=255)
                               atttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGt                   >  1:3448598‑M2/18‑50 (MQ=255)
                               atttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGt                   <  1:2539585‑M2/33‑1 (MQ=255)
                               atttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGt                   >  1:5719767‑M2/18‑50 (MQ=255)
                                tttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTg                  <  1:137702‑M2/34‑1 (MQ=255)
                                tttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTg                  >  1:5253689‑M2/17‑50 (MQ=255)
                                  tcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTgg                 <  1:6712824‑M2/35‑1 (MQ=255)
                                   catcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCg               >  1:1765964‑M2/14‑50 (MQ=255)
                                   catcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCg               >  1:4013972‑M2/14‑50 (MQ=255)
                                    atcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCgg              <  1:8305180‑M2/38‑1 (MQ=255)
                                    atcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCg               <  1:2647260‑M2/37‑1 (MQ=255)
                                    atcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCg               <  1:8292999‑M2/37‑1 (MQ=255)
                                     tcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGa             <  1:8101836‑M2/39‑1 (MQ=255)
                                       aactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGatt           <  1:5335248‑M2/41‑1 (MQ=255)
                                       aactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGatt           <  1:7778209‑M2/41‑1 (MQ=255)
                                       aactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGatt           <  1:1137992‑M2/41‑1 (MQ=255)
                                       aactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGatt           <  1:4281507‑M2/41‑1 (MQ=255)
                                       aactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGat            <  1:1226952‑M2/40‑1 (MQ=255)
                                       aactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGat            >  1:5231632‑M2/10‑49 (MQ=255)
                                         ctgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGattat         >  1:1475770‑M2/8‑50 (MQ=255)
                                         ctgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGattat         <  1:2348516‑M2/43‑1 (MQ=255)
                                         ctgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGattat         >  1:5381146‑M2/8‑50 (MQ=255)
                                         ctgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGattat         >  1:507888‑M2/8‑50 (MQ=255)
                                          tgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGattatt        >  1:5565632‑M2/7‑50 (MQ=255)
                                            aatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGattatt        <  1:443663‑M2/44‑1 (MQ=255)
                                            aatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGattatt        <  1:1643091‑M2/44‑1 (MQ=255)
                                            aatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATtgt      <  1:3463736‑M2/46‑1 (MQ=255)
                                            aatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATtgt      >  1:5067641‑M2/5‑50 (MQ=255)
                                            aatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATtg       <  1:1341950‑M2/45‑1 (MQ=255)
                                             atgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATtgtg     <  1:5142218‑M2/47‑1 (MQ=255)
                                             atgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATtgt      >  1:1285578‑M2/4‑49 (MQ=255)
                                             atgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATtgt      >  1:6244422‑M2/4‑49 (MQ=255)
                                               gGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGAc   >  1:1223225‑M2/2‑50 (MQ=255)
                                                gatgCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGAc   <  1:3359677/49‑1 (MQ=255)
                                                gatgCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGAc   <  1:5072154/49‑1 (MQ=255)
                                                gatgCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACg  <  1:3824430/50‑1 (MQ=255)
                                                |                                                 
CTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACG  >  NC_000913/4325118‑4325215

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