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 37 [0] [0] 36 [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)36 (0.960) 31/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

CTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGA  >  NC_000913/4325118‑4325213
                                                |                                               
acgggacctgcctgatacccgccgcccagcgatttcatcaactgaatgGa                                                <  1:2632415‑M2/2‑1 (MQ=255)
 cgggacctgcctgatacccgccgcccagcgatttcatcaactgaatgGAt                                               <  1:2503994‑M2/3‑1 (MQ=255)
  gggacctgcctgatacccgccgcccagcgatttcatcaactgaatgGATg                                              >  1:2549629‑M2/47‑50 (MQ=255)
   ggacctgcctgatacccgccgcccagcgatttcatcaactgaatgGATg                                              >  1:2719901‑M2/46‑49 (MQ=255)
    gacctgcctgatacccgccgcccagcgatttcatcaactgaatgGATg                                              >  1:1447678‑M2/45‑48 (MQ=255)
    gacctgcctgatacccgccgcccagcgatttcatcaactgaatgGATGCa                                            <  1:1558374‑M2/6‑1 (MQ=255)
       ctgcctgatacccgccgcccagcgatttcatcaactgaatgGATGCAGcc                                         >  1:1364972‑M2/42‑50 (MQ=255)
         gcctgatacccgccgcccagcgatttcatcaactgaatgGATGCAGCCTc                                       >  1:2771359‑M2/40‑50 (MQ=255)
          cctgatacccgccgcccagcgatttcatcaactgaatgGATGCAGCCTcc                                      <  1:1160786‑M2/12‑1 (MQ=255)
           ctgatacccgccgcccagcgatttcatcaactgaatgGATGCAGCCTccc                                     >  1:835729‑M2/38‑50 (MQ=255)
              atacccgccgcccagcgatttcatcaactgaatgGATGCAGCCTCCCt                                    <  1:361495‑M2/14‑1 (MQ=255)
                acccgccgcccagcgatttcatcaactgaatgGATGCAGCCTCCCTaaa                                 >  1:350569‑M2/33‑49 (MQ=255)
                 cccgccgcccagcgatttcatcaactgaatgGATGCAGCCTCCCTaaaa                                <  1:2712954‑M2/18‑1 (MQ=255)
                  ccgccgcccagcgatttcatcaactgaatgGATGCAGCCTCCCTaaaaaa                              <  1:2760778‑M2/20‑1 (MQ=255)
                       gcccagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGc                         >  1:1542493‑M2/26‑50 (MQ=255)
                          cagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTc                       >  1:1619469‑M2/23‑49 (MQ=255)
                           agcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCa                     >  1:2600887‑M2/22‑50 (MQ=255)
                             cgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATc                   <  1:1092774‑M2/31‑1 (MQ=255)
                             cgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATc                   <  1:669321‑M2/31‑1 (MQ=255)
                              gatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCg                  <  1:517935‑M2/32‑1 (MQ=255)
                              gatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCg                  <  1:1859771‑M2/32‑1 (MQ=255)
                                 ttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTgg               >  1:2563661‑M2/16‑50 (MQ=255)
                                 ttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTgg               >  1:2354283‑M2/16‑50 (MQ=255)
                                   catcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGc              <  1:1128493‑M2/36‑1 (MQ=255)
                                   catcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCg             >  1:1576591‑M2/14‑50 (MQ=255)
                                       aactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGat          >  1:741317‑M2/10‑49 (MQ=255)
                                       aactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGat          <  1:2621058‑M2/40‑1 (MQ=255)
                                        actgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGatta        >  1:2731000‑M2/9‑50 (MQ=255)
                                         ctgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGattat       >  1:77431‑M2/8‑50 (MQ=255)
                                         ctgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGattat       <  1:1426412‑M2/43‑1 (MQ=255)
                                          tgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGattatt      >  1:101102‑M2/7‑50 (MQ=255)
                                          tgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGattatt      >  1:761596‑M2/7‑50 (MQ=255)
                                            aatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATtgt    >  1:472197‑M2/5‑50 (MQ=255)
                                             atgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATtgt    >  1:2366262‑M2/4‑49 (MQ=255)
                                              tgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGa  >  1:542977‑M2/3‑50 (MQ=255)
                                              tgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGa  >  1:8535‑M2/3‑50 (MQ=255)
                                                |                                               
CTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGA  >  NC_000913/4325118‑4325213

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