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 58 [0] [0] 58 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)57 (0.780) 40/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

TGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTA  >  NC_000913/508823‑508921
                                                 |                                                 
tGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGt                                                   <  1:998694/50‑1 (MQ=255)
  cGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTg                                                  >  1:432680‑M1/1‑48 (MQ=255)
  cGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTg                                                  <  1:290853‑M1/49‑2 (MQ=255)
  cGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTg                                                  <  1:593635‑M1/49‑2 (MQ=255)
     cGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTggaac                                              >  1:5287413‑M1/1‑45 (MQ=255)
      gATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTggaacg                                             >  1:4803933‑M1/1‑44 (MQ=255)
      gATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTggaacg                                             >  1:3527776‑M1/1‑44 (MQ=255)
           cAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTggaacgcgcaa                                        <  1:1408171‑M1/50‑12 (MQ=255)
           cAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTggaacgcgcaa                                        >  1:2700851‑M1/1‑39 (MQ=255)
            aaaTGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTggaacgcgcaat                                       <  1:4908385‑M1/50‑13 (MQ=255)
            aaaTGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTggaacgcgcaat                                       <  1:4280990‑M1/50‑13 (MQ=255)
             aaTGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTggaacgcgcaat                                       >  1:1719992‑M1/1‑37 (MQ=255)
              aTGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTggaacgcgcaatga                                     <  1:1957377‑M1/50‑15 (MQ=255)
              aTGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTggaacgcgcaatga                                     <  1:2952981‑M1/50‑15 (MQ=255)
               tGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTggaacgcgcaatgac                                    <  1:2493988‑M1/50‑16 (MQ=255)
               tGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTggaacgcgcaatgac                                    <  1:250488‑M1/50‑16 (MQ=255)
               tGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTggaacgcgcaatgac                                    <  1:2586891‑M1/50‑16 (MQ=255)
                gcgcCACCCTAAAGCAGCGCATCCGCAATGATGTggaacgcgcaatgac                                    >  1:4548318‑M1/1‑34 (MQ=255)
                gcgcCACCCTAAAGCAGCGCATCCGCAATGATGTggaacgcgcaatgac                                    >  1:379784‑M1/1‑34 (MQ=255)
                gcgcCACCCTAAAGCAGCGCATCCGCAATGATGTggaacgcgcaatgac                                    >  1:1892552‑M1/1‑34 (MQ=255)
                 cgcCACCCTAAAGCAGCGCATCCGCAATGATGTggaacgcgcaatgacga                                  >  1:2103241‑M1/1‑33 (MQ=255)
                   ccACCCTAAAGCAGCGCATCCGCAATGATGTggaacgcgcaatgacgatc                                <  1:700639‑M1/50‑20 (MQ=255)
                    cACCCTAAAGCAGCGCATCCGCAATGATGTggaacgcgcaatgacgatca                               <  1:4503423‑M1/50‑21 (MQ=255)
                     aCCCTAAAGCAGCGCATCCGCAATGATGTggaacgcgcaatgacgatcat                              <  1:2929676‑M1/50‑22 (MQ=255)
                     aCCCTAAAGCAGCGCATCCGCAATGATGTggaacgcgcaatgacgatcat                              >  1:4305315‑M1/1‑29 (MQ=255)
                      cccTAAAGCAGCGCATCCGCAATGATGTggaacgcgcaatgacgatcatc                             >  1:565793‑M1/1‑28 (MQ=255)
                      cccTAAAGCAGCGCATCCGCAATGATGTggaacgcgcaatgacgatcatc                             <  1:434471‑M1/50‑23 (MQ=255)
                       ccTAAAGCAGCGCATCCGCAATGATGTggaacgcgcaatgacgatcatc                             >  1:1624864‑M1/1‑27 (MQ=255)
                        cTAAAGCAGCGCATCCGCAATGATGTggaacgcgcaatgacgatcatcca                           <  1:4630150‑M1/50‑25 (MQ=255)
                           aaGCAGCGCATCCGCAATGATGTggaacgcgcaatgacgatcatccatcc                        >  1:585147‑M1/1‑23 (MQ=255)
                           aaGCAGCGCATCCGCAATGATGTggaacgcgcaatgacgatcatccatcc                        <  1:1871491‑M1/50‑28 (MQ=255)
                           aaGCAGCGCATCCGCAATGATGTggaacgcgcaatgacgatcatccatc                         <  1:1269381‑M1/49‑27 (MQ=255)
                           aaGCAGCGCATCCGCAATGATGTggaacgcgcaatgacgatcatccatc                         <  1:4710433‑M1/49‑27 (MQ=255)
                            agcagcGCATCCGCAATGATGTggaacgcgcaatgacgatcatccatcct                       >  1:672701‑M1/1‑22 (MQ=255)
                             gcagcGCATCCGCAATGATGTggaacgcgcaatgacgatcatccatcctc                      <  1:1701226‑M1/50‑30 (MQ=255)
                             gcagcGCATCCGCAATGATGTggaacgcgcaatgacgatcatccatcctc                      <  1:3431284‑M1/50‑30 (MQ=255)
                              cagcGCATCCGCAATGATGTggaacgcgcaatgacgatcatccatcctct                     <  1:1458619‑M1/50‑31 (MQ=255)
                              cagcGCATCCGCAATGATGTggaacgcgcaatgacgatcatccatcctct                     <  1:5216683‑M1/50‑31 (MQ=255)
                                 cgcATCCGCAATGATGTggaacgcgcaatgacgatcatccatcctctgct                  >  1:5455747‑M1/1‑17 (MQ=255)
                                  gcATCCGCAATGATGTggaacgcgcaatgacgatcatccatcctctgctt                 >  1:3130571‑M1/1‑16 (MQ=255)
                                  gcATCCGCAATGATGTggaacgcgcaatgacgatcatccatcctctgctt                 <  1:4606217‑M1/50‑35 (MQ=255)
                                  gcATCCGCAATGATGTggaacgcgcaatgacgatcatccatcctctgctt                 <  1:4280795‑M1/50‑35 (MQ=255)
                                   cATCCGCAATGATGTggaacgcgcaatgacgatcatccatcctctgcttg                >  1:4831781‑M1/1‑15 (MQ=255)
                                    aTCCGCAATGATGTggaacgcgcaatgacgatcatccatcctctgcttgc               <  1:3899152‑M1/50‑37 (MQ=255)
                                    aTCCGCAATGATGTggaacgcgcaatgacgatcatccatcctctgcttgc               >  1:732553‑M1/1‑14 (MQ=255)
                                    aTCCGCAATGATGTggaacgcgcaatgacgatcatccatcctctgcttgc               >  1:3207746‑M1/1‑14 (MQ=255)
                                    aTCCGCAATGATGTggaacgcgcaatgacgatcatccatcctctgcttgc               >  1:1602480‑M1/1‑14 (MQ=255)
                                       cGCAATGATGTggaacgcgcaatgacgatcatccatcctctgcttgccag            >  1:1237753‑M1/1‑11 (MQ=255)
                                       cGCAATGATGTggaacgcgcaatgacgatcatccatcctctgcttgccag            <  1:2026207‑M1/50‑40 (MQ=255)
                                          aaTGATGTggaacgcgcaatgacgatcatccatcctctgcttgccagtag         <  1:4962306‑M1/50‑43 (MQ=255)
                                           atgatgTggaacgcgcaatgacgatcatccatcctctgcttgccagtagt        <  1:3709277‑M1/50‑44 (MQ=255)
                                           atgatgTggaacgcgcaatgacgatcatccatcctctgcttgccagtagt        <  1:4648155‑M1/50‑44 (MQ=255)
                                           atgatgTggaacgcgcaatgacgatcatccatcctctgcttgccagtagt        <  1:2277636‑M1/50‑44 (MQ=255)
                                             gatgTggaacgcgcaatgacgatcatccatcctctgcttgccagtagtag      >  1:3045030‑M1/1‑5 (MQ=255)
                                             gatgTggaacgcgcaatgacgatcatccatcctctgcttgccagtagtag      >  1:3806035‑M1/1‑5 (MQ=255)
                                              atgTggaacgcgcaatgacgatcatccatcctctgcttgccagtagtagc     <  1:96304‑M1/50‑47 (MQ=255)
                                               tgtggaacgcgcaatgacgatcatccatcctctgcttgccagtagtagcg    <  1:1492133‑M1/50‑48 (MQ=255)
                                                 tggaacgcgcaatgacgatcatccatcctctgcttgccagtagtagcgca  >  1:2234010‑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 < 33 ≤ 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.