Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 3,761,993 Δ8,250 bp rhsAyibV rhsA, yibA, yibJ, yibG, yibS, yibV

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 3761993 3770242 8250 117 [0] [0] 117 rhsA–yibV rhsA,yibA,yibJ,yibG,yibS,yibV

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 37619920 (0.000)115 (1.200) 64/96 0.0 100% intergenic (‑37/‑191) yibF/rhsA glutathione transferase‑like protein YibF/rhs element protein RhsA
?NC_000913 3770243 = 0 (0.000)coding (1137/1137 nt) yibH inner membrane protein YibH

Only 100 of 117 total aligned reads displayed.

CGACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGATGGTTTGTTACTTCCAA  >  NC_000913/3761943‑3762041
                                                 |                                                 
cgacgaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGC                                                   <  1:5071313‑M1/50‑1 (MQ=255)
 gacgaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCt                                                  >  1:5168727‑M1/1‑49 (MQ=255)
 gacgaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCt                                                  <  1:6792840‑M1/50‑2 (MQ=255)
  acgaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCt                                                  >  1:3737538‑M1/1‑48 (MQ=255)
  acgaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCt                                                  >  1:479303‑M1/1‑48 (MQ=255)
    gaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaa                                               <  1:3727799‑M1/50‑5 (MQ=255)
    gaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaa                                               >  1:6217464‑M1/1‑46 (MQ=255)
    gaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaa                                               >  1:4882094‑M1/1‑46 (MQ=255)
     aGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaat                                              <  1:520780‑M1/50‑6 (MQ=255)
     aGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaat                                              >  1:3632827‑M1/1‑45 (MQ=255)
     aGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaat                                              <  1:5979618‑M1/50‑6 (MQ=255)
     aGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaat                                              <  1:4906573‑M1/50‑6 (MQ=255)
      gTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatg                                             <  1:4236270‑M1/50‑7 (MQ=255)
      gTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatg                                             <  1:3821856‑M1/50‑7 (MQ=255)
         tCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgat                                           <  1:3998888‑M1/49‑9 (MQ=255)
          cATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatc                                          >  1:5977379‑M1/1‑40 (MQ=255)
           aTGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatcc                                         <  1:4106‑M1/49‑11 (MQ=255)
           aTGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatcc                                         >  1:6971346‑M1/1‑39 (MQ=255)
            tGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaa                                       <  1:7049733‑M1/50‑13 (MQ=255)
            tGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaa                                       <  1:6753086‑M1/50‑13 (MQ=255)
              ccGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaat                                     <  1:6048320‑M1/50‑15 (MQ=255)
              ccGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaat                                     <  1:6225895‑M1/50‑15 (MQ=255)
              ccGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaa                                      >  1:2147091‑M1/1‑36 (MQ=255)
              ccGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatacaaa                                      <  1:1307465‑M1/49‑14 (MQ=255)
                gAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataa                                   <  1:6815182‑M1/50‑17 (MQ=255)
                gAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataa                                   <  1:1673962‑M1/50‑17 (MQ=255)
                 aGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaata                                    >  1:3750126‑M1/1‑33 (MQ=255)
                  gTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaag                                 <  1:2430095‑M1/50‑19 (MQ=255)
                  gTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaag                                 <  1:6364522‑M1/50‑19 (MQ=255)
                   tCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaaga                                <  1:2512351‑M1/50‑20 (MQ=255)
                   tCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaata                                    >  1:475307‑M1/1‑31 (MQ=255)
                   tCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaata                                    >  1:4769265‑M1/1‑31 (MQ=255)
                    ccTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaaga                                >  1:3085040‑M1/1‑30 (MQ=255)
                     cTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagata                              <  1:691702‑M1/50‑22 (MQ=255)
                     cTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagata                              <  1:6961691‑M1/50‑22 (MQ=255)
                     cTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagat                               <  1:222547‑M1/49‑21 (MQ=255)
                       ttGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataat                            >  1:1151528‑M1/1‑27 (MQ=255)
                       ttGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataat                            >  1:1805946‑M1/1‑27 (MQ=255)
                        tgtgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataat                            <  1:7088509‑M1/49‑24 (MQ=255)
                        tgtgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataat                            <  1:3509341‑M1/49‑24 (MQ=255)
                        tgtgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataat                            <  1:4443871‑M1/49‑24 (MQ=255)
                        tgtgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataa                             <  1:3237360‑M1/48‑23 (MQ=255)
                        tgtgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataa                             <  1:6100198‑M1/48‑23 (MQ=255)
                         gtgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatg                           >  1:2533235‑M1/1‑25 (MQ=255)
                          tgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgca                         <  1:288594‑M1/50‑27 (MQ=255)
                          tgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgca                         >  1:5860306‑M1/1‑24 (MQ=255)
                          tgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgca                         <  1:3480895‑M1/50‑27 (MQ=255)
                          tgCGAGGAAAAATATCAGTAAGGCtcaatgatccaaataaagataatgca                         <  1:5645494‑M1/50‑27 (MQ=255)
                            cGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatc                       >  1:6455548‑M1/1‑22 (MQ=255)
                            cGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatc                       >  1:6845187‑M1/1‑22 (MQ=255)
                            cGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatc                       <  1:2938780‑M1/50‑29 (MQ=255)
                            cGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatc                       <  1:1356804‑M1/50‑29 (MQ=255)
                            cGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcat                        <  1:4621956‑M1/49‑28 (MQ=255)
                            cGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcat                        <  1:6855591‑M1/49‑28 (MQ=255)
                            cGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcat                        <  1:3512560‑M1/49‑28 (MQ=255)
                             gAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatcc                      <  1:1432166‑M1/50‑30 (MQ=255)
                             gAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatcc                      <  1:4979230‑M1/50‑30 (MQ=255)
                             gAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatc                       >  1:5836094‑M1/1‑21 (MQ=255)
                              aGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatcca                     <  1:3684199‑M1/50‑31 (MQ=255)
                              aGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatcca                     <  1:4662206‑M1/50‑31 (MQ=255)
                              aGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatc                       >  1:3263693‑M1/1‑20 (MQ=255)
                               ggAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatccag                    >  1:783155‑M1/1‑19 (MQ=255)
                               ggAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatccag                    <  1:3337261‑M1/50‑32 (MQ=255)
                               ggAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatccag                    <  1:6317455‑M1/50‑32 (MQ=255)
                               ggAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatccag                    <  1:1945036‑M1/50‑32 (MQ=255)
                                gAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatccagc                   <  1:3734480‑M1/50‑33 (MQ=255)
                                gAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatccagc                   >  1:6471440‑M1/1‑18 (MQ=255)
                                 aaaaaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagct                  <  1:871160‑M1/50‑34 (MQ=255)
                                 aaaaaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagct                  >  1:2721667‑M1/1‑17 (MQ=255)
                                  aaaaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagctg                 <  1:626348‑M1/50‑35 (MQ=255)
                                   aaaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagctgg                <  1:5564696‑M1/50‑36 (MQ=255)
                                   aaaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagctgg                >  1:606976‑M1/1‑15 (MQ=255)
                                    aaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagctggt               >  1:5556105‑M1/1‑14 (MQ=255)
                                     atatCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtc              <  1:924267‑M1/50‑38 (MQ=255)
                                     atatCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtc              <  1:1597913‑M1/50‑38 (MQ=255)
                                     atatCAGTATGGCtcaatgatccaaataaagataaggcatccagctggtc              <  1:6845291‑M1/50‑38 (MQ=255)
                                      tatCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtca             >  1:4541283‑M1/1‑12 (MQ=255)
                                        tCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtcatt           <  1:3147211‑M1/50‑41 (MQ=255)
                                        tCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtcat            <  1:2218969‑M1/49‑40 (MQ=255)
                                        tCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtcat            <  1:1457855‑M1/49‑40 (MQ=255)
                                         cAGTATGGCtcaatgatccaaataaagataatgcatccagctggtcattc          >  1:2586725‑M1/1‑9 (MQ=255)
                                          aGTATGGCtcaatgatccaaataaagataatgcatccagctggtcattct         <  1:6242477‑M1/50‑43 (MQ=255)
                                          aGTATGGCtcaatgatccaaataaagataatgcatccagctggtcattct         >  1:5202602‑M1/1‑8 (MQ=255)
                                          aGTATGGCtcaatgatccaaataaagataatgcatccagctggtcattct         >  1:3771214‑M1/1‑8 (MQ=255)
                                           gTATGGCtcaatgatccaaataaagataatgcatccagctggtcattctt        <  1:5045426‑M1/50‑44 (MQ=255)
                                           gTATGGCtcaatgatccaaataaagataatgcatccagctggtcattctt        <  1:1580619‑M1/50‑44 (MQ=255)
                                           gTATGGCtcaatgatccaaataaagataatgcatccagctggtcattctt        <  1:4768396‑M1/50‑44 (MQ=255)
                                           gTATGGCtcaatgatccaaataaagataatgcatccagctggtcattct         >  1:4510045‑M1/1‑7 (MQ=255)
                                           gTATGGCtcaatgatccaaataaagataatgcatccagctggtcattct         >  1:7057501‑M1/1‑7 (MQ=255)
                                           gTATGGCtcaatgatccaaataaagataatgcatccagctggtcattct         >  1:2366382‑M1/1‑7 (MQ=255)
                                           gTATGGCtcaatgatccaaataaagataatgcatccagctggtcattct         >  1:2239964‑M1/1‑7 (MQ=255)
                                            tATGGCtcaatgatccaaataaagataatgcatccagctggtcattctta       <  1:5765537‑M1/50‑45 (MQ=255)
                                            tATGGCtcaatgatccaaataaagataatgcatccagctggtcattctta       <  1:6409345‑M1/50‑45 (MQ=255)
                                              tGGCtcaatgatccaaataaagataatgcatccagctggtcattcttagc     <  1:7434576‑M1/50‑47 (MQ=255)
                                              tGGCtcaatgatccaaataaagataatgcatccagctggtcattcttag      <  1:3638083‑M1/49‑46 (MQ=255)
                                              tGGCtcaatgatccaaataaagataatgcatccagctggtcattcttag      >  1:2644187‑M1/1‑4 (MQ=255)
                                                gCtcaatgatccaaataaagataatgcatccagctggtcattcttagcag   <  1:1618935‑M1/50‑49 (MQ=255)
                                                gCtcaatgatccaaataaagataatgcatccagctggtcattcttagcag   <  1:100403‑M1/50‑49 (MQ=255)
                                                gCtcaatgatccaaataaagataatgcatccagctggtcattcttagca    <  1:7415628‑M1/49‑48 (MQ=255)
                                                 ctcaatgatccaaataaagataatgcatccagctggtcattcttagcagc  <  1:5229765‑M1/50‑50 (MQ=255)
                                                 |                                                 
CGACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGATGGTTTGTTACTTCCAA  >  NC_000913/3761943‑3762041

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