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 61 [0] [0] 61 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)60 (1.030) 45/98 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

TGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTA  >  NC_000913/508823‑508921
                                                 |                                                 
tGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGt                                                   >  1:1871882/1‑50 (MQ=255)
  cGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTg                                                  >  1:1992011‑M1/1‑48 (MQ=255)
   gTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTgga                                                <  1:601950‑M1/50‑4 (MQ=255)
   gTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTgg                                                 >  1:749811‑M1/1‑47 (MQ=255)
   gTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTgg                                                 >  1:1756866‑M1/1‑47 (MQ=255)
   gTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTgg                                                 >  1:1380638‑M1/1‑47 (MQ=255)
       aTGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTggaacgc                                            >  1:2566392‑M1/1‑43 (MQ=255)
         gCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTggaacgcg                                           >  1:757357‑M1/1‑41 (MQ=255)
           cAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTggaacgcgcaa                                        <  1:414851‑M1/50‑12 (MQ=255)
           cAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTggaacgcgcaa                                        <  1:3877658‑M1/50‑12 (MQ=255)
           cAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTggaacgcgcaa                                        >  1:2426448‑M1/1‑39 (MQ=255)
            aaaTGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTggaacgcgcaat                                       >  1:1552239‑M1/1‑38 (MQ=255)
            aaaTGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTggaacgcgcaat                                       <  1:836169‑M1/50‑13 (MQ=255)
             aaTGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTggaacgcgcaat                                       >  1:610899‑M1/1‑37 (MQ=255)
              aTGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTggaacgcgcaatga                                     >  1:3879888‑M1/1‑36 (MQ=255)
                 cgcCACCCTAAAGCAGCGCATCCGCAATGATGTggaacgcgcaatgacga                                  <  1:2008418‑M1/50‑18 (MQ=255)
                 cgcCACCCTAAAGCAGCGCATCCGCAATGATGTggaacgcgcaatgacg                                   <  1:2724613‑M1/49‑17 (MQ=255)
                 cgcCACCCTAAAGCAGCGCATCCGCAATGATGTggaacgcgcaatgacg                                   <  1:1381870‑M1/49‑17 (MQ=255)
                    cACCCTAAAGCAGCGCATCCGCAATGATGTggaacgcgcaatgacgatca                               >  1:3200510‑M1/1‑30 (MQ=255)
                    cACCCTAAAGCAGCGCATCCGCAATGATGTggaacgcgcaatgacgatca                               <  1:1612894‑M1/50‑21 (MQ=255)
                      cccTAAAGCAGCGCATCCGCAATGATGTggaacgcgcaatgacgatcatc                             >  1:2546667‑M1/1‑28 (MQ=255)
                      cccTAAAGCAGCGCATCCGCAATGATGTggaacgcgcaatgacgatcatc                             <  1:2662844‑M1/50‑23 (MQ=255)
                       ccTAAAGCAGCGCATCCGCAATGATGTggaacgcgcaatgacgatcatcc                            <  1:3081447‑M1/50‑24 (MQ=255)
                        cTAAAGCAGCGCATCCGCAATGATGTggaacgcgcaatgacgatcatcca                           <  1:1898802‑M1/50‑25 (MQ=255)
                         tAAAGCAGCGCATCCGCAATGATGTggaacgcgcaatgacgatcatccat                          >  1:3330607‑M1/1‑25 (MQ=255)
                           aaGCAGCGCATCCGCAATGATGTggaacgcgcaatgacgatcatccatcc                        <  1:1065551‑M1/50‑28 (MQ=255)
                           aaGCAGCGCATCCGCAATGATGTggaacgcgcaatgacgatcatccatc                         <  1:1592889‑M1/49‑27 (MQ=255)
                           aaGCAGCGCATCCGCAATGATGTggaacgcgcaatgacgatcatccatc                         <  1:1497784‑M1/49‑27 (MQ=255)
                            agcagcGCATCCGCAATGATGTggaacgcgcaatgacgatcatccatcct                       <  1:4191167‑M1/50‑29 (MQ=255)
                            agcagcGCATCCGCAATGATGTggaacgcgcaatgacgatcatccatcct                       <  1:274136‑M1/50‑29 (MQ=255)
                            agcagcGCATCCGCAATGATGTggaacgcgcaatgacgatcatccatcct                       >  1:1526428‑M1/1‑22 (MQ=255)
                             gcagcGCATCCGCAATGATGTggaacgcgcaatgacgatcatccatcctc                      <  1:443934‑M1/50‑30 (MQ=255)
                              cagcGCATCCGCAATGATGTggaacgcgcaatgacgatcatccatcctct                     <  1:3453225‑M1/50‑31 (MQ=255)
                              cagcGCATCCGCAATGATGTggaacgcgcaatgacgatcatccatcctct                     >  1:3164145‑M1/1‑20 (MQ=255)
                              cagcGCATCCGCAATGATGTggaacgcgcaatgacgatcatccatcctct                     <  1:368731‑M1/50‑31 (MQ=255)
                               agcGCATCCGCAATGATGTggaacgcgcaatgacgatcatccatcctctg                    <  1:1964255‑M1/50‑32 (MQ=255)
                                gcgcATCCGCAATGATGTggaacgcgcaatgacgatcatccatcctctgc                   >  1:4328917‑M1/1‑18 (MQ=255)
                                gcgcATCCGCAATGATGTggaacgcgcaatgacgatcatccatcctctgc                   <  1:4070006‑M1/50‑33 (MQ=255)
                                  gcATCCGCAATGATGTggaacgcgcaatgacgatcatccatcctctgctt                 <  1:1967067‑M1/50‑35 (MQ=255)
                                  gcATCCGCAATGATGTggaacgcgcaatgacgatcatccatcctctgctt                 <  1:247160‑M1/50‑35 (MQ=255)
                                      ccGCAATGATGTggaacgcgcaatgacgatcatccatcctctgcttgc               >  1:2559387‑M1/1‑12 (MQ=255)
                                      ccGCAATGATGTggaacgcgcaatgacgatcatccatcctctgcttgc               >  1:4264668‑M1/1‑12 (MQ=255)
                                       cGCAATGATGTggaacgcgcaatgacgatcatccatcctctgcttgccag            <  1:2165372‑M1/50‑40 (MQ=255)
                                       cGCAATGATGTggaacgcgcaatgacgatcatccatcctctgcttgccag            >  1:3887576‑M1/1‑11 (MQ=255)
                                        gCAATGATGTggaacgcgcaatgacgatcatccatcctctgcttgccagt           <  1:548645‑M1/50‑41 (MQ=255)
                                        gCAATGATGTggaacgcgcaatgacgatcatccatcctctgcttgccagt           <  1:3087798‑M1/50‑41 (MQ=255)
                                        gCAATGATGTggaacgcgcaatgacgatcatccatcctctgcttgccag            >  1:4120732‑M1/1‑10 (MQ=255)
                                        gCAATGATGTggaacgcgcaatgacgatcatccatcctctgcttgccag            >  1:2948719‑M1/1‑10 (MQ=255)
                                        gCAATGATGTggaacgcgcaatgacgatcatccatcctctgcttgccag            >  1:4338812‑M1/1‑10 (MQ=255)
                                         cAATGATGTggaacgcgcaatgacgatcatccatcctctgcttgccagta          <  1:3503682‑M1/50‑42 (MQ=255)
                                          aaTGATGTggaacgcgcaatgacgatcatccatcctctgcttgccagtag         <  1:967082‑M1/50‑43 (MQ=255)
                                          aaTGATGTggaacgcgcaatgacgatcatccatcctctgcttgccagtag         >  1:3230575‑M1/1‑8 (MQ=255)
                                            tgatgTggaacgcgcaatgacgatcatccatcctctgcttgccagtagta       <  1:3508115‑M1/50‑45 (MQ=255)
                                            tgatgTggaacgcgcaatgacgatcatccatcctctgcttgccagtagta       <  1:204631‑M1/50‑45 (MQ=255)
                                            tgatgTggaacgcgcaatgacgatcatccatcctctgcttgccagtag         >  1:4352950‑M1/1‑6 (MQ=255)
                                            tgatgTggaacgcgcaatgacgatcatccatcctctgcttgccagtag         >  1:2933167‑M1/1‑6 (MQ=255)
                                              atgTggaacgcgcaatgacgatcatccatcctctgcttgccagtagtagc     >  1:4138823‑M1/1‑4 (MQ=255)
                                              atgTggaacgcgcaatgacgatcatccatcctctgcttgccagtagtagc     <  1:3607960‑M1/50‑47 (MQ=255)
                                              atgTggaacgcgcaatgacgatcatccatcctctgcttgccagtagtagc     >  1:1275215‑M1/1‑4 (MQ=255)
                                              atgTggaacgcgcaatgacgatcatccatcctctgcttgccagtagtag      <  1:600748‑M1/49‑46 (MQ=255)
                                                 tggaacgcgcaatgacgatcatccatcctctgcttgccagtagtagcgca  >  1:483445‑M1/1‑1 (MQ=255)
                                                 |                                                 
TGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTA  >  NC_000913/508823‑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.