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

CGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTA  >  NC_000913/508828‑508921
                                            |                                                 
cGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTggaac                                              >  1:2270031‑M1/1‑45 (MQ=255)
cGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTggaac                                              >  1:157500‑M1/1‑45 (MQ=255)
 gATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTggaac                                              <  1:997365‑M1/49‑6 (MQ=255)
 gATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTggaac                                              <  1:2429641‑M1/49‑6 (MQ=255)
     ccAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTggaacgcgca                                         <  1:379418‑M1/50‑11 (MQ=255)
      cAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTggaacgcgcaa                                        <  1:1840012‑M1/50‑12 (MQ=255)
        aaTGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTggaacgcgcaatg                                      <  1:1783659‑M1/50‑14 (MQ=255)
              ccACCCTAAAGCAGCGCATCCGCAATGATGTggaacgcgcaatgacgatc                                <  1:1569270‑M1/50‑20 (MQ=255)
                aCCCTAAAGCAGCGCATCCGCAATGATGTggaacgcgcaatgacgatcat                              <  1:2172124‑M1/50‑22 (MQ=255)
                  ccTAAAGCAGCGCATCCGCAATGATGTggaacgcgcaatgacgatcatc                             >  1:1235300‑M1/1‑27 (MQ=255)
                   cTAAAGCAGCGCATCCGCAATGATGTggaacgcgcaatgacgatcatcca                           >  1:1661416‑M1/1‑26 (MQ=255)
                    tAAAGCAGCGCATCCGCAATGATGTggaacgcgcaatgacgatcatccat                          >  1:1245067‑M1/1‑25 (MQ=255)
                     aaaGCAGCGCATCCGCAATGATGTggaacgcgcaatgacgatcatccatc                         >  1:1255009‑M1/1‑24 (MQ=255)
                      aaGCAGCGCATCCGCAATGATGTggaacgcgcaatgacgatcatccatcc                        <  1:2621758‑M1/50‑28 (MQ=255)
                      aaGCAGCGCATCCGCAATGATGTggaacgcgcaatgacgatcatccat                          <  1:2675310‑M1/48‑26 (MQ=255)
                        gcagcGCATCCGCAATGATGTggaacgcgcaatgacgatcatccatcctc                      >  1:1719931‑M1/1‑21 (MQ=255)
                           gcgcATCCGCAATGATGTggaacgcgcaatgacgatcatccatcctctgc                   <  1:2035675‑M1/50‑33 (MQ=255)
                                tCCGCAATGATGTggaacgcgcaatgacgatcatccatcctctgcttgcc              <  1:387844‑M1/50‑38 (MQ=255)
                                 ccGCAATGATGTggaacgcgcaatgacgatcatccatcctctgcttgc               >  1:1797184‑M1/1‑12 (MQ=255)
                                   gCAATGATGTggaacgcgcaatgacgatcatccatcctctgcttgccag            >  1:1015743‑M1/1‑10 (MQ=255)
                                    cAATGATGTggaacgcgcaatgacgatcatccatcctctgcttgccagta          <  1:2630846‑M1/50‑42 (MQ=255)
                                      atgatgTggaacgcgcaatgacgatcatccatcctctgcttgccagtagt        <  1:971036‑M1/50‑44 (MQ=255)
                                      atgatgTggaacgcgcaatgacgatcatccatcctctgcttgccagtag         >  1:242808‑M1/1‑7 (MQ=255)
                                        gatgTggaacgcgcaatgacgatcatccatcctctgcttgccagtagtag      >  1:1531257‑M1/1‑5 (MQ=255)
                                          tgtggaacgcgcaatgacgatcatccatcctctgcttgccagtagtagcg    <  1:2109134‑M1/50‑48 (MQ=255)
                                           gtggaacgcgcaatgacgatcatccatcctctgcttgccagtagtagcg    >  1:838483‑M1/1‑2 (MQ=255)
                                            tggaacgcgcaatgacgatcatccatcctctgcttgccagtagtagcgca  >  1:2503130‑M1/1‑1 (MQ=255)
                                            |                                                 
CGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTA  >  NC_000913/508828‑508921

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: 
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.