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 96 [0] [0] 98 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)94 (1.320) 58/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

GACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGATGGTTTGTTACTTCCAA  >  NC_000913/3761944‑3762041
                                                |                                                 
gacgaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCt                                                  >  1:5476755‑M1/1‑49 (MQ=255)
   gaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaa                                               <  1:1662026‑M1/50‑5 (MQ=255)
   gaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaa                                               >  1:3250882‑M1/1‑46 (MQ=255)
    aGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaat                                              <  1:4120087‑M1/50‑6 (MQ=255)
     gTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatg                                             <  1:2709018‑M1/50‑7 (MQ=255)
     gTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatg                                             <  1:2328672‑M1/50‑7 (MQ=255)
     gTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatg                                             <  1:2455561‑M1/50‑7 (MQ=255)
     gTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatg                                             <  1:26403‑M1/50‑7 (MQ=255)
      tttCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatga                                            <  1:5285418‑M1/50‑8 (MQ=255)
      tttCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatga                                            >  1:2964926‑M1/1‑43 (MQ=255)
       ttCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgat                                           >  1:5380723‑M1/1‑42 (MQ=255)
         cATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatc                                          <  1:1278967‑M1/49‑10 (MQ=255)
         cATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatc                                          <  1:2457317‑M1/49‑10 (MQ=255)
          aTGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatcca                                        <  1:5300866‑M1/50‑12 (MQ=255)
          aTGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatcc                                         >  1:2264775‑M1/1‑39 (MQ=255)
           tGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaa                                       <  1:1680043‑M1/50‑13 (MQ=255)
             ccGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaa                                      <  1:134610‑M1/49‑14 (MQ=255)
              cGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaata                                    <  1:5352077‑M1/50‑16 (MQ=255)
              cGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaata                                    <  1:85100‑M1/50‑16 (MQ=255)
               gAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataa                                   <  1:3255327‑M1/50‑17 (MQ=255)
               gAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataa                                   <  1:411863‑M1/50‑17 (MQ=255)
               gAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataa                                   <  1:3106866‑M1/50‑17 (MQ=255)
               gAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaata                                    >  1:5138547‑M1/1‑34 (MQ=255)
                aGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaata                                    >  1:3946921‑M1/1‑33 (MQ=255)
                 gTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaag                                 <  1:3369940‑M1/50‑19 (MQ=255)
                  tCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaaga                                <  1:2968490‑M1/50‑20 (MQ=255)
                   ccTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaaga                                >  1:4320129‑M1/1‑30 (MQ=255)
                   ccTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaaga                                >  1:1419333‑M1/1‑30 (MQ=255)
                   ccTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaaga                                >  1:5015326‑M1/1‑30 (MQ=255)
                    cTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagata                              <  1:324605‑M1/50‑22 (MQ=255)
                    cTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagata                              <  1:3489618‑M1/50‑22 (MQ=255)
                    cTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagat                               <  1:40693‑M1/49‑21 (MQ=255)
                    cTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagat                               <  1:1684021‑M1/49‑21 (MQ=255)
                     tttGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataa                             >  1:5262626‑M1/1‑28 (MQ=255)
                      ttGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataat                            >  1:3892378‑M1/1‑27 (MQ=255)
                      ttGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataat                            >  1:2174927‑M1/1‑27 (MQ=255)
                       tgtgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataat                            <  1:2811152‑M1/49‑24 (MQ=255)
                       tgtgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataat                            <  1:3711570‑M1/49‑24 (MQ=255)
                       tgtgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataa                             <  1:93256‑M1/48‑23 (MQ=255)
                        gtgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatg                           <  1:3580348‑M1/49‑25 (MQ=255)
                         tgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgc                          <  1:908201‑M1/49‑26 (MQ=255)
                           cGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcat                        <  1:4280339‑M1/49‑28 (MQ=255)
                             aGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatcca                     <  1:1536495‑M1/50‑31 (MQ=255)
                             aGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatcca                     <  1:1555015‑M1/50‑31 (MQ=255)
                             aGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatc                       >  1:4947949‑M1/1‑20 (MQ=255)
                             aGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatc                       >  1:3993409‑M1/1‑20 (MQ=255)
                             aGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatc                       >  1:4710224‑M1/1‑20 (MQ=255)
                              ggAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatccag                    <  1:4329487‑M1/50‑32 (MQ=255)
                              ggAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatccag                    <  1:1343463‑M1/50‑32 (MQ=255)
                               gAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatccagc                   <  1:5144846‑M1/50‑33 (MQ=255)
                                aaaaaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagct                  <  1:3443853‑M1/50‑34 (MQ=255)
                                aaaaaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagct                  >  1:2349027‑M1/1‑17 (MQ=255)
                                aaaaaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagct                  <  1:605536‑M1/50‑34 (MQ=255)
                                aaaaaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagct                  <  1:3000138‑M1/50‑34 (MQ=255)
                                aaaaaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagct                  >  1:242870‑M1/1‑17 (MQ=255)
                                 aaaaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagctg                 >  1:4386258‑M1/1‑16 (MQ=255)
                                 aaaaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagctg                 <  1:4453817‑M1/50‑35 (MQ=255)
                                 aaaaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagctg                 >  1:948878‑M1/1‑16 (MQ=255)
                                 aaaaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagctg                 <  1:2621854‑M1/50‑35 (MQ=255)
                                  aaaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagctgg                <  1:2229185‑M1/50‑36 (MQ=255)
                                   aaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagctggt               >  1:547179‑M1/1‑14 (MQ=255)
                                   aaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagctggt               >  1:3689988‑M1/1‑14 (MQ=255)
                                   aaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagctggt               >  1:2184329‑M1/1‑14 (MQ=255)
                                    atatCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtc              <  1:5201367‑M1/50‑38 (MQ=255)
                                    atatCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtc              <  1:330755‑M1/50‑38 (MQ=255)
                                    atatCAGTATGGCtcaatgatccaaataaagataatgcatccagctggt               >  1:3389243‑M1/1‑13 (MQ=255)
                                     tatCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtca             >  1:2899281‑M1/1‑12 (MQ=255)
                                     tatCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtca             <  1:27785‑M1/50‑39 (MQ=255)
                                     tatCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtca             <  1:1164804‑M1/50‑39 (MQ=255)
                                     tatCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtca             <  1:4071267‑M1/50‑39 (MQ=255)
                                     tatCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtca             <  1:4370670‑M1/50‑39 (MQ=255)
                                       tCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtcatt           <  1:4496635‑M1/50‑41 (MQ=255)
                                       tCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtcatt           <  1:828121‑M1/50‑41 (MQ=255)
                                       tCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtcat            <  1:659807‑M1/49‑40 (MQ=255)
                                        cAGTATGGCtcaatgatccaaataaagataatgcatccagctggtcattc          >  1:700658‑M1/1‑9 (MQ=255)
                                        cAGTATGGCtcaatgatccaaataaagataatgcatccagctggtcattc          >  1:3096432‑M1/1‑9 (MQ=255)
                                         aGTATGGCtcaatgatccaaataaagataatgcatccagctggtcattct         >  1:4294824‑M1/1‑8 (MQ=255)
                                         aGTATGGCtcaatgatccaaataaagataatgcatccagctggtcattct         <  1:4307573‑M1/50‑43 (MQ=255)
                                         aGTATGGCtcaatgatccaaataaagataatgcatccagctggtcattct         <  1:402871‑M1/50‑43 (MQ=255)
                                         aGTATGGCtcaatgatccaaataaagataatgcatccagctggtcattc          <  1:5252125‑M1/49‑42 (MQ=255)
                                         aGTATGGCtcaatgatccaaataaagataatgcatccagctggtcattc          <  1:2282967‑M1/49‑42 (MQ=255)
                                          gTATGGCtcaatgatccaaataaagataatgcatccagctggtcattctt        <  1:3492382‑M1/50‑44 (MQ=255)
                                          gTATGGCtcaatgatccaaataaagataatgcatccagctggtcattctt        <  1:1999424‑M1/50‑44 (MQ=255)
                                           tATGGCtcaatgatccaaataaagataatgcatccagctggtcattctta       <  1:2471765‑M1/50‑45 (MQ=255)
                                           tATGGCtcaatgatccaaataaagataatgcatccagctggtcattctta       <  1:933326‑M1/50‑45 (MQ=255)
                                           tATGGCtcaatgatccaaataaagataatgcatccagctggtcattct         >  1:1106302‑M1/1‑6 (MQ=255)
                                            aTGGCtcaatgatccaaataaagataatgcatccagctggtcattcttag      >  1:2304911‑M1/1‑5 (MQ=255)
                                            aTGGCtcaatgatccaaataaagataatgcatccagctggtcattcttag      >  1:1170249‑M1/1‑5 (MQ=255)
                                             tGGCtcaatgatccaaataaagataatgcatccagctggtcattcttagc     <  1:3138323‑M1/50‑47 (MQ=255)
                                             tGGCtcaatgatccaaataaagataatgcatccagctggtcattcttag      <  1:1236708‑M1/49‑46 (MQ=255)
                                             tGGCtcaatgatccaaataaagataatgcatccagctggtcattcttag      <  1:5165094‑M1/49‑46 (MQ=255)
                                             tGGCtcaatgatccaaataaagataatgcatccagctggtcattcttag      >  1:115538‑M1/1‑4 (MQ=255)
                                               gCtcaatgatccaaataaagataatgcatccagctggtcattcttagcag   <  1:3660562‑M1/50‑49 (MQ=255)
                                               gCtcaatgatccaaataaagataatgcatccagctggtcattcttagca    <  1:4127625‑M1/49‑48 (MQ=255)
                                                ctcaatgatccaaataaagataatgcatccagctggtcattcttagcagc  <  1:208407‑M1/50‑50 (MQ=255)
                                                ctcaatgatccaaataaagataatgcatccagctggtcattcttagcag   >  1:2874393‑M1/1‑1 (MQ=255)
                                                |                                                 
GACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGATGGTTTGTTACTTCCAA  >  NC_000913/3761944‑3762041

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: