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 30 [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)27 (0.890) 23/98 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

GCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGA  >  NC_000913/4325117‑4325213
                                                 |                                               
gacgggacctgcctgatacccgccgcccagcgatttcatcaactgaatgg                                                 <  1:699651‑M2/1‑1 (MQ=255)
gacgggacctgcctgatacccgccgcccagcgatttcatcaactgaatgg                                                 <  1:1540610‑M2/1‑1 (MQ=255)
  cgggacctgcctgatacccgccgcccagcgatttcatcaactgaatgGAt                                               <  1:582453‑M2/3‑1 (MQ=255)
  cgggacctgcctgatacccgccgcccagcgatttcatcaactgaatgGAt                                               <  1:247631‑M2/3‑1 (MQ=255)
     gacctgcctgatacccgccgcccagcgatttcatcaactgaatgTATGCa                                            <  1:1890714‑M2/6‑1 (MQ=255)
     gacctgcctgatacccgccgcccagcgatttcatcaactgaatgGATGCa                                            <  1:2373631‑M2/6‑1 (MQ=255)
           cctgatacccgccgcccagcgatttcatcaactgaatgGATGCAGCCTcc                                      >  1:965567‑M2/39‑50 (MQ=255)
              gatacccgccgcccagcgatttcatcaactgaatgGATGCAGCCTCCCt                                    >  1:2331556‑M2/36‑49 (MQ=255)
               atacccgccgcccagcgatttcatcaactgaatgGATGCAGCCTCCCt                                    <  1:2082540‑M2/14‑1 (MQ=255)
               atacccgccgcccagcgatttcatcaactgaatgGATGCAGCCTCCCTaa                                  <  1:2497814‑M2/16‑1 (MQ=255)
                   ccgccgcccagcgatttcatcaactgaatgGATGCAGCCTCCCTaaaaaa                              <  1:325859‑M2/20‑1 (MQ=255)
                     gccgcccagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGt                            <  1:2115024‑M2/22‑1 (MQ=255)
                         cccagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCt                        <  1:2703519‑M2/26‑1 (MQ=255)
                          ccagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTc                       <  1:119492‑M2/27‑1 (MQ=255)
                          ccagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTc                       <  1:185086‑M2/27‑1 (MQ=255)
                           cagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTcc                      <  1:2612626‑M2/28‑1 (MQ=255)
                           cagcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTc                       >  1:1974966‑M2/23‑49 (MQ=255)
                            agcgatttcatcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCa                     <  1:455661‑M2/29‑1 (MQ=255)
                                     atcaactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCg             >  1:629949‑M2/13‑49 (MQ=255)
                                       caactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGat          >  1:1926229‑M2/11‑50 (MQ=255)
                                        aactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGat          >  1:1911256‑M2/10‑49 (MQ=255)
                                        aactgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGat          <  1:566615‑M2/40‑1 (MQ=255)
                                         actgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGatta        <  1:1202090‑M2/42‑1 (MQ=255)
                                          ctgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGattat       >  1:821506‑M2/8‑50 (MQ=255)
                                           tgaatgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGattatt      >  1:1731698‑M2/7‑50 (MQ=255)
                                              atgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATtgtg   <  1:2283439‑M2/47‑1 (MQ=255)
                                               tgGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGa  >  1:578957‑M2/3‑50 (MQ=255)
                                                 |                                               
GCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGA  >  NC_000913/4325117‑4325213

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