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 34 [0] [0] 33 [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)31 (1.060) 26/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

GATGTGATTAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGC  >  NC_000913/4596965‑4597062
                                                 |                                                
tttagctatctgtataacctgaatgttagtactcattcttcctggtagta                                                  <  1:1285326‑M2/1‑1 (MQ=255)
     ctatctgtataacctgaatgttagtactcattcttcctggtagtAAGGc                                              >  1:1261282‑M2/45‑49 (MQ=255)
          tgtataacctgaatgttagtactcattcttcctggtagtAAGGCGAGat                                         >  1:2035560‑M2/40‑49 (MQ=255)
            tataacctgaatgttagtactcattcttcctggtagtAAGGCGAGattat                                      >  1:750454‑M2/38‑50 (MQ=255)
            tataacctgaatgttagtactcattcttcctggtagtAAGGCGAGattat                                      >  1:291748‑M2/38‑50 (MQ=255)
            tataacctgaatgttagtactcattcttcctggtagtAAGGCGAGatta                                       <  1:984962‑M2/12‑1 (MQ=255)
            tataacctgaatgttagtactcattcttcctggtagtAAGGCGAGat                                         <  1:1028474‑M2/10‑1 (MQ=255)
              taacctgaatgttagtactcattcttcctggtagtAAGGCGAGattatt                                     <  1:226440‑M2/14‑1 (MQ=255)
                 cctgaatgttagtactcattcttcctggtagtAAGGCGAGATTATTAAAg                                 <  1:657650‑M2/18‑1 (MQ=255)
                  ctgaatgttagtactcattcttcctggtagtAAGGCGAGATTATTAAAGt                                <  1:398327‑M2/19‑1 (MQ=255)
                   tgaatgttagtactcattcttcctggtagtAAGGCGAGATTATTAAAGtt                               >  1:1921462‑M2/31‑50 (MQ=255)
                     aatgttagtactcattcttcctggtagtAAGGCGAGATTATTAAAGTTg                              <  1:841406‑M2/21‑1 (MQ=255)
                     aatgttagtactcattcttcctggtagtAAGGCGAGATTATTAAAGTTGc                             >  1:1278019‑M2/29‑50 (MQ=255)
                     aatgttagtactcattcttcctggtagtAAGGCGAGATTATTAAAGTTGc                             >  1:723326‑M2/29‑50 (MQ=255)
                      atgttagtactcattcttcctggtagtAAGGCGAGATTATTAAAGTTGc                             >  1:1972760‑M2/28‑49 (MQ=255)
                         ttagtactcattcttcctggtagtAAGGCGAGATTATTAAAGTTGCCAt                          <  1:940057‑M2/25‑1 (MQ=255)
                         ttagtactcattcttcctggtagtAAGGCGAGATTATTAAAGTTGCCAt                          <  1:1630561‑M2/25‑1 (MQ=255)
                         ttagtactcattcttcctggtagtAAGGCGAGATTATTAAAGTTGCCATg                         >  1:2033095‑M2/25‑50 (MQ=255)
                            gtactcattcttcctggtagtAAGGCGAGATTATTAAAGTTGCCATGCa                       <  1:1468067‑M2/28‑1 (MQ=255)
                            gtactcattcttcctggtagtAAGGCGAGATTATTAAAGTTGCCATGCAg                      <  1:321554‑M2/29‑1 (MQ=255)
                               ctcattcttcctggtagtAAGGCGAGATTATTAAAGTTGCCATGCAGCGt                   >  1:1281035‑M2/19‑50 (MQ=255)
                                tcattcttcctggtagtAAGGCGAGATTATTAAAGTTGCCATGCAGCGTc                  >  1:1630957‑M2/18‑50 (MQ=255)
                                  attcttcctggtagtAAGGCGAGATTATTAAAGTTGCCATGCAGCGTcc                 <  1:1850667‑M2/34‑1 (MQ=255)
                                   ttcttcctggtagtAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCgg               >  1:1082802‑M2/15‑50 (MQ=255)
                                      ttcctggtagtAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCgggg             >  1:902518‑M2/12‑49 (MQ=255)
                                        cctggtagtAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGaa           <  1:1972730‑M2/40‑1 (MQ=255)
                                            gtagtAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAgtgt       <  1:428435‑M2/44‑1 (MQ=255)
                                            gtagtAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGtt      >  1:735533‑M2/6‑50 (MQ=255)
                                            gtagtAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGtt      <  1:1611936‑M2/45‑1 (MQ=255)
                                            gtagtAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGtt      <  1:1943192‑M2/45‑1 (MQ=255)
                                               gtAAGGCGAGATTATTAAAGTTGCCATGCAGGGTCCGGGGAAGTGttgcg   >  1:1210396‑M2/3‑48 (MQ=255)
                                                tAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGc  >  1:1155278‑M2/2‑50 (MQ=255)
                                                tAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGc  <  1:2071333‑M2/49‑1 (MQ=255)
                                                 |                                                
GATGTGATTAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGC  >  NC_000913/4596965‑4597062

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑