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 97 [0] [0] 90 [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)88 (0.770) 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

TGATTGAAGGTTTTGGCGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGAC  >  NC_000913/4325119‑4325214
                                               |                                                
cgggacctgcctgatacccgccgcccagcgatttcatcaactgaatgGAt                                                <  1:7274003‑M2/3‑1 (MQ=255)
cgggacctgcctgatacccgccgcccagcgatttcatcaactgaatgGAt                                                <  1:5355726‑M2/3‑1 (MQ=255)
cgggacctgcctgatacccgccgcccagcgatttcatcaactgaatgGAt                                                >  1:5206388‑M2/48‑50 (MQ=255)
cgggacctgcctgatacccgccgcccagcgatttcatcaactgaatgGAt                                                <  1:4972689‑M2/3‑1 (MQ=255)
cgggacctgcctgatacccgccgcccagcgatttcatcaactgaatgGAt                                                <  1:335734‑M2/3‑1 (MQ=255)
 tggacctgcctgatacccgccgcccagcgatttcatcaactgaatgGATg                                               >  1:5272392‑M2/47‑50 (MQ=255)
 gggacctgcctgatacccgccgcccagcgatttcatcaactgaatgGATg                                               >  1:958884‑M2/47‑50 (MQ=255)
 gggacctgcctgatacccgccgcccagcgatttcatcaactgaatgGATg                                               <  1:6941357‑M2/4‑1 (MQ=255)
 gggacctgcctgatacccgccgcccagcgatttcatcaactgaatgGATg                                               <  1:4722441‑M2/4‑1 (MQ=255)
 gggacctgcctgatacccgccgcccagcgatttcatcaactgaatgGATg                                               >  1:1699092‑M2/47‑50 (MQ=255)
   gacctgcctgatacccgccgcccagcgatttcatcaactgaatgGATg                                               >  1:26018‑M2/45‑48 (MQ=255)
   gacctgcctgatacccgccgcccagcgatttcatcaactgaatgGATGCa                                             <  1:4238178‑M2/6‑1 (MQ=255)
    acctgcctgatacccgccgcccagcgatttcatcaactgaatgGATGCAg                                            >  1:2029733‑M2/44‑50 (MQ=255)
    acctgcctgatacccgccgcccagcgatttcatcaactgaatgGATGCAg                                            >  1:3751760‑M2/44‑50 (MQ=255)
    acctgcctgatacccgccgcccagcgatttcatcaactgaatgGATGCAg                                            <  1:6209291‑M2/7‑1 (MQ=255)
     cctgcctgatacccgccgcccagcgatttcatcaactgaatgGATGCAGc                                           <  1:3085588‑M2/8‑1 (MQ=255)
      ctgcctgatacccgccgcccagcgatttcatcaactgaatgGATGCAGcc                                          <  1:6660078‑M2/9‑1 (MQ=255)
      ctgcctgatacccgccgcccagcgatttcatcaactgaatgGATGCAGcc                                          <  1:8333786‑M2/9‑1 (MQ=255)
       tgcctgatacccgccgcccagcgatttcatcaactgaatgGATGCAGCCt                                         >  1:7635255‑M2/41‑50 (MQ=255)
       tgcctgatacccgccgcccagcgatttcatcaactgaatgGATGCAGCCt                                         >  1:7753509‑M2/41‑50 (MQ=255)
        gcctgatacccgccgcccagcgatttcatcaactgaatgGATGCAGCCTc                                        >  1:2612064‑M2/40‑50 (MQ=255)
        gcctgatacccgccgcccagcgatttcatcaactgaatgGATGCAGCCTc                                        >  1:4096470‑M2/40‑50 (MQ=255)
         cctgatacccgccgcccagcgatttcatcaactgaatgGATGCAGCCt                                         <  1:5323668‑M2/10‑1 (MQ=255)
         cctgatacccgccgcccagcgatttcatcaactgaatgGATGCAGCCTcc                                       >  1:917186‑M2/39‑50 (MQ=255)
         cctgatacccgccgcccagcgatttcatcaactgaatgGATGCAGCCTcc                                       <  1:4107534‑M2/12‑1 (MQ=255)
         cctgatacccgccgcccagcgatttcatcaactgaatgGATGCAGCCTcc                                       >  1:1535080‑M2/39‑50 (MQ=255)
          ctgatacccgccgcccagcgatttcatcaactgaatgGATGCAGCCTccc                                      >  1:3616214‑M2/38‑50 (MQ=255)
           tgatacccgccgcccagcgatttcatcaactgaatgGATGCAGCCTccc                                      >  1:2353245‑M2/37‑49 (MQ=255)
            gatacccgccgcccagcgatttcatcaactgaatgGATGCAGCCTCCCt                                     >  1:6743030‑M2/36‑49 (MQ=255)
             atacccgccgcccagcgatttcatcaactgaatgGATGCAGCCTCCCTa                                    <  1:8211700‑M2/15‑1 (MQ=255)
             atacccgccgcccagcgatttcatcaactgaatgGATGCAGCCTCCCTa                                    <  1:1725988‑M2/15‑1 (MQ=255)
              tacccgccgcccagcgatttcatcaactgaatgGATGCAGCCTCCCTaa                                   >  1:6355830‑M2/34‑49 (MQ=255)
               acccgccgcccagcgatttcatcaactgaatgGATGCAGCCTCCCTaaa                                  >  1:7778433‑M2/33‑49 (MQ=255)
               acccgccgcccagcgatttcatcaactgaatgGATGCAGCCTCCCTaaa                                  >  1:1841860‑M2/33‑49 (MQ=255)
                cccgccgcccagcgatttcatcaactgaatgGATGCAGCCTCCCTaaaa                                 <  1:3620890‑M2/18‑1 (MQ=255)
                cccgccgcccagcgatttcatcaactgaatgGATGCAGCCTCCCTaaaa                                 <  1:6458479‑M2/18‑1 (MQ=255)
                 ccgccgcccagcgatttcatcaactgaatgGATGCAGCCTCCCTaaaa                                 >  1:2496308‑M2/31‑48 (MQ=255)
                   gccgcccagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGt                             <  1:4722240‑M2/22‑1 (MQ=255)
                    ccgcccagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTg                            >  1:808957‑M2/28‑50 (MQ=255)
                      gcccagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGc                          <  1:3444525‑M2/25‑1 (MQ=255)
                       cccagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCt                         <  1:1799647‑M2/26‑1 (MQ=255)
                       cccagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCt                         <  1:5798997‑M2/26‑1 (MQ=255)
                       cccagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCt                         >  1:8252586‑M2/25‑50 (MQ=255)
                        ccagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTc                        >  1:4452476‑M2/24‑50 (MQ=255)
                         cagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTcc                       <  1:1234961‑M2/28‑1 (MQ=255)
                          agcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCa                      <  1:2351574‑M2/29‑1 (MQ=255)
                          agcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCa                      >  1:8193340‑M2/22‑50 (MQ=255)
                           gcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCAt                     <  1:3464173‑M2/30‑1 (MQ=255)
                            cgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATc                    <  1:1537519‑M2/31‑1 (MQ=255)
                            cgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATc                    <  1:5161068‑M2/31‑1 (MQ=255)
                            cgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATc                    >  1:6916506‑M2/20‑50 (MQ=255)
                            cgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATc                    <  1:2763056‑M2/31‑1 (MQ=255)
                            cgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATc                    <  1:4515738‑M2/31‑1 (MQ=255)
                            cgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATc                    >  1:2747563‑M2/20‑50 (MQ=255)
                            cgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATc                    <  1:7616754‑M2/31‑1 (MQ=255)
                             gatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCg                   <  1:4249604‑M2/32‑1 (MQ=255)
                             gatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCg                   >  1:6254331‑M2/19‑50 (MQ=255)
                              atttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGt                  <  1:2904724‑M2/33‑1 (MQ=255)
                              atttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGt                  >  1:6935464‑M2/18‑50 (MQ=255)
                               tttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTg                 <  1:390958‑M2/34‑1 (MQ=255)
                               tttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTg                 <  1:7024752‑M2/34‑1 (MQ=255)
                               tttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTg                 <  1:6086898‑M2/34‑1 (MQ=255)
                                ttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTgg                >  1:2139351‑M2/16‑50 (MQ=255)
                                 tcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTgg                <  1:4360132‑M2/35‑1 (MQ=255)
                                 tcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTgg                <  1:7526678‑M2/35‑1 (MQ=255)
                                 tcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTgg                <  1:5432666‑M2/35‑1 (MQ=255)
                                  catcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGc               <  1:8177281‑M2/36‑1 (MQ=255)
                                  catcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGc               <  1:8216055‑M2/36‑1 (MQ=255)
                                  catcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCg              >  1:3978167‑M2/14‑50 (MQ=255)
                                  catcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCg              >  1:6723813‑M2/14‑50 (MQ=255)
                                   atcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCgg             <  1:8503572‑M2/38‑1 (MQ=255)
                                   atcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCg              <  1:2358890‑M2/37‑1 (MQ=255)
                                      aactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGatt          <  1:6529555‑M2/41‑1 (MQ=255)
                                      aactgaatgGATGCAGCCTCCCTAAAAAAGAGGCTCCATCGTGGCGGat           >  1:1152736‑M2/10‑49 (MQ=255)
                                       actgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGatta         <  1:6484914‑M2/42‑1 (MQ=255)
                                       actgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGatta         <  1:7954045‑M2/42‑1 (MQ=255)
                                        ctgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGattat        <  1:3281629‑M2/43‑1 (MQ=255)
                                        ctgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGattat        <  1:6214985‑M2/43‑1 (MQ=255)
                                         tgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGattatt       >  1:3211727‑M2/7‑50 (MQ=255)
                                         tgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGattatt       >  1:1661658‑M2/7‑50 (MQ=255)
                                          gaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATtg      >  1:7804555‑M2/6‑50 (MQ=255)
                                           aatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGattatt       <  1:8212848‑M2/44‑1 (MQ=255)
                                           aatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGattatt       <  1:2749079‑M2/44‑1 (MQ=255)
                                           aatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGattatt       <  1:6643445‑M2/44‑1 (MQ=255)
                                           aatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATtgt     <  1:3371098‑M2/46‑1 (MQ=255)
                                           aatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATtg      <  1:6463130‑M2/45‑1 (MQ=255)
                                             tgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGa   <  1:5269021‑M2/48‑1 (MQ=255)
                                              gGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGAc  >  1:8301666‑M2/2‑50 (MQ=255)
                                               gatgCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGAc  >  1:5779267/1‑49 (MQ=255)
                                               gatgCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGAc  >  1:1140452/1‑49 (MQ=255)
                                               |                                                
TGATTGAAGGTTTTGGCGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGAC  >  NC_000913/4325119‑4325214

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: