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 90 [0] [0] 84 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)77 (0.850) 51/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

CATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGAATCGATGATTGGCATTTTTCACAGCATTCTAACAATCACGTTTCATCG  >  NC_000913/2558649‑2558747
                                                 |                                                 
cATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAG                                                   <  1:4289535‑M1/50‑1 (MQ=255)
cATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAG                                                   >  1:5730042‑M1/1‑50 (MQ=255)
cATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAG                                                   <  1:4905965‑M1/50‑1 (MQ=255)
 aTTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGg                                                  <  1:2241856‑M1/50‑2 (MQ=255)
 aTTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGg                                                  <  1:5492047‑M1/50‑2 (MQ=255)
 aTTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAG                                                   >  1:3013436‑M1/1‑49 (MQ=255)
 aTTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAG                                                   >  1:8127742‑M1/1‑49 (MQ=255)
 aTTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAG                                                   >  1:5097333‑M1/1‑49 (MQ=255)
  ttGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGg                                                  >  1:4144296‑M1/1‑48 (MQ=255)
     ccGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaat                                              >  1:5339907‑M1/1‑45 (MQ=255)
     ccGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaa                                               <  1:840167‑M1/49‑5 (MQ=255)
     ccGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaa                                               <  1:1180276‑M1/49‑5 (MQ=255)
      cGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatg                                             >  1:6162721‑M1/1‑44 (MQ=255)
      cGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatg                                             <  1:3988557‑M1/50‑7 (MQ=255)
      cGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatg                                             <  1:1171609‑M1/50‑7 (MQ=255)
       gAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgc                                            >  1:977655‑M1/1‑43 (MQ=255)
       gAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgc                                            >  1:3095663‑M1/1‑43 (MQ=255)
        aaCAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgc                                            >  1:4973013‑M1/1‑42 (MQ=255)
         aCAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcga                                          <  1:79006‑M1/50‑10 (MQ=255)
         aCAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgc                                            >  1:4127778‑M1/1‑41 (MQ=255)
          cAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgag                                         <  1:2755701‑M1/50‑11 (MQ=255)
          cAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgag                                         >  1:7924052‑M1/1‑40 (MQ=255)
           aaTGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagt                                        >  1:8043770‑M1/1‑39 (MQ=255)
           aaTGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagt                                        <  1:3906348‑M1/50‑12 (MQ=255)
             tgtgGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtga                                      >  1:6668536‑M1/1‑37 (MQ=255)
              gtgGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgat                                     >  1:5007202‑M1/1‑36 (MQ=255)
               tgGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgat                                     <  1:1660347‑M1/49‑15 (MQ=255)
                ggTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgattt                                   >  1:7424777‑M1/1‑34 (MQ=255)
                 gTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttc                                  <  1:4714310‑M1/50‑18 (MQ=255)
                 gTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttc                                  <  1:3506336‑M1/50‑18 (MQ=255)
                 gTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgattt                                   <  1:8251394‑M1/49‑17 (MQ=255)
                  tCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttca                                 >  1:287114‑M1/1‑32 (MQ=255)
                  tCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttca                                 <  1:5484736‑M1/50‑19 (MQ=255)
                  tCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttca                                 >  1:4655156‑M1/1‑32 (MQ=255)
                  tCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgaatgatttca                                 >  1:1571917‑M1/1‑32 (MQ=255)
                   cTTTAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcac                                >  1:7891097‑M1/1‑31 (MQ=255)
                    tttAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcacc                               >  1:390893‑M1/1‑30 (MQ=255)
                    tttAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcacc                               >  1:3789710‑M1/1‑30 (MQ=255)
                      tAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgt                             >  1:130749‑M1/1‑28 (MQ=255)
                      tAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccg                              <  1:3117554‑M1/49‑22 (MQ=255)
                      tAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcacc                               <  1:7365430‑M1/48‑21 (MQ=255)
                      tAGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcacc                               <  1:4393011‑M1/48‑21 (MQ=255)
                       aGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgt                             <  1:4609347‑M1/49‑23 (MQ=255)
                       aGTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgt                             <  1:2212678‑M1/49‑23 (MQ=255)
                        gTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtc                            <  1:7437792‑M1/49‑24 (MQ=255)
                        gTTTCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtc                            <  1:6864009‑M1/49‑24 (MQ=255)
                         tttCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcac                          >  1:5875467‑M1/1‑25 (MQ=255)
                          ttCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcacc                         >  1:4261479‑M1/1‑24 (MQ=255)
                          ttCATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcacc                         >  1:7833198‑M1/1‑24 (MQ=255)
                            cATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccgg                       >  1:1512941‑M1/1‑22 (MQ=255)
                            cATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccgg                       >  1:6949760‑M1/1‑22 (MQ=255)
                            cATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccgg                       >  1:6582022‑M1/1‑22 (MQ=255)
                            cATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccgg                       >  1:32268‑M1/1‑22 (MQ=255)
                            cATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccg                        <  1:7059394‑M1/49‑28 (MQ=255)
                            cATAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccg                        <  1:1629066‑M1/49‑28 (MQ=255)
                             aTAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggc                      <  1:4156849‑M1/50‑30 (MQ=255)
                             aTAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggc                      <  1:4763923‑M1/50‑30 (MQ=255)
                             aTAAGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccgg                       >  1:6315331‑M1/1‑21 (MQ=255)
                               aaGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggca                     >  1:5724307‑M1/1‑19 (MQ=255)
                               aaGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggca                     >  1:4945950‑M1/1‑19 (MQ=255)
                                aGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcaca                   <  1:844077‑M1/50‑33 (MQ=255)
                                aGTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcac                    >  1:6468377‑M1/1‑18 (MQ=255)
                                 gTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaa                  <  1:2707048‑M1/50‑34 (MQ=255)
                                 gTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaa                  >  1:4730399‑M1/1‑17 (MQ=255)
                                 gTCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaa                  >  1:5613396‑M1/1‑17 (MQ=255)
                                  tCGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaac                 >  1:5312028‑M1/1‑16 (MQ=255)
                                   cGTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaacc                >  1:3400684‑M1/1‑15 (MQ=255)
                                    gTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaacc                >  1:7176304‑M1/1‑14 (MQ=255)
                                    gTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaacc                <  1:2845379‑M1/49‑36 (MQ=255)
                                    gTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaacc                >  1:17719‑M1/1‑14 (MQ=255)
                                    gTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaacc                <  1:7556227‑M1/49‑36 (MQ=255)
                                    gTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaacc                <  1:1475030‑M1/49‑36 (MQ=255)
                                    gTTCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaacc                >  1:5995026‑M1/1‑14 (MQ=255)
                                     ttCCCTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcgcaaccga              <  1:1288473‑M1/50‑38 (MQ=255)
                                       cccTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgat             <  1:7118872‑M1/49‑39 (MQ=255)
                                        ccTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatcc           <  1:7057021‑M1/50‑41 (MQ=255)
                                        ccTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatcc           >  1:7459965‑M1/1‑10 (MQ=255)
                                        ccTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatcc           <  1:1567132‑M1/50‑41 (MQ=255)
                                        ccTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatcc           >  1:2455903‑M1/1‑10 (MQ=255)
                                        ccTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatc            <  1:3892603‑M1/49‑40 (MQ=255)
                                        ccTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatc            <  1:4849109‑M1/49‑40 (MQ=255)
                                         cTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatccg          <  1:1527877‑M1/50‑42 (MQ=255)
                                         cTCAGGAAGgaaatgcgagtgatttcaccgtcaccggcacaaccgatccg          >  1:3405875‑M1/1‑9 (MQ=255)
                                            aggaaggaaatgcgagtgatttcaccgtcaccggcacaaccgatccgcc        >  1:3182906‑M1/1‑6 (MQ=255)
                                            aggaaggaaatgcgagtgatttcaccgtcaccggcacaaccgatccgc         <  1:4042549‑M1/48‑43 (MQ=255)
                                              gaaggaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaa      >  1:1189779‑M1/1‑4 (MQ=255)
                                               aaggaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaa     >  1:2893070‑M1/1‑3 (MQ=255)
                                               aaggaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaa     >  1:2976244‑M1/1‑3 (MQ=255)
                                                aggaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaa     >  1:6916583‑M1/1‑2 (MQ=255)
                                                 ggaaatgcgagtgatttcaccgtcaccggcacaaccgatccgccaaaaag  <  1:2095915‑M1/50‑50 (MQ=255)
                                                 |                                                 
CATTGCCGAACAATGTGGTCTTTAGTTTCATAAGTCGTTCCCTCAGGAAGAATCGATGATTGGCATTTTTCACAGCATTCTAACAATCACGTTTCATCG  >  NC_000913/2558649‑2558747

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: