Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 4,555,491 Δ41,523 bp [yjiC]lgoD 37 genes

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 4555491 4597013 41523 58 [0] [0] 56 [yjiC]–lgoD [yjiC],ytiC,ytiD,idlP,iraD,hypT,iadA,yjiG,yjiH,kptA,yjiJ,yjiK,ytiA,yjiL,yjiM,yjiN,mdtM,rpnD,yjiR,yjiS,yjiT,yjiV,mcrC,mcrB,symE,symR,hsdS,hsdM,hsdR,mrr,yjiA,yjiX,btsT,tsr,lgoT,lgoR,lgoD

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 45554900 (0.000)53 (0.930) 41/96 0.0 100% coding (831/831 nt) yjiC uncharacterized protein YjiC
?NC_000913 4597014 = 0 (0.000)intergenic (+2/+136) lgoD/opgB L‑galactonate oxidoreductase/phosphoglycerol transferase I

GTGATTAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCG  >  NC_000913/4596968‑4597063
                                              |                                                 
agctatctgtataacctgaatgttagtactcattcttcctggtagtAAgg                                                >  1:1824331‑M2/47‑50 (MQ=255)
agctatctgtataacctgaatgttagtactcattcttcctggtagtAAg                                                 <  1:4327886‑M2/3‑1 (MQ=255)
 gctatctgtataacctgaatgttagtactcattcttcctggtagtAAgg                                                <  1:1831716‑M2/4‑1 (MQ=255)
 gctatctgtataacctgaatgttagtactcattcttcctggtagtAAgg                                                <  1:4429810‑M2/4‑1 (MQ=255)
  ctatctgtataacctgaatgttagtactcattcttcctggtagtAAGGc                                               >  1:341059‑M2/45‑49 (MQ=255)
     tctgtataacctgaatgttagtactcattcttcctggtagtAAGGc                                               >  1:3076843‑M2/42‑46 (MQ=255)
     tctgtataacctgaatgttagtactcattcttcctggtagtAAGGCgag                                            <  1:792920‑M2/8‑1 (MQ=255)
      ctgtataacctgaatgttagtactcattcttcctggtagtAAGGCGAGat                                          >  1:2867484‑M2/41‑50 (MQ=255)
       tgtataacctgaatgttagtactcattcttcctggtagtAAGGCGAGat                                          >  1:3240316‑M2/40‑49 (MQ=255)
       tgtataacctgaatgttagtactcattcttcctggtagtAAGGCGAGat                                          >  1:2758937‑M2/40‑49 (MQ=255)
         tataacctgaatgttagtactcattcttcctggtagtATGGCgaga                                           <  1:3125057‑M2/9‑1 (MQ=255)
         tataacctgaatgttagtactcattcttcctggtagtAAGGCgaga                                           <  1:2415651‑M2/9‑1 (MQ=255)
          ataacctgaatgttagtactcattcttcctggtagtAAGGCGAGattat                                       >  1:169932‑M2/37‑49 (MQ=255)
           taacctgaatgttagtactcattcttcctggtagtAAGGCGAGattatt                                      <  1:3777388‑M2/14‑1 (MQ=255)
             acctgaatgttagtactcattcttcctggtagtAAGGCGAGATTATTaa                                    <  1:2185047‑M2/16‑1 (MQ=255)
               ctgaatgttagtactcattcttcctggtagtAAGGCGAGATTATTAAAGt                                 >  1:1088428‑M2/32‑50 (MQ=255)
                tgaatgttagtactcattcttcctggtagtAAGGCGAGATTATTAAAGtt                                >  1:698588‑M2/31‑50 (MQ=255)
                 gaatgttagtactcattcttcctggtagtAAGGCGAGATTATTAAAGTTg                               >  1:161682‑M2/30‑50 (MQ=255)
                  aatgttagtactcattcttcctggtagtAAGGCGAGATTATTAAAGTTg                               <  1:2829865‑M2/21‑1 (MQ=255)
                   atgttagtactcattcttcctggtagtAAGGCGAGATTATTAAAGTTGc                              >  1:1133797‑M2/28‑49 (MQ=255)
                    tgttagtactcattcttcctggtagtAAGGCGAGATTATTAAAGTTGCCa                            >  1:4297706‑M2/27‑50 (MQ=255)
                     gttagtactcattcttcctggtagtAAGGCGAGATTATTAAAGTTGCCAt                           >  1:2692976‑M2/26‑50 (MQ=255)
                     gttagtactcattcttcctggtagtAAGGCGAGATTATTAAAGTTGCCAt                           >  1:4372880‑M2/26‑50 (MQ=255)
                      ttagtactcattcttcctggtagtAAGGCGAGATTATTAAAGTTGCCAt                           <  1:3681732‑M2/25‑1 (MQ=255)
                      ttagtactcattcttcctggtagtAAGGCGAGATTATTAAAGTTGCCAt                           <  1:942997‑M2/25‑1 (MQ=255)
                      ttagtactcattcttcctggtagtAAGGCGAGATTATTAAAGTTGCCAt                           <  1:3276358‑M2/25‑1 (MQ=255)
                      ttagtactcattcttcctggtagtAAGGCGAGATTATTAAAGTTGCCATg                          <  1:387307‑M2/26‑1 (MQ=255)
                        agtactcattcttcctggtagtAAGGCGAGATTATTAAAGTTGCCATg                          >  1:3351663‑M2/23‑48 (MQ=255)
                        agtactcattcttcctggtagtAAGGCGAGATTATTAAAGTTGCCATGc                         <  1:1101516‑M2/27‑1 (MQ=255)
                         gtactcattcttcctggtagtAAGGCGAGATTATTAAAGTTGCCATGCa                        <  1:3656799‑M2/28‑1 (MQ=255)
                         gtactcattcttcctggtagtAAGGCGAGATTATTAAAGTTGCCATGCa                        <  1:2461221‑M2/28‑1 (MQ=255)
                         gtactcattcttcctggtagtAAGGCGAGATTATTAAAGTTGCCATGCAg                       >  1:587540‑M2/22‑50 (MQ=255)
                          tactcattcttcctggtagtAAGGCGAGATTATTAAAGTTGCCATGCAGc                      <  1:3363844‑M2/30‑1 (MQ=255)
                          tactcattcttcctggtagtAAGGCGAGATTATTAAAGTTGCCATGCAGc                      <  1:2401961‑M2/30‑1 (MQ=255)
                          tactcattcttcctggtagtAAGGCGAGATTATTAAAGTTGCCATGCAGc                      <  1:3884060‑M2/30‑1 (MQ=255)
                           actcattcttcctggtagtAAGGCGAGATTATTAAAGTTGCCATGCAGCg                     >  1:3890202‑M2/20‑50 (MQ=255)
                            ctcattcttcctggtagtAAGGCGAGATTATTAAAGTTGCCATGCAGCg                     <  1:2621722‑M2/31‑1 (MQ=255)
                            ctcattcttcctggtagtAAGGCGAGATTATTAAAGTTGCCATGCAGCGt                    <  1:921667‑M2/32‑1 (MQ=255)
                            ctcattcttcctggtagtAAGGCGAGATTATTAAAGTTGCCATGCAGCGt                    <  1:833345‑M2/32‑1 (MQ=255)
                              cattcttcctggtagtAAGGCGAGATTATTAAAGTTGCCATGCAGCGTcc                  >  1:757030‑M2/17‑50 (MQ=255)
                                ttcttcctggtagtAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCgg                <  1:1281699‑M2/36‑1 (MQ=255)
                                ttcttcctggtagtAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCgg                <  1:4191028‑M2/36‑1 (MQ=255)
                                  cttcctggtagtAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCggg               >  1:1525887‑M2/13‑49 (MQ=255)
                                    tcctggtagtAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCgggg              >  1:517743‑M2/11‑48 (MQ=255)
                                     cctggtagtAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCgggg              <  1:786424‑M2/38‑1 (MQ=255)
                                     cctggtagtAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGa             <  1:3981782‑M2/39‑1 (MQ=255)
                                     cctggtagtAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAg           >  1:2429356‑M2/10‑50 (MQ=255)
                                     cctggtagtAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAg           >  1:2911298‑M2/10‑50 (MQ=255)
                                      ctggtagtAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAgt          <  1:3418989‑M2/42‑1 (MQ=255)
                                       tggtagtAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAgtg         >  1:1959467‑M2/8‑50 (MQ=255)
                                         gtagtAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGtt       <  1:4102570‑M2/45‑1 (MQ=255)
                                          tagtAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTg      <  1:4049204‑M2/46‑1 (MQ=255)
                                           agtAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTgg     >  1:3901378‑M2/4‑50 (MQ=255)
                                             tAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGc   <  1:1217912‑M2/49‑1 (MQ=255)
                                              aaGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGgcg  >  1:1832013/1‑50 (MQ=255)
                                              aaGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGgcg  >  1:4322485/1‑50 (MQ=255)
                                              |                                                 
GTGATTAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCG  >  NC_000913/4596968‑4597063

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: