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 3584133–3583427 1222–1928 203 [0] [64] 68 [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)181 (1.300) 98/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

GCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATAT  >  NC_000913/3584054‑3584233
                                                                                |                                                                                                   
gctgCATATCCTGAATCTGAGGCAGCACCTGGCACGTCTGGGACGGAAGTCGCTGTCGTGCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCa                                                                                   <  1:778778/96‑1 (MQ=1)
                                               aGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAt                                   >  1:1225753/1‑100 (MQ=11)
                                               aGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAt                                   >  2:1386631/1‑100 (MQ=11)
                                                gTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAtt                                  >  1:539169/1‑100 (MQ=31)
                                                 tCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTa                                 <  1:3236929/100‑1 (MQ=32)
                                                 tCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAAAACTATCAATAAGTTGGAGACATTACCGACATTa                                 >  1:986717/1‑100 (MQ=12)
                                                  cGCTGTCGTTCTCAAAATCGGTGGCGCTGCATGACAAAGTCATCGGGCATTATCTGAAAATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTaa                                >  1:1141264/1‑100 (MQ=16)
                                                    cTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTaaaa                              >  2:1006616/1‑100 (MQ=35)
                                                    cTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTaaaa                              >  1:231866/1‑100 (MQ=35)
                                                    cTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTaaaa                              <  1:3378145/100‑1 (MQ=35)
                                                     tGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAg                             >  1:3812634/1‑100 (MQ=38)
                                                        cGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATAAATAAGTTGGAGTCATTACCGACATTAAAAGTaa                          >  1:3137770/1‑100 (MQ=255)
                                                         gTTCTCAAAATCGGTTCAGCTGCATGACAAATTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGTGTCATTACCGACATTAAAAGTa                           <  2:1646804/98‑1 (MQ=14)
                                                         gTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAAc                         <  2:1488788/100‑1 (MQ=255)
                                                         gTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAAc                         >  1:1029095/1‑100 (MQ=255)
                                                          ttCTCAAAATCGGTGGAGCTGCATGACAAAGTCTTCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATGAAAAGTAACt                        >  2:2564381/1‑100 (MQ=255)
                                                            ctcAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACttt                      >  2:1254148/1‑100 (MQ=255)
                                                            ctcAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACttt                      >  2:88766/1‑100 (MQ=255)
                                                            ctcAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAAAACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACttt                      >  1:971314/1‑100 (MQ=255)
                                                             tcAAAATGGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAGGTTGGAGTCATTACCGACATTAAAAGTAACtttt                     >  2:2047719/1‑100 (MQ=255)
                                                             tcAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACtttt                     >  1:469920/1‑100 (MQ=255)
                                                             tcAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACtttt                     >  1:579092/1‑100 (MQ=255)
                                                             tcAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACtttt                     >  2:1074664/1‑100 (MQ=255)
                                                             tcAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACtttt                     >  2:2770021/1‑100 (MQ=255)
                                                             tcAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACtttt                     >  2:725747/1‑100 (MQ=255)
                                                              cAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACttttt                    >  1:3050824/1‑100 (MQ=255)
                                                              cAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACttttt                    >  1:989280/1‑100 (MQ=255)
                                                              cAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACttttt                    >  1:2770735/1‑100 (MQ=255)
                                                              cAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACttttt                    <  1:1780030/100‑1 (MQ=255)
                                                                aaaTCGGTGGAGCTGTTTGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTg                  <  1:1880764/100‑1 (MQ=255)
                                                                aaaTCGGTGGAGCTGCATTACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTg                  <  1:3590389/100‑1 (MQ=255)
                                                                aaaTCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTg                  <  1:1076429/100‑1 (MQ=255)
                                                                aaaTCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTg                  <  2:3433395/100‑1 (MQ=255)
                                                                aaaTCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTg                  <  1:645062/100‑1 (MQ=255)
                                                                aaaTCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTg                  <  1:1340775/100‑1 (MQ=255)
                                                                aaaTCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTg                  >  1:303800/1‑100 (MQ=255)
                                                                aaaTCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTg                  <  1:1881348/100‑1 (MQ=255)
                                                                aaaTCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGATGGAGTCATTACCGACATTAAAAGTAACTTTTTTg                  <  2:2317611/100‑1 (MQ=255)
                                                                 aaTCGGTGGATCTGCATGAAAAAGTTATTGGGCATTATCTGAATATAAAACACCATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACtttttt                   <  2:484104/98‑1 (MQ=17)
                                                                 aaTCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGa                 >  2:1994641/1‑100 (MQ=255)
                                                                    cGGTGGAGCTGCATTACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAAt              >  2:3476849/1‑100 (MQ=255)
                                                                    cGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAATTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAAt              <  1:631188/100‑1 (MQ=255)
                                                                    cGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAAt              <  1:1730349/100‑1 (MQ=255)
                                                                     ggTGGAGCTGCATGACAAAGTCATCTGGGATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTTCCGACATTAAAAGTAACTTTTTTGAAAtt             <  2:3461862/100‑1 (MQ=255)
                                                                     ggTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAAtt             >  2:610591/1‑100 (MQ=255)
                                                                     ggTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAAtt             >  1:3220775/1‑100 (MQ=255)
                                                                     ggTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAAtt             >  1:1381355/1‑100 (MQ=255)
                                                                     ggTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAATCATTACCGAAATTAAAAGTAACTTTTTTGAAAt              >  1:1908806/1‑99 (MQ=255)
                                                                     ggTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGAAGTCATTACCGACATTAAAAGTAACTTTTTTGAAAtt             <  2:2083249/100‑1 (MQ=255)
                                                                      gTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTa            <  2:945160/100‑1 (MQ=255)
                                                                       tGGAGCTGCATGACAAAGTCATCGGGCATTATCTGATCATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTag           <  1:1269687/100‑1 (MQ=255)
                                                                       tGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTag           <  1:2980661/100‑1 (MQ=255)
                                                                       tGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTGGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTa            >  2:3893012/1‑99 (MQ=255)
                                                                        ggAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATAAATAAATTAGAGTAATTAAAGAAATTAAAAGTAAATTTTTTGAAATTaga          >  1:1480176/1‑100 (MQ=17)
                                                                          aGCTGCATTACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTaga          <  2:122077/98‑1 (MQ=255)
                                                                           gCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGat       >  1:3263/1‑100 (MQ=255)
                                                                             tGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGataa     >  2:3213959/1‑100 (MQ=255)
                                                                             tGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGataa     >  1:2779551/1‑100 (MQ=255)
                                                                             tGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTAACGACATTAAAAGTAACTTTTTTGAAATTAGAGataa     >  1:23662/1‑100 (MQ=255)
                                                                             tGCATGACAAAGTCATCGGGCATTATCTGAACATAAAAAACTATCAATAAGTTGTAGTCATTACTGAAATTAAAAGTAACTTTTTTGAAATTAGAGataa     >  2:3697541/1‑100 (MQ=21)
                                                                              gCATGACAAAGTCATCGGGCGTTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGataat    >  2:2382518/1‑100 (MQ=255)
                                                                              gCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGataat    >  2:3439787/1‑100 (MQ=255)
                                                                              gCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGAAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGataat    >  1:1100020/1‑100 (MQ=255)
                                                                              gCATGACAAAGTCATCGGGCATGATAACAACATAAAAAACTATCAATAAGATAGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTATAGataa     >  2:1291264/1‑99 (MQ=18)
                                                                                aTGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAatat  >  2:3389383/1‑100 (MQ=255)
                                                                                aTGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAatat  >  2:3029843/1‑100 (MQ=255)
                                                                                aTGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAatat  >  1:2036176/1‑100 (MQ=255)
                                                                                aTGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAatat  >  1:2384581/1‑100 (MQ=255)
                                                                                |                                                                                                   
GCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATAT  >  NC_000913/3584054‑3584233

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: