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 77 [0] [0] 75 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)63 (0.950) 49/82 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

AGCACCCAGCCACCCAGCGCCACCGCTGGCGAATACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAG  >  NC_000913/2565440‑2565537
                                                 |                                                
attgccgaacaatgtggtctttagtttcataagtcgttccctcaggaagg                                                  <  1:4108756‑M2/1‑1 (MQ=255)
 ttgccgaacaatgtggtctttagtttcataagtcgttccctcaggaagg                                                  >  1:958506‑M2/49‑49 (MQ=255)
   gccgaacaatgtggtctttagtttcataagtcgttccctcaggaaggaa                                                <  1:3603634‑M2/3‑1 (MQ=255)
     cgaacaatgtggtctttagtttcataagtcgttccctcaggaagGAAATg                                             >  1:4690223‑M2/45‑50 (MQ=255)
     cgaacaatgtggtctttagtttcataagtcgttccctcaggaagGAAATg                                             <  1:4682121‑M2/6‑1 (MQ=255)
      gaacaatgtggtctttagtttcataagtcgttccctcaggaagGAAATGc                                            <  1:2276‑M2/7‑1 (MQ=255)
      gaacaatgtggtctttagtttcataagtcgttccctcaggaagGAAATGc                                            <  1:5054598‑M2/7‑1 (MQ=255)
      gaacaatgtggtctttagtttcataagtcgttccctcaggaagGAAATGc                                            >  1:415989‑M2/44‑50 (MQ=255)
      gaacaatgtggtctttagtttcataagtcgttccctcaggaagGAAATGc                                            >  1:393737‑M2/44‑50 (MQ=255)
       aacaatgtggtctttagtttcataagtcgttccctcaggaagGAAATGc                                            >  1:3046309‑M2/43‑49 (MQ=255)
       aacaatgtggtctttagtttcataagtcgttccctcaggaagGAAATGc                                            >  1:3551289‑M2/43‑49 (MQ=255)
        acaatgtggtctttagtttcataagtcgttccctcaggaagGAAATGc                                            >  1:109816‑M2/42‑48 (MQ=255)
        acaatgtggtctttagtttcataagtcgttccctcaggaagGAAATGCGa                                          <  1:1005955‑M2/9‑1 (MQ=255)
         caatgtggtctttagtttcataagtcgttccctcaggaagGAAATGCGAg                                         >  1:180020‑M2/41‑50 (MQ=255)
          aatgtggtctttagtttcataagtcgttccctcaggaagGAAATGCGAGt                                        >  1:807446‑M2/40‑50 (MQ=255)
          aatgtggtctttagtttcataagtcgttccctcaggaagGAAATGCGAGt                                        >  1:1083924‑M2/40‑50 (MQ=255)
          aatgtggtctttagtttcataagtcgttccctcaggaagGAAATGCGAGt                                        <  1:2473033‑M2/11‑1 (MQ=255)
           atgtggtctttagtttcataagtcgttccctcaggaagGAAATGCGAGTg                                       <  1:3057417‑M2/12‑1 (MQ=255)
            tgtggtctttagtttcataagtcgttccctcaggaagGAAATGCGAGTGa                                      >  1:1829197‑M2/38‑50 (MQ=255)
             gtggtctttagtttcataagtcgttccctcaggaagGAAATGCGAGTGAt                                     >  1:2273441‑M2/37‑50 (MQ=255)
              tggtctttagtttcataagtcgttccctcaggaagGAAATGCGAGTGAtt                                    >  1:1545931‑M2/36‑50 (MQ=255)
               ggtctttagtttcataagtcgttccctcaggaagGAAATGCGAGTGAttt                                   >  1:896514‑M2/35‑50 (MQ=255)
                gtctttagtttcataagtcgttccctcaggaagGAAATGCGAGTGAttt                                   <  1:5682776‑M2/16‑1 (MQ=255)
                gtctttagtttcataagtcgttccctcaggaagGAAATGCGAGTGAttt                                   >  1:5768149‑M2/34‑49 (MQ=255)
                gtctttagtttcataagtcgttccctcaggaagGAAATGCGAGTGAttt                                   >  1:962278‑M2/34‑49 (MQ=255)
                gtctttagtttcataagtcgttccctcaggaagGAAATGCGAGTGATTtc                                  <  1:1164755‑M2/17‑1 (MQ=255)
                 tctttagtttcataagtcgttccctcaggaagGAAATGCGAGTGATTtca                                 <  1:2223458‑M2/18‑1 (MQ=255)
                 tctttagtttcataagtcgttccctcaggaagGAAATGCGAGTGATTtca                                 <  1:3856642‑M2/18‑1 (MQ=255)
                 tctttagtttcataagtcgttccctcaggaagGAAATGCGAGTGATTtca                                 >  1:1050934‑M2/33‑50 (MQ=255)
                  ctttagtttcataagtcgttccctcaggaagGAAATGCGAGTGATTtcac                                >  1:1674844‑M2/32‑50 (MQ=255)
                   tttagtttcataagtcgttccctcaggaagGAAATGCGAGTGATTtcacc                               >  1:2087054‑M2/31‑50 (MQ=255)
                     tagtttcataagtcgttccctcaggaagGAAATGCGAGTGATTtcaccgt                             >  1:359631‑M2/29‑50 (MQ=255)
                     tagtttcataagtcgttccctcaggaagGAAATGCGAGTGATTtcaccg                              <  1:3437415‑M2/21‑1 (MQ=255)
                     tagtttcataagtcgttccctcaggaagGAAATGCGAGTGATTtcacc                               <  1:5974393‑M2/20‑1 (MQ=255)
                     tagtttcataagtcgttccctcaggaagGAAATGCGAGTGATTtcacc                               <  1:750168‑M2/20‑1 (MQ=255)
                     tagtttcataagtcgttccctcaggaagGAAATGCGAGTGATTtcacc                               <  1:159888‑M2/20‑1 (MQ=255)
                     tagtttcataagtcgttccctcaggaagGAAATGCGAGTGATTtcac                                <  1:810492‑M2/19‑1 (MQ=255)
                      agtttcataagtcgttccctcaggaagGAAATGCGAGTGATTtcaccgt                             <  1:3395078‑M2/22‑1 (MQ=255)
                       gtttcataagtcgttccctcaggaagGAAATGCGAGTGATTtcaccgtca                           <  1:2716893‑M2/24‑1 (MQ=255)
                       gtttcataagtcgttccctcaggaagGAAATGCGAGTGATTtcaccgtca                           <  1:5896324‑M2/24‑1 (MQ=255)
                       gtttcataagtcgttccctcaggaagGAAATGCGAGTGATTtcaccgtc                            <  1:180098‑M2/23‑1 (MQ=255)
                        tttcataagtcgttccctcaggaagGAAATGCGAGTGATTtcaccgtcac                          >  1:4307704‑M2/26‑50 (MQ=255)
                         ttcataagtcgttccctcaggaagGAAATGCGAGTGATTtcaccgtcacc                         >  1:3160671‑M2/25‑50 (MQ=255)
                         ttcataagtcgttccctcaggaagGAAATGCGAGTGATTtcaccgtcacc                         >  1:1249130‑M2/25‑50 (MQ=255)
                          tcataagtcgttccctcaggaagGAAATGCGAGTGATTtcaccgtcacc                         <  1:4236663‑M2/26‑1 (MQ=255)
                           cataagtcgttccctcaggaagGAAATGCGAGTGATTtcaccgtcaccg                        <  1:265654‑M2/27‑1 (MQ=255)
                           cataagtcgttccctcaggaagGAAATGCGAGTGATTtcaccgtcaccg                        <  1:605479‑M2/27‑1 (MQ=255)
                           cataagtcgttccctcaggaagGAAATGCGAGTGATTTCACCGTCACCgg                       >  1:1106712‑M2/23‑50 (MQ=255)
                            ataagtcgttccctcaggaagGAAATGCGAGTGATTTCACCGTCACCgg                       >  1:1050302‑M2/22‑49 (MQ=255)
                               agtcgttccctcaggaagGAAATGCGAGTGATTTCACCGTCACCGGcaca                   <  1:3819691‑M2/32‑1 (MQ=255)
                               agtcgttccctcaggaagGAAATGCGAGTGATTTCACCGTCACCGGcaca                   <  1:3848056‑M2/32‑1 (MQ=255)
                               agtcgttccctcaggaagGAAATGCGAGTGATTTCACCGTCACCGGcac                    >  1:4534120‑M2/19‑49 (MQ=255)
                                gtcgttccctcaggaagGAAATGCGAGTGATTTCACCGTCACCGGCACaa                  >  1:4868972‑M2/18‑50 (MQ=255)
                                gtcgttccctcaggaagGAAATGCGAGTGATTTCACCGTCACCGGCACaa                  <  1:2863707‑M2/33‑1 (MQ=255)
                                gtcgttccctcaggaagGAAATGCGAGTGATTTCACCGTCACCGGCACaa                  >  1:4375577‑M2/18‑50 (MQ=255)
                                 tcgttccctcaggaagGAAATGCGAGTGATTTCACCGTCACCGGCACAAc                 >  1:5035680‑M2/17‑50 (MQ=255)
                                   gttccctcaggaagGAAATGCGAGTGATTTCACCGTCACCGGCACAAcc                >  1:2526042‑M2/15‑49 (MQ=255)
                                   gttccctcaggaagGAAATGCGAGTGATTTCACCGTCACCGGCACAAcc                <  1:2700191‑M2/35‑1 (MQ=255)
                                   gttccctcaggaagGAAATGCGAGTGATTTCACCGTCACCGGCACAAcc                <  1:2106426‑M2/35‑1 (MQ=255)
                                    ttccctcaggaagGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGa              >  1:1025346‑M2/14‑50 (MQ=255)
                                      ccctcaggaagGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATc            >  1:3017490‑M2/12‑50 (MQ=255)
                                       cctcaggaagGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATcc           <  1:1289520‑M2/40‑1 (MQ=255)
                                        ctcaggaagGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCg          >  1:2343675‑M2/10‑50 (MQ=255)
                                         tcaggaagGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGc         >  1:3391896‑M2/9‑50 (MQ=255)
                                         tcaggaagGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGc         >  1:4857057‑M2/9‑50 (MQ=255)
                                         tcaggaagGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGc         >  1:314998‑M2/9‑50 (MQ=255)
                                          caggaagGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGc         >  1:4283915‑M2/8‑49 (MQ=255)
                                           aggaagGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGcc        >  1:5643257‑M2/7‑49 (MQ=255)
                                            ggaagGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCa       >  1:5765436‑M2/6‑49 (MQ=255)
                                              aagGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCaaa     >  1:2139227‑M2/4‑49 (MQ=255)
                                               agGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCaaaaa   <  1:5566202‑M2/48‑1 (MQ=255)
                                               agGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCaaa     >  1:2649683‑M2/3‑48 (MQ=255)
                                                gGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCaaa     >  1:5954864‑M2/2‑47 (MQ=255)
                                                gGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAg  <  1:1805574‑M2/49‑1 (MQ=255)
                                                 gaaATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCaaaaa   >  1:5006679/1‑48 (MQ=255)
                                                 |                                                
AGCACCCAGCCACCCAGCGCCACCGCTGGCGAATACCGCATTCAGGAAGGAAATGCGAGTGATTTCACCGTCACCGGCACAACCGATCCGCCAAAAAG  >  NC_000913/2565440‑2565537

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