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 149 [0] [0] 147 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)146 (1.380) 76/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 149 total aligned reads displayed.

CGACGAGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCCTGGTGATGGGAAAAAGGGCAGAAAATGTTGATGGTTTGTTACTTCCA  >  NC_000913/3761943‑3762040
                                                 |                                                
cgacgaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGC                                                  >  1:4072167‑M1/1‑50 (MQ=255)
cgacgaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGC                                                  <  1:4501748‑M1/50‑1 (MQ=255)
   cgaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtc                                                >  1:4220306‑M1/1‑47 (MQ=255)
    gaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaa                                              >  1:3718296‑M1/1‑46 (MQ=255)
    gaGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaa                                              <  1:252053‑M1/50‑5 (MQ=255)
     aGTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaat                                             <  1:184490‑M1/50‑6 (MQ=255)
      gTTTCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaat                                             >  1:4981911‑M1/1‑44 (MQ=255)
         tCATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgat                                          <  1:6320545‑M1/49‑9 (MQ=255)
          cATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatc                                         <  1:6148299‑M1/49‑10 (MQ=255)
          cATGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatc                                         <  1:3278855‑M1/49‑10 (MQ=255)
           aTGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatcc                                        <  1:1397519‑M1/49‑11 (MQ=255)
           aTGCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatcc                                        >  1:7506017‑M1/1‑39 (MQ=255)
             gCCGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaa                                     >  1:3397418‑M1/1‑37 (MQ=255)
              ccGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaat                                    <  1:4657572‑M1/50‑15 (MQ=255)
               cGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaata                                   <  1:3069510‑M1/50‑16 (MQ=255)
               cGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaata                                   <  1:1234757‑M1/50‑16 (MQ=255)
               cGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaata                                   <  1:3865869‑M1/50‑16 (MQ=255)
               cGAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaat                                    >  1:278431‑M1/1‑35 (MQ=255)
                gAGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataa                                  <  1:7943849‑M1/50‑17 (MQ=255)
                 aGTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaata                                   >  1:4920163‑M1/1‑33 (MQ=255)
                  gTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaag                                <  1:5592592‑M1/50‑19 (MQ=255)
                  gTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaag                                <  1:2369035‑M1/50‑19 (MQ=255)
                  gTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaag                                <  1:1848893‑M1/50‑19 (MQ=255)
                  gTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaag                                <  1:2719628‑M1/50‑19 (MQ=255)
                  gTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaata                                   >  1:1203109‑M1/1‑32 (MQ=255)
                  gTCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaata                                   >  1:7269093‑M1/1‑32 (MQ=255)
                   tCCTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaata                                   >  1:5937694‑M1/1‑31 (MQ=255)
                    ccTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaaga                               >  1:1646496‑M1/1‑30 (MQ=255)
                    ccTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaaga                               >  1:6822773‑M1/1‑30 (MQ=255)
                    ccTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaaga                               >  1:829955‑M1/1‑30 (MQ=255)
                     cTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagata                             <  1:1477723‑M1/50‑22 (MQ=255)
                     cTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagata                             <  1:3491092‑M1/50‑22 (MQ=255)
                     cTTTGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagat                              <  1:7725447‑M1/49‑21 (MQ=255)
                      tttGTGCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataa                            >  1:708286‑M1/1‑28 (MQ=255)
                        tgtgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatg                          >  1:1999503‑M1/1‑26 (MQ=255)
                        tgtgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataat                           <  1:3190881‑M1/49‑24 (MQ=255)
                        tgtgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataa                            <  1:337242‑M1/48‑23 (MQ=255)
                         gtgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatg                          >  1:1599065‑M1/1‑25 (MQ=255)
                         gtgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatg                          <  1:5445071‑M1/49‑25 (MQ=255)
                          tgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgca                        <  1:6729543‑M1/50‑27 (MQ=255)
                          tgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgca                        <  1:7178691‑M1/50‑27 (MQ=255)
                          tgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgca                        >  1:5402268‑M1/1‑24 (MQ=255)
                          tgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgca                        >  1:8191503‑M1/1‑24 (MQ=255)
                          tgCGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgca                        <  1:2364912‑M1/50‑27 (MQ=255)
                            cGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatc                      >  1:1957408‑M1/1‑22 (MQ=255)
                            cGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatc                      <  1:1073990‑M1/50‑29 (MQ=255)
                            cGAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatc                      <  1:4222296‑M1/50‑29 (MQ=255)
                             gAGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatcc                     <  1:5584800‑M1/50‑30 (MQ=255)
                              aGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatcca                    <  1:3515964‑M1/50‑31 (MQ=255)
                              aGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatcca                    <  1:6198722‑M1/50‑31 (MQ=255)
                              aGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatcca                    <  1:3534997‑M1/50‑31 (MQ=255)
                              aGGAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatcca                    <  1:5069517‑M1/50‑31 (MQ=255)
                                gAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatccagc                  <  1:1093269‑M1/50‑33 (MQ=255)
                                gAAAAATATCAGTATGGCtcaatgatccaaataaagataatgcatccagc                  <  1:4699981‑M1/50‑33 (MQ=255)
                                 aaaaaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagct                 >  1:1700793‑M1/1‑17 (MQ=255)
                                 aaaaaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagct                 <  1:6306697‑M1/50‑34 (MQ=255)
                                 aaaaaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagct                 <  1:7025153‑M1/50‑34 (MQ=255)
                                 aaaaaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagct                 <  1:2600353‑M1/50‑34 (MQ=255)
                                 aaaaaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagct                 <  1:7501473‑M1/50‑34 (MQ=255)
                                 aaaaaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagct                 >  1:1109191‑M1/1‑17 (MQ=255)
                                 aaaaaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagct                 <  1:784667‑M1/50‑34 (MQ=255)
                                 aaaaaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagct                 >  1:7880181‑M1/1‑17 (MQ=255)
                                  aaaaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagctg                <  1:6422368‑M1/50‑35 (MQ=255)
                                  aaaaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagctg                <  1:2421493‑M1/50‑35 (MQ=255)
                                  aaaaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagctg                <  1:5635437‑M1/50‑35 (MQ=255)
                                  aaaaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagctg                >  1:2890447‑M1/1‑16 (MQ=255)
                                  aaaaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagctg                >  1:3854145‑M1/1‑16 (MQ=255)
                                  aaaaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagctg                <  1:8263753‑M1/50‑35 (MQ=255)
                                   aaaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagctgg               <  1:892409‑M1/50‑36 (MQ=255)
                                    aaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagctggt              <  1:5048281‑M1/50‑37 (MQ=255)
                                    aaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagctggt              <  1:1949090‑M1/50‑37 (MQ=255)
                                    aaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagctggt              <  1:84325‑M1/50‑37 (MQ=255)
                                    aaTATCAGTATGGCtcaatgatccaaataaagataatgcatccagctggt              <  1:1999784‑M1/50‑37 (MQ=255)
                                     atatCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtc             <  1:1314408‑M1/50‑38 (MQ=255)
                                      tatCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtca            <  1:603085‑M1/50‑39 (MQ=255)
                                        tCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtcatt          >  1:1446610‑M1/1‑10 (MQ=255)
                                        tCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtcatt          <  1:1274553‑M1/50‑41 (MQ=255)
                                        tCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtcatt          <  1:5825396‑M1/50‑41 (MQ=255)
                                        tCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtcatt          <  1:3369483‑M1/50‑41 (MQ=255)
                                        tCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtcat           <  1:5260586‑M1/49‑40 (MQ=255)
                                        tCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtcat           <  1:3432825‑M1/49‑40 (MQ=255)
                                        tCAGTATGGCtcaatgatccaaataaagataatgcatccagctggtcat           <  1:4707007‑M1/49‑40 (MQ=255)
                                         cAGTATGGCtcaatgatccaaataaagataatgcatccagctggtcattc         >  1:2155578‑M1/1‑9 (MQ=255)
                                          aGTATGGCtcaatgatccaaataaagataatgcatccagctggtcattct        >  1:1516196‑M1/1‑8 (MQ=255)
                                          aGTATGGCtcaatgatccaaataaagataatgcatccagctggtcattct        <  1:7176822‑M1/50‑43 (MQ=255)
                                           gTATGGCtcaatgatccaaataaagataatgcatccagctggtcattctt       <  1:2436325‑M1/50‑44 (MQ=255)
                                            tATGGCtcaatgatccaaataaagataatgcatccagctggtcattctta      <  1:4199519‑M1/50‑45 (MQ=255)
                                            tATGGCtcaatgatccaaataaagataatgcatccagctggtcattct        >  1:7708076‑M1/1‑6 (MQ=255)
                                             aTGGCtcaatgatccaaataaagataatgcatccagctggtcattcttag     >  1:1670515‑M1/1‑5 (MQ=255)
                                             aTGGCtcaatgatccaaataaagataatgcatccagctggtcattcttag     >  1:4589521‑M1/1‑5 (MQ=255)
                                              tGGCtcaatgatccaaataaagataatgcatccagctggtcattcttagc    <  1:2729826‑M1/50‑47 (MQ=255)
                                              tGGCtcaatgatccaaataaagataatgcatccagctggtcattcttag     >  1:6835619‑M1/1‑4 (MQ=255)
                                              tGGCtcaatgatccaaataaagataatgcatccagctggtcattcttag     <  1:447734‑M1/49‑46 (MQ=255)
                                              tGGCtcaatgatccaaataaagataatgcatccagctggtcattcttag     <  1:7300375‑M1/49‑46 (MQ=255)
                                                gCtcaatgatccaaataaagataatgcatccagctggtcattcttagcag  <  1:4221878‑M1/50‑49 (MQ=255)
                                                gCtcaatgatccaaataaagataatgcatccagctggtcattcttagcag  <  1:5353436‑M1/50‑49 (MQ=255)
                                                gCtcaatgatccaaataaagataatgcatccagctggtcattcttagca   <  1:8169273‑M1/49‑48 (MQ=255)
                                                gCtcaatgatccaaataaagataatgcatccagctggtcattcttagca   <  1:2834090‑M1/49‑48 (MQ=255)
                                                gCtcaatgatccaaataaagataatgcatccagctggtcattcttagca   <  1:2148764‑M1/49‑48 (MQ=255)
                                                gCtcaatgatccaaataaagataatgcatccagctggtcattcttagca   <  1:905710‑M1/49‑48 (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: