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 3584146–3583427 1222–1941 113 [3] [52] 53 [yhhZ]–[insB1] [yhhZ], yrhA, insA, [insB1]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 35822053 (0.030)99 (1.080) 70/178 NT 97.1% 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

Only 100 of 113 total aligned reads displayed.

CAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCT  >  NC_000913/3582106‑3582304
                                                                                                   |                                                                                                   
cAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATtc                                                                                                     >  1:1625280/1‑100 (MQ=255)
  attattTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGTATTACTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCgg                                                                                                   <  2:2247892‑M1/100‑3 (MQ=255)
    tatttatGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtg                                                                                                 <  2:1792961‑M1/100‑5 (MQ=255)
       ttatGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatg                                                                                              >  2:888732‑M1/1‑93 (MQ=255)
        tatGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgc                                                                                             >  1:1589024‑M1/1‑92 (MQ=255)
         atGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgct                                                                                            >  2:2515209‑M1/1‑91 (MQ=255)
          tGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCtgttatgctg                                                                                           >  1:1695402‑M1/1‑90 (MQ=255)
          tGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctg                                                                                           <  1:2175292‑M1/100‑11 (MQ=255)
          tGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctg                                                                                           >  2:1610876‑M1/1‑90 (MQ=255)
             gATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgcca                                                                                        >  1:2335285‑M1/1‑87 (MQ=255)
             gATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgcca                                                                                        >  2:125274‑M1/1‑87 (MQ=255)
              aTTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaa                                                                                       >  1:1326166‑M1/1‑86 (MQ=255)
               tttGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaac                                                                                      >  1:1615874‑M1/1‑85 (MQ=255)
               tttGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaac                                                                                      <  2:1326173‑M1/100‑16 (MQ=255)
                  gACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaactta                                                                                   <  2:215251‑M1/100‑19 (MQ=255)
                  gACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaactta                                                                                   <  2:215346‑M1/100‑19 (MQ=255)
                   aCTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttac                                                                                  <  2:2324973‑M1/100‑20 (MQ=255)
                    cTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTTTTTTTTTTTTggtgatgctgccaacttact                                                                                 >  1:29952‑M1/1‑80 (MQ=255)
                    cTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttact                                                                                 <  2:1526250‑M1/100‑21 (MQ=255)
                    cTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttact                                                                                 <  2:2191775‑M1/100‑21 (MQ=255)
                    cTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttact                                                                                 >  1:1528267‑M1/1‑80 (MQ=255)
                      tCTATCGAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactga                                                                               <  1:2442666‑M1/100‑23 (MQ=255)
                      tCTATCGAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttact                                                                                 >  1:977615‑M1/1‑78 (MQ=255)
                         aTCGAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgattt                                                                            <  1:116393‑M1/100‑26 (MQ=255)
                           cGAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttag                                                                          >  2:2324858‑M1/1‑73 (MQ=255)
                            gAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagt                                                                         >  1:2501505‑M1/1‑72 (MQ=255)
                             aataaataGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtg                                                                        <  2:1615883‑M1/100‑30 (MQ=255)
                             aataaataGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtg                                                                        <  1:2100052‑M1/100‑30 (MQ=255)
                               taaataGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgta                                                                      <  1:543487‑M1/100‑32 (MQ=255)
                                  ataGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatga                                                                   <  1:1041233‑M1/100‑35 (MQ=255)
                                   taGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGATATTCggtgatgatgccaacttacttatttagtgtatgat                                                                  >  1:137725‑M1/1‑65 (MQ=255)
                                      aTTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggt                                                               <  1:970620‑M1/100‑39 (MQ=255)
                                            tAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgttttt                                                         <  2:572112‑M1/100‑45 (MQ=255)
                                              gATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttga                                                       >  2:1446717‑M1/1‑54 (MQ=255)
                                               tctggAAAAGTATTACTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgag                                                      <  1:1942510‑M1/98‑48 (MQ=255)
                                               aTTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTAtt                                                      >  1:2587168/1‑100 (MQ=255)
                                                ttGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgagg                                                     <  2:2074972‑M1/100‑49 (MQ=255)
                                                ttGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgagg                                                     <  2:1470540‑M1/100‑49 (MQ=255)
                                                  ggAAAAGTATTACTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgagaga                                                   >  2:1942520‑M1/1‑50 (MQ=255)
                                                   gAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgc                                                  <  1:125273‑M1/100‑52 (MQ=255)
                                                   gAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgc                                                  >  2:1236346‑M1/1‑49 (MQ=255)
                                                   gAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgc                                                  >  2:1041566‑M1/1‑49 (MQ=255)
                                                      aaGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctcc                                               <  1:690828‑M1/100‑55 (MQ=255)
                                                      aaGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctcc                                               <  1:699456‑M1/100‑55 (MQ=255)
                                                           ttattatATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtgg                                          >  2:1636141‑M1/1‑41 (MQ=255)
                                                           ttattatATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtgg                                          <  2:1665363‑M1/100‑60 (MQ=255)
                                                            tattatATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggc                                         <  2:581882‑M1/100‑61 (MQ=255)
                                                             attatATACAACTCAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggct                                        <  1:363313‑M1/100‑62 (MQ=255)
                                                             attatATACAACTAAGAGGCGCTTTTTTATTTGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggct                                        <  1:990566‑M1/100‑62 (MQ=255)
                                                             attatATACAACTAAGAGGCGCTTTTTTATCGTCTATTCcgtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggct                                        <  1:1816588‑M1/100‑62 (MQ=255)
                                                             attatATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggct                                        <  1:216668‑M1/100‑62 (MQ=255)
                                                             attatATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggct                                        >  2:1258563‑M1/1‑39 (MQ=255)
                                                              ttatATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggctt                                       >  1:952147‑M1/1‑38 (MQ=255)
                                                               tatataCAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttc                                      <  1:18668‑M1/100‑64 (MQ=255)
                                                                 tataCAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctg                                    <  2:534382‑M1/100‑66 (MQ=255)
                                                                  ataCAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgt                                   <  2:1636239‑M1/100‑67 (MQ=255)
                                                                   taCAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtt                                  <  2:267786‑M1/100‑68 (MQ=255)
                                                                   taCAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtt                                  >  2:137442‑M1/1‑33 (MQ=255)
                                                                   taCAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtt                                  >  2:2019712‑M1/1‑33 (MQ=255)
                                                                    aCAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgttt                                 <  1:1138120‑M1/100‑69 (MQ=255)
                                                                     cAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtttatgatggtgtttttgaggtgctccagtggcttctgtttc                                >  1:1232364‑M1/1‑31 (MQ=255)
                                                                     cAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttc                                <  2:758670‑M1/100‑70 (MQ=255)
                                                                     cAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttc                                >  2:2579754‑M1/1‑31 (MQ=255)
                                                                       aCTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttcta                              <  1:561545‑M1/100‑72 (MQ=255)
                                                                       aCTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttcta                              <  2:2104968‑M1/100‑72 (MQ=255)
                                                                       aCTAAGAGGCGCTTTTTTATCGGCTATTCggggatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttct                               >  2:2555613‑M1/1‑29 (MQ=255)
                                                                        cTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctat                             >  1:1150613‑M1/1‑28 (MQ=255)
                                                                        cTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctat                             >  1:849456‑M1/1‑28 (MQ=255)
                                                                         tAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatc                            <  1:263210‑M1/100‑74 (MQ=255)
                                                                          aaGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatca                           <  1:1610867‑M1/100‑75 (MQ=255)
                                                                          aaGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatca                           >  1:818667‑M1/1‑26 (MQ=255)
                                                                           agagGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcag                          <  1:403362‑M1/100‑76 (MQ=255)
                                                                           agagGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcag                          >  2:1103057‑M1/1‑25 (MQ=255)
                                                                            gagGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagc                         <  1:1103053‑M1/100‑77 (MQ=255)
                                                                            gagGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagc                         <  1:2497872‑M1/100‑77 (MQ=255)
                                                                             agGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagct                        <  2:1167754‑M1/100‑78 (MQ=255)
                                                                               gcgcTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTAt                      <  2:118765/100‑1 (MQ=255)
                                                                                cgcTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtc                     <  2:618970‑M1/100‑81 (MQ=255)
                                                                                cgcTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtc                     <  2:442254‑M1/100‑81 (MQ=255)
                                                                                  cTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccc                   >  2:204112‑M1/1‑18 (MQ=255)
                                                                                  cTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtcc                    <  1:939618‑M1/99‑82 (MQ=255)
                                                                                   ttttttATCGGCTATTCggtgatgctgccaaattactgatttagtgtatgatggtgcttttgaggtgctccagtggcttctgtttctatccgctgtccct                  >  1:2008418‑M1/1‑17 (MQ=255)
                                                                                     ttttATCGGCTATTTggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcc                <  1:1098569‑M1/100‑86 (MQ=255)
                                                                                      tttATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcct               >  1:2064088‑M1/1‑14 (MQ=255)
                                                                                         aTCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgtt            >  2:1058604‑M1/1‑11 (MQ=255)
                                                                                           cGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttca          >  2:684543‑M1/1‑9 (MQ=255)
                                                                                            ggCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcag         >  2:2103632‑M1/1‑8 (MQ=255)
                                                                                            ggCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcag         >  1:842139‑M1/1‑8 (MQ=255)
                                                                                            ggCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcag         >  2:291053‑M1/1‑8 (MQ=255)
                                                                                            ggCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcag         <  1:701343‑M1/100‑93 (MQ=255)
                                                                                            ggCTATTCggtgatgctgccaacttacagatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcag         >  1:858031‑M1/1‑8 (MQ=255)
                                                                                             gCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagc        <  2:93178‑M1/100‑94 (MQ=255)
                                                                                             gCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagc        <  1:21592‑M1/100‑94 (MQ=255)
                                                                                               tATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagct       <  1:209794‑M1/99‑95 (MQ=255)
                                                                                                aTTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctac     <  2:1976605‑M1/100‑97 (MQ=255)
                                                                                                aTTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctac     <  2:1540488‑M1/100‑97 (MQ=255)
                                                                                                 ttCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctact    <  2:417464‑M1/100‑98 (MQ=255)
                                                                                                 tgcggtgatgctgccaacttacggatttattgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctact    >  1:381110‑M1/3‑3 (MQ=255)
                                                                                                  tCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactg   >  1:820752‑M1/1‑2 (MQ=255)
                                                                                                   cggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactga  >  1:2215270‑M1/1‑1 (MQ=255)
                                                                                                   |                                                                                                   
CAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCT  >  NC_000913/3582106‑3582304

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: