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 39 [0] [0] 39 [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)37 (0.480) 23/96 0.9 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

GCGAGCGGCACGTAGTGTGGATGCCTTACACGCCGCATTTAATCAATAACCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTATTCC  >  NC_000913/2791202‑2791300
                                                 |                                                 
acgaaaagctctattgtttacgttgggctaaacgcagggagactccccac                                                   <  1:2043793‑M2/1‑1 (MQ=255)
  gaaaagctctattgtttacgttgggctaaacgcagggagactccccaCCt                                                 <  1:463649‑M2/3‑1 (MQ=255)
    aaagctctattgtttacgttgggctaaacgcagggagactccccaCCttt                                               <  1:1948158‑M2/5‑1 (MQ=255)
     aagctctattgtttacgttgggctaaacgcagggagactccccaCCTTTg                                              <  1:3144922‑M2/6‑1 (MQ=255)
      agctctattgtttacgttgggctaaacgcagggagactccccaCCTTTGa                                             <  1:1912124‑M2/7‑1 (MQ=255)
       gctctattgtttacgttgggctaaacgcagggagactccccaCCTTTGaa                                            <  1:2843086‑M2/8‑1 (MQ=255)
       gctctattgtttacgttgggctaaacgcagggagactccccaCCTTTGaa                                            <  1:2138591‑M2/8‑1 (MQ=255)
             ttgtttacgttgggctaaacgcagggagactccccaCCTTTGAAAACAg                                       <  1:4146227‑M2/13‑1 (MQ=255)
             ttgtttacgttgggctaaacgcagggagactccccaCCTTTGAAAACAg                                       <  1:1633771‑M2/13‑1 (MQ=255)
                tttacgttgggctaaacgcagggagactccccaCCTTTGAAAACAGGAt                                    <  1:5302527‑M2/16‑1 (MQ=255)
                tttacgttgggctaaacgcagggagactccccaCCTTTGAAAACAGGAt                                    <  1:5586390‑M2/16‑1 (MQ=255)
                  tacgttgggctaaacgcagggagactccccaCCTTTGAAAACAGGATGt                                  <  1:320780‑M2/18‑1 (MQ=255)
                   acgttgggctaaacgcagggagactccccaCCTTTGAAAACAGGATGTa                                 <  1:1698181‑M2/19‑1 (MQ=255)
                   acgttgggctaaacgcagggagactccccaCCTTTGAAAACAGGATGTa                                 <  1:247029‑M2/19‑1 (MQ=255)
                   acgttgggctaaacgcagggagactccccaCCTTTGAAAACAGGATGTa                                 <  1:1636020‑M2/19‑1 (MQ=255)
                   acgttgggctaaacgcagggagactccccaCCTTTGAAAACAGGATGTa                                 <  1:4194364‑M2/19‑1 (MQ=255)
                    cgttgggctaaacgcagggagactccccaCCTTTGAAAACAGGATGTAGc                               <  1:5351138‑M2/21‑1 (MQ=255)
                     gttgggctaaacgcagggagactccccaCCTTTGAAAACAGGATGTAGc                               >  1:5690871‑M2/29‑49 (MQ=255)
                         ggctaaacgcagggagactccccaCCTTTGAAAACAGGATGTAGCGATg                           <  1:1462038‑M2/25‑1 (MQ=255)
                         ggctaaacgcagggagactccccaCCTTTGAAAACAGGATGTAGCGATg                           <  1:2379174‑M2/25‑1 (MQ=255)
                         ggctaaacgcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGa                          <  1:317766‑M2/26‑1 (MQ=255)
                          gctaaacgcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGaa                         <  1:2506883‑M2/27‑1 (MQ=255)
                          gctaaacgcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGaa                         <  1:859155‑M2/27‑1 (MQ=255)
                              aacgcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAAc                       <  1:4019494‑M2/29‑1 (MQ=255)
                              aacgcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACt                      <  1:297013‑M2/30‑1 (MQ=255)
                               acgcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTa                    <  1:4667078‑M2/32‑1 (MQ=255)
                                 gcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAAc                  <  1:5514878‑M2/34‑1 (MQ=255)
                                 gcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAAc                  <  1:1595385‑M2/34‑1 (MQ=255)
                                   agggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGa                <  1:2122220‑M2/36‑1 (MQ=255)
                                      gagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAg             <  1:3724879‑M2/39‑1 (MQ=255)
                                       agactccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGt            <  1:3105705‑M2/40‑1 (MQ=255)
                                        gactccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTa           <  1:4116232‑M2/41‑1 (MQ=255)
                                          ctccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAAc         <  1:2214376‑M2/43‑1 (MQ=255)
                                          ctccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAAc         <  1:561443‑M2/43‑1 (MQ=255)
                                             cccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACtt       <  1:2769537‑M2/45‑1 (MQ=255)
                                               caCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTAtt    <  1:3768977‑M2/48‑1 (MQ=255)
                                               caCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTAtt    <  1:1598903‑M2/48‑1 (MQ=255)
                                                aCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTATTc   <  1:573633‑M2/49‑1 (MQ=255)
                                                 ccTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTATTcc  <  1:2749896/50‑1 (MQ=255)
                                                 |                                                 
GCGAGCGGCACGTAGTGTGGATGCCTTACACGCCGCATTTAATCAATAACCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTATTCC  >  NC_000913/2791202‑2791300

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: