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

GAAAAGCTCTATTGTTTACGTTGGGCCTAAACGCAGGGAGACTCCCCATGGCAAGAAAAACCAAGCCGTTAACTGATACGGAAATCAAAGCCGCC  >  NC_000913/2756112‑2756206
                                               |                                               
gAAAAGCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcc                                               >  1:1545440‑M1/1‑47 (MQ=255)
 aaaaGCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcct                                              <  1:1123471‑M1/49‑4 (MQ=255)
  aaaGCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcctt                                             >  1:489601‑M1/1‑45 (MQ=255)
   aaGCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAccttt                                            <  1:1428397‑M1/49‑6 (MQ=255)
    aGCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttg                                           <  1:2728860‑M1/49‑7 (MQ=255)
    aGCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttg                                           <  1:2837932‑M1/49‑7 (MQ=255)
      ctctATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaa                                         <  1:96056‑M1/49‑9 (MQ=255)
       tctATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaa                                         >  1:1166318‑M1/1‑40 (MQ=255)
       tctATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaa                                         >  1:1574137‑M1/1‑40 (MQ=255)
       tctATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaa                                         >  1:315216‑M1/1‑40 (MQ=255)
       tctATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaa                                         >  1:2428496‑M1/1‑40 (MQ=255)
       tctATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaa                                         >  1:679087‑M1/1‑40 (MQ=255)
         tATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaa                                        <  1:319417‑M1/47‑10 (MQ=255)
           ttGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacag                                    <  1:398348‑M1/49‑14 (MQ=255)
           ttGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaa                                       <  1:892805‑M1/46‑11 (MQ=255)
             gTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacagga                                  >  1:2714389‑M1/1‑34 (MQ=255)
             gTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacagg                                   <  1:1747073‑M1/48‑15 (MQ=255)
             gTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacagg                                   <  1:918422‑M1/48‑15 (MQ=255)
             gTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacagg                                   <  1:2440648‑M1/48‑15 (MQ=255)
             gTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacagg                                   <  1:1633945‑M1/48‑15 (MQ=255)
             gTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacagg                                   <  1:1264129‑M1/48‑15 (MQ=255)
              tttACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacaggat                                 <  1:528186‑M1/49‑17 (MQ=255)
              tttACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacaggat                                 <  1:1906189‑M1/49‑17 (MQ=255)
              tttACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacaggat                                 >  1:901551‑M1/1‑33 (MQ=255)
                tACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacaggatg                                >  1:2611210‑M1/1‑31 (MQ=255)
                  cGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacaggatgtag                             <  1:2032052‑M1/49‑21 (MQ=255)
                   gTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacaggatgtagc                            <  1:2073546‑M1/49‑22 (MQ=255)
                    ttGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacaggatgtagcg                           <  1:2568447‑M1/49‑23 (MQ=255)
                          cTAAACGCAGGGAGACTCCCCAcctttgaaaacaggatgtagcgatgaa                      <  1:1520463‑M1/49‑28 (MQ=255)
                           tAAACGCAGGGAGACTCCCCAcctttgaaaacaggatgtagcgatgaaa                     >  1:965704‑M1/1‑21 (MQ=255)
                             aaCGCAGGGAGACTCCCCAcctttgaaaacaggatgtagcgatgaaact                   <  1:1681240‑M1/49‑31 (MQ=255)
                             aaCGCAGGGAGACTCCCCAcctttgaaaacaggatgtagcgatgaaa                     <  1:2789126‑M1/47‑29 (MQ=255)
                              aCGCAGGGAGACTCCCCAcctttgaaaacaggatgtagcgatgaaactt                  <  1:947425‑M1/49‑32 (MQ=255)
                              aCGCAGGGAGACTCCCCAcctttgaaaacaggatgtagcgatgaaactt                  <  1:1638807‑M1/49‑32 (MQ=255)
                              aCGCAGGGAGACTCCCCAcctttgaaaacaggatgtagcgatgaaactt                  <  1:825162‑M1/49‑32 (MQ=255)
                                gCAGGGAGACTCCCCAcctttgaaaacaggatgtagcgatgaaacttaa                <  1:2408175‑M1/49‑34 (MQ=255)
                                gCAGGGAGACTCCCCAcctttgaaaacaggatgtagcgatgaaacttaa                <  1:2239337‑M1/49‑34 (MQ=255)
                                gCAGGGAGACTCCCCAcctttgaaaacaggatggagagatgaaaatta                 >  1:930540‑M1/1‑16 (MQ=255)
                                 cAGGGAGACTCCCCAcctttgaaaacaggatgtagcgatgaaacttaac               <  1:1123825‑M1/49‑35 (MQ=255)
                                 cAGGGAGACTCCCCAcctttgaaaacaggatgtagcgatgaaacttaac               >  1:89097‑M1/1‑15 (MQ=255)
                                  aGGGAGACTCCCCAcctttgaaaacaggatgtagcgatgaaacttaacg              <  1:2734477‑M1/49‑36 (MQ=255)
                                  aGGGAGACTCCCCAcctttgaaaacaggatgtagcgatgaaacttaac               >  1:2670003‑M1/1‑14 (MQ=255)
                                     gagaCTCCCCAcctttgaaaacaggatgtagcgatgaaacttaacgaca           <  1:1846650‑M1/49‑39 (MQ=255)
                                      agaCTCCCCAcctttgaaaacaggatgtagcgatgaaacttaacgacag          <  1:2154479‑M1/49‑40 (MQ=255)
                                      agaCTCCCCAcctttgaaaacaggatgtagcgatgaaacttaacgacag          <  1:1717512‑M1/49‑40 (MQ=255)
                                       gaCTCCCCAcctttgaaaacaggatgtagcgatgaaacttaacgacagt         <  1:2351603‑M1/49‑41 (MQ=255)
                                        aCTCCCCAcctttgaaaacaggatgtagcgatgaaacttaacgacagta        <  1:455554‑M1/49‑42 (MQ=255)
                                        aCTCCCCAcctttgaaaacaggatgtagcgatgaaacttaacgacagt         >  1:481101‑M1/1‑8 (MQ=255)
                                         cTCCCCAcctttgaaaacaggatgtagcgatgaaacttaacgacagtaa       <  1:304936‑M1/49‑43 (MQ=255)
                                          cccccacctttgaaaacaggatgtagcgatgaaacttaacgacagtaac      <  1:2460360‑M1/48‑44 (MQ=255)
                                            cccAcctttgaaaacaggatgtagcgatgaaacttaacgacagtaactt    <  1:2579694‑M1/49‑46 (MQ=255)
                                            cccAcctttgaaaacaggatgtagcgatgaaacttaacgacagtaact     <  1:859774‑M1/48‑45 (MQ=255)
                                            cccAcctttgaaaacaggatgtagcgatgaaacttaacgacagtaact     <  1:702873‑M1/48‑45 (MQ=255)
                                             ccAcctttgaaaacaggatgtagcgatgaaacttaacgacagtaactta   <  1:467951‑M1/49‑47 (MQ=255)
                                             ccAcctttgaaaacaggatgtagcgatgaaacttaacgacagtaactta   >  1:666566‑M1/1‑3 (MQ=255)
                                             ccAcctttgaaaacaggatgtagcgatgaaacttaacgacagtaactta   <  1:474912‑M1/49‑47 (MQ=255)
                                              cAcctttgaaaacaggatgtagcgatgaaacttaacgacagtaacttat  >  1:580326‑M1/1‑2 (MQ=255)
                                              cAcctttgaaaacaggatgtagcgatgaaacttaacgacagtaacttat  <  1:1155346‑M1/49‑48 (MQ=255)
                                               |                                               
GAAAAGCTCTATTGTTTACGTTGGGCCTAAACGCAGGGAGACTCCCCATGGCAAGAAAAACCAAGCCGTTAACTGATACGGAAATCAAAGCCGCC  >  NC_000913/2756112‑2756206

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.