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 98 [0] [0] 101 [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)92 (1.040) 57/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 101 total aligned reads displayed.

CTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACG  >  NC_000913/4325118‑4325215
                                                |                                                 
acgggacctgcctgatacccgccgcccagcgatttcatcaactgaatgGa                                                  <  1:2999345‑M2/2‑1 (MQ=255)
acgggacctgcctgatacccgccgcccagcgatttcatcaactgaatgGa                                                  <  1:2571315‑M2/2‑1 (MQ=255)
acgggacctgcctgatacccgccgcccagcgatttcatcaactgaatgGa                                                  >  1:6401824‑M2/49‑50 (MQ=255)
acgggacctgcctgatacccgccgcccagcgatttcatcaactgaatgGa                                                  <  1:6162423‑M2/2‑1 (MQ=255)
 cgggacctgcctgatacccgccgcccagcgatttcatcaactgaatgGAt                                                 >  1:4874126‑M2/48‑50 (MQ=255)
 cgggacctgcctgatacccgccgcccagcgatttcatcaactgaatgGAt                                                 <  1:5602077‑M2/3‑1 (MQ=255)
 cgggacctgcctgatacccgccgcccagcgatttcatcaactgaatgGAt                                                 <  1:1092374‑M2/3‑1 (MQ=255)
 cgggacctgcctgatacccgccgcccagcgatttcatcaactgaatgGAt                                                 >  1:3515659‑M2/48‑50 (MQ=255)
  gggacctgcctgatacccgccgcccagcgatttcatcaactgaatgGATg                                                >  1:5200123‑M2/47‑50 (MQ=255)
  gggacctgcctgatacccgccgcccagcgatttcatcaactgaatgGATg                                                >  1:3952359‑M2/47‑50 (MQ=255)
  gggacctgcctgatacccgccgcccagcgatttcatcaactgaatgGATg                                                >  1:1284087‑M2/47‑50 (MQ=255)
  gggacctgcctgatacccgccgcccagcgatttcatcaactgaatgGATg                                                >  1:268173‑M2/47‑50 (MQ=255)
  gggacctgcctgatacccgccgcccagcgatttcatcaactgaatgGATg                                                >  1:1003943‑M2/47‑50 (MQ=255)
   ggacctgcctgatacccgccgcccagcgatttcatcaactgaatgGATg                                                >  1:6085369‑M2/46‑49 (MQ=255)
    gacctgcctgatacccgccgcccagcgatttcatcaactgaatgGATGCa                                              <  1:3894473‑M2/6‑1 (MQ=255)
    gacctgcctgatacccgccgcccagcgatttcatcaactgaatgGATGCa                                              <  1:5761577‑M2/6‑1 (MQ=255)
    gacctgcctgatacccgccgcccagcgatttcatcaactgaatgGATGCa                                              <  1:4790971‑M2/6‑1 (MQ=255)
     acctgcctgatacccgccgcccagcgatttcatcaactgaatgGATGCAg                                             <  1:1162922‑M2/7‑1 (MQ=255)
     acctgcctgatacccgccgcccagcgatttcatcaactgaatgGATGCAg                                             <  1:2670762‑M2/7‑1 (MQ=255)
      cctgcctgatacccgccgcccagcgatttcatcaactgaatgGATGCAGc                                            >  1:5358579‑M2/43‑50 (MQ=255)
      cctgcctgatacccgccgcccagcgatttcatcaactgaatgGATGCAGc                                            >  1:6414258‑M2/43‑50 (MQ=255)
      cctgcctgatacccgccgcccagcgatttcatcaactgaatgGATGCAGc                                            <  1:1557857‑M2/8‑1 (MQ=255)
       ctgcctgatacccgccgcccagcgatttcatcaactgaatgGATGCAGcc                                           >  1:2544371‑M2/42‑50 (MQ=255)
       ctgcctgatacccgccgcccagcgatttcatcaactgaatgGATGCAGcc                                           >  1:2795906‑M2/42‑50 (MQ=255)
        tgcctgatacccgccgcccagcgatttcatcaactgaatgGATGCAGCCt                                          >  1:5947143‑M2/41‑50 (MQ=255)
          cctgatacccgccgcccagcgatttcatcaactgaatgGATGCAGCCt                                          <  1:2333152‑M2/10‑1 (MQ=255)
          cctgatacccgccgcccagcgatttcatcaactgaatgGATGCAGCCTcc                                        <  1:3776695‑M2/12‑1 (MQ=255)
           ctgatacccgccgcccagcgatttcatcaactgaatgGATGCAGCCTccc                                       <  1:4053489‑M2/13‑1 (MQ=255)
           ctgatacccgccgcccagcgatttcatcaactgaatgGATGCAGCCTccc                                       <  1:5202298‑M2/13‑1 (MQ=255)
            tgatacccgccgcccagcgatttcatcaactgaatgGATGCAGCCTccc                                       >  1:3414778‑M2/37‑49 (MQ=255)
              atacccgccgcccagcgatttcatcaactgaatgGATGCAGCCTCCCt                                      <  1:1267594‑M2/14‑1 (MQ=255)
              atacccgccgcccagcgatttcatcaactgaatgGATGCAGCCTCCCt                                      <  1:6425530‑M2/14‑1 (MQ=255)
               tgcccgccgcccagcgatttcatcaactgaatgGATGCAGCCTCCCTaaa                                   <  1:2892120‑M2/17‑1 (MQ=255)
                acccgccgcccagcgatttcatcaactgaatgGATGCAGCCTCCCTaaa                                   >  1:6556901‑M2/33‑49 (MQ=255)
                 cccgccgcccagcgatttcatcaactgaatgGATGCAGCCTCCCTaaaa                                  >  1:1309535‑M2/32‑49 (MQ=255)
                   cgccgcccagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAg                               >  1:3495689‑M2/30‑50 (MQ=255)
                   cgccgcccagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAg                               >  1:5366426‑M2/30‑50 (MQ=255)
                   cgccgcccagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAg                               <  1:755375‑M2/21‑1 (MQ=255)
                    gccgcccagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGt                              >  1:1004815‑M2/29‑50 (MQ=255)
                     ccgcccagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTg                             <  1:813345‑M2/23‑1 (MQ=255)
                     ccgcccagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTg                             >  1:121284‑M2/28‑50 (MQ=255)
                      cgcccagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTgg                            >  1:3996741‑M2/27‑50 (MQ=255)
                      cgcccagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTgg                            <  1:6392150‑M2/24‑1 (MQ=255)
                       gcccagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGc                           >  1:5945847‑M2/26‑50 (MQ=255)
                       cgcccagcgattcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGc                           <  1:2728295‑M2/25‑1 (MQ=255)
                        cccagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCt                          <  1:1026987‑M2/26‑1 (MQ=255)
                         ccagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTc                         >  1:278300‑M2/24‑50 (MQ=255)
                         ccagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTc                         <  1:2498944‑M2/27‑1 (MQ=255)
                          cagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTcc                        <  1:975542‑M2/28‑1 (MQ=255)
                           agcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCa                       >  1:3263240‑M2/22‑50 (MQ=255)
                           agcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCa                       >  1:1447579‑M2/22‑50 (MQ=255)
                           agcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCa                       <  1:1633075‑M2/29‑1 (MQ=255)
                            gcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCAt                      <  1:1179171‑M2/30‑1 (MQ=255)
                            gcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCAt                      <  1:4769400‑M2/30‑1 (MQ=255)
                             cgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATc                     <  1:5741045‑M2/31‑1 (MQ=255)
                              gatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCg                    >  1:4128735‑M2/19‑50 (MQ=255)
                              gatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCg                    >  1:209540‑M2/19‑50 (MQ=255)
                              gatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCg                    >  1:4474289‑M2/19‑50 (MQ=255)
                               atttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGt                   >  1:1022975‑M2/18‑50 (MQ=255)
                               atttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGt                   <  1:2352375‑M2/33‑1 (MQ=255)
                               atttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGCGGCTCCATCGt                   >  1:1511445‑M2/18‑50 (MQ=255)
                               atttcagcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGt                   <  1:5870541‑M2/33‑1 (MQ=255)
                                tttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTg                  <  1:2822891‑M2/34‑1 (MQ=255)
                                tttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTg                  <  1:4217850‑M2/34‑1 (MQ=255)
                                 ttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTgg                 >  1:6726041‑M2/16‑50 (MQ=255)
                                  tcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTgg                 <  1:217‑M2/35‑1 (MQ=255)
                                  tcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGc                >  1:2260155‑M2/15‑50 (MQ=255)
                                   catcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGc                <  1:4547256‑M2/36‑1 (MQ=255)
                                    atcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCg               >  1:4700725‑M2/13‑49 (MQ=255)
                                    atcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCg               <  1:2560945‑M2/37‑1 (MQ=255)
                                    atcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCg               >  1:866367‑M2/13‑49 (MQ=255)
                                    atcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCg               <  1:65212‑M2/37‑1 (MQ=255)
                                     tcaactgaatgGATGCAGCCTCCCTAAAAAAGAGGCTCCATCGTGGCGGa             >  1:5369139‑M2/12‑50 (MQ=255)
                                      caactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGat            >  1:4570810‑M2/11‑50 (MQ=255)
                                         ctgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGattat         <  1:2716038‑M2/43‑1 (MQ=255)
                                         ctgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGattat         <  1:2633611‑M2/43‑1 (MQ=255)
                                         ctgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGattat         <  1:4853031‑M2/43‑1 (MQ=255)
                                         ctgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGattat         <  1:2556494‑M2/43‑1 (MQ=255)
                                           gaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATtg       >  1:3375386‑M2/6‑50 (MQ=255)
                                            aatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGattatt        <  1:5844582‑M2/44‑1 (MQ=255)
                                            aatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGattatt        <  1:1554100‑M2/44‑1 (MQ=255)
                                            aatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATtgt      <  1:3951099‑M2/46‑1 (MQ=255)
                                            aatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATtgt      <  1:5726519‑M2/46‑1 (MQ=255)
                                            aatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATtg       <  1:1078216‑M2/45‑1 (MQ=255)
                                            aatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATtg       <  1:5089717‑M2/45‑1 (MQ=255)
                                             atgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATtgtg     <  1:1856546‑M2/47‑1 (MQ=255)
                                             atgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATtgtg     <  1:4184128‑M2/47‑1 (MQ=255)
                                             atgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATtgtg     <  1:96246‑M2/47‑1 (MQ=255)
                                             atgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATtgt      >  1:3717297‑M2/4‑49 (MQ=255)
                                              tgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGa    >  1:1335328‑M2/3‑50 (MQ=255)
                                               gGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGAc   >  1:1096160‑M2/2‑50 (MQ=255)
                                               gGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGAc   >  1:2352214‑M2/2‑50 (MQ=255)
                                                gatgCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGAc   >  1:1609641/1‑49 (MQ=255)
                                                gatgCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGAc   <  1:2443089/49‑1 (MQ=255)
                                                gatgCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGAc   <  1:3072903/49‑1 (MQ=255)
                                                gatgCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGAc   >  1:289394/1‑49 (MQ=255)
                                                gatgCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGAc   <  1:973385/49‑1 (MQ=255)
                                                gatgCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGAc   >  1:2837863/1‑49 (MQ=255)
                                                gatgCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACg  <  1:5593538/50‑1 (MQ=255)
                                                gatgCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACg  <  1:3094249/50‑1 (MQ=255)
                                                |                                                 
CTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACG  >  NC_000913/4325118‑4325215

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