Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 2,756,160 Δ35,091 bp [intA][lhgO] 33 genes

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 2756160 2791250 35091 72 [0] [1] 72 [intA]–[lhgO] [intA],yfjH,alpA,yfjI,yfjJ,abpB,abpA,yfjM,rnlA,rnlB,yfjP,yfjQ,yfjR,ypjK,yfjS,yfjT,yfjU,yfjV,yfjW,ypjI,yfjX,yfjY,ypjJ,yfjZ,ypjF,ypjA,pinH,ypjB,ypjC,ileY,ygaQ,csiD,[lhgO]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 27561590 (0.000)69 (0.980) 34/66 0.1 99.3% coding (1/1242 nt) intA CP4‑57 prophage; integrase
?NC_000913 2791251 = 1 (0.010)intergenic (+1/‑22) lhgO/gabD L‑2‑hydroxyglutarate oxidase/NADP(+)‑dependent succinate‑semialdehyde dehydrogenase

TGTTTACGTTGGGCCTAAACGCAGGGAGACTCCCCATGGCAAGAAAAACCAAGCCGTTAACTGATACGGA  >  NC_000913/2756124‑2756193
                                   |                                  
tGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCA                                    <  1:4517438‑M1/35‑1 (MQ=255)
tGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCA                                    <  1:5038894‑M1/35‑1 (MQ=255)
 gTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAc                                   <  1:5524957‑M1/35‑2 (MQ=255)
 gTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAc                                   <  1:5188291‑M1/35‑2 (MQ=255)
 gTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAc                                   <  1:7889841‑M1/35‑2 (MQ=255)
  tttACGTTGGGCTAAACGCAGTGAGACTCCCCAcc                                  >  1:2450369‑M1/1‑33 (MQ=255)
  tttACGTTGGGCTAAACGCAGGGAGACTCCCCAcc                                  >  1:6803093‑M1/1‑33 (MQ=255)
  tttACGTTGGGCTAAACGCAGGGAGACTCCCCAcc                                  <  1:3506016‑M1/35‑3 (MQ=255)
  tttACGTTGGGCTAAACGCAGGGAGACTCCCCAcc                                  <  1:3587960‑M1/35‑3 (MQ=255)
  tttACGTTGGGCTAAACGCAGGGAGACTCCCCAcc                                  <  1:5900939‑M1/35‑3 (MQ=255)
   ttACGTTGGGCTAAACGCAGGGAGACTCCCCAcct                                 >  1:1519758‑M1/1‑32 (MQ=255)
    tACGTTGGGCTAAACGCAGGGAGACTCCCCAcctt                                >  1:7361378‑M1/1‑31 (MQ=255)
    tACGTTGGGCTAAACGCAGGGAGACTCCCCAcctt                                <  1:2095215‑M1/35‑5 (MQ=255)
    tACGTTGGGCTAAACGCAGGGAGACTCCCCAcctt                                >  1:5312253‑M1/1‑31 (MQ=255)
     aCGTTGGGCTAAACGCAGGGAGACTCCCCAccttt                               <  1:4042105‑M1/35‑6 (MQ=255)
     aCGTTGGGCTAAACGCAGGGAGACTCCCCAccttt                               <  1:7061829‑M1/35‑6 (MQ=255)
     aCGTTGGGCTAAACGCAGGGAGACTCCCCAccttt                               >  1:5313097‑M1/1‑30 (MQ=255)
     aCGTTGGGCTAAACGCAGGGAGACTCCCCAccttt                               <  1:3465330‑M1/35‑6 (MQ=255)
      cGTTGGGCTAAACGCAGGGAGACTCCCCAcctttg                              <  1:1814632‑M1/35‑7 (MQ=255)
      cGTTGGGCTAAACGCAGGGAGACTCCCCAcctttg                              <  1:1952363‑M1/35‑7 (MQ=255)
      cGTTGGGCTAAACGCAGGGAGACTCCCCAcctttg                              <  1:5641922‑M1/35‑7 (MQ=255)
      cGTTGGGCTAAACGCAGGGAGACTCCCCAcctttg                              <  1:2397150‑M1/35‑7 (MQ=255)
      cGTTGGGCTAAACGCAGGGAGACTCCCCAcctttg                              <  1:6471714‑M1/35‑7 (MQ=255)
      cGTTGGGCTAAACGCAGGGAGACTCCCCAcctttg                              >  1:3314158‑M1/1‑29 (MQ=255)
              cTAAACGCAGGGAGACTCCCCAcctttgaaaacag                       >  1:6815828‑M1/1‑22 (MQ=255)
               tAAACGCAGGTAGACTCCCCAcctttgaaaacagg                      >  1:3686511‑M1/1‑21 (MQ=255)
               tAAACGCAGGGAGACTCCCCAcctttgaaaacagg                      >  1:3658313‑M1/1‑21 (MQ=255)
               tAAACGCAGGGAGACTCCCCAcctttgaaaacagg                      >  1:4248668‑M1/1‑21 (MQ=255)
                aaaCGCAGGGAGACTCCCCAcctttgaaaacagga                     >  1:3862996‑M1/1‑20 (MQ=255)
                aaaCGCAGGGAGACTCCCCAcctttgaaaacagga                     <  1:2119998‑M1/35‑16 (MQ=255)
                 aaCGCAGGGAGACTCCCCAcctttgaaaacaggat                    >  1:3725313‑M1/1‑19 (MQ=255)
                 aaCGCAGGGAGACTCCCCAcctttgaaaacaggat                    >  1:553712‑M1/1‑19 (MQ=255)
                 aaCGCAGGGAGACTCCCCAcctttgaaaacaggat                    >  1:1494537‑M1/1‑19 (MQ=255)
                 aaCGCAGGGAGACTCCCCAcctttgaaaacaggat                    >  1:1677454‑M1/1‑19 (MQ=255)
                 aaCGCAGGGAGACTCCCCAcctttgaaaacaggat                    >  1:5021940‑M1/1‑19 (MQ=255)
                  aCGCAGGGAGACTCCCCAcctttgaaaacaggatg                   <  1:1254242‑M1/35‑18 (MQ=255)
                  aCGCAGGGAGACTCCCCAcctttgaaaacaggatg                   <  1:6810276‑M1/35‑18 (MQ=255)
                  aCGCAGGGAGACTCCCCAcctttgaaaacaggatg                   <  1:6579162‑M1/35‑18 (MQ=255)
                  aCGCAGGGAGACTCCCCAcctttgaaaacaggatg                   <  1:4112864‑M1/35‑18 (MQ=255)
                  aCGCAGGGAGACTCCCCAcctttgaaaacaggatg                   <  1:1955722‑M1/35‑18 (MQ=255)
                   cGCAGGGAGACTCCCCAcctttgaaaacaggatgt                  <  1:1570188‑M1/35‑19 (MQ=255)
                   cGCAGGGAGACTCCCCAcctttgaaaacaggatgt                  <  1:5801186‑M1/35‑19 (MQ=255)
                    gCAGGGAGACTCCCCAcctttgaaaacaggatgta                 <  1:1524199‑M1/35‑20 (MQ=255)
                     cAGGGAGACTCCCCAcctttgaaaacaggatgtag                <  1:7201638‑M1/35‑21 (MQ=255)
                     cAGGGAGACTCCCCAcctttgaaaacaggatgtag                <  1:2342731‑M1/35‑21 (MQ=255)
                      aGGGAGACTCCCCAcctttgaaaacaggatgtagc               >  1:4342151‑M1/1‑14 (MQ=255)
                      aGGGAGACTCCCCAcctttgaaaacaggatgtagc               <  1:1531246‑M1/35‑22 (MQ=255)
                       gggAGACTCCCCAcctttgaaaacaggatgtagcg              >  1:6599516‑M1/1‑13 (MQ=255)
                        ggAGACTCCCCAcctttgaaaacaggatgtagcga             <  1:4839565‑M1/35‑24 (MQ=255)
                         gagaCTCCCCAcctttgaaaacaggatgtagcgat            <  1:6959982‑M1/35‑25 (MQ=255)
                         gagaCTCCCCAcctttgaaaacaggatgtagcgat            >  1:6539738‑M1/1‑11 (MQ=255)
                         gagaCTCCCCAcctttgaaaacaggatgtagcgat            >  1:2504858‑M1/1‑11 (MQ=255)
                         gagaCTCCCCAcctttgaaaacaggatgtagcgat            >  1:5189504‑M1/1‑11 (MQ=255)
                         gagaCTCCCCAcctttgaaaacaggatgtagcgat            >  1:101633‑M1/1‑11 (MQ=255)
                          agaCTCCCCAcctttgaaaacaggatgtagcgatg           >  1:1638073‑M1/1‑10 (MQ=255)
                          agaCTCCCCAcctttgaaaacaggatgtagcgatg           >  1:3055868‑M1/1‑10 (MQ=255)
                          agaCTCCCCAcctttgaaaacaggatgtagcgatg           >  1:4854284‑M1/1‑10 (MQ=255)
                           gaCTCCCCAcctttgaaaacaggatgtagcgatga          >  1:5545890‑M1/1‑9 (MQ=255)
                           gaCTCCCCAcctttgaaaacaggatgtagcgatga          <  1:2460155‑M1/35‑27 (MQ=255)
                           gaCTCCCCAcctttgaaaacaggatgtagcgatga          <  1:7151548‑M1/35‑27 (MQ=255)
                           gaCTCCCCAcctttgaaaacaggatgtagcgatga          <  1:760649‑M1/35‑27 (MQ=255)
                           gaCTCCCCAcctttgaaaacaggatgtagcgatga          <  1:1433863‑M1/35‑27 (MQ=255)
                            aCTCCCCAcctttgaaaacaggatgtagcgatgaa         <  1:7270399‑M1/35‑28 (MQ=255)
                             cTCCCCAcctttgaaaacaggatgtagcgatgaaa        >  1:4492330‑M1/1‑7 (MQ=255)
                             cTCCCCAcctttgaaaacaggatgtagcgatgaaa        >  1:5003033‑M1/1‑7 (MQ=255)
                              tCCCCAcctttgaaaacaggatgtagcgatgaaac       <  1:6450007‑M1/35‑30 (MQ=255)
                               ccccAcctttgaaaacaggatgtagcgatgaaact      <  1:632328‑M1/35‑31 (MQ=255)
                               ccccAcctttgaaaacaggatgtagcgatgaaact      <  1:2866672‑M1/35‑31 (MQ=255)
                                cccAcctttgaaaacaggatgtagcgatgaaactt     <  1:3911672‑M1/35‑32 (MQ=255)
                                  cAcctttgaaaacaggatgtagcgatgaaacttaa   <  1:4467106‑M1/35‑34 (MQ=255)
                                  cAcctttgaaaacaggatgtagcgatgaaacttaa   <  1:839791‑M1/35‑34 (MQ=255)
                                   acctttgaaaacaggatgtagcgatgaaacttaac  <  1:6368444‑M1/35‑35 (MQ=255)
                                   |                                  
TGTTTACGTTGGGCCTAAACGCAGGGAGACTCCCCATGGCAAGAAAAACCAAGCCGTTAACTGATACGGA  >  NC_000913/2756124‑2756193

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