Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 2,558,699 Δ6,790 bp intZ[eutA] intZ, yffL, yffM, yffN, yffO, yffP, yffQ, yffR, yffS, [eutA]

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 2558699 2565488 6790 68 [0] [0] 68 intZ–[eutA] intZ,yffL,yffM,yffN,yffO,yffP,yffQ,yffR,yffS,[eutA]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 25586980 (0.000)50 (0.900) 29/52 0.0 100% intergenic (‑19/‑160) eutB/intZ ethanolamine ammonia‑lyase subunit alpha/CPZ‑55 prophage; putative phage integrase IntZ
?NC_000913 2565489 = 0 (0.000)coding (1396/1404 nt) eutA reactivating factor for ethanolamine ammonia‑lyase EutBC

AGCGCCACCGCTGGCGAATACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACC  >  NC_000913/2565455‑2565523
                                  |                                  
ggtctttagtttcataagtcgttccctcaggaagg                                    <  1:1605960‑M2/1‑1 (MQ=255)
ggtctttagtttcataagtcgttccctcaggaagg                                    <  1:3408691‑M2/1‑1 (MQ=255)
 gtctttagtttcataagtcgttccctcaggaagga                                   >  1:6653517‑M2/34‑35 (MQ=255)
 gtctttagtttcataagtcgttccctcaggaagga                                   >  1:5280303‑M2/34‑35 (MQ=255)
 gtctttagttgcataagtcgttccctcaggaagga                                   <  1:7136229‑M2/2‑1 (MQ=255)
   ctttagtttcataagtcgttccctcaggaagGaaa                                 <  1:1828535‑M2/4‑1 (MQ=255)
    tttagtttcataagtcgttccctcaggaagGAAAt                                >  1:3662334‑M2/31‑35 (MQ=255)
      tagtttcataagtcgttccctcaggaagGAAATGc                              >  1:7622232‑M2/29‑35 (MQ=255)
      tagtttcataagtcgttccctcaggaagGAAATGc                              <  1:5379779‑M2/7‑1 (MQ=255)
      tagtttcataagtcgttccctcaggaagGAAATGc                              >  1:4759701‑M2/29‑35 (MQ=255)
      tagtttcataagtcgttccctcaggaagGAAATGc                              >  1:3165566‑M2/29‑35 (MQ=255)
       agtttcataagtcgttccctcaggaagGAAATGCg                             >  1:5374317‑M2/28‑35 (MQ=255)
        gtttcataagtcgttccctcaggaagGAAATGCGa                            >  1:78094‑M2/27‑35 (MQ=255)
        gtttcataagtcgttccctcaggaagGAAATGCGa                            >  1:48797‑M2/27‑35 (MQ=255)
        gtttcataagtcgttccctcaggaagGAAATGCGa                            >  1:305201‑M2/27‑35 (MQ=255)
        gtttcataagtcgttccctcaggaagGAAATGCGa                            >  1:2822707‑M2/27‑35 (MQ=255)
         tttcataagtcgttccctcaggaagGAAATGCGAg                           <  1:5774757‑M2/10‑1 (MQ=255)
         tttcataagtcgttccctcaggaagGAAATGCGAg                           >  1:3888368‑M2/26‑35 (MQ=255)
          ttcataagtcgttccctcaggaagGAAATGCGAGt                          <  1:5079985‑M2/11‑1 (MQ=255)
          ttcataagtcgttccctcaggaagGAAATGCGAGt                          <  1:6996141‑M2/11‑1 (MQ=255)
           tcataagtcgttccctcaggaagGAAATGCGAGTg                         >  1:6599761‑M2/24‑35 (MQ=255)
            cataagtcgttccctcaggaagGAAATGCGAGTGa                        >  1:718134‑M2/23‑35 (MQ=255)
            cataagtcgttccctcaggaagGAAATGCGAGTGa                        >  1:7161038‑M2/23‑35 (MQ=255)
            cataagtcgttccctcaggaagGAAATGCGAGTGa                        >  1:1538603‑M2/23‑35 (MQ=255)
            cataagtcgttccctcaggaagGAAATGCGAGTGa                        >  1:3727566‑M2/23‑35 (MQ=255)
            cataagtcgttccctcaggaagGAAATGCGAGTGa                        >  1:6530784‑M2/23‑35 (MQ=255)
             ataagtcgttccctcaggaagGAAATGCGAGTGAt                       <  1:5712979‑M2/14‑1 (MQ=255)
             ataagtcgttccctcaggaagGAAATGCGAGTGAt                       >  1:6220289‑M2/22‑35 (MQ=255)
              taagtcgttccctcaggaagGAAATGCGAGTGAtt                      >  1:7547320‑M2/21‑35 (MQ=255)
              taagtcgttccctcaggaagGAAATGCGAGTGAtt                      >  1:5519552‑M2/21‑35 (MQ=255)
              taagtcgttccctcaggaagGAAATGCGAGTGAtt                      <  1:3052376‑M2/15‑1 (MQ=255)
               aagtcgttccctcaggaagGAAATGCGAGTGAttt                     <  1:3822816‑M2/16‑1 (MQ=255)
               aagtcgttccctcaggaagGAAATGCGAGTGAttt                     <  1:1444283‑M2/16‑1 (MQ=255)
               aagtcgttccctcaggaagGAAATGCGAGTGAttt                     <  1:2041532‑M2/16‑1 (MQ=255)
               aagtcgttccctcaggaagGAAATGCGAGTGAttt                     >  1:5514449‑M2/20‑35 (MQ=255)
               aagtcgttccctcaggaagGAAATGCGAGTGAttt                     >  1:200748‑M2/20‑35 (MQ=255)
                agtcgttccctcaggaagGAAATGCGAGTGATTtc                    >  1:6188939‑M2/19‑35 (MQ=255)
                 gtcgttccctcaggaagGAAATGCGAGTGATTtca                   <  1:3042867‑M2/18‑1 (MQ=255)
                 gtcgttccctcaggaagGAAATGCGAGTGATTtca                   <  1:1784330‑M2/18‑1 (MQ=255)
                 gtcgttccctcaggaagGAAATGCGAGTGATTtca                   <  1:4338371‑M2/18‑1 (MQ=255)
                  tcgttccctcaggaagGAAATGCGAGTGATTtcac                  >  1:1194935‑M2/17‑35 (MQ=255)
                  tcgttccctcaggaagGAAATGCGAGTGATTtcac                  <  1:2815594‑M2/19‑1 (MQ=255)
                  tcgttccctcaggaagGAAATGCGAGTGATTtcac                  <  1:2513426‑M2/19‑1 (MQ=255)
                  tcgttccctcaggaagGAAATGCGAGTGATTtcac                  <  1:4356522‑M2/19‑1 (MQ=255)
                    gttccctcaggaagGAAATGCGAGTGATTtcaccg                <  1:1324314‑M2/21‑1 (MQ=255)
                     ttccctcaggaagGAAATGCGAGTGATttcaccgc               >  1:2804091‑M2/14‑34 (MQ=255)
                     ttccctcaggaagGAAATGCGAGTGATTtcaccgt               <  1:4557687‑M2/22‑1 (MQ=255)
                     ttccctcaggaagGAAATGCGAGTGATTtcaccgt               >  1:340540‑M2/14‑35 (MQ=255)
                       ccctcaggaagGAAATGCGAGTGATTtcaccgtca             >  1:1579104‑M2/12‑35 (MQ=255)
                        cctcaggaagGAAATGCGAGTGATTtcaccgtcac            <  1:6269760‑M2/25‑1 (MQ=255)
                          tcaggaagGAAATGCGAGTGATTtcaccgtcaccg          >  1:2151343‑M2/9‑35 (MQ=255)
                          tcaggaagGAAATGCGAGTGATTtcaccgtcaccg          >  1:2387113‑M2/9‑35 (MQ=255)
                           caggaagGAAATGCGAGTGATTTCACCGTCACCgg         <  1:7215175‑M2/28‑1 (MQ=255)
                            aggaagGAAATGCGAGTGATTTCACCGTCACCGGc        >  1:4642258‑M2/7‑35 (MQ=255)
                             ggaagGAAATGCGAGTGATTTCACCGTCACCGGca       >  1:4817878‑M2/6‑35 (MQ=255)
                             ggaagGAAATGCGAGTGATTTCACCGTCACCGGca       <  1:752545‑M2/30‑1 (MQ=255)
                              gaagGAAATGCGAGTGATTTCACCGTCACCGGcac      >  1:6208858‑M2/5‑35 (MQ=255)
                              gaagGAAATGCGAGTGATTTCACCGTCACCGGcac      >  1:1851375‑M2/5‑35 (MQ=255)
                              gaagGAAATGCGAGTGATTTCACCGTCACCGGcac      >  1:5480219‑M2/5‑35 (MQ=255)
                              gaagGAAATGCGAGTGATTTCACCGTCACCGGcac      >  1:6803683‑M2/5‑35 (MQ=255)
                              gaagGAAATGCGAGTGATTTCACCGTCACCGGcac      >  1:4331919‑M2/5‑35 (MQ=255)
                              gaagGAAATGCGAGTGATTTCACCGTCACCGGcac      <  1:6213086‑M2/31‑1 (MQ=255)
                               aagGAAATGCGAGTGATTTCACCGTCACCGGcaca     >  1:6808745‑M2/4‑35 (MQ=255)
                               aagGAAATGCGAGTGATTTCACCGTCACCGGcaca     >  1:7238788‑M2/4‑35 (MQ=255)
                                 gGAAATGCGAGTGATTTCACCGTCACCGGCACAAc   >  1:4413614‑M2/2‑35 (MQ=255)
                                 gGAAATGCGAGTGATTTCACCGTCACCGGCACAAc   >  1:568499‑M2/2‑35 (MQ=255)
                                 gGAAATGCGAGTGATTTCACCGTCACCGGCACAAc   >  1:700714‑M2/2‑35 (MQ=255)
                                  gaaATGCGAGTGATTTCACCGTCACCGGCACAAcc  <  1:4265247/35‑1 (MQ=255)
                                  |                                  
AGCGCCACCGCTGGCGAATACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACC  >  NC_000913/2565455‑2565523

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ 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.