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 86 [0] [0] 89 [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)82 (0.860) 52/96 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

GCGAGCGGCACGTAGTGTGGATGCCTTACACGCCGCATTTAATCAATAACCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTATTCC  >  NC_000913/2791202‑2791300
                                                 |                                                 
acgaaaagctctattgtttacgttgggctaaacgcagggagactccccac                                                   <  1:354608‑M2/1‑1 (MQ=255)
acgaaaagctctattgtttacgttgggctaaacgcagggagactccccac                                                   >  1:5572700‑M2/50‑50 (MQ=255)
acgaaaagctctattgtttacgttgggctaaacgcagggagactccccac                                                   >  1:2407317‑M2/50‑50 (MQ=255)
 cgaaaagctctattgtttacgttgggctaaacgcagggagactccccacc                                                  <  1:4775738‑M2/2‑1 (MQ=255)
  gaaaagctctattgtttacgttgggctaaacgcagggagactccccaCCt                                                 >  1:6321282‑M2/48‑50 (MQ=255)
  gaaaagctctattgtttacgttgggctaaacgcagggagactccccaCCt                                                 >  1:7257527‑M2/48‑50 (MQ=255)
   aaaagctctattgtttacgttgggctaaacgcagggagactccccaCCtt                                                >  1:4202561‑M2/47‑50 (MQ=255)
    aaagctctattgtttacgttgggctaaacgcagggagactccccaCCttt                                               <  1:261280‑M2/5‑1 (MQ=255)
    aaagctctattgtttacgttgggctaaacgcagggagactccccaCCttt                                               <  1:2633908‑M2/5‑1 (MQ=255)
    aaagctctattgtttacgttgggctaaacgcagggagactccccaCCttt                                               >  1:5492408‑M2/46‑50 (MQ=255)
    aaagctctattgtttacgttgggctaaacgcagggagactccccaCCttt                                               <  1:619703‑M2/5‑1 (MQ=255)
    aaagctctattgattacgttgggctaaacgcagggagactccccaCCttt                                               <  1:7307705‑M2/5‑1 (MQ=255)
      agctctattgtttacgttgggctaaacgcagggagactccccaCCTTTGa                                             <  1:6685315‑M2/7‑1 (MQ=255)
       gctctattgtttacgttgggctaaacgcagggagactccccaCCTTTGaa                                            <  1:6745377‑M2/8‑1 (MQ=255)
           tattgtttacgttgggctaaacgcagggagactccccaCCTTTGaaaa                                          <  1:4821009‑M2/10‑1 (MQ=255)
           tattgtttacgttgggctaaacgcagggagactccccaCCTTTGaaaa                                          <  1:2596026‑M2/10‑1 (MQ=255)
            attgtttacgttgggctaaacgcagggagactccccaCCTTTGAAAACAg                                       >  1:5461796‑M2/38‑50 (MQ=255)
             ttgtttacgttgggctaaacgcagggagactccccaCCTTTGAAAAc                                         <  1:4455397‑M2/11‑1 (MQ=255)
             ttgtttacgttgggctaaacgcagggagactccccaCCTTTGAAAACa                                        <  1:2623662‑M2/12‑1 (MQ=255)
             ttgtttacgttgggctaaacgcagggagactccccaCCTTTGAAAACa                                        <  1:1309721‑M2/12‑1 (MQ=255)
             ttgtttacgttgggctaaacgcagggagactccccaCCTTTGAAAACAgg                                      <  1:6008071‑M2/14‑1 (MQ=255)
             ttgtttacgttgggctaaacgcagggagactccccaCCTTTGAAAACAg                                       <  1:3141240‑M2/13‑1 (MQ=255)
             ttgtttacgttgggctaaacgcagggagactccccaCCTTTGAAAACAg                                       >  1:3369407‑M2/37‑49 (MQ=255)
             ttgtttacgttgggctaaacgcagggagactccccaCCTTTGAAAACAg                                       >  1:6940793‑M2/37‑49 (MQ=255)
              tgtttacgttgggctaaacgcagggagactccccaCCTTTGAAAACAGGa                                     >  1:3390696‑M2/36‑50 (MQ=255)
               gtttacgttgggctaaacgcagggagactccccaCCTTTGAAAACAGGa                                     <  1:2135677‑M2/15‑1 (MQ=255)
                tttacgttgggctaaacgcagggagactccccaCCTTTGAAAACAGGAt                                    <  1:1215201‑M2/16‑1 (MQ=255)
                  tacgttgggctaaacgcagggagactccccaCCTTTGAAAACAGGATg                                   >  1:6436045‑M2/32‑48 (MQ=255)
                  tacgttgggctaaacgcagggagactccccaCCTTTGAAAACAGGATg                                   >  1:3715142‑M2/32‑48 (MQ=255)
                   acgttgggctaaacgcagggagactccccaCCTTTGAAAACAGGATGTa                                 <  1:126735‑M2/19‑1 (MQ=255)
                   acgttgggctaaacgcagggagactccccaCCTTTGAAAACAGGATGTa                                 <  1:5348202‑M2/19‑1 (MQ=255)
                   acgttgggctaaacgcagggagactccccaCCTTTGAAAACAGGATGTa                                 <  1:7363130‑M2/19‑1 (MQ=255)
                    cgttgggctaaacgcagggagactccccaCCTTTGAAAACAGGATGTAGc                               <  1:1165538‑M2/21‑1 (MQ=255)
                    cgttgggctaaacgcagggagactccccaCCTTTGAAAACAGGATGTAGc                               >  1:6137964‑M2/30‑50 (MQ=255)
                    cgttgggctaaacgcagggagactccccaCCTTTGAAAACAGGATGTAGc                               <  1:6311017‑M2/21‑1 (MQ=255)
                     gttgggctaaacgcagggagactccccaCCTTTGAAAACAGGATGTAGc                               >  1:4814675‑M2/29‑49 (MQ=255)
                     gttgggctaaacgcagggagactccccaCCTTTGAAAACAGGATGTAGCg                              <  1:6986476‑M2/22‑1 (MQ=255)
                     gttgggctaaacgcagggagactccccaCCTTTGAAAACAGGATGTAGCg                              <  1:4899791‑M2/22‑1 (MQ=255)
                         ggctaaacgcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGa                          <  1:6334048‑M2/26‑1 (MQ=255)
                            taaacgcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAAc                       >  1:4384087‑M2/22‑50 (MQ=255)
                            taaacgcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAAc                       >  1:4445865‑M2/22‑50 (MQ=255)
                            taaacgcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAAc                       >  1:3723719‑M2/22‑50 (MQ=255)
                            taaacgcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAAc                       >  1:4007271‑M2/22‑50 (MQ=255)
                            taaacgcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAAc                       >  1:7328264‑M2/22‑50 (MQ=255)
                             aaacgcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACt                      <  1:921487‑M2/30‑1 (MQ=255)
                             aaacgcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACt                      >  1:7214993‑M2/21‑50 (MQ=255)
                              aacgcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAAc                       <  1:7172379‑M2/29‑1 (MQ=255)
                              aacgcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACt                      <  1:4217418‑M2/30‑1 (MQ=255)
                               acgcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACtt                     >  1:7347252‑M2/19‑49 (MQ=255)
                               acgcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTa                    <  1:489780‑M2/32‑1 (MQ=255)
                                cgcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTaa                   <  1:1746488‑M2/33‑1 (MQ=255)
                                cgcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTaa                   >  1:6563961‑M2/18‑50 (MQ=255)
                                 gcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAAc                  <  1:90672‑M2/34‑1 (MQ=255)
                                 gcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAAc                  <  1:6864259‑M2/34‑1 (MQ=255)
                                 gcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAAc                  <  1:5713120‑M2/34‑1 (MQ=255)
                                  cagggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACg                 <  1:4220950‑M2/35‑1 (MQ=255)
                                  cagggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACg                 <  1:3936209‑M2/35‑1 (MQ=255)
                                   agggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGa                >  1:84967‑M2/15‑50 (MQ=255)
                                   agggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGa                >  1:4530343‑M2/15‑50 (MQ=255)
                                    gggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGa                >  1:7029471‑M2/14‑49 (MQ=255)
                                    gggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGAc               <  1:5217734‑M2/37‑1 (MQ=255)
                                      gagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAg             <  1:781460‑M2/39‑1 (MQ=255)
                                      gagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAg             >  1:2454699‑M2/12‑50 (MQ=255)
                                      gagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAg             >  1:6970198‑M2/12‑50 (MQ=255)
                                      gagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAg             <  1:1679448‑M2/39‑1 (MQ=255)
                                       agactccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAg             >  1:1057107‑M2/11‑49 (MQ=255)
                                       agactccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGt            <  1:2772830‑M2/40‑1 (MQ=255)
                                       agactccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGt            <  1:4824799‑M2/40‑1 (MQ=255)
                                         actccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTaa          >  1:4365742‑M2/9‑50 (MQ=255)
                                         actccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTaa          <  1:7129222‑M2/42‑1 (MQ=255)
                                          ctccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAAc         <  1:3348877‑M2/43‑1 (MQ=255)
                                          ctccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAAc         <  1:5382161‑M2/43‑1 (MQ=255)
                                          ctccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAAc         >  1:5621‑M2/8‑50 (MQ=255)
                                            ccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACt        >  1:3075870‑M2/6‑49 (MQ=255)
                                             cccgCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACtt       <  1:3540737‑M2/45‑1 (MQ=255)
                                             cccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACtt       <  1:3200572‑M2/45‑1 (MQ=255)
                                             cccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACtt       <  1:5307069‑M2/45‑1 (MQ=255)
                                             cccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACt        <  1:1437001‑M2/44‑1 (MQ=255)
                                             cccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTa      <  1:6871694‑M2/46‑1 (MQ=255)
                                             cccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTa      <  1:317769‑M2/46‑1 (MQ=255)
                                             cccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTa      >  1:6321449‑M2/5‑50 (MQ=255)
                                              ccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTAt     <  1:2802580‑M2/47‑1 (MQ=255)
                                                aCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTATTc   <  1:157686‑M2/49‑1 (MQ=255)
                                                aCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTATTc   <  1:5559905‑M2/49‑1 (MQ=255)
                                                aCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTATTc   >  1:5908839‑M2/2‑50 (MQ=255)
                                                aCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTATTc   <  1:404904‑M2/49‑1 (MQ=255)
                                                 ccTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTAtttc  >  1:3611888/1‑48 (MQ=255)
                                                 ccTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTATTcc  <  1:6354218/50‑1 (MQ=255)
                                                 ccTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTATTcc  >  1:1664133/1‑50 (MQ=255)
                                                 |                                                 
GCGAGCGGCACGTAGTGTGGATGCCTTACACGCCGCATTTAATCAATAACCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTATTCC  >  NC_000913/2791202‑2791300

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: