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 173 [0] [0] 177 [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)167 (1.140) 74/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

Only 100 of 173 total aligned reads displayed.

AACGAAAAGCTCTATTGTTTACGTTGGGCCTAAACGCAGGGAGACTCCCCATGGCAAGAAAAACCAAGCCGTTAACTGATACGGAAATCAAAGCCGCCAA  >  NC_000913/2756109‑2756208
                                                  |                                                 
aaCGAAAAGCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCA                                                   <  2:3004416‑M1/50‑1 (MQ=255)
aaCGAAAAGCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCA                                                   >  1:1805293‑M1/1‑50 (MQ=255)
 aCGAAAAGCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAc                                                  <  2:3576193‑M1/50‑2 (MQ=255)
  cGAAAAGCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcc                                                 <  1:3878455‑M1/50‑3 (MQ=255)
   gAAAAGCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcct                                                <  2:1788101‑M1/50‑4 (MQ=255)
   gAAAAGCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcct                                                >  2:3943938‑M1/1‑47 (MQ=255)
   gAAAAGCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcct                                                <  1:2398821‑M1/50‑4 (MQ=255)
   gAAAAGCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcct                                                >  2:155441‑M1/1‑47 (MQ=255)
   gAAAAGCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcct                                                >  2:1532038‑M1/1‑47 (MQ=255)
   gAAAAGCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcct                                                <  1:5696415‑M1/50‑4 (MQ=255)
   gAAAAGCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcct                                                <  1:866646‑M1/50‑4 (MQ=255)
   gAAAAGCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcct                                                <  2:106403‑M1/50‑4 (MQ=255)
   gAAAAGCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcct                                                >  2:964954‑M1/1‑47 (MQ=255)
    aaaaGCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcctt                                               >  2:1322085‑M1/1‑46 (MQ=255)
     aaaGCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAccttt                                              >  1:4849388‑M1/1‑45 (MQ=255)
      aaGCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAccttt                                              >  1:4305721‑M1/1‑44 (MQ=255)
      aaGCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAccttt                                              >  2:4934875‑M1/1‑44 (MQ=255)
       aactctatTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttga                                            <  1:3539834‑M1/48‑8 (MQ=255)
       aGCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttga                                            <  1:4603030‑M1/50‑8 (MQ=255)
       aGCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttga                                            <  2:2692759‑M1/50‑8 (MQ=255)
        gCTCTATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaa                                           <  1:4940396‑M1/50‑9 (MQ=255)
         ctctATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaa                                          <  2:2027045‑M1/50‑10 (MQ=255)
         ctctATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaa                                          <  2:4415343‑M1/50‑10 (MQ=255)
         ctctATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaa                                           >  2:5054846‑M1/1‑41 (MQ=255)
            tATTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaa                                         >  2:3418418‑M1/1‑38 (MQ=255)
             aTTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacag                                      >  1:3481207‑M1/1‑37 (MQ=255)
             aTTGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacag                                      >  2:3784857‑M1/1‑37 (MQ=255)
              ttGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacagg                                     <  1:1325064‑M1/50‑15 (MQ=255)
              ttGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacagg                                     <  1:3142116‑M1/50‑15 (MQ=255)
              ttGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacag                                      <  1:5378532‑M1/49‑14 (MQ=255)
              ttGTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaac                                        <  2:2760668‑M1/47‑12 (MQ=255)
                gTTTACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacaggat                                   >  2:1282981‑M1/1‑34 (MQ=255)
                 tttACGTTGGGCTAACCGCAGGGAGATTCCCCAcctttgaaaacaggat                                   <  1:1872424‑M1/49‑17 (MQ=255)
                 tttACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacaggatg                                  <  2:4396924‑M1/50‑18 (MQ=255)
                 tttACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacaggatg                                  <  2:1662442‑M1/50‑18 (MQ=255)
                 tttACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacaggat                                   <  2:654687‑M1/49‑17 (MQ=255)
                  ttACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacaggatg                                  >  2:5005234‑M1/1‑32 (MQ=255)
                   tACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacaggatgt                                 <  2:4027171‑M1/49‑19 (MQ=255)
                   tACGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacaggatg                                  >  1:599956‑M1/1‑31 (MQ=255)
                    aCGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacaggatgta                                <  2:5016630‑M1/49‑20 (MQ=255)
                    aCGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacaggatgta                                <  1:2041178‑M1/49‑20 (MQ=255)
                     cGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacaggatgtagc                              <  2:3426184‑M1/50‑22 (MQ=255)
                     cGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacaggatgtagc                              >  2:4001660‑M1/1‑29 (MQ=255)
                     cGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacaggatgtagc                              <  2:5401559‑M1/50‑22 (MQ=255)
                     cGTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacaggatgtag                               <  2:5471346‑M1/49‑21 (MQ=255)
                      gTTGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacaggatgtagcg                             <  2:4454403‑M1/50‑23 (MQ=255)
                       ttGGGCTAAACGCAGGGAGACTCCCCAcctttgaaaacaggatgtagcga                            >  1:3263996‑M1/1‑27 (MQ=255)
                           ggctaAACGCAGGGAGACTCCCCAcctttgaaaacaggatgtagcgatga                         <  1:1127931‑M1/48‑27 (MQ=255)
                           ggctaAACGCAGGGAGACTCCCCAcctttgaaaacaggatgtagcgatga                         <  2:1987585‑M1/48‑27 (MQ=255)
                           ggctaAACGCAGGGAGACTCCCCAcctttgaaaacaggatgtagcgatg                          >  2:4747443‑M1/3‑24 (MQ=255)
                           ggctaAACGCAGGGAGACTCCCCAcctttgaaaacaggatgtagcgatg                          <  1:3824928‑M1/47‑26 (MQ=255)
                            gctAAACGCAGGGAGACTCCCCAcctttgaaaacaggatgtagcgatga                         >  1:4446673‑M1/2‑23 (MQ=255)
                             cTAAACGCAGGGAGACTCCCCAcctttgaaaacaggatgtagcgatgaaa                       >  1:2616019‑M1/1‑22 (MQ=255)
                              tAAACGCAGGGAGACTCCCCAcctttgaaaacaggatgtagcgatgaaac                      >  2:2008071‑M1/1‑21 (MQ=255)
                               aaaCGCAGGGAGACTCCCCAcctttgaaaacaggatgtagcgatgaaact                     >  1:3352469‑M1/1‑20 (MQ=255)
                                aaCGCAGGGAGACTCCCCAcctttgaaaacaggatgtagcgatgaaactt                    <  1:1115726‑M1/50‑32 (MQ=255)
                                aaCGCAGGGAGACTCCCCAcctttgaaaacaggatgtagcgatgaaactt                    <  1:3444647‑M1/50‑32 (MQ=255)
                                aaCGCAGGGAGACTCCCCAcctttgaaaacaggatgtagcgatgaaactt                    <  2:3238567‑M1/50‑32 (MQ=255)
                                aaCGCAGGGAGACTCCCCAcctttgaaaacaggatgtagcgatgaaact                     <  1:498178‑M1/49‑31 (MQ=255)
                                aaCGCAGGGAGACTCCCCAcctttgaaaacaggatgtagcgatgaaac                      <  2:62900‑M1/48‑30 (MQ=255)
                                aaCGCAGGGAGACTCCCCAcctttgaaaacaggatgtagcgatgaaac                      <  1:198502‑M1/48‑30 (MQ=255)
                                aaCGCAGGGAGACTCCCCAcctttgaaaacaggatgtagcgatgaaac                      <  1:5443297‑M1/48‑30 (MQ=255)
                                 aCGCAGGGAGACTCCCCGcctttgaaaacaggatgtagcgatgaaactta                   <  2:379473‑M1/50‑33 (MQ=255)
                                 aCGCAGGGAGACTCCCCAcctttgaaaacaggatgtagcgatgaaactta                   <  2:1177035‑M1/50‑33 (MQ=255)
                                 aCGCAGGGAGACTCCCCAcctttgaaaacaggatgtagcgatgaaactt                    >  1:5408804‑M1/1‑18 (MQ=255)
                                  cGCAGGGAGACTCCCCAcctttgaaaacaggatgtagcgattaaacttaa                  <  1:28621‑M1/50‑34 (MQ=255)
                                  cGCAGGGAGACTCCCCAcctttgaaaacaggatgtagcgatgaaacttaa                  <  2:112873‑M1/50‑34 (MQ=255)
                                  cGCAGGGAGACTCCCCAcctttgaaaacaggatgtagcgatgaaacttaa                  >  2:852957‑M1/1‑17 (MQ=255)
                                   gCAGGGAGACTCCCCAcctttgaaaacaggatgtagcgatgaaacttaac                 <  1:1504466‑M1/50‑35 (MQ=255)
                                   gCAGGGAGACTCCCCAcctttgaaaacaggatgtagcgatgaaacttaac                 >  1:268198‑M1/1‑16 (MQ=255)
                                    cAGGGAGACTCCCCAcctttgaaaacaggatgtagcgatgaaacttaacg                <  2:705342‑M1/50‑36 (MQ=255)
                                    cAGGGAGACTCCCCAcctttgaaaacaggatgtagcgatgaaacttaacg                <  2:4486591‑M1/50‑36 (MQ=255)
                                     aGGGAGACTCCCCAcctttgaaaacaggatgtagcgatgaaacttaacga               <  1:1612646‑M1/50‑37 (MQ=255)
                                     aGGGAGACTCCCCAcctttgaaaacaggatgtagcgatgaaacttaacga               <  1:3559295‑M1/50‑37 (MQ=255)
                                      gggAGACTCCCCAcctttgaaaacaggatgtagcgatgaaacttaacgac              <  1:5560589‑M1/50‑38 (MQ=255)
                                      gggAGACTCCCCAcctttgaaaacaggatgtagcgatgaaacttaacga               >  1:2496158‑M1/1‑13 (MQ=255)
                                       ggAGACTCCCCAcctttgaaaacaggatgtagcgatgaaacttaacga               >  2:4201273‑M1/1‑12 (MQ=255)
                                        gagaCTCCCCAcctttgaaagcaggatgtagcgatgaaacttaacgacag            <  2:267339‑M1/50‑40 (MQ=255)
                                        gagaCTCCCCAcctttgaaaacaggatgtagcgatgaaacttaacgacag            >  2:4351859‑M1/1‑11 (MQ=255)
                                        gagaCTCCCCAcctttgaaaacaggatgtagcgatgaaacttaacgacag            <  2:1634015‑M1/50‑40 (MQ=255)
                                          gaCTCCCCAcctttgaaaacaggatgtagcgatgaaacttaacgacagta          <  2:3162690‑M1/50‑42 (MQ=255)
                                          gaCTCCCCAcctttgaaaacaggatgtagcgatgaaacttaacgacagt           >  1:2779426‑M1/1‑9 (MQ=255)
                                           aCTCCCCAcctttgaaaacaggatgtagcgatgaaacttaacgacagtaa         >  1:2475520‑M1/1‑8 (MQ=255)
                                           aCTCCCCAcctttgaaaacaggatgtagcgatgaaacttaacgacagtaa         >  1:3998027‑M1/1‑8 (MQ=255)
                                           aCTCCCCAcctttgaaaacaggatgtagcgatgaaacttaacgacagtaa         <  1:1670689‑M1/50‑43 (MQ=255)
                                           aCTCCCCAcctttgaaaacaggatgtagcgatgaaacttaacgacagtaa         <  1:4775705‑M1/50‑43 (MQ=255)
                                           aCTCCCCAcctttgaaaacaggatgtagcgatgaaacttaacgacagtaa         <  1:3432077‑M1/50‑43 (MQ=255)
                                            cTCCCCAcctttgaaaacaggatgtagcgatgaaacttaacgacagtaac        <  1:4618708‑M1/50‑44 (MQ=255)
                                              ccccAcctttgaaaacaggatgtagcgatgaaacttaacgacagtaact       >  1:3770679‑M1/1‑5 (MQ=255)
                                               cccAcctttgaaaacaggatgtagcgatgaaacttaacgacagtaactt      <  1:351606‑M1/49‑46 (MQ=255)
                                               cccAcctttgaaaacaggatgtagcgatgaaacttaacgacagtaactt      <  2:5458163‑M1/49‑46 (MQ=255)
                                               cccAcctttgaaaacaggatgtagcgatgaaacttaacgacagtaactt      <  2:4314678‑M1/49‑46 (MQ=255)
                                                ccAcctttgaaaacaggatgtagcgatgaaacttaacgacagtaacttat    <  2:2810029‑M1/50‑48 (MQ=255)
                                                ccAcctttgaaaacaggatgtagcgatgaaacttaacgacagtaacttat    >  2:3663142‑M1/1‑3 (MQ=255)
                                                ccAcctttgaaaacaggatgtagcgatgaaacttaacgacagtaacttat    <  1:3948290‑M1/50‑48 (MQ=255)
                                                 cAcctttgaaaacaggatgtagcgatgaaacttaacgacagtaacttatt   <  1:5128531‑M1/50‑49 (MQ=255)
                                                 cAcctttgaaaacaggatgtagcgatgaaacttaacgacagtaacttatt   <  2:110119‑M1/50‑49 (MQ=255)
                                                 cAcctttgaaaacaggatgtagcgatgaaacttaacgacagtaacttatt   <  2:2544114‑M1/50‑49 (MQ=255)
                                                 cAcctttgaaaacaggatgtagcgatgaaacttaacgacagtaacttatt   <  2:89490‑M1/50‑49 (MQ=255)
                                                  acctttgaaaacaggatgtagcgatgaaacttaacgacagtaacttattc  <  2:215168‑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 < 20 ≤ ATCG/ATCG < 33 ≤ 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.