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 133 [0] [0] 131 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)129 (1.490) 74/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 133 total aligned reads displayed.

CGACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGATGGTTTGTTACTTCCAA  >  NC_000913/3761943‑3762041
                                                 |                                                 
cgacgaGTTTCATGCCGAGTCCTTTGTGCGATGAAAAATATCAGTATGGC                                                   >  1:2838897‑M1/1‑50 (MQ=255)
cgacgaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGC                                                   >  1:4172647‑M1/1‑50 (MQ=255)
cgacgaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGC                                                   <  1:1094868‑M1/50‑1 (MQ=255)
 gacgaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCt                                                  >  1:80115‑M1/1‑49 (MQ=255)
 gacgaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCt                                                  <  1:3794707‑M1/50‑2 (MQ=255)
   cgaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtca                                                <  1:252609‑M1/50‑4 (MQ=255)
   cgaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtca                                                <  1:657008‑M1/50‑4 (MQ=255)
    gaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaa                                               <  1:3961173‑M1/50‑5 (MQ=255)
     aGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaat                                              >  1:6345408‑M1/1‑45 (MQ=255)
     aGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaat                                              >  1:2062711‑M1/1‑45 (MQ=255)
      gTTTCATACCGAGTCCTTTGTCCGATGATAAATATCAGTATGGCtcaat                                              <  1:1154482‑M1/49‑6 (MQ=255)
       tttCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatga                                            >  1:123231‑M1/1‑43 (MQ=255)
         tCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatc                                          >  1:1353551‑M1/1‑41 (MQ=255)
         tCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgat                                           <  1:3274147‑M1/49‑9 (MQ=255)
         tCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgat                                           <  1:2937149‑M1/49‑9 (MQ=255)
          cATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatc                                          <  1:5931733‑M1/49‑10 (MQ=255)
           aTGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatcca                                        <  1:2340534‑M1/50‑12 (MQ=255)
           aTGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatcca                                        <  1:640190‑M1/50‑12 (MQ=255)
           aTGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatcc                                         <  1:4624594‑M1/49‑11 (MQ=255)
            tGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaa                                       <  1:2976401‑M1/50‑13 (MQ=255)
            tGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaa                                       <  1:3047677‑M1/50‑13 (MQ=255)
             gCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaa                                      >  1:2122957‑M1/1‑37 (MQ=255)
              ccGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaat                                     <  1:4035232‑M1/50‑15 (MQ=255)
              ccGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaa                                      <  1:5795408‑M1/49‑14 (MQ=255)
              ccGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaa                                      <  1:4429806‑M1/49‑14 (MQ=255)
              ccGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaa                                      >  1:4665805‑M1/1‑36 (MQ=255)
               cGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaata                                    <  1:5519135‑M1/50‑16 (MQ=255)
               cGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaat                                     >  1:6246046‑M1/1‑35 (MQ=255)
                gAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaata                                    >  1:4372432‑M1/1‑34 (MQ=255)
                gAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaata                                    >  1:5220763‑M1/1‑34 (MQ=255)
                 aGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaa                                  <  1:377993‑M1/50‑18 (MQ=255)
                  gTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaag                                 <  1:1770590‑M1/50‑19 (MQ=255)
                    ccTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaaga                                >  1:2777887‑M1/1‑30 (MQ=255)
                     cTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagata                              <  1:3529891‑M1/50‑22 (MQ=255)
                     cTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagat                               >  1:3560025‑M1/1‑29 (MQ=255)
                     cTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagat                               <  1:5614622‑M1/49‑21 (MQ=255)
                      tttGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataa                             >  1:2128358‑M1/1‑28 (MQ=255)
                       ttGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataat                            >  1:2730336‑M1/1‑27 (MQ=255)
                        tgtgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatg                           >  1:1326364‑M1/1‑26 (MQ=255)
                        tgtgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataat                            <  1:301234‑M1/49‑24 (MQ=255)
                        tgtgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataat                            <  1:5375380‑M1/49‑24 (MQ=255)
                        tgtgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataa                             <  1:1681698‑M1/48‑23 (MQ=255)
                         gtgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatg                           <  1:2769707‑M1/49‑25 (MQ=255)
                         gtgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatg                           <  1:3805973‑M1/49‑25 (MQ=255)
                         gtgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatg                           <  1:250114‑M1/49‑25 (MQ=255)
                         gtgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatg                           <  1:6084363‑M1/49‑25 (MQ=255)
                         gtgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatg                           >  1:1483531‑M1/1‑25 (MQ=255)
                          tgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgca                         >  1:4379544‑M1/1‑24 (MQ=255)
                          tgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgca                         >  1:6621074‑M1/1‑24 (MQ=255)
                            cGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatc                       <  1:1526710‑M1/50‑29 (MQ=255)
                            cGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcat                        <  1:5310864‑M1/49‑28 (MQ=255)
                            cGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcat                        <  1:1104481‑M1/49‑28 (MQ=255)
                             gAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatcc                      <  1:3274943‑M1/50‑30 (MQ=255)
                             gAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatcc                      <  1:185677‑M1/50‑30 (MQ=255)
                             gAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatcc                      <  1:4044667‑M1/50‑30 (MQ=255)
                             gAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatc                       >  1:1707108‑M1/1‑21 (MQ=255)
                              aGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatcca                     <  1:2641075‑M1/50‑31 (MQ=255)
                              aGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatc                       >  1:922183‑M1/1‑20 (MQ=255)
                              aGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatc                       >  1:64635‑M1/1‑20 (MQ=255)
                               ggAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatccag                    >  1:3554523‑M1/1‑19 (MQ=255)
                               ggAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatccag                    <  1:5505086‑M1/50‑32 (MQ=255)
                               ggAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatccag                    <  1:2137637‑M1/50‑32 (MQ=255)
                               ggAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatccag                    <  1:1927399‑M1/50‑32 (MQ=255)
                               ggAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatccag                    >  1:3979506‑M1/1‑19 (MQ=255)
                                gAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatccagc                   >  1:4875943‑M1/1‑18 (MQ=255)
                                gAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatccagc                   >  1:2864263‑M1/1‑18 (MQ=255)
                                 aaaaaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagct                  <  1:3731195‑M1/50‑34 (MQ=255)
                                  aaaaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagctg                 <  1:735364‑M1/50‑35 (MQ=255)
                                  aaaaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagctg                 <  1:4524709‑M1/50‑35 (MQ=255)
                                  aaaaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagctg                 <  1:2142536‑M1/50‑35 (MQ=255)
                                  aaaaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagctg                 <  1:2296146‑M1/50‑35 (MQ=255)
                                  aaaaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagctg                 >  1:2345129‑M1/1‑16 (MQ=255)
                                    aaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagctggt               <  1:890207‑M1/50‑37 (MQ=255)
                                    aaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagctggt               <  1:31830‑M1/50‑37 (MQ=255)
                                     atatCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtc              <  1:6019706‑M1/50‑38 (MQ=255)
                                     atatCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtc              <  1:5184397‑M1/50‑38 (MQ=255)
                                     atatCAGTATGGCtcaatgatccaaataaagataatgcatccagctggt               >  1:1357325‑M1/1‑13 (MQ=255)
                                      tatCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtca             >  1:3393269‑M1/1‑12 (MQ=255)
                                      tatCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtca             >  1:4183943‑M1/1‑12 (MQ=255)
                                      tatCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtca             <  1:3033805‑M1/50‑39 (MQ=255)
                                       atCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtcat            >  1:3486383‑M1/1‑11 (MQ=255)
                                        tCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtcatt           <  1:5360427‑M1/50‑41 (MQ=255)
                                        tCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtcat            <  1:3567200‑M1/49‑40 (MQ=255)
                                          aGTATGGCtcaatgatccaaataaagataatgcatccagctggtcattct         <  1:1923498‑M1/50‑43 (MQ=255)
                                          aGTATGGCtcaatgatccaaataaagataatgcatccagctggtcattct         >  1:6248123‑M1/1‑8 (MQ=255)
                                          aGTATGGCtcaatgatccaaataaagataatgcatccagctggtcattct         >  1:519208‑M1/1‑8 (MQ=255)
                                          aGTATGGCtcaatgatccaaataaagataatgcatccagctggtcattc          <  1:6188207‑M1/49‑42 (MQ=255)
                                          aGTATGGCtcaatgatccaaataaagataatgcatccagctggtcattc          <  1:312956‑M1/49‑42 (MQ=255)
                                           gTATGGCtcaatgatccaaataaagataatgcatccagctggtcattctt        <  1:4907101‑M1/50‑44 (MQ=255)
                                           gTATGGCtcaatgatccaaataaagataatgcatccagctggtcattctt        <  1:4467628‑M1/50‑44 (MQ=255)
                                           gTATGGCtcaatgatccaaataaagataatgcatccagctggtcattct         >  1:4836843‑M1/1‑7 (MQ=255)
                                            tATGGCtcaatgatccaaataaagataatgcatccagctggtcattct         >  1:3499135‑M1/1‑6 (MQ=255)
                                            tATGGCtcaatgatccaaataaagataatgcatccagctggtcattct         >  1:3713214‑M1/1‑6 (MQ=255)
                                            tATGGCtcaatgatccaaataaagacaatgcatccagctggtcattctta       <  1:6420558‑M1/50‑45 (MQ=255)
                                              tGGCtcaatgatccaaataaagataatgcatccagctggtcattcttag      >  1:5844960‑M1/1‑4 (MQ=255)
                                               ggCtcaatgatccaaataaagataatgcatccagctggtcattcttag      >  1:3976914‑M1/1‑3 (MQ=255)
                                               ggCtcaatgatccaaataaagataatgcatccagctggtcattcttag      >  1:3694025‑M1/1‑3 (MQ=255)
                                                gCtcaatgatccaaataaagataatgcatccagctggtcattcttagcag   >  1:2137213‑M1/1‑2 (MQ=255)
                                                gCtcaatgatccaaataaagataatgcatccagctggtcattcttagcag   >  1:4744429‑M1/1‑2 (MQ=255)
                                                 ctcaatgatccaaataaagataatgcatccagctggtcattcttagcagc  <  1:1710484‑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 < 30 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: