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 26 [0] [0] 27 [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)26 (0.970) 19/94 0.2 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

AGCGGCACGTAGTGTGGATGCCTTACACGCCGCATTTAATCAATAACCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTATTC  >  NC_000913/2791205‑2791299
                                              |                                                
aaaagctctattgtttacgttgggctaaacgcagggagactccccaCCt                                                >  1:460619‑M2/47‑49 (MQ=255)
  aagctctattgtttacgttgggctaaacgcagggagactccccaCCttt                                              <  1:1324050‑M2/5‑1 (MQ=255)
  aagctctattgtttacgttgggctaaacgcagggagactccccaCCttt                                              <  1:1735543‑M2/5‑1 (MQ=255)
  aagctctattgtttacgttgggctaaacgcagggagactccccaCCttt                                              >  1:1426100‑M2/45‑49 (MQ=255)
    gctctattgtttacgttgggctaaacgcagggagactccccaCCTTTg                                             >  1:1753671‑M2/43‑48 (MQ=255)
      tctattgtttacgttgggctaaacgcagggagactccccaCCTTTGaa                                           >  1:802037‑M2/41‑48 (MQ=255)
      tctattgtttacgttgggctaaacgcagggagactccccaCCTTTGaa                                           >  1:1157767‑M2/41‑48 (MQ=255)
             tttacgttgggctaaacgcagggagactccccaCCTTTGAAAACAGGAt                                   <  1:1426755‑M2/16‑1 (MQ=255)
             tttacgttgggctaaacgcagggagactccccaCCTTTGAAAACAGGAt                                   <  1:704744‑M2/16‑1 (MQ=255)
               tacgttgggctaaacgcagggagactccccaCCTTTGAAAACAGGATg                                  <  1:703357‑M2/17‑1 (MQ=255)
               tacgttgggctaaacgcagggagactccccaCCTTTGAAAACAGGATg                                  <  1:58650‑M2/17‑1 (MQ=255)
                 cgttgggctaaacgcagggagactccccaCCTTTGAAAACAGGATGTa                                <  1:603987‑M2/19‑1 (MQ=255)
                   ttgggctaaacgcagggagactccccaCCTTTGAAAACAGGATGTAGCg                             <  1:1294233‑M2/22‑1 (MQ=255)
                   ttgggctaaacgcagggagactccccaCCTTTGAAAACAGGATGTAGCg                             <  1:1317715‑M2/22‑1 (MQ=255)
                        ctaaacgcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGa                         >  1:1598230‑M2/23‑48 (MQ=255)
                           aacgcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACt                     >  1:26028‑M2/20‑49 (MQ=255)
                           aacgcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACt                     >  1:1702935‑M2/20‑49 (MQ=255)
                             cgcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACtt                    >  1:234411‑M2/18‑48 (MQ=255)
                                agggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACg                <  1:582622‑M2/35‑1 (MQ=255)
                                    agactccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAg            <  1:309804‑M2/39‑1 (MQ=255)
                                    agactccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAg            <  1:1311861‑M2/39‑1 (MQ=255)
                                     gactccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGt           <  1:1944449‑M2/40‑1 (MQ=255)
                                      actccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTa          <  1:1449563‑M2/41‑1 (MQ=255)
                                          cccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACt       <  1:1810009‑M2/44‑1 (MQ=255)
                                           ccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTa     >  1:1462581‑M2/4‑49 (MQ=255)
                                           ccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTa     <  1:1417498‑M2/46‑1 (MQ=255)
                                              ccTTTGAAAACAGGATGTAGCGATGAAACTTAACGATAGTAACTTATTc  <  1:1402581/49‑1 (MQ=255)
                                              |                                                
AGCGGCACGTAGTGTGGATGCCTTACACGCCGCATTTAATCAATAACCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTATTC  >  NC_000913/2791205‑2791299

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: