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 27 [0] [0] 27 [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)25 (0.920) 21/96 0.1 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

GACGGGACCTGCCTGATACCCGCCGCCCAGCGATTTCATCAACTGAATGCTTTGGATCACCCGGCGGCTGTCCAGCATCAGTAATGACATCTCTT  >  NC_000913/4299579‑4299673
                                                |                                              
gacgGGACCTGCCTGATACCCGCCGCCCAGCGATTTCATCAACTGAATg                                                >  1:1593816/1‑49 (MQ=255)
gacgGGACCTGCCTGATACCCGCCGCCCAGCGATTTCATCAACTGAATg                                                >  1:209566/1‑49 (MQ=255)
   gggACCTGCCTGATACCCGCCGCCCAGCGATTTCATCAACTGAATGgat                                             >  1:1145440‑M1/1‑46 (MQ=255)
   gggACCTGCCTGATACCCGCCGCCCAGCGATTTCATCAACTGAATGgat                                             >  1:1577995‑M1/1‑46 (MQ=255)
   gggACCTGCCTGATACCCGCCGCCCAGCGATTTCATCAACTGAATGgat                                             >  1:360268‑M1/1‑46 (MQ=255)
    ggACCTGCCTGATACCCGCCGCCCAGCGATTTCATCAACTGAATGgatg                                            >  1:178277‑M1/1‑45 (MQ=255)
         tgcctgATACCCGCCGCCCAGCGATTTCATCAACTGAATGgatgcagcc                                       >  1:1059448‑M1/1‑40 (MQ=255)
          gcctgATACCCGCCGCCCAGCGATTTCATCAACTGAATGgatgcagcct                                      >  1:497519‑M1/1‑39 (MQ=255)
           cctgATACCCGCCGCCCAGCGATTTCATCAACTGAATGgatgcagcc                                       <  1:175975‑M1/47‑10 (MQ=255)
               aTACCCGCCGCCCAGCGATTTCATCAACTGAATGgatgcagcctccct                                  <  1:738903‑M1/48‑15 (MQ=255)
                tACCCGCCGCCCAGCGATTTCATCAACTGAATGgatgcagcctccctaa                                <  1:476274‑M1/49‑17 (MQ=255)
                  cccGCCGCCCAGCGATTTCATCAACTGAATGgatgcagcctccctaaaa                              <  1:89955‑M1/49‑19 (MQ=255)
                  cccGCCGCCCAGCGATTTCATCAACTGAATGgatgcagcctccctaaa                               <  1:534708‑M1/48‑18 (MQ=255)
                    cgccgcCCAGCGATTTCATCAACTGAATGgatgcagcctccctaaaaaa                            <  1:1210717‑M1/49‑21 (MQ=255)
                         cccAGCGATTTCATCAACTGAATGgatgcagcctccctaaaaaagtggc                       >  1:1946479‑M1/1‑24 (MQ=255)
                          ccAGCGATTTCATCAACTGAATGgatgcagcctccctaaaaaagtggct                      >  1:1005970‑M1/1‑23 (MQ=255)
                              cGATTTCATCAACTGAATGgatgcagcctccctaaaaaagtggctccat                  <  1:1751686‑M1/49‑31 (MQ=255)
                               gATTTCATCAACTGAATGgatgcagcctccctaaaaaagtggctccatc                 >  1:606862‑M1/1‑18 (MQ=255)
                                aTTTCATCAACTGAATGgatgcagcctccctaaaaaagtggctccatcg                >  1:622174‑M1/1‑17 (MQ=255)
                                  ttCATCAACTGAATGgatgcagcctccctaaaaaagtggctccatcgtg              >  1:20266‑M1/1‑15 (MQ=255)
                                      tcaACTGAATGgatgcagcctccctaaaaaagtggctccatcgtggcg           >  1:2057993‑M1/1‑11 (MQ=255)
                                         aCTGAATGgatgcagcctccctaaaaaagtggctccatcgtggcggatt       <  1:1378244‑M1/49‑42 (MQ=255)
                                          cTGAATGgatgcagcctccctaaaaaagtggctccatcgtggcggatta      <  1:1993572‑M1/49‑43 (MQ=255)
                                          cTGAATGgatgcagcctccctaaaaaagtggctccatcgtggcggatta      <  1:1924298‑M1/49‑43 (MQ=255)
                                          cTGAATGgatgcagcctccctaaaaaagtggctccatcgtggcggatta      <  1:522984‑M1/49‑43 (MQ=255)
                                             aaTGgatgcagcctccctaaaaaagtggctccatcgtggcggattattg   >  1:459884‑M1/1‑4 (MQ=255)
                                              aTGgatgcagcctccctaaaaaagtggctccatcgtggcggattattgt  <  1:282680‑M1/49‑47 (MQ=255)
                                                |                                              
GACGGGACCTGCCTGATACCCGCCGCCCAGCGATTTCATCAACTGAATGCTTTGGATCACCCGGCGGCTGTCCAGCATCAGTAATGACATCTCTT  >  NC_000913/4299579‑4299673

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.