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 39 [0] [0] 39 [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)37 (0.480) 23/96 0.9 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

AACGAAAAGCTCTATTGTTTACGTTGGGCCTAAACGCAGGGAGACTCCCCATGGCAAGAAAAACCAAGCCGTTAACTGATACGGAAATCAAAGCCGCCAA  >  NC_000913/2756109‑2756208
                                                  |                                                 
aaCGAAAAGCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCA                                                   <  1:1275569‑M1/50‑1 (MQ=255)
 aCGAAAAGCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAc                                                  <  1:2043793‑M1/50‑2 (MQ=255)
   gAAAAGCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcct                                                <  1:463649‑M1/50‑4 (MQ=255)
     aaaGCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAccttt                                              <  1:1948158‑M1/50‑6 (MQ=255)
      aaGCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttg                                             <  1:3144922‑M1/50‑7 (MQ=255)
       aGCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttga                                            <  1:1912124‑M1/50‑8 (MQ=255)
        gCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaa                                           <  1:2843086‑M1/50‑9 (MQ=255)
        gCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaa                                           <  1:2138591‑M1/50‑9 (MQ=255)
              ttGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacag                                      <  1:4146227‑M1/49‑14 (MQ=255)
              ttGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacag                                      <  1:1633771‑M1/49‑14 (MQ=255)
                 tttACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacaggat                                   <  1:5586390‑M1/49‑17 (MQ=255)
                 tttACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacaggat                                   <  1:5302527‑M1/49‑17 (MQ=255)
                   tACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacaggatgt                                 <  1:320780‑M1/49‑19 (MQ=255)
                    aCGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacaggatgta                                <  1:247029‑M1/49‑20 (MQ=255)
                    aCGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacaggatgta                                <  1:1698181‑M1/49‑20 (MQ=255)
                    aCGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacaggatgta                                <  1:1636020‑M1/49‑20 (MQ=255)
                    aCGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacaggatgta                                <  1:4194364‑M1/49‑20 (MQ=255)
                     cGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacaggatgtagc                              <  1:5351138‑M1/50‑22 (MQ=255)
                      gTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacaggatgtagc                              >  1:5690871‑M1/1‑28 (MQ=255)
                           ggctaAACGCAGGGAGACTCCCCAcctttgaaaacaggatgtagcgatga                         <  1:317766‑M1/48‑27 (MQ=255)
                           ggctaAACGCAGGGAGACTCCCCAcctttgaaaacaggatgtagcgatg                          <  1:2379174‑M1/47‑26 (MQ=255)
                           ggctaAACGCAGGGAGACTCCCCAcctttgaaaacaggatgtagcgatg                          <  1:1462038‑M1/47‑26 (MQ=255)
                            gctAAACGCAGGGAGACTCCCCAcctttgaaaacaggatgtagcgatgaa                        <  1:859155‑M1/49‑28 (MQ=255)
                            gctAAACGCAGGGAGACTCCCCAcctttgaaaacaggatgtagcgatgaa                        <  1:2506883‑M1/49‑28 (MQ=255)
                                aaCGCAGGGAGACTCCCCAcctttgaaaacaggatgtagcgatgaaact                     <  1:297013‑M1/49‑31 (MQ=255)
                                aaCGCAGGGAGACTCCCCAcctttgaaaacaggatgtagcgatgaaac                      <  1:4019494‑M1/48‑30 (MQ=255)
                                 aCGCAGGGAGACTCCCCAcctttgaaaacaggatgtagcgatgaaactta                   <  1:4667078‑M1/50‑33 (MQ=255)
                                   gCAGGGAGACTCCCCAcctttgaaaacaggatgtagcgatgaaacttaac                 <  1:5514878‑M1/50‑35 (MQ=255)
                                   gCAGGGAGACTCCCCAcctttgaaaacaggatgtagcgatgaaacttaac                 <  1:1595385‑M1/50‑35 (MQ=255)
                                     aGGGAGACTCCCCAcctttgaaaacaggatgtagcgatgaaacttaacga               <  1:2122220‑M1/50‑37 (MQ=255)
                                        gagaCTCCCCAcctttgaaaacaggatgtagcgatgaaacttaacgacag            <  1:3724879‑M1/50‑40 (MQ=255)
                                         agaCTCCCCAcctttgaaaacaggatgtagcgatgaaacttaacgacagt           <  1:3105705‑M1/50‑41 (MQ=255)
                                          gaCTCCCCAcctttgaaaacaggatgtagcgatgaaacttaacgacagta          <  1:4116232‑M1/50‑42 (MQ=255)
                                            cTCCCCAcctttgaaaacaggatgtagcgatgaaacttaacgacagtaac        <  1:2214376‑M1/50‑44 (MQ=255)
                                            cTCCCCAcctttgaaaacaggatgtagcgatgaaacttaacgacagtaac        <  1:561443‑M1/50‑44 (MQ=255)
                                               cccAcctttgaaaacaggatgtagcgatgaaacttaacgacagtaactt      <  1:2769537‑M1/49‑46 (MQ=255)
                                                 cAcctttgaaaacaggatgtagcgatgaaacttaacgacagtaacttatt   <  1:3768977‑M1/50‑49 (MQ=255)
                                                 cAcctttgaaaacaggatgtagcgatgaaacttaacgacagtaacttatt   <  1:1598903‑M1/50‑49 (MQ=255)
                                                  acctttgaaaacaggatgtagcgatgaaacttaacgacagtaacttattc  <  1:573633‑M1/50‑50 (MQ=255)
                                                  |                                                 
AACGAAAAGCTCTATTGTTTACGTTGGGCCTAAACGCAGGGAGACTCCCCATGGCAAGAAAAACCAAGCCGTTAACTGATACGGAAATCAAAGCCGCCAA  >  NC_000913/2756109‑2756208

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: