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

TTTGAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCAC  >  NC_000913/549569‑549666
                                                |                                                 
cgtcgatgccaaatgcgccaccctaaagcagcgcatccgcaatgatgtg                                                   >  1:432680‑M2/49‑49 (MQ=255)
cgtcgatgccaaatgcgccaccctaaagcagcgcatccgcaatgatgtg                                                   <  1:593635‑M2/1‑1 (MQ=255)
cgtcgatgccaaatgcgccaccctaaagcagcgcatccgcaatgatgtg                                                   <  1:290853‑M2/1‑1 (MQ=255)
   cgatgccaaatgcgccaccctaaagcagcgcatccgcaatgatgtGGAAc                                               >  1:5287413‑M2/46‑50 (MQ=255)
    gatgccaaatgcgccaccctaaagcagcgcatccgcaatgatgtGGAAcg                                              >  1:4803933‑M2/45‑50 (MQ=255)
    gatgccaaatgcgccaccctaaagcagcgcatccgcaatgatgtGGAAcg                                              >  1:3527776‑M2/45‑50 (MQ=255)
         caaatgcgccaccctaaagcagcgcatccgcaatgatgtGGAACGCGCaa                                         <  1:1408171‑M2/11‑1 (MQ=255)
         caaatgcgccaccctaaagcagcgcatccgcaatgatgtGGAACGCGCaa                                         >  1:2700851‑M2/40‑50 (MQ=255)
          aaatgcgccaccctaaagcagcgcatccgcaatgatgtGGAACGCGCAAt                                        <  1:4908385‑M2/12‑1 (MQ=255)
          aaatgcgccaccctaaagcagcgcatccgcaatgatgtGGAACGCGCAAt                                        <  1:4280990‑M2/12‑1 (MQ=255)
           aatgcgccaccctaaagcagcgcatccgcaatgatgtGGAACGCGCAAt                                        >  1:1719992‑M2/38‑49 (MQ=255)
            atgcgccaccctaaagcagcgcatccgcaatgatgtGGAACGCGCAATGa                                      <  1:1957377‑M2/14‑1 (MQ=255)
            atgcgccaccctaaagcagcgcatccgcaatgatgtGGAACGCGCAATGa                                      <  1:2952981‑M2/14‑1 (MQ=255)
             tgcgccaccctaaagcagcgcatccgcaatgatgtGGAACGCGCAATGAc                                     <  1:2586891‑M2/15‑1 (MQ=255)
             tgcgccaccctaaagcagcgcatccgcaatgatgtGGAACGCGCAATGAc                                     <  1:250488‑M2/15‑1 (MQ=255)
             tgcgccaccctaaagcagcgcatccgcaatgatgtGGAACGCGCAATGAc                                     <  1:2493988‑M2/15‑1 (MQ=255)
              gcgccaccctaaagcagcgcatccgcaatgatgtGGAACGCGCAATGAc                                     >  1:379784‑M2/35‑49 (MQ=255)
              gcgccaccctaaagcagcgcatccgcaatgatgtGGAACGCGCAATGAc                                     >  1:1892552‑M2/35‑49 (MQ=255)
              gcgccaccctaaagcagcgcatccgcaatgatgtGGAACGCGCAATGAc                                     >  1:4548318‑M2/35‑49 (MQ=255)
               cgccaccctaaagcagcgcatccgcaatgatgtGGAACGCGCAATGACGa                                   >  1:2103241‑M2/34‑50 (MQ=255)
                 ccaccctaaagcagcgcatccgcaatgatgtGGAACGCGCAATGACGatc                                 <  1:700639‑M2/19‑1 (MQ=255)
                  caccctaaagcagcgcatccgcaatgatgtGGAACGCGCAATGACGatca                                <  1:4503423‑M2/20‑1 (MQ=255)
                   accctaaagcagcgcatccgcaatgatgtGGAACGCGCAATGACGatcat                               <  1:2929676‑M2/21‑1 (MQ=255)
                   accctaaagcagcgcatccgcaatgatgtGGAACGCGCAATGACGatcat                               >  1:4305315‑M2/30‑50 (MQ=255)
                    ccctaaagcagcgcatccgcaatgatgtGGAACGCGCAATGACGatcatc                              <  1:434471‑M2/22‑1 (MQ=255)
                    ccctaaagcagcgcatccgcaatgatgtGGAACGCGCAATGACGatcatc                              >  1:565793‑M2/29‑50 (MQ=255)
                     cctaaagcagcgcatccgcaatgatgtGGAACGCGCAATGACGatcatc                              >  1:1624864‑M2/28‑49 (MQ=255)
                      ctaaagcagcgcatccgcaatgatgtGGAACGCGCAATGACGATcatcca                            <  1:4630150‑M2/24‑1 (MQ=255)
                         aagcagcgcatccgcaatgatgtGGAACGCGCAATGACGATcatccatc                          <  1:4710433‑M2/26‑1 (MQ=255)
                         aagcagcgcatccgcaatgatgtGGAACGCGCAATGACGATcatccatc                          <  1:1269381‑M2/26‑1 (MQ=255)
                         aagcagcgcatccgcaatgatgtGGAACGCGCAATGACGATCatccatcc                         >  1:585147‑M2/24‑50 (MQ=255)
                         aagcagcgcatccgcaatgatgtGGAACGCGCAATGACGATCatccatcc                         <  1:1871491‑M2/27‑1 (MQ=255)
                          agcagcgcatccgcaatgatgtGGAACGCGCAATGACGATCATCCATCct                        >  1:672701‑M2/23‑50 (MQ=255)
                           gcagcgcatccgcaatgatgtGGAACGCGCAATGACGATCATCCATCctc                       <  1:3431284‑M2/29‑1 (MQ=255)
                           gcagcgcatccgcaatgatgtGGAACGCGCAATGACGATCATCCATCctc                       <  1:1701226‑M2/29‑1 (MQ=255)
                            cagcgcatccgcaatgatgtGGAACGCGCAATGACGATCATCCATCctct                      <  1:5216683‑M2/30‑1 (MQ=255)
                            cagcgcatccgcaatgatgtGGAACGCGCAATGACGATCATCCATCctct                      <  1:1458619‑M2/30‑1 (MQ=255)
                               cgcatccgcaatgatgtGGAACGCGCAATGACGATCATCCATCCTCTGCt                   >  1:5455747‑M2/18‑50 (MQ=255)
                                gcatccgcaatgatgtGGAACGCGCAATGACGATCATCCATCCTCTGCtt                  <  1:4606217‑M2/34‑1 (MQ=255)
                                gcatccgcaatgatgtGGAACGCGCAATGACGATCATCCATCCTCTGCtt                  >  1:3130571‑M2/17‑50 (MQ=255)
                                gcatccgcaatgatgtGGAACGCGCAATGACGATCATCCATCCTCTGCtt                  <  1:4280795‑M2/34‑1 (MQ=255)
                                 catccgcaatgatgtGGAACGCGCAATGACGATCATCCATCCTCTGCTTg                 >  1:4831781‑M2/16‑50 (MQ=255)
                                  atccgcaatgatgtGGAACGCGCAATGACGATCATCCATCCTCTGCTTGc                >  1:3207746‑M2/15‑50 (MQ=255)
                                  atccgcaatgatgtGGAACGCGCAATGACGATCATCCATCCTCTGCTTGc                >  1:1602480‑M2/15‑50 (MQ=255)
                                  atccgcaatgatgtGGAACGCGCAATGACGATCATCCATCCTCTGCTTGc                >  1:732553‑M2/15‑50 (MQ=255)
                                  atccgcaatgatgtGGAACGCGCAATGACGATCATCCATCCTCTGCTTGc                <  1:3899152‑M2/36‑1 (MQ=255)
                                     cgcaatgatgtGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCag             <  1:2026207‑M2/39‑1 (MQ=255)
                                     cgcaatgatgtGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCag             >  1:1237753‑M2/12‑50 (MQ=255)
                                        aatgatgtGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCagtag          <  1:4962306‑M2/42‑1 (MQ=255)
                                         atgatgtGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCagtagt         <  1:3709277‑M2/43‑1 (MQ=255)
                                         atgatgtGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCagtagt         <  1:4648155‑M2/43‑1 (MQ=255)
                                         atgatgtGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCagtagt         <  1:2277636‑M2/43‑1 (MQ=255)
                                           gatgtGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGtagtag       >  1:3806035‑M2/6‑50 (MQ=255)
                                           gatgtGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGtagtag       >  1:3045030‑M2/6‑50 (MQ=255)
                                            atgtGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAgc      <  1:96304‑M2/46‑1 (MQ=255)
                                             tgtGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAgcg     <  1:1492133‑M2/47‑1 (MQ=255)
                                               tGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCa   >  1:2234010‑M2/2‑50 (MQ=255)
                                                ggAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCAc  >  1:2582110/1‑50 (MQ=255)
                                                |                                                 
TTTGAACAGGCGCTGGAAGCACTGGCTGAAAGCATGGGTATTGGTTGAGGAACGCGCAATGACGATCATCCATCCTCTGCTTGCCAGTAGTAGCGCAC  >  NC_000913/549569‑549666

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: