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 3584153–3583427 1222–1948 189 [1] [93] 96 [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.180) 92/178 NT 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

AGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCA  >  NC_000913/3584101‑3584253
                                                     |                                                                                                   
aGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAt                                                       <  2:2979597/100‑1 (MQ=11)
aGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAt                                                       <  1:2726966/100‑1 (MQ=11)
aGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAt                                                       >  2:3496773/1‑100 (MQ=11)
   cGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGCGTCATTACCGACATTaa                                                    >  1:1417979/1‑100 (MQ=18)
   cGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTaa                                                    >  1:3381458/1‑100 (MQ=33)
   cGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTaa                                                    >  1:960454/1‑100 (MQ=33)
   cGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCCGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTaa                                                    <  1:165066/100‑1 (MQ=18)
    gCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATCAAACACTATCAATAAGTTGGAGTCATTACCGACATTaaa                                                   >  1:3020679/1‑100 (MQ=25)
    gCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTaaa                                                   >  1:978138/1‑100 (MQ=35)
     cTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCTGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTaaa                                                   <  1:3682403/99‑1 (MQ=25)
     cTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTaaaa                                                  <  2:1300701/100‑1 (MQ=35)
     cTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTaaaa                                                  <  2:3619893/100‑1 (MQ=35)
     cTGTCGTTCTCAAAATCGGTGAAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGAAGTCATTACCGACATTaaaa                                                  <  1:840739/100‑1 (MQ=255)
        tCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAg                                                 <  2:1597337/98‑1 (MQ=38)
         cGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTaa                                              <  1:3052297/100‑1 (MQ=40)
         cGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTaa                                              >  1:2675634/1‑100 (MQ=40)
         cGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTaa                                              <  1:1021242/100‑1 (MQ=40)
         cGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTaa                                              >  1:1914608/1‑100 (MQ=40)
          gTTCTCAAAATCGGTGGAGCTGCATGACAAATTCATCTGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAAc                                             <  1:2557664/100‑1 (MQ=255)
          gTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAAc                                             >  2:2979711/1‑100 (MQ=255)
            tctcAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACtt                                           >  1:550057/1‑100 (MQ=255)
            tctcAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACtt                                           >  2:2045606/1‑100 (MQ=255)
             ctcAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACttt                                          >  1:1405538/1‑100 (MQ=255)
              tcAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTAGCGACATTAAAAGTAACtttt                                         <  1:2323376/100‑1 (MQ=255)
              tcAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACtttt                                         >  2:3384922/1‑100 (MQ=255)
               cAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACttttt                                        <  2:3520035/100‑1 (MQ=255)
               cAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACttttt                                        >  2:1182877/1‑100 (MQ=255)
               cAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACttttt                                        >  2:762519/1‑100 (MQ=255)
               cAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACttttt                                        >  2:1452209/1‑100 (MQ=255)
               cAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACttttt                                        >  1:471509/1‑100 (MQ=255)
                aaaaTCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACtttttt                                       <  2:1687003/100‑1 (MQ=255)
                 aaaTCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTg                                      <  2:340010/100‑1 (MQ=255)
                 aaaTCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTg                                      <  2:3793901/100‑1 (MQ=255)
                 aaaTCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTg                                      <  1:727165/100‑1 (MQ=255)
                  aaTCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGa                                     >  2:2534332/1‑100 (MQ=255)
                  aaTCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGa                                     >  2:397315/1‑100 (MQ=255)
                    tCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGaaa                                   >  1:3082288/1‑100 (MQ=255)
                    tCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGaaa                                   <  2:1471005/100‑1 (MQ=255)
                      ggTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTGTCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAAtt                                 >  1:1369508/1‑100 (MQ=255)
                      ggTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATATGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAAt                                  <  1:3093819/99‑1 (MQ=255)
                      ggTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAAtt                                 >  1:663328/1‑100 (MQ=255)
                      ggTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAAtt                                 <  2:788244/100‑1 (MQ=255)
                      ggTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAAtt                                 <  2:2511702/100‑1 (MQ=255)
                      ggTGGAGCCGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTAGCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACtttttt                                       <  1:3476608/94‑1 (MQ=255)
                       gTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTa                                >  2:3455435/1‑100 (MQ=255)
                       gTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTa                                >  1:3707416/1‑100 (MQ=255)
                       gTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTa                                <  1:4025/100‑1 (MQ=255)
                       gTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTa                                >  1:23219/1‑100 (MQ=255)
                        tGGAGCTACATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTag                               <  2:1645865/100‑1 (MQ=255)
                         ggAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTTGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTaga                              >  1:1176191/1‑100 (MQ=255)
                         ggAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTaga                              >  1:3149649/1‑100 (MQ=255)
                         ggAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTaga                              >  1:2910685/1‑100 (MQ=255)
                          gAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTagag                             <  1:1549236/100‑1 (MQ=255)
                          gAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTagag                             <  1:1559835/100‑1 (MQ=255)
                           aGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAgaga                            <  2:1409200/100‑1 (MQ=255)
                           aGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAgaga                            >  2:3048641/1‑100 (MQ=255)
                           aGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAgaga                            >  2:3264641/1‑100 (MQ=255)
                            gCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGat                           <  2:2642735/100‑1 (MQ=255)
                             cTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGata                          <  2:609990/100‑1 (MQ=255)
                             cTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGata                          >  1:1247646/1‑100 (MQ=255)
                              tGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGataa                         >  2:2990307/1‑100 (MQ=255)
                              tGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGataa                         >  2:3085732/1‑100 (MQ=255)
                              tGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGataa                         >  1:1025895/1‑100 (MQ=255)
                               gCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGataat                        >  1:1325957/1‑100 (MQ=255)
                               gCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGataat                        >  2:1435230/1‑100 (MQ=255)
                               gCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTAACAATAAGTTGGAGACATTACCGACATAGACAGAAACTTTATAGaaagtgg                               >  1:3238660/1‑89 (MQ=11)
                                cATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGataata                       >  1:1680913/1‑100 (MQ=255)
                                 aTGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAatat                      <  2:3418228/100‑1 (MQ=255)
                                  tGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATAtt                     >  2:2091958/1‑100 (MQ=255)
                                   gACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTg                    <  2:3675511/100‑1 (MQ=255)
                                    aCAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTgg                   >  1:3087644/1‑100 (MQ=255)
                                     cAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAatat                      >  1:2106875/1‑96 (MQ=255)
                                     cAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTggn                  <  2:1120739/100‑2 (MQ=255)
                                     cAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGc                  <  1:2038523/100‑1 (MQ=255)
                                     cAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGc                  >  2:3641269/1‑100 (MQ=255)
                                     cAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGc                  <  1:2138085/100‑1 (MQ=255)
                                     cAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGc                  >  1:3454249/1‑100 (MQ=255)
                                      aaaGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGca                 <  2:373346/100‑1 (MQ=255)
                                      aaaGTCATCGGGCATTATCTGAACAGAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGca                 <  1:1297308/100‑1 (MQ=255)
                                        aGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGcaca               <  2:456466/100‑1 (MQ=255)
                                              cGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAat         <  1:276180/100‑1 (MQ=255)
                                               gggCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAata        >  2:809363/1‑100 (MQ=255)
                                               gggCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAata        <  2:323598/100‑1 (MQ=255)
                                               gggCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAata        >  1:683283/1‑100 (MQ=255)
                                               gggCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAata        <  2:1452561/100‑1 (MQ=255)
                                               gggCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAat         <  1:964748/99‑1 (MQ=255)
                                               gggCATTATCTGAACATAAAACACTAACAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAat         <  1:396435/99‑1 (MQ=255)
                                                ggCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTTTTTTTTTTTTTAATATTGGCACAGAAAatgt       >  1:2292769/1‑98 (MQ=17)
                                                 gCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAatata      <  2:1194409/100‑1 (MQ=255)
                                                 gCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAatat       <  1:680362/99‑1 (MQ=255)
                                                  cATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAatatat     >  2:2487403/1‑100 (MQ=255)
                                                   aTTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATAtt    <  2:1836799/100‑1 (MQ=255)
                                                    ttATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTc   <  1:931302/100‑1 (MQ=255)
                                                     tATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCa  >  1:3454845/1‑100 (MQ=255)
                                                     tATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCa  >  2:3257817/1‑100 (MQ=255)
                                                     tATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCa  <  2:1859755/100‑1 (MQ=255)
                                                     |                                                                                                   
AGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCA  >  NC_000913/3584101‑3584253

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