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 26 [0] [0] 26 [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)26 (0.370) 20/96 1.4 100% coding (1/1242 nt) intA CP4‑57 prophage; integrase
?NC_000913 2791251 = 0 (0.000)intergenic (+1/‑22) lhgO/gabD L‑2‑hydroxyglutarate oxidase/NADP(+)‑dependent succinate‑semialdehyde dehydrogenase

GGCACGTAGTGTGGATGCCTTACACGCCGCATTTAATCAATAACCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTATT  >  NC_000913/2791208‑2791298
                                           |                                               
agctctattgtttacgttgggctaaacgcagggagactccccaCCTTTGa                                           <  1:855300‑M2/7‑1 (MQ=255)
 gctctattgtttacgttgggctaaacgcagggagactccccaCCTTTGaa                                          <  1:1183860‑M2/8‑1 (MQ=255)
       ttgtttacgttgggctaaacgcagggagactccccaCCTTTGAAAAc                                       <  1:2587862‑M2/11‑1 (MQ=255)
       ttgtttacgttgggctaaacgcagggagactccccaCCTTTGAAAACAgg                                    <  1:4581397‑M2/14‑1 (MQ=255)
       ttgtttacgttgggctaaacgcagggagactccccaCCTTTGAAAACAg                                     <  1:2668513‑M2/13‑1 (MQ=255)
          tttacgttgggctaaacgcagggagactccccaCCTTTGAAAACAGGATg                                 <  1:849462‑M2/17‑1 (MQ=255)
          tttacgttgggctaaacgcagggagactccccaCCTTTGAAAACAGGATg                                 <  1:841190‑M2/17‑1 (MQ=255)
            tacgttgggctaaacgcagggagactccccaCCTTTGAAAACAGGATGt                                <  1:154241‑M2/18‑1 (MQ=255)
              cgttgggctaaacgcagggagactccccaCCTTTGAAAACAGGATGTAg                              <  1:3320317‑M2/20‑1 (MQ=255)
              cgttgggctaaacgcagggagactccccaCCTTTGAAAACAGGATGTAGc                             <  1:216595‑M2/21‑1 (MQ=255)
                   ggctaaacgcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGa                        <  1:1585491‑M2/26‑1 (MQ=255)
                        aacgcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAAc                     <  1:3232455‑M2/29‑1 (MQ=255)
                        aacgcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACtt                   <  1:4783764‑M2/31‑1 (MQ=255)
                         acgcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTa                  <  1:3894576‑M2/32‑1 (MQ=255)
                               ggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACa            <  1:1319308‑M2/38‑1 (MQ=255)
                               ggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACa            <  1:164997‑M2/38‑1 (MQ=255)
                                 agactccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGt          <  1:1795611‑M2/40‑1 (MQ=255)
                                  gactccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTa         <  1:1875764‑M2/41‑1 (MQ=255)
                                  gactccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTa         <  1:3631790‑M2/41‑1 (MQ=255)
                                   actccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTaa        <  1:10661‑M2/42‑1 (MQ=255)
                                   actccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTaa        <  1:5234137‑M2/42‑1 (MQ=255)
                                     tccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACt      >  1:3019043‑M2/7‑50 (MQ=255)
                                       cccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACtt     <  1:535135‑M2/45‑1 (MQ=255)
                                        ccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTAt   <  1:4820805‑M2/47‑1 (MQ=255)
                                         caCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTAtt  <  1:500931‑M2/48‑1 (MQ=255)
                                         caCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTAtt  <  1:1986567‑M2/48‑1 (MQ=255)
                                           |                                               
GGCACGTAGTGTGGATGCCTTACACGCCGCATTTAATCAATAACCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTATT  >  NC_000913/2791208‑2791298

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