Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 508,873 Δ40,744 bp copA[ylbE] 37 genes

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 508873 549616 40744 75 [0] [0] 75 copA–[ylbE] copA,copA,copA,glsA,ybaT,cueR,ybbJ,qmcA,fetA,fetB,cnoX,ybbO,tesA,ybbA,ybbP,rhsD,ybbC,ylbH,ybbD,ylbI,ylbG,selU,allS,allA,allR,gcl,hyi,glxR,ybbW,allB,ybbY,glxK,allE,allC,allD,fdrA,[ylbE]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 5088720 (0.000)73 (1.010) 39/68 0.0 100% intergenic (+313/+3) ybaQ/copA putative DNA‑binding transcriptional regulator YbaQ/Cu(+) exporting P‑type ATPase
?NC_000913 549617 = 0 (0.000)intergenic (+1/‑10) ylbE/ylbF DUF1116 domain‑containing protein YlbE/DUF2877 domain‑containing protein YlbF

GGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTG  >  NC_000913/549583‑549651
                                  |                                  
gcgccaccctaaagcagcgcatccgcaatgatgtg                                    <  1:3868785‑M2/1‑1 (MQ=255)
gcgccaccctaaagcagcgcatccgcaatgatgtg                                    <  1:7478691‑M2/1‑1 (MQ=255)
gcgccaccctaaagcagcgcatccgcaatgatgtg                                    <  1:1492636‑M2/1‑1 (MQ=255)
gcgccaccctaaagcagcgcatccgcaatgatgtg                                    <  1:2119802‑M2/1‑1 (MQ=255)
 cgccaccctaaagcagcgcatccgcaatgatgtgg                                   <  1:7220350‑M2/2‑1 (MQ=255)
 cgccaccctaaagcagcgcatccgcaatgatgtgg                                   <  1:5446512‑M2/2‑1 (MQ=255)
 cgccaccctaaagcagcgcatccgcaatgatgtgg                                   <  1:6949428‑M2/2‑1 (MQ=255)
 cgccacccgaaagcagcgcatccgcaatgatgtgg                                   <  1:6068953‑M2/2‑1 (MQ=255)
  gccaccctaaagcagcgcatccgcaatgatgtGGa                                  <  1:6511802‑M2/3‑1 (MQ=255)
   ccaccctaaagcagcgcatccgcaatgatgtGGaa                                 >  1:7559087‑M2/32‑35 (MQ=255)
      ccctaaagcagcgcatccgcaatgatgtGGAAcgc                              >  1:1500211‑M2/29‑35 (MQ=255)
      ccctaaagcagcgcatccgcaatgatgtGGAAcgc                              >  1:153015‑M2/29‑35 (MQ=255)
       cctaaagcagcgcatccgcaatgatgtGGAAcgcg                             >  1:5639526‑M2/28‑35 (MQ=255)
        ctaaagcagcgcatccgcaatgatgtGGAACgcgc                            >  1:7341803‑M2/27‑35 (MQ=255)
        ctaaagcagcgcatccgcaatgatgtGGAACgcgc                            >  1:7485703‑M2/27‑35 (MQ=255)
         taaagcagcgcatccgcaatgatgtGGAACGCGCa                           <  1:1147861‑M2/10‑1 (MQ=255)
          aaagcagcgcatccgcaatgatgtGGAACGCGCaa                          <  1:7465615‑M2/11‑1 (MQ=255)
          aaagcagcgcatccgcaatgatgtGGAACGCGCaa                          <  1:6323903‑M2/11‑1 (MQ=255)
           aagcagcgcatccgcaatgatgtGGAACGCGCAAt                         <  1:1990480‑M2/12‑1 (MQ=255)
           aagcagcgcatccgcaatgatgtGGAACGCGCAAt                         <  1:5036049‑M2/12‑1 (MQ=255)
           aagcagcgcatccgcaatgatgtGGAACGCGCAAt                         >  1:187292‑M2/24‑35 (MQ=255)
           aagcagcgcatccgcaatgatgtGGAACGCGCAAt                         <  1:4564518‑M2/12‑1 (MQ=255)
            agcagcgcatccgcaatgatgtGGAACGCGCAATg                        >  1:661883‑M2/23‑35 (MQ=255)
            agcagcgcatccgcaatgatgtGGAACGCGCAATg                        >  1:6250829‑M2/23‑35 (MQ=255)
              cagcgcatccgcaatgatgtGGAACGCGCAATGAc                      <  1:5147374‑M2/15‑1 (MQ=255)
               agcgcatccgcaatgatgtGGAACGCGCAATGACg                     <  1:5381856‑M2/16‑1 (MQ=255)
                gcgcatccgcaatgatgtGGAACGCGCAATGACGa                    <  1:2092827‑M2/17‑1 (MQ=255)
                gcgcatccgcaatgatgtGGAACGCGCAATGACGa                    <  1:4891059‑M2/17‑1 (MQ=255)
                 cgcatccgcaatgatgtGGAACGCGCAATGACGat                   <  1:3778200‑M2/18‑1 (MQ=255)
                 cgcatccgcaatgatgtGGAACGCGCAATGACGat                   <  1:2046794‑M2/18‑1 (MQ=255)
                  gcatccgcaatgatgtGGAACGCGCAATGACGatc                  <  1:3166892‑M2/19‑1 (MQ=255)
                   catccgcaatgatgtGGAACGCGCAATGACGatca                 <  1:456797‑M2/20‑1 (MQ=255)
                    ttccgcaatgatgtGGAACGCGCAATGACGatcat                <  1:1481041‑M2/21‑1 (MQ=255)
                     tccgcaatgatgtGGAACGCGCAATGACGatcatc               <  1:5533712‑M2/22‑1 (MQ=255)
                     tccgcaatgatgtGGAACGCGCAATGACGatcatc               <  1:7513366‑M2/22‑1 (MQ=255)
                     tccgcaatgatgtGGAACGCGCAATGACGatcatc               >  1:7494645‑M2/14‑35 (MQ=255)
                       cgcaatgatgtGGAACGCGCAATGACGATcatcca             >  1:7673111‑M2/12‑35 (MQ=255)
                       cgcaatgatgtGGAACGCGCAATGACGATcatcca             <  1:1138155‑M2/24‑1 (MQ=255)
                       cgcaatgatgtGGAACGCGCAATGACGATcatcca             >  1:1622811‑M2/12‑35 (MQ=255)
                       cgcaatgatgtGGAACGCGCAATGACGATcatcca             <  1:5497760‑M2/24‑1 (MQ=255)
                        gcaatgatgtGGAACGCGCAATGACGATcatccat            <  1:7312134‑M2/25‑1 (MQ=255)
                         caatgatgtGGAACGCGCAATGACGATcatccatc           >  1:3163823‑M2/10‑35 (MQ=255)
                         caatgatgtGGAACGCGCAATGACGATcatccatc           <  1:7358401‑M2/26‑1 (MQ=255)
                         caatgatgtGGAACGCGCAATGACGATcatccatc           <  1:1237211‑M2/26‑1 (MQ=255)
                         caatgatgtGGAACGCGCAATGACGATcatccatc           <  1:5615734‑M2/26‑1 (MQ=255)
                          aatgatgtGGAACGCGCAATGACGATCatccatcc          <  1:6126698‑M2/27‑1 (MQ=255)
                          aatgatgtGGAACGCGCAATGACGATCatccatcc          >  1:5562007‑M2/9‑35 (MQ=255)
                          aatgatgtGGAACGCGCAATGACGATCatccatcc          <  1:7763515‑M2/27‑1 (MQ=255)
                          aatgatgtGGAACGCGCAATGACGATCatccatcc          >  1:1947017‑M2/9‑35 (MQ=255)
                          aatgatgtGGAACGCGCAATGACGATCatccatcc          >  1:4779970‑M2/9‑35 (MQ=255)
                          aatgatgtGGAACGCGCAATGACGATCatccatcc          <  1:1050208‑M2/27‑1 (MQ=255)
                          aatgatgtGGAACGCGCAATGACGATCatccatcc          <  1:442719‑M2/27‑1 (MQ=255)
                           atgatgtGGAACGCGCAATGNCGATCATCCATCct         <  1:4090069‑M2/28‑1 (MQ=255)
                           atgatgtGGAACGCGCAATGACGATCATCCATCct         >  1:5395768‑M2/8‑35 (MQ=255)
                           atgatgtGGAACGCGCAATGACGATCATCCATCct         <  1:1246670‑M2/28‑1 (MQ=255)
                           atgatgtGGAACGCGCAATGACGATCATCCATCct         <  1:7357001‑M2/28‑1 (MQ=255)
                            tgatgtGGAACGCGCAATGACGATCATCCATCctc        <  1:4728296‑M2/29‑1 (MQ=255)
                            tgatgtGGAACGCGCAATGACGATCATCCATCctc        <  1:1807430‑M2/29‑1 (MQ=255)
                             gatgtGGAACGCGCAATGACGATCATCCATCctct       <  1:4384339‑M2/30‑1 (MQ=255)
                             gatgtGGAACGCGCAATGACGATCATCCATCctct       >  1:3753390‑M2/6‑35 (MQ=255)
                             gatgtGGAACGCGCAATGACGATCATCCATCctct       <  1:2741460‑M2/30‑1 (MQ=255)
                              atgtGGAACGCGCAATGACGATCATCCATCCTCTg      <  1:1123473‑M2/31‑1 (MQ=255)
                              atgtGGAACGCGCAATGACGATCATCCATCCTCTg      >  1:5390371‑M2/5‑35 (MQ=255)
                              atgtGGAACGCGCAATGACGATCATCCATCCTCTg      >  1:2292546‑M2/5‑35 (MQ=255)
                              atgtGGAACGCGCAATGACGATCATCCATCCTCTg      >  1:6362094‑M2/5‑35 (MQ=255)
                               tgtGGAACGCGCAATGACGATCATCCATCCTCTGc     <  1:7608867‑M2/32‑1 (MQ=255)
                               tgtGGAACGCGCAATGACGATCATCCATCCTCTGc     >  1:2325630‑M2/4‑35 (MQ=255)
                                gtGGAACGCGCAATGACGATCATCCATCCTCTGCt    >  1:2545591‑M2/3‑35 (MQ=255)
                                gtGGAACGCGCAATGACGATCATCCATCCTCTGCt    <  1:3857709‑M2/33‑1 (MQ=255)
                                gtGGAACGCGCAATGACGATCATCCATCCTCTGCt    <  1:5779668‑M2/33‑1 (MQ=255)
                                 tGGAACGCGCAATGACGATCATCCATCCTCTGCtt   >  1:1953449‑M2/2‑35 (MQ=255)
                                 tGGAACGCGCAATGACGATCATCCATCCTCTGCtt   >  1:1853728‑M2/2‑35 (MQ=255)
                                 tGGAACGCGCAATGACGATCATCCATCCTCTGCtt   >  1:4097488‑M2/2‑35 (MQ=255)
                                  ggAACGCGCAATGACGATCATCCATCCTCTGCTTg  <  1:6734038/35‑1 (MQ=255)
                                  ggAACGCGCAATGACGATCATCCATCCTCTGCTTg  <  1:1238224/35‑1 (MQ=255)
                                  |                                  
GGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTG  >  NC_000913/549583‑549651

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