Predicted mutation
evidence seq id position mutation freq annotation gene description
MC JC NC_000913 3,582,206 Δ1,222 bp 100% IS1‑mediated [yhhZ]yrhA [yhhZ], yrhA

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 3582206 3584138–3583427 1222–1933 146 [2] [70] 74 [yhhZ]–[insB1] [yhhZ], yrhA, insA, [insB1]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 35822052 (0.020)127 (1.060) 82/178 NT 98.4% coding (343/1179 nt) yhhZ putative Hcp1 family polymorphic toxin protein; putative colicin‑like DNase/tRNase activity
?NC_000913 3583428 = NA (NA)noncoding (1/768 nt) IS1 repeat region

CGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCA  >  NC_000913/3584104‑3584238
                                   |                                                                                                   
cGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAGAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTaa                                     <  1:2280956/100‑1 (MQ=18)
 gCTGTCGTTCTCAAAATCGGTGGAGCTGCATTACAAAGTCATCGGGCATTATCTGCACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTaaa                                    <  2:2773342/100‑1 (MQ=255)
 gCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTaaa                                    <  1:879742/100‑1 (MQ=35)
 gCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTaaa                                    <  2:3152986/100‑1 (MQ=35)
  cTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTaaaa                                   <  1:2962358/100‑1 (MQ=35)
  cTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTaaaa                                   <  1:157032/100‑1 (MQ=35)
  cTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTaaaa                                   <  1:2155790/100‑1 (MQ=35)
  cTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTaaaa                                   >  2:3135754/1‑100 (MQ=35)
   tGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAg                                  <  2:1112615/100‑1 (MQ=38)
   tGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAg                                  <  1:2035428/100‑1 (MQ=38)
    gTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGt                                 >  2:430955/1‑100 (MQ=38)
     tCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTa                                <  1:1849948/100‑1 (MQ=40)
      cGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTaa                               >  1:1078952/1‑100 (MQ=40)
      cGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTaa                               <  1:436625/100‑1 (MQ=40)
       gTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAAc                              >  2:370970/1‑100 (MQ=255)
        ttCTCAAAATCTTTGTAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAAc                              <  1:1882479/99‑1 (MQ=255)
         tctcAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACtt                            >  1:1827770/1‑100 (MQ=255)
         tctcAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACtt                            >  2:534147/1‑100 (MQ=255)
           tcAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCCATAAGTTGGAGTCATTACCGACATTAAAAGTAACttt                           <  1:432501/99‑1 (MQ=255)
           tcAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACtttt                          >  2:2033369/1‑100 (MQ=255)
           tcAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACtttt                          >  2:3233972/1‑100 (MQ=255)
           tcAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACtttt                          >  1:1494089/1‑100 (MQ=255)
           tcAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACtttt                          >  1:279640/1‑100 (MQ=255)
           tcAAAATCGGTGGAGCTGCATGACAAAGTCAGCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATCAAAAGTAACtt                            >  1:1479435/1‑98 (MQ=255)
            cAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACttttt                         <  2:962440/100‑1 (MQ=255)
             aaaaTCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACtttttt                        <  1:1490073/100‑1 (MQ=255)
             aaaaTCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACtttttt                        >  1:2509831/1‑100 (MQ=255)
              aaaTCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTg                       <  2:1327001/100‑1 (MQ=255)
              aaaTCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTg                       >  2:621273/1‑100 (MQ=255)
              aaaTCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTg                       <  1:327984/100‑1 (MQ=255)
              aaaTCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTg                       <  2:1828762/100‑1 (MQ=255)
               aaTCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGa                      <  1:526065/100‑1 (MQ=255)
                 tCGTTGGAAATGCATGACAAAGTCATAGGGCATTATCTGAATATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGaaa                    <  2:765620/100‑1 (MQ=17)
                  cGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAAt                   >  2:405388/1‑100 (MQ=255)
                   gttGGAGCTGCATGACAAAGTCATCGGGCATTATCTTAACATAAATCACTATCAATAAGTTGTAGTCATTACCGACATTAAAAGTAACTTTTTTGaaa                    <  2:765581/96‑1 (MQ=17)
                   ggTGGAGTTTCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAAtt                  <  1:68349/100‑1 (MQ=255)
                   ggTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACttttttttattt                  >  2:2524977/1‑94 (MQ=255)
                   ggTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAAtt                  >  1:689271/1‑100 (MQ=255)
                   ggTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTGTTTTGAAAtt                  >  2:2136762/1‑100 (MQ=255)
                    gTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTa                 <  1:3026132/100‑1 (MQ=255)
                     tGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTag                <  2:1870438/100‑1 (MQ=255)
                     tGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTag                <  2:1048643/100‑1 (MQ=255)
                      ggAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTaga               <  2:1048360/100‑1 (MQ=255)
                      ggAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTaga               <  1:1268513/100‑1 (MQ=255)
                       gAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTagag              <  2:940864/100‑1 (MQ=255)
                       gAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTagag              <  1:1017550/100‑1 (MQ=255)
                       gAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTagag              <  1:2280249/100‑1 (MQ=255)
                        aGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAgaga             <  1:783488/100‑1 (MQ=255)
                        aGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTATTTTGAAATTAgaga             >  1:3206919/1‑100 (MQ=255)
                         gCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGat            >  1:3217722/1‑100 (MQ=255)
                         gCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGat            <  1:714506/100‑1 (MQ=255)
                         gCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGat            <  1:743576/100‑1 (MQ=255)
                          cTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGata           <  2:2566098/100‑1 (MQ=255)
                          cTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGata           <  1:2443408/100‑1 (MQ=255)
                          cTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGat            <  1:1439877/99‑1 (MQ=255)
                           tGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGataa          >  2:2313696/1‑100 (MQ=255)
                           tGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGataa          >  2:195045/1‑100 (MQ=255)
                           tGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGataa          <  1:3234000/100‑1 (MQ=255)
                            gCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGataat         >  1:2207982/1‑100 (MQ=255)
                            gCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGataat         >  2:2794500/1‑100 (MQ=255)
                            gCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGataat         >  2:364168/1‑100 (MQ=255)
                            gCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGataa          <  2:955676/99‑1 (MQ=255)
                              aTGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAatat       >  1:3168877/1‑100 (MQ=255)
                               tGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATAtt      >  1:2030685/1‑100 (MQ=255)
                                 aCAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTgg    <  2:1941448/100‑1 (MQ=255)
                                 aCAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTgg    >  1:1578189/1‑100 (MQ=255)
                                 aCAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTgg    >  1:181753/1‑100 (MQ=255)
                                 aCAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTgg    <  2:2967027/100‑1 (MQ=255)
                                  cAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGc   >  1:2099168/1‑100 (MQ=255)
                                  cAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGc   >  1:3179365/1‑100 (MQ=255)
                                   aaaGTCATCGGGCATTATCTGAACATAAAACACTGTCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGc   <  2:1468251/99‑1 (MQ=255)
                                   aaaGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTgg  <  2:1468206/99‑2 (MQ=255)
                                   aaaGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGca  <  2:571688/100‑1 (MQ=255)
                                   aaaGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGca  >  2:1230885/1‑100 (MQ=255)
                                   |                                                                                                   
CGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCA  >  NC_000913/3584104‑3584238

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