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 177 total aligned reads displayed.

CGAGCGGCACGTAGTGTGGATGCCTTACACGCCGCATTTAATCAATAACCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTATTCC  >  NC_000913/2791203‑2791300
                                                |                                                 
cgaaaagctctattgtttacgttgggctaaacgcagggagactccccacc                                                  >  1:2020708‑M2/49‑50 (MQ=255)
 gaaaagctctattgtttacgttgggctaaacgcagggagactccccaCCt                                                 >  2:964954‑M2/48‑50 (MQ=255)
 gaaaagctctattgtttacgttgggctaaacgcagggagactccccaCCt                                                 <  1:5696415‑M2/3‑1 (MQ=255)
 gaaaagctctattgtttacgttgggctaaacgcagggagactccccaCCt                                                 >  2:1532038‑M2/48‑50 (MQ=255)
 gaaaagctctattgtttacgttgggctaaacgcagggagactccccaCCt                                                 >  2:155441‑M2/48‑50 (MQ=255)
 gaaaagctctattgtttacgttgggctaaacgcagggagactccccaCCt                                                 >  2:4102198‑M2/48‑50 (MQ=255)
 gaaaagctctattgtttacgttgggctaaacgcagggagactccccaCCt                                                 >  2:3943938‑M2/48‑50 (MQ=255)
 gaaaagctctattgtttacgttgggctaaacgcagggagactccccaCCt                                                 <  2:1788101‑M2/3‑1 (MQ=255)
 gaaaagctctattgtttacgttgggctaaacgcagggagactccccaCCt                                                 <  1:2398821‑M2/3‑1 (MQ=255)
 gaaaagctctattgtttacgttgggctaaacgcagggagactccccaCCt                                                 >  2:3324433‑M2/48‑50 (MQ=255)
  aaaagctctattgtttacgttgggctaaacgcagggagactccccaCCtt                                                >  2:1322085‑M2/47‑50 (MQ=255)
    aagctctattgtttacgttgggctaaacgcagggagactccccaCCttt                                               >  1:4305721‑M2/45‑49 (MQ=255)
    aagctctattgtttacgttgggctaaacgcagggagactccccaCCTTTg                                              <  1:3377177‑M2/6‑1 (MQ=255)
    aagctctattgtttacgttgggctaaacgcagggagactccccaCCTTTg                                              <  1:4803839‑M2/6‑1 (MQ=255)
     agctctattgtttacgttgggctaaacgcagggagactccccaCCTTTGa                                             <  2:2692759‑M2/7‑1 (MQ=255)
     agctctattgtttacgttgggctaaacgcagggagactccacaCCTTTg                                              >  2:4236844‑M2/44‑49 (MQ=255)
       ctctattgtttacgttgggctaaacgcagggagactccccaCCTTTGaa                                            >  1:1374739‑M2/42‑49 (MQ=255)
       ctctattgtttacgttgggctaaacgcagggagactccccaCCTTTGaa                                            >  2:5054846‑M2/42‑49 (MQ=255)
         ctattgtttacgttgggctaaacgcagggagactccccaCCTTTGaaaa                                          >  2:4968159‑M2/40‑49 (MQ=255)
          tattgtttacgttgggctaaacgcagggagactccccaCCTTTGaaaa                                          >  2:3418418‑M2/39‑48 (MQ=255)
          tattgtttacgttgggctaaacgcagggagactccccaCCTTTGaaaa                                          >  2:1094098‑M2/39‑48 (MQ=255)
           attgtttacgttgggctaaacgcagggagactccccaCCTTTGAAAACAg                                       >  1:3481207‑M2/38‑50 (MQ=255)
            ttgtttacgttgggctaaacgcagggagactccccaCCTTTGAAAACAg                                       <  1:5378532‑M2/13‑1 (MQ=255)
            ttgtttacgttgggctaaacgcagggagactccccaCCTTTGAAAACAg                                       <  2:1243084‑M2/13‑1 (MQ=255)
             tgtttacgttgggctaaacgcagggagactccccaCCTTTGAAAACAGGa                                     >  1:4912317‑M2/36‑50 (MQ=255)
               tttacgttgggctaaccgcagggagattccccaCCTTTGAAAACAGGAt                                    <  1:1872424‑M2/16‑1 (MQ=255)
               tttacgttgggctaaacgcagggagactccccaCCTTTGAAAACAGGAt                                    <  2:654687‑M2/16‑1 (MQ=255)
               tttacgttgggctaaacgcagggagactccccaCCTTTGAAAACAGGATg                                   >  1:3942319‑M2/34‑50 (MQ=255)
               tttacgttgggctaaacgcagggagactccccaCCTTTGAAAACAGGATg                                   <  2:4396924‑M2/17‑1 (MQ=255)
                 tacgttgggctaaacgcagggagactccccaCCTTTGAAAACAGGATGt                                  <  1:4358845‑M2/18‑1 (MQ=255)
                 tacgttgggctaaacgcagggagactccccaCCTTTGAAAACAGGATGt                                  <  2:4027171‑M2/18‑1 (MQ=255)
                  acgttgggctaaacgcagggagactccccaCCTTTGAAAACAGGATGTa                                 <  1:2041178‑M2/19‑1 (MQ=255)
                  acgttgggctaaacgcagggagactccccaCCTTTGAAAACAGGATGTa                                 <  2:5016630‑M2/19‑1 (MQ=255)
                   cgttgggctaaacgcagggagactccccaCCTTTGAAAACAGGATGTAg                                <  1:75545‑M2/20‑1 (MQ=255)
                   cgttgggctaaacgcagggagactccccaCCTTTGAAAACAGGATGTAg                                <  2:5471346‑M2/20‑1 (MQ=255)
                   cgttgggctaaacgcagggagactccccaCCTTTGAAAACAGGATGTAGc                               <  1:4581826‑M2/21‑1 (MQ=255)
                   cgttgggctaaacgcagggagactccccaCCTTTGAAAACAGGATGTAGc                               <  2:5401559‑M2/21‑1 (MQ=255)
                   cgttgggctaaacgcagggagactccccaCCTTTGAAAACAGGATGTAGc                               >  2:4001660‑M2/30‑50 (MQ=255)
                    gttgggctaaacgcagggagactccccaCCTTTGAAAACAGGATGTAGc                               >  2:1167418‑M2/29‑49 (MQ=255)
                    gttgggctaaacgcagggagactccccaCCTTTGAAAACAGGATGTAGCg                              <  2:2112509‑M2/22‑1 (MQ=255)
                    gttgggctaaacgcagggagactccccaCCTTTGAAAACAGGATGTAGCg                              <  2:1691535‑M2/22‑1 (MQ=255)
                     ttgggctaaacgcagggagactccccaCCTTTGAAAACAGGATGTAGCGa                             >  1:3263996‑M2/28‑50 (MQ=255)
                     ttgggctaaacgcagggagactccccaCCTTTGAAAACAGGATGTAGCGa                             <  2:3231876‑M2/23‑1 (MQ=255)
                        ggctaaacgcagggagactccccaCCTTTGAAAACAGGATGTAGCGATg                           >  2:4469758‑M2/25‑49 (MQ=255)
                        ggctaaacgcagggagactccccaCCTTTGAAAACAGGATGTAGCGATg                           >  2:4747443‑M2/25‑49 (MQ=255)
                        ggctaaacgcagggagactccccaCCTTTGAAAACAGGATGTAGCGATg                           <  1:3824928‑M2/25‑1 (MQ=255)
                        ggctaaacgcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGa                          <  2:1987585‑M2/26‑1 (MQ=255)
                        ggctaaacgcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGa                          <  1:1127931‑M2/26‑1 (MQ=255)
                         gctaaacgcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGaa                         <  1:2818538‑M2/27‑1 (MQ=255)
                          ctaaacgcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGaaa                        <  2:2806419‑M2/28‑1 (MQ=255)
                           taaacgcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAAc                       >  2:2008071‑M2/22‑50 (MQ=255)
                            aaacgcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACt                      >  1:3521058‑M2/21‑50 (MQ=255)
                             aacgcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAAc                       <  2:62900‑M2/29‑1 (MQ=255)
                             aacgcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAAc                       <  2:2409626‑M2/29‑1 (MQ=255)
                             aacgcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAAc                       <  1:5443297‑M2/29‑1 (MQ=255)
                             aacgcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAAc                       <  1:3438983‑M2/29‑1 (MQ=255)
                             aacgcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACtt                     <  1:1115726‑M2/31‑1 (MQ=255)
                             aacgcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACt                      <  2:373994‑M2/30‑1 (MQ=255)
                             aacgcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACt                      <  1:498178‑M2/30‑1 (MQ=255)
                             aacgcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACt                      <  1:3554347‑M2/30‑1 (MQ=255)
                              acgcagggagactccccgCCTTTGAAAACAGGATGTAGCGATGAAACTTa                    <  2:379473‑M2/32‑1 (MQ=255)
                              acgcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACtt                     >  1:5408804‑M2/19‑49 (MQ=255)
                              acgcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTa                    <  1:926528‑M2/32‑1 (MQ=255)
                               cgcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTaa                   >  2:852957‑M2/18‑50 (MQ=255)
                               cgcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTaa                   >  1:3468026‑M2/18‑50 (MQ=255)
                               cgcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTaa                   >  2:3593373‑M2/18‑50 (MQ=255)
                               cgcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTaa                   >  1:1663457‑M2/18‑50 (MQ=255)
                                gcagggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAAc                  >  1:268198‑M2/17‑50 (MQ=255)
                                 cagggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAAc                  >  2:3005222‑M2/16‑49 (MQ=255)
                                 cagggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACg                 <  1:2563487‑M2/35‑1 (MQ=255)
                                 cagggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACg                 <  2:5418252‑M2/35‑1 (MQ=255)
                                  agggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGa                <  1:1612646‑M2/36‑1 (MQ=255)
                                   gggagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGAc               <  1:5560589‑M2/37‑1 (MQ=255)
                                     gagactccccaCCTTTGAAAGCAGGATGTAGCGATGAAACTTAACGACAg             <  2:267339‑M2/39‑1 (MQ=255)
                                     gagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAg             <  2:1634015‑M2/39‑1 (MQ=255)
                                     gagactccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAg             <  1:3088969‑M2/39‑1 (MQ=255)
                                      agactccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGt            <  2:47651‑M2/40‑1 (MQ=255)
                                       gactccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTa           <  1:3230735‑M2/41‑1 (MQ=255)
                                        actccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTaa          >  1:3998027‑M2/9‑50 (MQ=255)
                                        actccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTaa          >  1:2475520‑M2/9‑50 (MQ=255)
                                        actccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTaa          <  1:1747539‑M2/42‑1 (MQ=255)
                                        actccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTaa          <  1:4775705‑M2/42‑1 (MQ=255)
                                         ctccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAAc         <  1:4618708‑M2/43‑1 (MQ=255)
                                          tccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACt        >  1:2628297‑M2/7‑50 (MQ=255)
                                           ccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACt        >  1:3770679‑M2/6‑49 (MQ=255)
                                           ccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACt        >  2:4408512‑M2/6‑49 (MQ=255)
                                           ccccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGCAACt        >  2:3488390‑M2/6‑49 (MQ=255)
                                            cccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACt        <  2:4314984‑M2/44‑1 (MQ=255)
                                            cccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTa      <  2:681537‑M2/46‑1 (MQ=255)
                                             ccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTAt     <  2:4453371‑M2/47‑1 (MQ=255)
                                             ccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTAt     <  2:2810029‑M2/47‑1 (MQ=255)
                                             ccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTAt     <  1:3829807‑M2/47‑1 (MQ=255)
                                             ccaCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTAt     >  2:3663142‑M2/4‑50 (MQ=255)
                                              caCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTAtt    <  1:5128531‑M2/48‑1 (MQ=255)
                                              caCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTAtt    <  2:110119‑M2/48‑1 (MQ=255)
                                              caCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTAtt    <  2:89490‑M2/48‑1 (MQ=255)
                                               aCCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTATTc   <  2:215168‑M2/49‑1 (MQ=255)
                                                ccTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTATTcc  <  1:1287979/50‑1 (MQ=255)
                                                ccTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTATTcc  >  1:5700857/1‑50 (MQ=255)
                                                ccTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTATTcc  <  2:1531653/50‑1 (MQ=255)
                                                |                                                 
CGAGCGGCACGTAGTGTGGATGCCTTACACGCCGCATTTAATCAATAACCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTATTCC  >  NC_000913/2791203‑2791300

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ 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.