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

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 3582206 3584130–3583428 1223–1925 166 [1] [56] 61 [yhhZ]–[insB1] [yhhZ], yrhA, insA, [insB1]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 35822051 (0.010)164 (1.350) 91/198 0.1 99.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

AGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAA  >  NC_000913/3584101‑3584230
                              |                                                                                                   
aGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAt                                >  1:1401928/1‑100 (MQ=11)
aGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAt                                >  2:289047/1‑100 (MQ=11)
 gTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAtt                               >  1:2975459/1‑100 (MQ=31)
   cGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTaa                             >  1:1104301/1‑100 (MQ=33)
   cGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTaa                             <  1:134177/100‑1 (MQ=33)
    gCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTaaa                            <  2:2089567/100‑1 (MQ=35)
      tGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAg                          <  2:1241764/100‑1 (MQ=38)
       gTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGt                         >  2:303624/1‑100 (MQ=38)
         cGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTaa                       >  2:2172340/1‑100 (MQ=40)
         cGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTaa                       >  2:2757716/1‑100 (MQ=40)
          gTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAAc                      >  2:366606/1‑100 (MQ=255)
          gTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAAc                      <  1:274671/100‑1 (MQ=255)
             ctcAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACttt                   >  1:10658/1‑100 (MQ=255)
             ctcAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACttt                   <  2:596656/100‑1 (MQ=255)
             ctcAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACttt                   >  1:2846843/1‑100 (MQ=255)
              tcAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACtttt                  >  2:2554229/1‑100 (MQ=255)
               cAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATGAGTTGGAGTCATTACCGACATTAAAAGTAACttttt                 >  2:2148623/1‑100 (MQ=255)
               cAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACttttt                 >  2:2205798/1‑100 (MQ=255)
               cAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACttttt                 >  2:2151388/1‑100 (MQ=255)
               cAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACttttt                 >  1:2914009/1‑100 (MQ=255)
               cAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACttttt                 >  2:1363923/1‑100 (MQ=255)
               cAAAATAGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACttttt                 >  1:1785373/1‑100 (MQ=255)
                acaaacagTGGAGCTGCCTGACACCGTCATCGGGCATTATCTGAAAATAAAACACTATCAATATTTTGGAGTCATTACCGACATTAAAAGTAACtttttt                <  1:563344/93‑1 (MQ=17)
                aaaaTCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACtttttt                <  1:1381949/100‑1 (MQ=255)
                aaaaTCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACtttttt                <  2:1088211/100‑1 (MQ=255)
                 aaaTCGGTGGAGCTGCATGACAGAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTg               <  1:1930327/100‑1 (MQ=255)
                 aaaTCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTg               <  1:1930294/100‑1 (MQ=255)
                  tctCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGa              <  1:624486/98‑1 (MQ=255)
                  aaTCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTgn              <  1:1953967/100‑2 (MQ=255)
                  aaTCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGa              >  2:1873415/1‑100 (MQ=255)
                   aTCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGaa             >  2:935984/1‑100 (MQ=255)
                    tCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGaaa            <  2:970500/100‑1 (MQ=255)
                    tCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGaaa            <  1:366605/100‑1 (MQ=255)
                     cGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGaaga           >  2:624487/1‑98 (MQ=255)
                      ggTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATTAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAAtt          <  1:640484/100‑1 (MQ=255)
                      ggTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAAtt          >  1:1653824/1‑100 (MQ=255)
                      ggTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAAtt          >  2:3058502/1‑100 (MQ=255)
                       tccgatcTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTa         <  2:756989/94‑1 (MQ=255)
                       gTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTa         >  2:2242606/1‑100 (MQ=255)
                        tGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTag        <  1:2533024/100‑1 (MQ=255)
                        tGGAGCTGCATGAAAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTag        <  2:1982690/100‑1 (MQ=255)
                         ggAGCTGCATGACAAAGTCATCTGGCATTATCTGAACTTAAAACACTCTCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTaga       <  1:2067132/100‑1 (MQ=255)
                         ggAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTaga       <  1:1350570/100‑1 (MQ=255)
                         ggAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTaga       >  2:202973/1‑100 (MQ=255)
                         ggAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTaga       <  1:819507/100‑1 (MQ=255)
                         ggAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTaga       >  2:3049452/1‑100 (MQ=255)
                          gAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTagag      <  2:1420303/100‑1 (MQ=255)
                          gAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTagag      >  2:132986/1‑100 (MQ=255)
                          gAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTagag      <  1:2037995/100‑1 (MQ=255)
                           aGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTagag      >  1:1145394/1‑99 (MQ=255)
                           aGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAgaga     >  2:2244771/1‑100 (MQ=255)
                            gCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGat    <  1:1535637/100‑1 (MQ=255)
                             cTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTGTCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGata   <  1:1956459/100‑1 (MQ=255)
                             cTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGata   >  1:1628018/1‑100 (MQ=255)
                             cTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGata   <  2:391973/100‑1 (MQ=255)
                             cTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGata   >  1:514141/1‑100 (MQ=255)
                              tGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGataa  >  1:2060459/1‑100 (MQ=255)
                              tGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGataa  >  1:2525432/1‑100 (MQ=255)
                              tGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGataa  >  2:175607/1‑100 (MQ=255)
                              tGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGataa  >  1:793339/1‑100 (MQ=255)
                              tGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGata   >  1:1923131/1‑99 (MQ=255)
                              |                                                                                                   
AGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAA  >  NC_000913/3584101‑3584230

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