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 131 [0] [0] 129 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.190) 75/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 131 total aligned reads displayed.

CGACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGATGGTTTGTTACTTCCA  >  NC_000913/3761943‑3762040
                                                 |                                                
cgacgaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGC                                                  <  1:3478690‑M1/50‑1 (MQ=255)
cgacgaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGC                                                  >  1:1531594‑M1/1‑50 (MQ=255)
 gacgaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCt                                                 >  1:3350118‑M1/1‑49 (MQ=255)
 gacgaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCt                                                 >  1:8274092‑M1/1‑49 (MQ=255)
 gacgaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCt                                                 <  1:3195338‑M1/50‑2 (MQ=255)
 gacgaGTTTCATGCCGAGTCCTATGTGCGAGGAAAAATATCAGTATGGCt                                                 <  1:5408466‑M1/50‑2 (MQ=255)
  acgaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCt                                                 >  1:4425952‑M1/1‑48 (MQ=255)
   cgaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtca                                               <  1:4164096‑M1/50‑4 (MQ=255)
    gaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaa                                              >  1:913214‑M1/1‑46 (MQ=255)
    gaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaa                                              >  1:7416797‑M1/1‑46 (MQ=255)
     aGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCttaat                                             <  1:1520202‑M1/50‑6 (MQ=255)
      gTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaat                                             >  1:5754773‑M1/1‑44 (MQ=255)
       tttCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatga                                           >  1:2747581‑M1/1‑43 (MQ=255)
        ttCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgat                                          >  1:6102757‑M1/1‑42 (MQ=255)
         tCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatc                                         >  1:1182045‑M1/1‑41 (MQ=255)
          cATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatc                                         >  1:2999728‑M1/1‑40 (MQ=255)
          cATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatc                                         <  1:3806554‑M1/49‑10 (MQ=255)
          cATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatc                                         <  1:8107515‑M1/49‑10 (MQ=255)
           aTGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatcca                                       <  1:2830649‑M1/50‑12 (MQ=255)
           aTGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatcca                                       <  1:2206114‑M1/50‑12 (MQ=255)
            tGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaa                                      <  1:2595421‑M1/50‑13 (MQ=255)
             gCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaa                                     >  1:3308099‑M1/1‑37 (MQ=255)
              ccGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaa                                     <  1:2517257‑M1/49‑14 (MQ=255)
              ccGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaa                                     <  1:1061901‑M1/49‑14 (MQ=255)
               cGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaata                                   <  1:3571832‑M1/50‑16 (MQ=255)
               cGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaata                                   <  1:5929819‑M1/50‑16 (MQ=255)
               cGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaat                                    >  1:2317649‑M1/1‑35 (MQ=255)
                gAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataa                                  <  1:7744882‑M1/50‑17 (MQ=255)
                gAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaata                                   >  1:3536750‑M1/1‑34 (MQ=255)
                 aGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaa                                 <  1:6922896‑M1/50‑18 (MQ=255)
                 aGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaa                                 <  1:756277‑M1/50‑18 (MQ=255)
                 aGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaata                                   >  1:4106645‑M1/1‑33 (MQ=255)
                  gTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaag                                <  1:422213‑M1/50‑19 (MQ=255)
                  gTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaata                                   >  1:180658‑M1/1‑32 (MQ=255)
                   tCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaaga                               <  1:2249521‑M1/50‑20 (MQ=255)
                    ccTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaaga                               >  1:6983705‑M1/1‑30 (MQ=255)
                     cTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagata                             <  1:5027221‑M1/50‑22 (MQ=255)
                     cTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagata                             <  1:1578056‑M1/50‑22 (MQ=255)
                     cTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagat                              <  1:5236434‑M1/49‑21 (MQ=255)
                     cTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagat                              >  1:7021288‑M1/1‑29 (MQ=255)
                     cTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagat                              >  1:3193972‑M1/1‑29 (MQ=255)
                      tttGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataa                            >  1:191794‑M1/1‑28 (MQ=255)
                       ttGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataat                           >  1:4244258‑M1/1‑27 (MQ=255)
                       ttGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataat                           >  1:7187586‑M1/1‑27 (MQ=255)
                        tgtgCGAGGAAAAATATCAGTATGGCtcaattatccaaataaagataat                           <  1:4952047‑M1/49‑24 (MQ=255)
                        tgtgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataa                            <  1:8307865‑M1/48‑23 (MQ=255)
                         gtgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatg                          >  1:478319‑M1/1‑25 (MQ=255)
                          tgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgca                        <  1:2137344‑M1/50‑27 (MQ=255)
                          tgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgca                        >  1:6884826‑M1/1‑24 (MQ=255)
                          tgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgca                        <  1:1998011‑M1/50‑27 (MQ=255)
                          tgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgca                        <  1:7546628‑M1/50‑27 (MQ=255)
                          tgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgc                         <  1:6695776‑M1/49‑26 (MQ=255)
                          tgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgc                         <  1:2219000‑M1/49‑26 (MQ=255)
                           gCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcat                       >  1:7876047‑M1/1‑23 (MQ=255)
                            cGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatc                      >  1:5935566‑M1/1‑22 (MQ=255)
                            cGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcat                       <  1:181070‑M1/49‑28 (MQ=255)
                            cGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcat                       <  1:5397533‑M1/49‑28 (MQ=255)
                             gAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatcc                     <  1:3888426‑M1/50‑30 (MQ=255)
                              aGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatcca                    <  1:3382752‑M1/50‑31 (MQ=255)
                              aGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatcca                    <  1:6955305‑M1/50‑31 (MQ=255)
                              aGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatc                      >  1:641656‑M1/1‑20 (MQ=255)
                              aGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatc                      >  1:7962041‑M1/1‑20 (MQ=255)
                              aGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatc                      >  1:1143603‑M1/1‑20 (MQ=255)
                              aGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatc                      >  1:4607239‑M1/1‑20 (MQ=255)
                               ggAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatccag                   <  1:6810967‑M1/50‑32 (MQ=255)
                                gAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatccagc                  <  1:4659793‑M1/50‑33 (MQ=255)
                                gAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatccagc                  <  1:7212703‑M1/50‑33 (MQ=255)
                                 aaaaaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagct                 >  1:2482223‑M1/1‑17 (MQ=255)
                                 aaaaaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagct                 >  1:5939771‑M1/1‑17 (MQ=255)
                                 aaaaaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagct                 >  1:6084940‑M1/1‑17 (MQ=255)
                                  aaaaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagctg                >  1:3183711‑M1/1‑16 (MQ=255)
                                  aaaaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagctg                <  1:1844164‑M1/50‑35 (MQ=255)
                                  aaaaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagctg                <  1:1690530‑M1/50‑35 (MQ=255)
                                    aaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagctggt              >  1:252696‑M1/1‑14 (MQ=255)
                                    aaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagctggt              <  1:4034633‑M1/50‑37 (MQ=255)
                                     atatCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtc             <  1:8209379‑M1/50‑38 (MQ=255)
                                     atatCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtc             <  1:5740246‑M1/50‑38 (MQ=255)
                                     atatCAGTATGGCtcaatgatccaaataaagataatgcatccagctggt              >  1:3041309‑M1/1‑13 (MQ=255)
                                      tatCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtca            <  1:1334884‑M1/50‑39 (MQ=255)
                                       atCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtcat           >  1:5824168‑M1/1‑11 (MQ=255)
                                        tCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtcatt          <  1:1349629‑M1/50‑41 (MQ=255)
                                        tCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtcatt          <  1:2401589‑M1/50‑41 (MQ=255)
                                        tCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtcatt          <  1:2583564‑M1/50‑41 (MQ=255)
                                         cAGTATGGCtcaatgatccaaataaagataatgcatccagctggtcattc         >  1:2708963‑M1/1‑9 (MQ=255)
                                          aGTATGGCtcaatgatccaaataaagataatgcatccagctggtcattct        <  1:6275413‑M1/50‑43 (MQ=255)
                                          aGTATGGCtcaatgatccaaataaagataatgcatccagctggtcattct        <  1:8058153‑M1/50‑43 (MQ=255)
                                          aGTATGGCtcaatgatccaaataaagataatgcatccagctggtcattc         <  1:3857431‑M1/49‑42 (MQ=255)
                                          aGTATGGCtcaatgatccaaataaagataatgcatccagctggtcattc         <  1:1364038‑M1/49‑42 (MQ=255)
                                           gTATGGCtcaatgatccaaataaagataatgcatccagctggtcattctt       <  1:1775910‑M1/50‑44 (MQ=255)
                                            tATGGCtcaatgatccaaataaagataatgcatccagctggtcattctta      <  1:3949770‑M1/50‑45 (MQ=255)
                                            tATGGCtcaatgatccaaataaagataatgcatccagctggtcattct        >  1:3636207‑M1/1‑6 (MQ=255)
                                              tGGCtcaatgatccaaataaagataatgcatccagctggtcattcttag     >  1:770903‑M1/1‑4 (MQ=255)
                                              tGGCtcaatgatccaaataaagataatgcatccagctggtcattcttag     >  1:1764068‑M1/1‑4 (MQ=255)
                                              tGGCtcaatgatccaaataaagataatgcatccagctggtcattcttag     <  1:2107039‑M1/49‑46 (MQ=255)
                                              tGGCtcaatgatccaaataaagataatgcatccagctggtcattcttag     <  1:7403761‑M1/49‑46 (MQ=255)
                                                gCtcaatgatccaaataaagataatgcatccagctggtcattcttagcag  >  1:4967367‑M1/1‑2 (MQ=255)
                                                gCtcaatgatccaaataaagataatgcatccagctggtcattcttagcag  >  1:6315096‑M1/1‑2 (MQ=255)
                                                gCtcaatgatccaaataaagataatgcatccagctggtcattcttagcag  >  1:2207517‑M1/1‑2 (MQ=255)
                                                gCtcaatgatccaaataaagataatgcatccagctggtcattcttagcag  >  1:6605579‑M1/1‑2 (MQ=255)
                                                gCtcaatgatccaaataaagataatgcatccagctggtcattcttagcag  <  1:1177099‑M1/50‑49 (MQ=255)
                                                 |                                                
CGACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGATGGTTTGTTACTTCCA  >  NC_000913/3761943‑3762040

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