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 121 [0] [0] 119 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)116 (1.490) 67/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 121 total aligned reads displayed.

CGACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGATGGTTTGTTACTTCCA  >  NC_000913/3761943‑3762040
                                                 |                                                
cgacgaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGC                                                  <  1:3457679‑M1/50‑1 (MQ=255)
cgacgaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGC                                                  <  1:6002897‑M1/50‑1 (MQ=255)
 gacgaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCt                                                 <  1:3146554‑M1/50‑2 (MQ=255)
 gacgaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCt                                                 <  1:5370395‑M1/50‑2 (MQ=255)
 gacgaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCt                                                 <  1:5012352‑M1/50‑2 (MQ=255)
  acgaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCt                                                 >  1:4985468‑M1/1‑48 (MQ=255)
   cgaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtc                                                >  1:3090891‑M1/1‑47 (MQ=255)
    gaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaa                                              >  1:4447116‑M1/1‑46 (MQ=255)
     aGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaat                                             <  1:2703805‑M1/50‑6 (MQ=255)
      gTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatg                                            <  1:5306888‑M1/50‑7 (MQ=255)
       tttCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatga                                           <  1:4461958‑M1/50‑8 (MQ=255)
        ttCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgat                                          >  1:210389‑M1/1‑42 (MQ=255)
         tCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatc                                         >  1:5030762‑M1/1‑41 (MQ=255)
         tCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgat                                          <  1:862967‑M1/49‑9 (MQ=255)
         tCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgat                                          <  1:4279135‑M1/49‑9 (MQ=255)
         tCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgat                                          <  1:813221‑M1/49‑9 (MQ=255)
         tCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgat                                          <  1:540927‑M1/49‑9 (MQ=255)
           aTGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatcca                                       <  1:1625037‑M1/50‑12 (MQ=255)
           aTGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatcca                                       <  1:3657773‑M1/50‑12 (MQ=255)
           aTGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatcc                                        >  1:4747259‑M1/1‑39 (MQ=255)
              ccGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaat                                    <  1:340195‑M1/50‑15 (MQ=255)
              ccGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaat                                    <  1:5744146‑M1/50‑15 (MQ=255)
              ccGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaa                                     <  1:5269809‑M1/49‑14 (MQ=255)
              ccGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaa                                     <  1:5945031‑M1/49‑14 (MQ=255)
              ccGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaa                                     <  1:4926841‑M1/49‑14 (MQ=255)
                gAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataa                                  <  1:3460123‑M1/50‑17 (MQ=255)
                 aGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaa                                 <  1:1201450‑M1/50‑18 (MQ=255)
                 aGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaa                                 <  1:784181‑M1/50‑18 (MQ=255)
                 aGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaa                                 <  1:4498532‑M1/50‑18 (MQ=255)
                 aGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaata                                   >  1:2643506‑M1/1‑33 (MQ=255)
                 aGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaata                                   >  1:4264788‑M1/1‑33 (MQ=255)
                  gTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaag                                <  1:103525‑M1/50‑19 (MQ=255)
                   tCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaaga                               <  1:3247093‑M1/50‑20 (MQ=255)
                   tCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaaga                               <  1:613598‑M1/50‑20 (MQ=255)
                   tCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaata                                   >  1:4363845‑M1/1‑31 (MQ=255)
                    ccTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaaga                               >  1:4079228‑M1/1‑30 (MQ=255)
                     cTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagata                             <  1:1351888‑M1/50‑22 (MQ=255)
                     cTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagata                             <  1:3241521‑M1/50‑22 (MQ=255)
                     cTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagata                             <  1:1353279‑M1/50‑22 (MQ=255)
                     cTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagat                              >  1:5054220‑M1/1‑29 (MQ=255)
                     cTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagat                              <  1:1495683‑M1/49‑21 (MQ=255)
                        tgtgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatg                          >  1:930196‑M1/1‑26 (MQ=255)
                        tgtgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatg                          >  1:642090‑M1/1‑26 (MQ=255)
                        tgtgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataat                           <  1:3605055‑M1/49‑24 (MQ=255)
                        tgtgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataat                           <  1:4884119‑M1/49‑24 (MQ=255)
                        tgtgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataa                            <  1:3350145‑M1/48‑23 (MQ=255)
                         gtgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatg                          <  1:2919794‑M1/49‑25 (MQ=255)
                          tgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgca                        <  1:5219202‑M1/50‑27 (MQ=255)
                          tgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgca                        <  1:1226834‑M1/50‑27 (MQ=255)
                          tgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgca                        >  1:557404‑M1/1‑24 (MQ=255)
                          tgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgc                         <  1:1002423‑M1/49‑26 (MQ=255)
                          tgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgc                         <  1:2275106‑M1/49‑26 (MQ=255)
                          tgCGAGGAAAAATATCAGTATGGCtcaatgatcaaaataaagataatgc                         <  1:840129‑M1/49‑26 (MQ=255)
                            cGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatc                      <  1:2239701‑M1/50‑29 (MQ=255)
                            cGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcat                       <  1:1975115‑M1/49‑28 (MQ=255)
                             gAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatcc                     <  1:1149806‑M1/50‑30 (MQ=255)
                             gAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatcc                     <  1:3361131‑M1/50‑30 (MQ=255)
                             gAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatc                      >  1:5699796‑M1/1‑21 (MQ=255)
                             gAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatc                      >  1:1175000‑M1/1‑21 (MQ=255)
                              aGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatcca                    <  1:5613195‑M1/50‑31 (MQ=255)
                              aGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatcca                    <  1:5624366‑M1/50‑31 (MQ=255)
                              aGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatcca                    <  1:5961778‑M1/50‑31 (MQ=255)
                              aGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatcca                    <  1:952018‑M1/50‑31 (MQ=255)
                               ggAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatccag                   <  1:4831300‑M1/50‑32 (MQ=255)
                               ggAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatccag                   <  1:5976335‑M1/50‑32 (MQ=255)
                                gAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatccagc                  >  1:5092359‑M1/1‑18 (MQ=255)
                                gAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatccagc                  <  1:3899407‑M1/50‑33 (MQ=255)
                                gAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatccagc                  <  1:1269936‑M1/50‑33 (MQ=255)
                                 aaaaaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagct                 >  1:2715060‑M1/1‑17 (MQ=255)
                                  aaaaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagctg                >  1:2216275‑M1/1‑16 (MQ=255)
                                  aaaaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagctg                <  1:4667615‑M1/50‑35 (MQ=255)
                                   aaaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagctgg               >  1:5015869‑M1/1‑15 (MQ=255)
                                    aaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagctggt              >  1:2865241‑M1/1‑14 (MQ=255)
                                    aaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagctggt              <  1:4793418‑M1/50‑37 (MQ=255)
                                     atatCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtc             <  1:1497916‑M1/50‑38 (MQ=255)
                                     atatCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtc             <  1:673221‑M1/50‑38 (MQ=255)
                                     atatCAGTATGGCtcaatgatccaaataaagataatgcatccagctggt              >  1:768641‑M1/1‑13 (MQ=255)
                                       atCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtcat           >  1:191758‑M1/1‑11 (MQ=255)
                                        tCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtcatt          >  1:4549191‑M1/1‑10 (MQ=255)
                                        tCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtcatt          <  1:928323‑M1/50‑41 (MQ=255)
                                        tCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtcat           <  1:3468194‑M1/49‑40 (MQ=255)
                                          aGTATGGCtcaatgatccaaataaagataatgcatccagctggtcattct        >  1:2592987‑M1/1‑8 (MQ=255)
                                          aGTATGGCtcaatgatccaaataaagataatgcatccagctggtcattc         <  1:1169768‑M1/49‑42 (MQ=255)
                                          aGTATGGCtcaatgatccaaataaagataatgcatccagctggtcattc         <  1:2798222‑M1/49‑42 (MQ=255)
                                           gTATGGCtcaatgatccaaataaagataatgcatccagctggtcattctt       <  1:4306492‑M1/50‑44 (MQ=255)
                                           gTATGGCtcaatgatccaaataaagataatgcatccagctggtcattctt       <  1:3958968‑M1/50‑44 (MQ=255)
                                           gTATGGCtcaatgatccaaataaagataatgcatccagctggtcattctt       <  1:87473‑M1/50‑44 (MQ=255)
                                           gTATGGCtcaatgatccaaataaagataatgcatccagctggtcattct        >  1:2535258‑M1/1‑7 (MQ=255)
                                            tATGGCtcaatgatccaaataaagataatgcatccagctggtcattctta      <  1:5567229‑M1/50‑45 (MQ=255)
                                            tATGGCtcaatgatccaaataaagataatgcatccagctggtcattctta      <  1:5253277‑M1/50‑45 (MQ=255)
                                             aTGGCtcaatgatccaaataaagataatgcatccagctggtcattcttag     >  1:901053‑M1/1‑5 (MQ=255)
                                              tGGCtcaatgatccaaataaagataatgcatccagctggtcattcttag     >  1:588562‑M1/1‑4 (MQ=255)
                                               ggCtcaatgatccaaataaagataatgcatccagctggtcattcttag     >  1:2161368‑M1/1‑3 (MQ=255)
                                                gCtcaatgatccaaataaagataatgcatccagctggtcattcttagcag  <  1:3788623‑M1/50‑49 (MQ=255)
                                                gCtcaatgatccaaataaagataatgcatccagctggtcattcttagcag  >  1:5452859‑M1/1‑2 (MQ=255)
                                                gCtcaatgatccaaataaagataatgcatccagctggtcattcttagcag  >  1:2344290‑M1/1‑2 (MQ=255)
                                                gCtcaatgatccaaataaagataatgcatccagctggtcattcttagca   <  1:612445‑M1/49‑48 (MQ=255)
                                                gCtcaatgatccaaataaagataatgcatccagctggtcattcttagca   <  1:4515531‑M1/49‑48 (MQ=255)
                                                 ctcaatgatccaaataaagataatgcatccagctggtcattcttagcag  >  1:3581720‑M1/1‑1 (MQ=255)
                                                 ctcaatgatccaaataaagataatgcatccagctggtcattcttagcag  >  1:5469235‑M1/1‑1 (MQ=255)
                                                 |                                                
CGACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGATGGTTTGTTACTTCCA  >  NC_000913/3761943‑3762040

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: 
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.