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 69 [0] [0] 68 [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)65 (0.990) 45/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

GATGTGATTAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCG  >  NC_000913/4596965‑4597063
                                                 |                                                 
tttagctatctgtataacctgaatgttagtactcattcttcctggtagta                                                   <  1:3476677‑M2/1‑1 (MQ=255)
  tagctatctgtataacctgaatgttagtactcattcttcctggtagtAAg                                                 >  1:930632‑M2/48‑50 (MQ=255)
   agctatctgtataacctgaatgttagtactcattcttcctggtagtAAg                                                 <  1:3216563‑M2/3‑1 (MQ=255)
    gctatctgtataacctgaatgttagtactcattcttcctggtagtAAgg                                                <  1:2721519‑M2/4‑1 (MQ=255)
    gctatctgtataacctgaatgttagtactcattcttcctggtagtAAgg                                                <  1:848888‑M2/4‑1 (MQ=255)
     ctatctgtataacctgaatgttagtactcattcttcctggtagtAAGGCg                                              <  1:3522069‑M2/6‑1 (MQ=255)
      tatctgtataacctgaatgttagtactcattcttcctggtagtAAGGc                                               >  1:4350017‑M2/44‑48 (MQ=255)
       atctgtataacctgaatgttagtactcattcttcctggtagtAAGGc                                               >  1:1428696‑M2/43‑47 (MQ=255)
       atctgtataacctgaatgttagtactcattcttcctggtagtAAGGc                                               >  1:2659333‑M2/43‑47 (MQ=255)
        tctgtataacctgaatgttagtactcattcttcctggtagtAAGGc                                               >  1:3662656‑M2/42‑46 (MQ=255)
        tctgtataacctgaatgttagtactcattcttcctggtagtAAGGc                                               >  1:4624518‑M2/42‑46 (MQ=255)
        tctgtataacctgaatgttagtactcattcttcctggtagtAAGGc                                               >  1:2268852‑M2/42‑46 (MQ=255)
        tctgtataacctgaatgttagtactcattcttcctggtagtAAGGc                                               >  1:1199610‑M2/42‑46 (MQ=255)
        tctgtataacctgaatgttagtactcattcttcctggtagtAAGGCga                                             <  1:3002225‑M2/7‑1 (MQ=255)
         ctgtataacctgaatgttagtactcattcttcctggtagtAAGGCGAGat                                          >  1:398609‑M2/41‑50 (MQ=255)
          tgtataacctgaatgttagtactcattcttcctggtagtAAGGCGAGat                                          >  1:1665185‑M2/40‑49 (MQ=255)
           gtataacctgaatgttagtactcattcttcctggtagtAAGGCGAGatta                                        >  1:838843‑M2/39‑50 (MQ=255)
            tataacctgaatgttagtactcattcttcctggtagtAAGGCgaga                                           <  1:4353887‑M2/9‑1 (MQ=255)
             ataacctgaatgttagtactcattcttcctggtagtAAGGCGAGattat                                       >  1:3431974‑M2/37‑49 (MQ=255)
              taacctgaatgttagtactcattcttcctggtagtAAGGCGAGattatt                                      >  1:1608375‑M2/36‑49 (MQ=255)
              taacctgaatgttagtactcattcttcctggtagtAAGGCGAGAttatta                                     <  1:4814789‑M2/15‑1 (MQ=255)
               aacctgaatgttagtactcattcttcctggtagtAAGGCGAGAttatta                                     >  1:4862450‑M2/35‑49 (MQ=255)
                acctgaatgttagtactcattcttcctggtagtAAGGCGAGAttatta                                     >  1:4263368‑M2/34‑48 (MQ=255)
                acctgaatgttagtactcattcttcctggtagtAAGGCGAGATTATTaaa                                   <  1:2435823‑M2/17‑1 (MQ=255)
                acctgaatgttagtactcattcttcctggtagtAAGGCGAGATTATTaa                                    <  1:5007602‑M2/16‑1 (MQ=255)
                 cctgaatgttagtactcattcttcctggtagtAAGGCGAGATTATTAAAg                                  >  1:4906470‑M2/33‑50 (MQ=255)
                 cctgaatgttagtactcattcttcctggtagtAAGGCGAGATTATTAAAg                                  >  1:1693761‑M2/33‑50 (MQ=255)
                 cctgaatgttagtactcattcttcctggtagtAAGGCGAGATTATTAAAg                                  >  1:4396691‑M2/33‑50 (MQ=255)
                   tgaatgttagtactcattcttcctggtagtAAGGCGAGATTATTAAAGtt                                >  1:183322‑M2/31‑50 (MQ=255)
                   tgaatgttagtactcattcttcctggtagtAAGGCGAGATTATTAAAGtt                                >  1:1029776‑M2/31‑50 (MQ=255)
                     aatgttagtactcattcttcctggtagtAAGGCGAGATTATTAAAGTTg                               <  1:2290973‑M2/21‑1 (MQ=255)
                     aatgttagtactcattcttcctggtagtAAGGCGAGATTATTAAAGTTg                               <  1:3412423‑M2/21‑1 (MQ=255)
                     aatgttagtactcattcttcctggtagtAAGGCGAGATTATTAAAGTTGc                              >  1:2964053‑M2/29‑50 (MQ=255)
                      atgttagtactcattcttcctggtagtAAGGCGAGATTATTAAAGTTGcc                             <  1:3281214‑M2/23‑1 (MQ=255)
                      atgttagtactcattcttcctggtagtAAGGCGAGATTATTAAAGTTGcc                             <  1:3336892‑M2/23‑1 (MQ=255)
                       tgttagtactcattcttcctggtagtAAGGCGAGATTATTAAAGTTGCCa                            <  1:348527‑M2/24‑1 (MQ=255)
                       tgttagtactcattcttcctggtagtAAGGCGAGATTATTAAAGTTGCCa                            <  1:5134408‑M2/24‑1 (MQ=255)
                       tgttagtactcattcttcctggtagtAAGGCGAGATTATTAAAGTTGCCa                            <  1:3793052‑M2/24‑1 (MQ=255)
                         ttagtactcattcttcctggtagtAAGGCGAGATTATTAAAGTTGCCATg                          <  1:4242339‑M2/26‑1 (MQ=255)
                          tagtactcattcttcctggtagtAAGGCGAGATTATTAAAGTTGCCATg                          >  1:1987436‑M2/24‑49 (MQ=255)
                          tagtactcattcttcctggtagtAAGGCGAGATTATTAAAGTTGCCATg                          >  1:662857‑M2/24‑49 (MQ=255)
                            gtactcattcttcctggtagtAAGGCGAGATTATTAAAGTTGCCATGCa                        <  1:2449755‑M2/28‑1 (MQ=255)
                            gtactcattcttcctggtagtAAGGCGAGATTATTAAAGTTGCCATGCa                        <  1:2344051‑M2/28‑1 (MQ=255)
                             tactcattcttcctggtagtAAGGCGAGATTATTAAAGTTGCCATGCAGc                      <  1:4290637‑M2/30‑1 (MQ=255)
                              actcattcttcctggtagtAAGGCGAGATTATTAAAGTTGCCATGCAGCg                     >  1:2905131‑M2/20‑50 (MQ=255)
                              actcattcttcctggtagtAAGGCGAGATTATTAAAGTTGCCATGCAGCg                     >  1:1321888‑M2/20‑50 (MQ=255)
                               ctcattcttcctggtagtAAGGCGAGATTATTAAAGTTGCCATGCAGCGt                    >  1:751894‑M2/19‑50 (MQ=255)
                               ctcattcttcctggtagtAAGGCGAGATTATTAAAGTTGCCATGCAGCGt                    >  1:2283656‑M2/19‑50 (MQ=255)
                                 cattcttcctggtagtAAGGCGAGATTATTAAAGTTGCCATGCAGCGTcc                  >  1:2916521‑M2/17‑50 (MQ=255)
                                  attcttcctggtagtAAGGCGAGATTATTAAAGTTGCCATGCAGCGTcc                  <  1:2152017‑M2/34‑1 (MQ=255)
                                  attcttcctggtagtAAGGCGAGATTATTAAAGTTGCCATGCAGCGTcc                  <  1:5084346‑M2/34‑1 (MQ=255)
                                  attcttcctggtagtAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCg                 >  1:411937‑M2/16‑50 (MQ=255)
                                   ttcttcctggtagtAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCgg                >  1:4455025‑M2/15‑50 (MQ=255)
                                     cttcctggtagtAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCggg               >  1:4413213‑M2/13‑49 (MQ=255)
                                      ttcctggtagtAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCgggg              >  1:4005650‑M2/12‑49 (MQ=255)
                                      ttcctggtagtAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCgggg              >  1:3885144‑M2/12‑49 (MQ=255)
                                         ctggtagtAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAgt          <  1:203750‑M2/42‑1 (MQ=255)
                                           ggtagtAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAgtgt        >  1:3163163‑M2/7‑50 (MQ=255)
                                            gtagtAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAgtgt        <  1:1060867‑M2/44‑1 (MQ=255)
                                            gtagtAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGtt       <  1:3954469‑M2/45‑1 (MQ=255)
                                            gtagtAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGtt       <  1:2188239‑M2/45‑1 (MQ=255)
                                            gtagtAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGtt       >  1:253187‑M2/6‑50 (MQ=255)
                                             tagtAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTg      <  1:1863910‑M2/46‑1 (MQ=255)
                                               gtAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTggg    <  1:4934124‑M2/48‑1 (MQ=255)
                                               gtAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTggg    <  1:441033‑M2/48‑1 (MQ=255)
                                                tAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGc   <  1:3596754‑M2/49‑1 (MQ=255)
                                                 aaGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGgcg  >  1:2914619/1‑50 (MQ=255)
                                                 aaGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGgcg  >  1:1302772/1‑50 (MQ=255)
                                                 |                                                 
GATGTGATTAACAATCGTGAGTTAATTAAAGGCGTAATTACTTTCTGATAAGGCGAGATTATTAAAGTTGCCATGCAGCGTCCGGGGAAGTGTTGGGCG  >  NC_000913/4596965‑4597063

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: