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 27 [0] [0] 27 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)27 (0.890) 25/98 0.1 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

GAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCA  >  NC_000913/549572‑549665
                                             |                                                
cgatgccaaatgcgccaccctaaagcagcgcatccgcaatgatgtGGAAc                                              >  1:2270031‑M2/46‑50 (MQ=255)
cgatgccaaatgcgccaccctaaagcagcgcatccgcaatgatgtGGAAc                                              >  1:157500‑M2/46‑50 (MQ=255)
 gatgccaaatgcgccaccctaaagcagcgcatccgcaatgatgtGGAAc                                              <  1:997365‑M2/5‑1 (MQ=255)
 gatgccaaatgcgccaccctaaagcagcgcatccgcaatgatgtGGAAc                                              <  1:2429641‑M2/5‑1 (MQ=255)
     ccaaatgcgccaccctaaagcagcgcatccgcaatgatgtGGAACGCGCa                                         <  1:379418‑M2/10‑1 (MQ=255)
      caaatgcgccaccctaaagcagcgcatccgcaatgatgtGGAACGCGCaa                                        <  1:1840012‑M2/11‑1 (MQ=255)
        aatgcgccaccctaaagcagcgcatccgcaatgatgtGGAACGCGCAATg                                      <  1:1783659‑M2/13‑1 (MQ=255)
              ccaccctaaagcagcgcatccgcaatgatgtGGAACGCGCAATGACGatc                                <  1:1569270‑M2/19‑1 (MQ=255)
                accctaaagcagcgcatccgcaatgatgtGGAACGCGCAATGACGatcat                              <  1:2172124‑M2/21‑1 (MQ=255)
                  cctaaagcagcgcatccgcaatgatgtGGAACGCGCAATGACGatcatc                             >  1:1235300‑M2/28‑49 (MQ=255)
                   ctaaagcagcgcatccgcaatgatgtGGAACGCGCAATGACGATcatcca                           >  1:1661416‑M2/27‑50 (MQ=255)
                    taaagcagcgcatccgcaatgatgtGGAACGCGCAATGACGATcatccat                          >  1:1245067‑M2/26‑50 (MQ=255)
                     aaagcagcgcatccgcaatgatgtGGAACGCGCAATGACGATcatccatc                         >  1:1255009‑M2/25‑50 (MQ=255)
                      aagcagcgcatccgcaatgatgtGGAACGCGCAATGACGATcatccat                          <  1:2675310‑M2/25‑1 (MQ=255)
                      aagcagcgcatccgcaatgatgtGGAACGCGCAATGACGATCatccatcc                        <  1:2621758‑M2/27‑1 (MQ=255)
                        gcagcgcatccgcaatgatgtGGAACGCGCAATGACGATCATCCATCctc                      >  1:1719931‑M2/22‑50 (MQ=255)
                           gcgcatccgcaatgatgtGGAACGCGCAATGACGATCATCCATCCTCTGc                   <  1:2035675‑M2/32‑1 (MQ=255)
                                tccgcaatgatgtGGAACGCGCAATGACGATCATCCATCCTCTGCTTGcc              <  1:387844‑M2/37‑1 (MQ=255)
                                 ccgcaatgatgtGGAACGCGCAATGACGATCATCCATCCTCTGCTTGc               >  1:1797184‑M2/13‑48 (MQ=255)
                                   gcaatgatgtGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCag            >  1:1015743‑M2/11‑49 (MQ=255)
                                    caatgatgtGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCagta          <  1:2630846‑M2/41‑1 (MQ=255)
                                      atgatgtGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCagtagt        <  1:971036‑M2/43‑1 (MQ=255)
                                      atgatgtGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCagtag         >  1:242808‑M2/8‑49 (MQ=255)
                                        gatgtGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGtagtag      >  1:1531257‑M2/6‑50 (MQ=255)
                                          tgtGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAgcg    <  1:2109134‑M2/47‑1 (MQ=255)
                                           gtGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAgcg    >  1:838483‑M2/3‑49 (MQ=255)
                                            tGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCa  >  1:2503130‑M2/2‑50 (MQ=255)
                                             |                                                
GAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCA  >  NC_000913/549572‑549665

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