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 3584130–3583428 1223–1925 190 [2] [56] 57 [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.010)171 (1.250) 93/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

AGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAA  >  NC_000913/3584044‑3584230
                                                                                       |                                                                                                   
agtgAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAg                                                                                          <  2:1245528/96‑1 (MQ=1)
                                                         aGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAt                                >  1:2077665/1‑100 (MQ=11)
                                                          gTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAtt                               <  1:1631188/100‑1 (MQ=31)
                                                           tCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTa                              <  2:42728/100‑1 (MQ=17)
                                                           tCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTa                              <  2:42793/100‑1 (MQ=32)
                                                           tCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTa                              <  2:1377975/100‑1 (MQ=32)
                                                           tCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACAGACATTa                              <  1:1994364/100‑1 (MQ=17)
                                                            cGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTCTCAATAAGTTGGAGTCATTACCGACATTaa                             >  2:620338/1‑100 (MQ=18)
                                                            cGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTaa                             >  2:1593384/1‑100 (MQ=33)
                                                            cGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTaa                             >  2:1939249/1‑100 (MQ=33)
                                                            cGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTa                              <  2:1793552/99‑1 (MQ=32)
                                                             gCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTaaa                            <  1:1454131/100‑1 (MQ=35)
                                                             gCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTaaa                            <  1:21208/100‑1 (MQ=35)
                                                             gCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGGCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTaaa                            <  1:801660/100‑1 (MQ=25)
                                                                gTCGTTCTCAACATCGGTGGAGCGGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCGATAAGTTGGAGTCATTACCGACATCAAAAg                          <  2:758599/99‑1 (MQ=11)
                                                                 tCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCACTAAGTTGGAGTCATTACCGACATTAAAAGTa                        >  2:19906/1‑100 (MQ=255)
                                                                  cGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTaa                       >  1:2010482/1‑100 (MQ=40)
                                                                  cGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTaa                       >  1:2873446/1‑100 (MQ=40)
                                                                   gTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAAc                      >  1:1416612/1‑100 (MQ=255)
                                                                   gTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAAc                      <  1:1057375/100‑1 (MQ=255)
                                                                   gTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAAc                      <  1:3662413/100‑1 (MQ=255)
                                                                   gTTCTCAAAATCGGAGGAGCTGCATGACCAAGTCATCGGGCATTATCTGAACATAACACACCATCAATCAGTTAGAGTCATTACCGCCATTAAAAGTAAc                      >  2:806626/1‑100 (MQ=14)
                                                                    ttCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATCAGTTGGAGTCATTACCGACATTAAAAGTAACt                     >  1:647824/1‑100 (MQ=255)
                                                                    ttCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACt                     >  2:3505267/1‑100 (MQ=255)
                                                                    ttCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACt                     <  1:1951309/100‑1 (MQ=255)
                                                                     tcccAAAATCGGTGGAGCTGCATGCCAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAAc                      >  2:1124936/4‑98 (MQ=255)
                                                                      ctcAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGACCATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACttt                   >  1:162808/1‑100 (MQ=255)
                                                                       tcAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACtttt                  >  2:3181157/1‑100 (MQ=255)
                                                                       tcAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACtttt                  >  1:2385197/1‑100 (MQ=255)
                                                                        cAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACttttt                 >  1:3083784/1‑100 (MQ=255)
                                                                        cAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACttttt                 >  1:3543649/1‑100 (MQ=255)
                                                                        cAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACttttt                 >  2:1324123/1‑100 (MQ=255)
                                                                        cAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACttttt                 >  1:1142356/1‑100 (MQ=255)
                                                                         aaaaTCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACtttttt                >  2:2798759/1‑100 (MQ=255)
                                                                         aaaaTCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACtttttt                <  2:2628754/100‑1 (MQ=255)
                                                                         aaaaTCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACtttttt                <  2:1142299/100‑1 (MQ=255)
                                                                         aaaaTCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACtttttt                <  1:1897534/100‑1 (MQ=255)
                                                                          aaaTCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTg               <  1:386230/100‑1 (MQ=255)
                                                                          aaaTCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTg               <  2:918973/100‑1 (MQ=255)
                                                                          aaaTCCGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTg               <  2:987827/100‑1 (MQ=255)
                                                                           aaTCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGa              <  1:2197134/100‑1 (MQ=255)
                                                                             tCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGaaa            <  1:3739278/100‑1 (MQ=255)
                                                                             tCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGaaa            <  1:2100801/100‑1 (MQ=255)
                                                                              cGGTGGAGCTGCCTGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAAt           <  1:3691916/100‑1 (MQ=255)
                                                                              cGGTGGACGTCCTTGCCAAAGTCATCGGGCATTATTTGAACATAAAACACTATCCATAAATTGGAGTCATTACCGACATTCAAAGTAACTTTTTTCaact           >  2:22562/1‑98 (MQ=14)
                                                                               ggTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAAtt          >  1:1214091/1‑100 (MQ=255)
                                                                               ggTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAAtt          >  1:1680966/1‑100 (MQ=255)
                                                                               ggTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAAtt          <  1:1913368/100‑1 (MQ=255)
                                                                                gTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTa         >  1:3072589/1‑100 (MQ=255)
                                                                                 tGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTag        <  2:1540950/100‑1 (MQ=255)
                                                                                  ggAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACtttttttttttttt        >  2:3753309/1‑91 (MQ=255)
                                                                                  ggAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTaga       <  1:2507780/100‑1 (MQ=255)
                                                                                   gAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTagag      <  1:2762056/100‑1 (MQ=255)
                                                                                    aGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAgaga     >  1:1929186/1‑100 (MQ=255)
                                                                                      cTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGata   <  2:1014176/100‑1 (MQ=255)
                                                                                      cTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGata   >  1:1374003/1‑100 (MQ=255)
                                                                                       tGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGataa  >  2:1681709/1‑100 (MQ=255)
                                                                                       |                                                                                                   
AGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAA  >  NC_000913/3584044‑3584230

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: