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

TGTGGATGCCTTACACGCCGCATTTAATCAATAACCTTTGAAAACAGGATGTAGCGATGAAACTTAACG  >  NC_000913/2791217‑2791285
                                  |                                  
tgtgGATGCCTTACACGCCGCATTTAATCAATAAc                                    <  1:5805455/35‑1 (MQ=255)
gtttacgttgggctaaacgcagggagactccccac                                    <  1:7889841‑M2/1‑1 (MQ=255)
gtttacgttgggctaaacgcagggagactccccac                                    <  1:5188291‑M2/1‑1 (MQ=255)
gtttacgttgggctaaacgcagggagactccccac                                    <  1:5524957‑M2/1‑1 (MQ=255)
 tttacgttgggctaaacgcagtgagactccccacc                                   >  1:2450369‑M2/34‑35 (MQ=255)
 tttacgttgggctaaacgcagggagactccccacc                                   <  1:3587960‑M2/2‑1 (MQ=255)
 tttacgttgggctaaacgcagggagactccccacc                                   <  1:3506016‑M2/2‑1 (MQ=255)
 tttacgttgggctaaacgcagggagactccccacc                                   >  1:6803093‑M2/34‑35 (MQ=255)
 tttacgttgggctaaacgcagggagactccccacc                                   <  1:5900939‑M2/2‑1 (MQ=255)
  ttacgttgggctaaacgcagggagactccccaCCt                                  >  1:1519758‑M2/33‑35 (MQ=255)
   tacgttgggctaaacgcagggagactccccaCCtt                                 >  1:7361378‑M2/32‑35 (MQ=255)
   tacgttgggctaaacgcagggagactccccaCCtt                                 >  1:5312253‑M2/32‑35 (MQ=255)
   tacgttgggctaaacgcagggagactccccaCCtt                                 <  1:2095215‑M2/4‑1 (MQ=255)
    acgttgggctaaacgcagggagactccccaCCttt                                >  1:5313097‑M2/31‑35 (MQ=255)
    acgttgggctaaacgcagggagactccccaCCttt                                <  1:4042105‑M2/5‑1 (MQ=255)
    acgttgggctaaacgcagggagactccccaCCttt                                <  1:3465330‑M2/5‑1 (MQ=255)
    acgttgggctaaacgcagggagactccccaCCttt                                <  1:7061829‑M2/5‑1 (MQ=255)
     cgttgggctaaacgcagggagactccccaCCTTTg                               <  1:1952363‑M2/6‑1 (MQ=255)
     cgttgggctaaacgcagggagactccccaCCTTTg                               <  1:1814632‑M2/6‑1 (MQ=255)
     cgttgggctaaacgcagggagactccccaCCTTTg                               <  1:5641922‑M2/6‑1 (MQ=255)
     cgttgggctaaacgcagggagactccccaCCTTTg                               <  1:2397150‑M2/6‑1 (MQ=255)
     cgttgggctaaacgcagggagactccccaCCTTTg                               <  1:6471714‑M2/6‑1 (MQ=255)
     cgttgggctaaacgcagggagactccccaCCTTTg                               >  1:3314158‑M2/30‑35 (MQ=255)
            ctaaacgcagggagactccccaCCTTTGAAAACAg                        >  1:6815828‑M2/23‑35 (MQ=255)
             taaacgcaggtagactccccaCCTTTGAAAACAgg                       >  1:3686511‑M2/22‑35 (MQ=255)
             taaacgcagggagactccccaCCTTTGAAAACAgg                       >  1:3658313‑M2/22‑35 (MQ=255)
             taaacgcagggagactccccaCCTTTGAAAACAgg                       >  1:4248668‑M2/22‑35 (MQ=255)
              aaacgcagggagactccccaCCTTTGAAAACAGGa                      >  1:3862996‑M2/21‑35 (MQ=255)
              aaacgcagggagactccccaCCTTTGAAAACAGGa                      <  1:2119998‑M2/15‑1 (MQ=255)
               aacgcagggagactccccaCCTTTGAAAACAGGAt                     >  1:3725313‑M2/20‑35 (MQ=255)
               aacgcagggagactccccaCCTTTGAAAACAGGAt                     >  1:1677454‑M2/20‑35 (MQ=255)
               aacgcagggagactccccaCCTTTGAAAACAGGAt                     >  1:553712‑M2/20‑35 (MQ=255)
               aacgcagggagactccccaCCTTTGAAAACAGGAt                     >  1:5021940‑M2/20‑35 (MQ=255)
               aacgcagggagactccccaCCTTTGAAAACAGGAt                     >  1:1494537‑M2/20‑35 (MQ=255)
                acgcagggagactccccaCCTTTGAAAACAGGATg                    <  1:1254242‑M2/17‑1 (MQ=255)
                acgcagggagactccccaCCTTTGAAAACAGGATg                    <  1:6579162‑M2/17‑1 (MQ=255)
                acgcagggagactccccaCCTTTGAAAACAGGATg                    <  1:6810276‑M2/17‑1 (MQ=255)
                acgcagggagactccccaCCTTTGAAAACAGGATg                    <  1:1955722‑M2/17‑1 (MQ=255)
                acgcagggagactccccaCCTTTGAAAACAGGATg                    <  1:4112864‑M2/17‑1 (MQ=255)
                 cgcagggagactccccaCCTTTGAAAACAGGATGt                   <  1:5801186‑M2/18‑1 (MQ=255)
                 cgcagggagactccccaCCTTTGAAAACAGGATGt                   <  1:1570188‑M2/18‑1 (MQ=255)
                  gcagggagactccccaCCTTTGAAAACAGGATGTa                  <  1:1524199‑M2/19‑1 (MQ=255)
                   cagggagactccccaCCTTTGAAAACAGGATGTAg                 <  1:2342731‑M2/20‑1 (MQ=255)
                   cagggagactccccaCCTTTGAAAACAGGATGTAg                 <  1:7201638‑M2/20‑1 (MQ=255)
                    agggagactccccaCCTTTGAAAACAGGATGTAGc                >  1:4342151‑M2/15‑35 (MQ=255)
                    agggagactccccaCCTTTGAAAACAGGATGTAGc                <  1:1531246‑M2/21‑1 (MQ=255)
                     gggagactccccaCCTTTGAAAACAGGATGTAGCg               >  1:6599516‑M2/14‑35 (MQ=255)
                      ggagactccccaCCTTTGAAAACAGGATGTAGCGa              <  1:4839565‑M2/23‑1 (MQ=255)
                       gagactccccaCCTTTGAAAACAGGATGTAGCGAt             <  1:6959982‑M2/24‑1 (MQ=255)
                       gagactccccaCCTTTGAAAACAGGATGTAGCGAt             >  1:6539738‑M2/12‑35 (MQ=255)
                       gagactccccaCCTTTGAAAACAGGATGTAGCGAt             >  1:5189504‑M2/12‑35 (MQ=255)
                       gagactccccaCCTTTGAAAACAGGATGTAGCGAt             >  1:101633‑M2/12‑35 (MQ=255)
                       gagactccccaCCTTTGAAAACAGGATGTAGCGAt             >  1:2504858‑M2/12‑35 (MQ=255)
                        agactccccaCCTTTGAAAACAGGATGTAGCGATg            >  1:4854284‑M2/11‑35 (MQ=255)
                        agactccccaCCTTTGAAAACAGGATGTAGCGATg            >  1:1638073‑M2/11‑35 (MQ=255)
                        agactccccaCCTTTGAAAACAGGATGTAGCGATg            >  1:3055868‑M2/11‑35 (MQ=255)
                         gactccccaCCTTTGAAAACAGGATGTAGCGATGa           <  1:7151548‑M2/26‑1 (MQ=255)
                         gactccccaCCTTTGAAAACAGGATGTAGCGATGa           <  1:1433863‑M2/26‑1 (MQ=255)
                         gactccccaCCTTTGAAAACAGGATGTAGCGATGa           >  1:5545890‑M2/10‑35 (MQ=255)
                         gactccccaCCTTTGAAAACAGGATGTAGCGATGa           <  1:760649‑M2/26‑1 (MQ=255)
                         gactccccaCCTTTGAAAACAGGATGTAGCGATGa           <  1:2460155‑M2/26‑1 (MQ=255)
                          actccccaCCTTTGAAAACAGGATGTAGCGATGaa          <  1:7270399‑M2/27‑1 (MQ=255)
                           ctccccaCCTTTGAAAACAGGATGTAGCGATGaaa         >  1:5003033‑M2/8‑35 (MQ=255)
                           ctccccaCCTTTGAAAACAGGATGTAGCGATGaaa         >  1:4492330‑M2/8‑35 (MQ=255)
                            tccccaCCTTTGAAAACAGGATGTAGCGATGAAAc        <  1:6450007‑M2/29‑1 (MQ=255)
                             ccccaCCTTTGAAAACAGGATGTAGCGATGAAACt       <  1:632328‑M2/30‑1 (MQ=255)
                             ccccaCCTTTGAAAACAGGATGTAGCGATGAAACt       <  1:2866672‑M2/30‑1 (MQ=255)
                              cccaCCTTTGAAAACAGGATGTAGCGATGAAACtt      <  1:3911672‑M2/31‑1 (MQ=255)
                                caCCTTTGAAAACAGGATGTAGCGATGAAACTTaa    <  1:839791‑M2/33‑1 (MQ=255)
                                caCCTTTGAAAACAGGATGTAGCGATGAAACTTaa    <  1:4467106‑M2/33‑1 (MQ=255)
                                 aCCTTTGAAAACAGGATGTAGCGATGAAACTTAAc   <  1:6368444‑M2/34‑1 (MQ=255)
                                  ccTTTGAAAACAGGATGTAGCGATGAAACTTAACg  >  1:1876209/1‑35 (MQ=255)
                                  |                                  
TGTGGATGCCTTACACGCCGCATTTAATCAATAACCTTTGAAAACAGGATGTAGCGATGAAACTTAACG  >  NC_000913/2791217‑2791285

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.