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 3584135–3583427 1222–1930 196 [2] [53] 59 [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)171 (1.510) 103/178 NT 98.8% 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

CGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTG  >  NC_000913/3584104‑3584235
                                |                                                                                                   
cGCTGTTGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTaa                                  >  2:328312/1‑100 (MQ=18)
 gCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTaaa                                 <  1:1943138/100‑1 (MQ=35)
  cTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTaaaa                                <  2:1707343/100‑1 (MQ=35)
  cTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTaaa                                 <  1:1007795/99‑1 (MQ=35)
   tGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAg                               <  1:687079/100‑1 (MQ=38)
   tGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATAAATAAGTTGGAGTCATTACCGACATTAAAAg                               >  2:1780776/1‑100 (MQ=38)
    gTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGt                              >  2:2669281/1‑100 (MQ=38)
    gTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGt                              >  1:3016291/1‑100 (MQ=38)
     tCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTa                             <  2:1832894/100‑1 (MQ=40)
      cGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTaa                            >  1:2836815/1‑100 (MQ=40)
      cGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTaa                            >  1:2081603/1‑100 (MQ=40)
          ctcAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATATGTTGGAGTCATTACCGACATTAAAAGTAACttt                        <  1:1991604/100‑1 (MQ=255)
           tcAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACtttt                       >  2:2518246/1‑100 (MQ=255)
           tcAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACtttt                       >  2:2041613/1‑100 (MQ=255)
            cAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACttttt                      >  1:758589/1‑100 (MQ=255)
             aaaaTCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACtttttt                     <  1:1561834/100‑1 (MQ=255)
             aaaaTCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACtttttt                     <  1:1966004/100‑1 (MQ=255)
              aaaTCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTg                    <  1:110349/100‑1 (MQ=255)
               aaTCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGa                   >  1:363266/1‑100 (MQ=255)
               aaTCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGa                   >  1:326995/1‑100 (MQ=255)
                aTCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTg                    <  1:413254/98‑1 (MQ=255)
                 tCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGaaa                 <  1:2755944/100‑1 (MQ=255)
                   ggTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAAtt               >  2:1631738/1‑100 (MQ=255)
                   ggTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAAtt               >  1:1829927/1‑100 (MQ=255)
                   ggTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAAtt               <  1:1363812/100‑1 (MQ=255)
                     tGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTag             >  2:199618/1‑100 (MQ=255)
                     tGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTag             <  1:1736712/100‑1 (MQ=255)
                     tGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTag             <  1:787560/100‑1 (MQ=255)
                     tGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACAGAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTag             <  2:2962130/100‑1 (MQ=255)
                       gAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTagag           <  2:678404/100‑1 (MQ=255)
                       gAGCTGCATGACAAAGTCATCGGGCATGATCTGAACATAAAACACTATCAATTAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTagag           <  1:491623/100‑1 (MQ=255)
                         gCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGat         >  1:2938981/1‑100 (MQ=255)
                         gCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGat         >  1:800897/1‑100 (MQ=255)
                         gCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGat         >  1:1816832/1‑100 (MQ=255)
                         gCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTAGGCGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGat         >  1:131585/1‑100 (MQ=255)
                          cTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGata        <  1:1612613/100‑1 (MQ=255)
                          cTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGata        >  2:1020266/1‑100 (MQ=255)
                           tGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTCTCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGataa       <  1:2149011/100‑1 (MQ=255)
                           tGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGataa       >  2:1062571/1‑100 (MQ=255)
                           tGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGataa       <  1:740359/100‑1 (MQ=255)
                           tGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGataa       <  1:1526437/100‑1 (MQ=255)
                           tGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGataa       >  2:2280425/1‑100 (MQ=255)
                            gCATGACAAAGTCATCGGGCCTTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGata        >  1:1861698/1‑98 (MQ=255)
                            gCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGataat      >  2:369709/1‑100 (MQ=255)
                            gCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGataat      <  1:1430497/100‑1 (MQ=255)
                            gCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGataat      >  2:2411686/1‑100 (MQ=255)
                             cATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGataata     >  1:2304260/1‑100 (MQ=255)
                             cATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGataata     >  2:2614345/1‑100 (MQ=255)
                              aTGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAatat    >  1:882259/1‑100 (MQ=255)
                              aTGACAAAGTCATCGGGCATTATCTGAACATAAAACACAATCAATAAGTTGGAGTCATTACCGACATTACAAGTAACTATTTTGAAATTCGAGata        >  1:2068088/1‑96 (MQ=17)
                               tGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATAtt   >  2:1846339/1‑100 (MQ=255)
                               tGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATAtt   >  1:1859132/1‑100 (MQ=255)
                               tGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATAtt   >  2:603492/1‑100 (MQ=255)
                                gACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTg  >  2:739156/1‑100 (MQ=255)
                                gACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTg  <  2:2200844/100‑1 (MQ=255)
                                gACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTg  <  2:1756499/100‑1 (MQ=255)
                                gACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTg  <  2:154954/100‑1 (MQ=255)
                                gACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTg  >  1:879203/1‑100 (MQ=255)
                                gACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTg  <  1:2652777/100‑1 (MQ=255)
                                |                                                                                                   
CGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTG  >  NC_000913/3584104‑3584235

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: