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 131 total aligned reads displayed.

ATAATTAGTCCTTACAAGGAGCAAAGTTCCATAATCTCCCCCCTCCCCCTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCATTCTTAGCAGCA  >  NC_000913/3770194‑3770292
                                                 |                                                 
gacgagtttcatgccgagtcctttgtgcgaggaaaaatatcagtatggct                                                   <  1:3794707‑M2/1‑1 (MQ=255)
gacgagtttcatgccgagtcctttgtgcgaggaaaaatatcagtatggct                                                   >  1:80115‑M2/50‑50 (MQ=255)
  cgagtttcatgccgagtcctttgtgcgaggaaaaatatcagtatggcTCa                                                 <  1:657008‑M2/3‑1 (MQ=255)
   gagtttcatgccgagtcctttgtgcgaggaaaaatatcagtatggcTCaa                                                <  1:3961173‑M2/4‑1 (MQ=255)
    agtttcatgccgagtcctttgtgcgaggaaaaatatcagtatggcTCAAt                                               >  1:2062711‑M2/46‑50 (MQ=255)
    agtttcatgccgagtcctttgtgcgaggaaaaatatcagtatggcTCAAt                                               >  1:6345408‑M2/46‑50 (MQ=255)
     gtttcatgccgagtcctttgtgcgaggaaaaatatcagtatggcTCAAt                                               >  1:4174468‑M2/45‑49 (MQ=255)
     gtttcataccgagtcctttgtccgatgataaatatcagtatggcTCAAt                                               <  1:1154482‑M2/5‑1 (MQ=255)
      tttcatgccgagtcctttgtgcgaggaaaaatatcagtatggcTCAATGa                                             >  1:123231‑M2/44‑50 (MQ=255)
      tttcatgccgagtcctttgtgcgaggaaaaatatcagtatggcTCAATGa                                             <  1:5030833‑M2/7‑1 (MQ=255)
        tcatgccgagtcctttgtgcgaggaaaaatatcagtatggcTCAATGAt                                            <  1:3274147‑M2/8‑1 (MQ=255)
        tcatgccgagtcctttgtgcgaggaaaaatatcagtatggcTCAATGAt                                            <  1:2937149‑M2/8‑1 (MQ=255)
        tcatgccgagtcctttgtgcgaggaaaaatatcagtatggcTCAATGATc                                           >  1:1353551‑M2/42‑50 (MQ=255)
          atgccgagtcctttgtgcgaggaaaaatatcagtatggcTCAATGATcc                                          <  1:4624594‑M2/10‑1 (MQ=255)
          atgccgagtcctttgtgcgaggaaaaatatcagtatggcTCAATGATCCa                                         <  1:640190‑M2/11‑1 (MQ=255)
           tgccgagtcctttgtgcgaggaaaaatatcagtatggcTCAATGATCCaa                                        <  1:2976401‑M2/12‑1 (MQ=255)
           tgccgagtcctttgtgcgaggaaaaatatcagtatggcTCAATGATCCaa                                        <  1:3047677‑M2/12‑1 (MQ=255)
            gccgagtcctttgtgcgaggaaaaatatcagtatggcTCAATGATCCaaa                                       >  1:170355‑M2/38‑50 (MQ=255)
             ccgagtcctttgtgcgaggaaaaatatcagtatggcTCAATGATCCaaa                                       <  1:4429806‑M2/13‑1 (MQ=255)
             ccgagtcctttgtgcgaggaaaaatatcagtatggcTCAATGATCCaaa                                       <  1:5795408‑M2/13‑1 (MQ=255)
              cgagtcctttgtgcgaggaaaaatatcagtatggcTCAATGATCCAAAt                                      >  1:6246046‑M2/36‑49 (MQ=255)
              cgagtcctttgtgcgaggaaaaatatcagtatggcTCAATGATCCAAATa                                     <  1:5519135‑M2/15‑1 (MQ=255)
               gagtcctttgtgcgaggaaaaatatcagtatggcTCAATGATCCAAATaa                                    <  1:2686779‑M2/16‑1 (MQ=255)
                agtcctttgtgcgaggaaaaatatcagtatggcTCAATGATCCAAATaaa                                   <  1:377993‑M2/17‑1 (MQ=255)
                agtcctttgtgcgaggaaaaatatcagtatggcTCAATGATCCAAATaaa                                   <  1:5675615‑M2/17‑1 (MQ=255)
                 gtcctttgtgcgaggaaaaatatcagtatggcTCAATGATCCAAATAAAg                                  <  1:1770590‑M2/18‑1 (MQ=255)
                    ctttgtgcgaggaaaaatatcagtatggcTCAATGATCCAAATAAAGAt                                >  1:3560025‑M2/30‑49 (MQ=255)
                    ctttgtgcgaggaaaaatatcagtatggcTCAATGATCCAAATAAAGATa                               <  1:3529891‑M2/21‑1 (MQ=255)
                     tttgtgcgaggaaaaatatcagtatggcTCAATGATCCAAATAAAGATaa                              >  1:2128358‑M2/29‑50 (MQ=255)
                      ttgtgcgaggaaaaatatcagtatggcTCAATGATCCAAATAAAGATAAt                             >  1:2730336‑M2/28‑50 (MQ=255)
                       tgtgcgaggaaaaatatcagtatggcTCAATGATCCAAATAAAGATaa                              <  1:1681698‑M2/22‑1 (MQ=255)
                       tgtgcgaggaaaaatatcagtatggcTCAATGATCCAAATAAAGATAAt                             <  1:301234‑M2/23‑1 (MQ=255)
                       tgtgcgaggaaaaatatcagtatggcTCAATGATCCAAATAAAGATAATg                            >  1:1363088‑M2/27‑50 (MQ=255)
                        gtgcgaggaaaaatatcagtatggcTCAATGATCCAAATAAAGATAATg                            <  1:6084363‑M2/24‑1 (MQ=255)
                        gtgcgaggaaaaatatcagtatggcTCAATGATCCAAATAAAGATAATg                            <  1:2769707‑M2/24‑1 (MQ=255)
                        gtgcgaggaaaaatatcagtatggcTCAATGATCCAAATAAAGATAATg                            <  1:5972859‑M2/24‑1 (MQ=255)
                        gtgcgaggaaaaatatcagtatggcTCAATGATCCAAATAAAGATAATg                            <  1:3805973‑M2/24‑1 (MQ=255)
                        gtgcgaggaaaaatatcagtatggcTCAATGATCCAAATAAAGATAATg                            <  1:250114‑M2/24‑1 (MQ=255)
                         tgcgaggaaaaatatcagtatggcTCAATGATCCAAATAAAGATAATGCa                          >  1:4379544‑M2/25‑50 (MQ=255)
                         tgcgaggaaaaatatcagtatggcTCAATGATCCAAATAAAGATAATGCa                          >  1:6621074‑M2/25‑50 (MQ=255)
                           cgaggaaaaatatcagtatggcTCAATGATCCAAATAAAGATAATGCAt                         <  1:1104481‑M2/27‑1 (MQ=255)
                           cgaggaaaaatatcagtatggcTCAATGATCCAAATAAAGATAATGCAt                         <  1:5310864‑M2/27‑1 (MQ=255)
                           cgaggaaaaatatcagtatggcTCAATGATCCAAATAAAGATAATGCATc                        <  1:1526710‑M2/28‑1 (MQ=255)
                            gaggaaaaatatcagtatggcTCAATGATCCAAATAAAGATAATGCATcc                       <  1:4044667‑M2/29‑1 (MQ=255)
                            gaggaaaaatatcagtatggcTCAATGATCCAAATAAAGATAATGCATcc                       <  1:3274943‑M2/29‑1 (MQ=255)
                             aggaaaaatatcagtatggcTCAATGATCCAAATAAAGATAATGCATc                        >  1:64635‑M2/21‑48 (MQ=255)
                             aggaaaaatatcagtatggcTCAATGATCCAAATAAAGATAATGCATc                        >  1:922183‑M2/21‑48 (MQ=255)
                             aggaaaaatatcagtatggcTCAATGATCCAAATAAAGATAATGCATCCa                      <  1:2641075‑M2/30‑1 (MQ=255)
                              ggaaaaatatcagtatggcTCAATGATCCAAATAAAGATAATGCATCCAg                     >  1:3554523‑M2/20‑50 (MQ=255)
                              ggaaaaatatcagtatggcTCAATGATCCAAATAAAGATAATGCATCCAg                     >  1:2202672‑M2/20‑50 (MQ=255)
                              ggaaaaatatcagtatggcTCAATGATCCAAATAAAGATAATGCATCCAg                     >  1:3979506‑M2/20‑50 (MQ=255)
                              ggaaaaatatcagtatggcTCAATGATCCAAATAAAGATAATGCATCCAg                     <  1:5505086‑M2/31‑1 (MQ=255)
                              ggaaaaatatcagtatggcTCAATGATCCAAATAAAGATAATGCATCCAg                     <  1:1927399‑M2/31‑1 (MQ=255)
                               gaaaaatatcagtatggcTCAATGATCCAAATAAAGATAATGCATCCAGc                    >  1:1979909‑M2/19‑50 (MQ=255)
                               gaaaaatatcagtatggcTCAATGATCCAAATAAAGATAATGCATCCAGc                    >  1:2864263‑M2/19‑50 (MQ=255)
                                aaaaatatcagtatggcTCAATGATCCAAATAAAGATAATGCATCCAGCt                   <  1:3731195‑M2/33‑1 (MQ=255)
                                 aaaatatcagtatggcTCAATGATCCAAATAAAGATAATGCATCCAGCTg                  >  1:2345129‑M2/17‑50 (MQ=255)
                                 aaaatatcagtatggcTCAATGATCCAAATAAAGATAATGCATCCAGCTg                  <  1:735364‑M2/34‑1 (MQ=255)
                                 aaaatatcagtatggcTCAATGATCCAAATAAAGATAATGCATCCAGCTg                  <  1:2296146‑M2/34‑1 (MQ=255)
                                 aaaatatcagtatggcTCAATGATCCAAATAAAGATAATGCATCCAGCTg                  <  1:2142536‑M2/34‑1 (MQ=255)
                                 aaaatatcagtatggcTCAATGATCCAAATAAAGATAATGCATCCAGCTg                  <  1:5494381‑M2/34‑1 (MQ=255)
                                 aaaatatcagtatggcTCAATGATCCAAATAAAGATAATGCATCCAGCTg                  <  1:4524709‑M2/34‑1 (MQ=255)
                                   aatatcagtatggcTCAATGATCCAAATAAAGATAATGCATCCAGCTGGt                <  1:31830‑M2/36‑1 (MQ=255)
                                   aatatcagtatggcTCAATGATCCAAATAAAGATAATGCATCCAGCTGGt                <  1:890207‑M2/36‑1 (MQ=255)
                                    atatcagtatggcTCAATGATCCAAATAAAGATAATGCATCCAGCTGGt                >  1:1357325‑M2/14‑49 (MQ=255)
                                    atatcagtatggcTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTc               <  1:4819148‑M2/37‑1 (MQ=255)
                                    atatcagtatggcTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTc               <  1:6019706‑M2/37‑1 (MQ=255)
                                    atatcagtatggcTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTc               <  1:2136562‑M2/37‑1 (MQ=255)
                                    atatcagtatggcTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTc               <  1:5184397‑M2/37‑1 (MQ=255)
                                     tatcagtatggcTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCa              >  1:4183943‑M2/13‑50 (MQ=255)
                                     tatcagtatggcTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCa              >  1:2787812‑M2/13‑50 (MQ=255)
                                     tatcagtatggcTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCa              <  1:3033805‑M2/38‑1 (MQ=255)
                                      atcagtatggcTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCAt             >  1:3486383‑M2/12‑50 (MQ=255)
                                       tcagtatggcTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCAtt            <  1:1284418‑M2/40‑1 (MQ=255)
                                       tcagtatggcTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCAtt            <  1:5360427‑M2/40‑1 (MQ=255)
                                       tcagtatggcTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCAtt            >  1:2411729‑M2/11‑50 (MQ=255)
                                       tcagtatggcTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCAt             <  1:3567200‑M2/39‑1 (MQ=255)
                                         agtatggcTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCATTc           <  1:6188207‑M2/41‑1 (MQ=255)
                                         agtatggcTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCATTc           <  1:3984745‑M2/41‑1 (MQ=255)
                                         agtatggcTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCATTc           <  1:312956‑M2/41‑1 (MQ=255)
                                         agtatggcTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCATTCt          >  1:519208‑M2/9‑50 (MQ=255)
                                         agtatggcTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCATTCt          <  1:1923498‑M2/42‑1 (MQ=255)
                                         agtatggcTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCATTCt          >  1:6248123‑M2/9‑50 (MQ=255)
                                          gtatggcTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCATTCtt         <  1:4907101‑M2/43‑1 (MQ=255)
                                          gtatggcTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCATTCt          >  1:4565885‑M2/8‑49 (MQ=255)
                                          gtatggcTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCATTCt          >  1:4836843‑M2/8‑49 (MQ=255)
                                           tatggcTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCATTCt          >  1:438346‑M2/7‑48 (MQ=255)
                                           tatggcTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCATTCt          >  1:3713214‑M2/7‑48 (MQ=255)
                                           tatggcTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCATTCt          >  1:3499135‑M2/7‑48 (MQ=255)
                                           tatggcTCAATGATCCAAATAAAGACAATGCATCCAGCTGGTCATTCTTa        <  1:6420558‑M2/44‑1 (MQ=255)
                                             tggcTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCATTCTTag       <  1:5271749‑M2/45‑1 (MQ=255)
                                             tggcTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCATTCTTag       >  1:5844960‑M2/5‑49 (MQ=255)
                                              ggcTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCATTCTTag       >  1:3976914‑M2/4‑48 (MQ=255)
                                              ggcTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCATTCTTag       >  1:3694025‑M2/4‑48 (MQ=255)
                                               gcTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCATTCTTagcag    >  1:4744429‑M2/3‑50 (MQ=255)
                                               gcTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCATTCTTagcag    >  1:2137213‑M2/3‑50 (MQ=255)
                                               gcTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCATTCTTagca     <  1:3405474‑M2/47‑1 (MQ=255)
                                               gcTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCATTCTTagca     <  1:1833567‑M2/47‑1 (MQ=255)
                                                cTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCATTCTTagcagc   <  1:1710484‑M2/49‑1 (MQ=255)
                                                 tcAATGATCCAAATAAAGATAATGCATCCAGCTGGTCATTCTTAgcagca  >  1:3356079/1‑50 (MQ=255)
                                                 |                                                 
ATAATTAGTCCTTACAAGGAGCAAAGTTCCATAATCTCCCCCCTCCCCCTCAATGATCCAAATAAAGATAATGCATCCAGCTGGTCATTCTTAGCAGCA  >  NC_000913/3770194‑3770292

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