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 134 [0] [0] 132 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)127 (1.300) 70/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 134 total aligned reads displayed.

CGACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGATGGTTTGTTACTTCCAA  >  NC_000913/3761943‑3762041
                                                 |                                                 
cgacgaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGC                                                   >  1:4883855‑M1/1‑50 (MQ=255)
cgacgaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGC                                                   >  1:226380‑M1/1‑50 (MQ=255)
cgacgaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGC                                                   <  1:3813907‑M1/50‑1 (MQ=255)
 gacgaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCt                                                  <  1:4739456‑M1/50‑2 (MQ=255)
 gacgaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCt                                                  >  1:1342559‑M1/1‑49 (MQ=255)
 gacgaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCt                                                  <  1:2483734‑M1/50‑2 (MQ=255)
  acgaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtc                                                 <  1:5501144‑M1/50‑3 (MQ=255)
   cgaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtca                                                <  1:2531758‑M1/50‑4 (MQ=255)
   cgaGTTTCATGCCAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaa                                               <  1:1421250‑M1/50‑5 (MQ=255)
    gaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaa                                               >  1:1559333‑M1/1‑46 (MQ=255)
     aGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaat                                              <  1:1686372‑M1/50‑6 (MQ=255)
     aGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaat                                              <  1:6029742‑M1/50‑6 (MQ=255)
     aGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaat                                              <  1:5981026‑M1/50‑6 (MQ=255)
      gTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatg                                             <  1:6517729‑M1/50‑7 (MQ=255)
      gTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatg                                             <  1:1182130‑M1/50‑7 (MQ=255)
      gTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatg                                             <  1:6581510‑M1/50‑7 (MQ=255)
       tttCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatga                                            <  1:71358‑M1/50‑8 (MQ=255)
        ttCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgat                                           >  1:5749864‑M1/1‑42 (MQ=255)
         tCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgat                                           <  1:4252873‑M1/49‑9 (MQ=255)
          cATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatc                                          <  1:6359121‑M1/49‑10 (MQ=255)
          cATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatc                                          >  1:317774‑M1/1‑40 (MQ=255)
           aTGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatcca                                        <  1:4395314‑M1/50‑12 (MQ=255)
           aTGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatcca                                        <  1:267063‑M1/50‑12 (MQ=255)
            tGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCttaatgatccaa                                       <  1:5301122‑M1/50‑13 (MQ=255)
              ccGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgttccaaa                                      <  1:1232369‑M1/49‑14 (MQ=255)
              ccGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaat                                     <  1:2115326‑M1/50‑15 (MQ=255)
              ccGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaa                                      >  1:1089317‑M1/1‑36 (MQ=255)
              ccGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaa                                      <  1:3689200‑M1/49‑14 (MQ=255)
                 aGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaa                                  <  1:7189071‑M1/50‑18 (MQ=255)
                  gTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaata                                    >  1:318681‑M1/1‑32 (MQ=255)
                   tCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaaga                                <  1:2462356‑M1/50‑20 (MQ=255)
                   tCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaaga                                <  1:1094618‑M1/50‑20 (MQ=255)
                    ccTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaaga                                >  1:2725437‑M1/1‑30 (MQ=255)
                     cTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagata                              <  1:7008641‑M1/50‑22 (MQ=255)
                     cTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagat                               >  1:3481402‑M1/1‑29 (MQ=255)
                       ttGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataat                            >  1:4206755‑M1/1‑27 (MQ=255)
                       ttGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataat                            >  1:5818457‑M1/1‑27 (MQ=255)
                        tgtgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataat                            <  1:2411578‑M1/49‑24 (MQ=255)
                        tgtgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataa                             <  1:4852900‑M1/48‑23 (MQ=255)
                        tgtgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataa                             <  1:5993525‑M1/48‑23 (MQ=255)
                        tgtgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataa                             <  1:38809‑M1/48‑23 (MQ=255)
                         gtgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatg                           <  1:338258‑M1/49‑25 (MQ=255)
                          tgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgca                         <  1:4346215‑M1/50‑27 (MQ=255)
                          tgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgca                         >  1:1779180‑M1/1‑24 (MQ=255)
                          tgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgc                          <  1:6970546‑M1/49‑26 (MQ=255)
                          tgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgc                          <  1:6782974‑M1/49‑26 (MQ=255)
                          tgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgc                          <  1:3441224‑M1/49‑26 (MQ=255)
                            cGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatc                       >  1:1408514‑M1/1‑22 (MQ=255)
                            cGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcat                        <  1:403447‑M1/49‑28 (MQ=255)
                            cGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcat                        <  1:4220441‑M1/49‑28 (MQ=255)
                            cGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcat                        <  1:428155‑M1/49‑28 (MQ=255)
                            cGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcat                        <  1:336544‑M1/49‑28 (MQ=255)
                             gAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatc                       >  1:2253579‑M1/1‑21 (MQ=255)
                             gAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatc                       >  1:3526517‑M1/1‑21 (MQ=255)
                              aGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatc                       >  1:537725‑M1/1‑20 (MQ=255)
                               ggAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatccag                    <  1:6566394‑M1/50‑32 (MQ=255)
                                 aaaaaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagct                  <  1:2970218‑M1/50‑34 (MQ=255)
                                 aaaaaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagct                  <  1:4468139‑M1/50‑34 (MQ=255)
                                 aaaaaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagct                  >  1:4079664‑M1/1‑17 (MQ=255)
                                 aaaaaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagct                  <  1:4733012‑M1/50‑34 (MQ=255)
                                  aaaaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagctg                 >  1:4512829‑M1/1‑16 (MQ=255)
                                  aaaaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagctg                 >  1:5217822‑M1/1‑16 (MQ=255)
                                  aaaaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagctg                 >  1:3263135‑M1/1‑16 (MQ=255)
                                  aaaaTATCAGTATGGCtcaatgatccaaataaagataaagcatccagctg                 <  1:5950940‑M1/50‑35 (MQ=255)
                                   aaaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagctgg                >  1:2279417‑M1/1‑15 (MQ=255)
                                    aaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagctggt               >  1:6913145‑M1/1‑14 (MQ=255)
                                    aaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagctggt               >  1:7367273‑M1/1‑14 (MQ=255)
                                    aaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagctggt               >  1:6719842‑M1/1‑14 (MQ=255)
                                     atatCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtc              <  1:120600‑M1/50‑38 (MQ=255)
                                     atatCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtc              <  1:4132498‑M1/50‑38 (MQ=255)
                                      tatCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtca             <  1:6766181‑M1/50‑39 (MQ=255)
                                      tatCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtca             >  1:7243381‑M1/1‑12 (MQ=255)
                                       atCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtcat            >  1:5740065‑M1/1‑11 (MQ=255)
                                       atCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtcat            >  1:4981962‑M1/1‑11 (MQ=255)
                                        tCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtcatt           >  1:4681275‑M1/1‑10 (MQ=255)
                                        tCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtcatt           >  1:1009024‑M1/1‑10 (MQ=255)
                                        tCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtcatt           <  1:7148761‑M1/50‑41 (MQ=255)
                                        tCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtcat            <  1:786533‑M1/49‑40 (MQ=255)
                                        tCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtcat            <  1:6491276‑M1/49‑40 (MQ=255)
                                        tCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtcat            <  1:2290985‑M1/49‑40 (MQ=255)
                                        tCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtcat            <  1:5605233‑M1/49‑40 (MQ=255)
                                        tCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtcat            <  1:6030149‑M1/49‑40 (MQ=255)
                                         cAGTATGGCtcaatgatccaaataaagataatgcatccagctggtcattc          >  1:3738927‑M1/1‑9 (MQ=255)
                                         cAGTATGGCtcaatgatccaaataaagataatgcatccagctggtcattc          >  1:6057468‑M1/1‑9 (MQ=255)
                                          aGTATGGCtcaatgatccaaataaagataatgcatccagctggtcattct         <  1:6315801‑M1/50‑43 (MQ=255)
                                          aGTATGGCtcaatgatccaaataaagataatgcatccagctggtcattct         >  1:2689551‑M1/1‑8 (MQ=255)
                                          aGTATGGCtcaatgatccaaataaagataatgcatccagctggtcattc          <  1:2956824‑M1/49‑42 (MQ=255)
                                           gTATGGCtcaatgatccaaataaagataatgcatccagctggtcattctt        <  1:7110868‑M1/50‑44 (MQ=255)
                                           gTATGGCtcaatgatccaaataaagataatgcatccagctggtcattct         >  1:2058070‑M1/1‑7 (MQ=255)
                                            tATGGCtcaatgatccaaataaagataatgcatccagctggtcattctta       <  1:2593703‑M1/50‑45 (MQ=255)
                                            tATGGCtcaatgatccaaataaagataatgcatccagctggtcattctta       <  1:2525838‑M1/50‑45 (MQ=255)
                                            tATGGCtcaatgatccaaataaagataatgcatccagctggtcattctta       <  1:7606463‑M1/50‑45 (MQ=255)
                                            tATGGCtcaatgatccaaataaagataatgcatccagctggtcattct         >  1:7425783‑M1/1‑6 (MQ=255)
                                              tGGCtcaatgatccaaataaagataatgcatccagctggtcattcttagc     <  1:2689008‑M1/50‑47 (MQ=255)
                                              tGGCtcaatgatccaaataaagataatgcatccagctggtcattcttag      <  1:3849757‑M1/49‑46 (MQ=255)
                                              tGGCtcaatgatccaaataaagataatgcatccagctggtcattcttag      <  1:5650491‑M1/49‑46 (MQ=255)
                                                 ctcaatgatccaaataaagataatgcatccagctggtcattcttagcagc  <  1:3153504‑M1/50‑50 (MQ=255)
                                                 ctcaatgatccaaataaagataatgcatccagctggtcattcttagcagc  <  1:7376811‑M1/50‑50 (MQ=255)
                                                 ctcaatgatccaaataaagataatgcatccagctggtcattcttagcag   >  1:3709892‑M1/1‑1 (MQ=255)
                                                 ctcaatgatccaaataaagataatgcatccagctggtcattcttagcag   >  1:7290631‑M1/1‑1 (MQ=255)
                                                 |                                                 
CGACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGATGGTTTGTTACTTCCAA  >  NC_000913/3761943‑3762041

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