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 58 [0] [0] 59 [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)58 (1.600) 41/94 0.0 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

GAGCGGCACGTAGTGTGGATGCCTTACACGCCGCATTTAATCAATAACCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTATTC  >  NC_000913/2791204‑2791299
                                               |                                                
gaaaagctctattgtttacgttgggctaaacgcagggagactccccacc                                                 >  1:1545440‑M2/48‑49 (MQ=255)
 aaaagctctattgtttacgttgggctaaacgcagggagactccccaCCt                                                <  1:1123471‑M2/3‑1 (MQ=255)
  aaagctctattgtttacgttgggctaaacgcagggagactccccaCCtt                                               >  1:489601‑M2/46‑49 (MQ=255)
   aagctctattgtttacgttgggctaaacgcagggagactccccaCCttt                                              <  1:1428397‑M2/5‑1 (MQ=255)
    agctctattgtttacgttgggctaaacgcagggagactccccaCCTTTg                                             <  1:2837932‑M2/6‑1 (MQ=255)
    agctctattgtttacgttgggctaaacgcagggagactccccaCCTTTg                                             <  1:2728860‑M2/6‑1 (MQ=255)
      ctctattgtttacgttgggctaaacgcagggagactccccaCCTTTGaa                                           <  1:96056‑M2/8‑1 (MQ=255)
       tctattgtttacgttgggctaaacgcagggagactccccaCCTTTGaa                                           >  1:2428496‑M2/41‑48 (MQ=255)
       tctattgtttacgttgggctaaacgcagggagactccccaCCTTTGaa                                           >  1:1574137‑M2/41‑48 (MQ=255)
       tctattgtttacgttgggctaaacgcagggagactccccaCCTTTGaa                                           >  1:1166318‑M2/41‑48 (MQ=255)
       tctattgtttacgttgggctaaacgcagggagactccccaCCTTTGaa                                           >  1:315216‑M2/41‑48 (MQ=255)
       tctattgtttacgttgggctaaacgcagggagactccccaCCTTTGaa                                           >  1:679087‑M2/41‑48 (MQ=255)
         tattgtttacgttgggctaaacgcagggagactccccaCCTTTGaaa                                          <  1:319417‑M2/9‑1 (MQ=255)
           ttgtttacgttgggctaaacgcagggagactccccaCCTTTGaaaa                                         <  1:892805‑M2/10‑1 (MQ=255)
           ttgtttacgttgggctaaacgcagggagactccccaCCTTTGAAAACAg                                      <  1:398348‑M2/13‑1 (MQ=255)
             gtttacgttgggctaaacgcagggagactccccaCCTTTGAAAACAgg                                     <  1:1747073‑M2/14‑1 (MQ=255)
             gtttacgttgggctaaacgcagggagactccccaCCTTTGAAAACAgg                                     <  1:2440648‑M2/14‑1 (MQ=255)
             gtttacgttgggctaaacgcagggagactccccaCCTTTGAAAACAgg                                     <  1:918422‑M2/14‑1 (MQ=255)
             gtttacgttgggctaaacgcagggagactccccaCCTTTGAAAACAgg                                     <  1:1633945‑M2/14‑1 (MQ=255)
             gtttacgttgggctaaacgcagggagactccccaCCTTTGAAAACAgg                                     <  1:1264129‑M2/14‑1 (MQ=255)
             gtttacgttgggctaaacgcagggagactccccaCCTTTGAAAACAGGa                                    >  1:2714389‑M2/35‑49 (MQ=255)
              tttacgttgggctaaacgcagggagactccccaCCTTTGAAAACAGGAt                                   >  1:901551‑M2/34‑49 (MQ=255)
              tttacgttgggctaaacgcagggagactccccaCCTTTGAAAACAGGAt                                   <  1:528186‑M2/16‑1 (MQ=255)
              tttacgttgggctaaacgcagggagactccccaCCTTTGAAAACAGGAt                                   <  1:1906189‑M2/16‑1 (MQ=255)
                tacgttgggctaaacgcagggagactccccaCCTTTGAAAACAGGATg                                  >  1:2611210‑M2/32‑48 (MQ=255)
                  cgttgggctaaacgcagggagactccccaCCTTTGAAAACAGGATGTAg                               <  1:2032052‑M2/20‑1 (MQ=255)
                   gttgggctaaacgcagggagactccccaCCTTTGAAAACAGGATGTAGc                              <  1:2073546‑M2/21‑1 (MQ=255)
                    ttgggctaaacgcagggagactccccaCCTTTGAAAACAGGATGTAGCg                             <  1:2568447‑M2/22‑1 (MQ=255)
                         ctaaacgcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGaa                        <  1:1520463‑M2/27‑1 (MQ=255)
                          taaacgcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGaaa                       >  1:965704‑M2/22‑49 (MQ=255)
                            aacgcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGaaa                       <  1:2789126‑M2/28‑1 (MQ=255)
                            aacgcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACt                     <  1:1681240‑M2/30‑1 (MQ=255)
                             acgcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACtt                    <  1:1638807‑M2/31‑1 (MQ=255)
                             acgcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACtt                    <  1:825162‑M2/31‑1 (MQ=255)
                             acgcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACtt                    <  1:947425‑M2/31‑1 (MQ=255)
                               gcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTaa                  <  1:2239337‑M2/33‑1 (MQ=255)
                               gcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTaa                  <  1:2408175‑M2/33‑1 (MQ=255)
                               gcagggagactccccaCCTTTGAAAACAGGATGGAGAGATGAAAATTa                   >  1:930540‑M2/17‑48 (MQ=255)
                                cagggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAAc                 >  1:89097‑M2/16‑49 (MQ=255)
                                cagggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAAc                 <  1:1123825‑M2/34‑1 (MQ=255)
                                 agggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAAc                 >  1:2670003‑M2/15‑48 (MQ=255)
                                 agggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACg                <  1:2734477‑M2/35‑1 (MQ=255)
                                    gagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACa             <  1:1846650‑M2/38‑1 (MQ=255)
                                     agactccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAg            <  1:1717512‑M2/39‑1 (MQ=255)
                                     agactccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAg            <  1:2154479‑M2/39‑1 (MQ=255)
                                      gactccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGt           <  1:2351603‑M2/40‑1 (MQ=255)
                                       actccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGt           >  1:481101‑M2/9‑48 (MQ=255)
                                       actccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTa          <  1:455554‑M2/41‑1 (MQ=255)
                                        ctccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTaa         <  1:304936‑M2/42‑1 (MQ=255)
                                         cccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAAc        <  1:2460360‑M2/43‑1 (MQ=255)
                                           cccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACtt      <  1:2579694‑M2/45‑1 (MQ=255)
                                           cccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACt       <  1:702873‑M2/44‑1 (MQ=255)
                                           cccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACt       <  1:859774‑M2/44‑1 (MQ=255)
                                            ccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTa     <  1:474912‑M2/46‑1 (MQ=255)
                                            ccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTa     >  1:666566‑M2/4‑49 (MQ=255)
                                            ccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTa     <  1:467951‑M2/46‑1 (MQ=255)
                                             caCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTAt    >  1:580326‑M2/3‑49 (MQ=255)
                                             caCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTAt    <  1:1155346‑M2/47‑1 (MQ=255)
                                               ccTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTATTc  <  1:2187604/49‑1 (MQ=255)
                                               |                                                
GAGCGGCACGTAGTGTGGATGCCTTACACGCCGCATTTAATCAATAACCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTATTC  >  NC_000913/2791204‑2791299

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