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 3584144–3583428 1223–1939 120 [0] [47] 48 [yhhZ]–[insB1] [yhhZ], yrhA, insA, [insB1]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 35822050 (0.000)107 (1.250) 72/178 NT 100% 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

CGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAA  >  NC_000913/3584104‑3584244
                                         |                                                                                                   
cGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTaa                                           >  2:1127580/1‑100 (MQ=33)
 gCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTaaa                                          >  1:2089280/1‑100 (MQ=35)
 gCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTaaa                                          >  2:59286/1‑100 (MQ=35)
 gCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTaaa                                          >  1:2238861/1‑100 (MQ=35)
    ggcGTTCTCCAAATGGGGGGAGCTGCATGACACAGTCAGCGGGAATTATCTGAACATAAAACACTATCAATAAGTTGGGGTCCTCACCGACATTAAAAg                                        <  2:1941225/97‑1 (MQ=11)
      cgtggtcAAAATCGGTGGAGCTGCATGACAAAGTTATCGGGCATTATCTGAACATAAAAGACTATAAATAAGGTAGAGTCATTACCGAAATTAAAAGTaa                                     >  2:2184537/6‑100 (MQ=14)
       gTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAAc                                    >  1:239658/1‑100 (MQ=255)
        ttCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACt                                   >  2:2198831/1‑100 (MQ=255)
         tctcAAAATCGGTGGAGCTGCATGACAACGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGACGGAGTCATTACCGACATTAAAAGCAAc                                    >  2:29382/1‑98 (MQ=14)
          ctcAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACttt                                 >  1:647688/1‑100 (MQ=255)
           tcAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACtttt                                >  2:672631/1‑100 (MQ=255)
           tcAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACtttt                                >  2:1211310/1‑100 (MQ=255)
            cAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCACTAAGTTGGAGTCATTACCGACATTAAAAGTAACttttt                               >  2:1488762/1‑100 (MQ=255)
            cAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACttttt                               >  2:1917837/1‑100 (MQ=255)
             aaaaTCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACtttttt                              <  2:1982282/100‑1 (MQ=255)
              aaaTCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTg                             <  1:57640/100‑1 (MQ=255)
                aTCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGaa                           <  1:731025/100‑1 (MQ=255)
                    gTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTa                       >  1:958204/1‑100 (MQ=255)
                    gTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTa                       >  1:352250/1‑100 (MQ=255)
                    gTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAACTa                       >  1:2185526/1‑100 (MQ=255)
                     ttgAGCTGAATGACAAAGTCATCGGGCACAATCGGAAAATAAAATACTATCAATAAATTGGAGAAATTAACAACTATAAAAGtccgtgt                                 >  2:1405580/3‑83 (MQ=2)
                     tGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTag                      <  1:154976/100‑1 (MQ=255)
                     acaaGCTGCATTACACACTCATCCCGCATTACCAAACCATTAAACACAATCCATAAGTTGCAATCATTACCTACATTGAAAGTAACTTTTTTGAAATTac                      >  2:614833/4‑99 (MQ=255)
                       gAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAGCACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTaga                     <  2:1313784/99‑1 (MQ=255)
                       gAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTa                       <  2:1938661/97‑1 (MQ=255)
                        aGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAgaga                   >  2:1964997/1‑100 (MQ=255)
                         gCTGCATGACAAAGTCATCGGGCATTATCTAAAGATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGCAATTAgaga                   >  2:1569286/1‑99 (MQ=255)
                           tGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGataa                >  1:2293670/1‑100 (MQ=255)
                           tGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGataa                >  2:2097002/1‑100 (MQ=255)
                           tGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGataa                >  1:623176/1‑100 (MQ=255)
                            gCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGACATTATAGata                 >  2:770224/1‑98 (MQ=255)
                            gCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGataat               >  1:1988702/1‑100 (MQ=255)
                            gCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGataat               >  1:1712284/1‑100 (MQ=255)
                            gCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGataat               >  1:530884/1‑100 (MQ=255)
                             cATGACAAAGTAATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGATTCATTACCGACATTAAAAGTAACTTTTTTGAAATTagag                    <  2:2019565/94‑1 (MQ=255)
                                gACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTg           <  1:842780/100‑1 (MQ=255)
                                 aCAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTgg          <  2:1767663/100‑1 (MQ=255)
                                 aCAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTgg          >  1:507877/1‑100 (MQ=255)
                                  cAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGc         <  2:593096/100‑1 (MQ=255)
                                  cAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGc         <  2:801944/100‑1 (MQ=255)
                                  cAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGc         >  2:1266850/1‑100 (MQ=255)
                                   aaaGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGca        <  2:4957/100‑1 (MQ=255)
                                     aGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGcaca      >  1:1329904/1‑100 (MQ=255)
                                      gTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAg     >  1:10017/1‑100 (MQ=255)
                                       tCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGa    >  2:1981257/1‑100 (MQ=255)
                                        cATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGaa   >  2:2035042/1‑100 (MQ=255)
                                        cATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGaa   >  1:896839/1‑100 (MQ=255)
                                         aTCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGaaa  <  1:495456/100‑1 (MQ=255)
                                         |                                                                                                   
CGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAA  >  NC_000913/3584104‑3584244

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