Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 3,582,206 Δ1,222 bp IS1‑mediated [yhhZ]yrhA [yhhZ], yrhA

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 3582206 3584102–3583429 1224–1897 69 [0] [20] 25 [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)69 (1.030) 44/298 0.2 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

TAAAAGCACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAATTCAAA  >  NC_000913/3582063‑3582353
                                                                                                                                              |                                                                                                                                                    
aaaaagCACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgc                                                                                                                                              <  1:668427‑M1/150‑9 (MQ=255)
   aaGCACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctg                                                                                                                                            >  1:220702‑M1/1‑140 (MQ=255)
                 ttattaATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgattt                                                                                                                             <  1:976014‑M1/151‑26 (MQ=255)
                 ttattaATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgattt                                                                                                                             <  1:26471‑M1/151‑26 (MQ=255)
                  tattaATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGATTTTTTATCGGCTATTCggtgatgctgccaacttactgattta                                                                                                                            >  1:731005‑M1/1‑125 (MQ=255)
                     taATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtg                                                                                                                         <  2:813340‑M1/151‑30 (MQ=255)
                           cATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgat                                                                                                                   <  1:878544‑M1/151‑36 (MQ=255)
                            aTTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatg                                                                                                                  >  2:285747‑M1/1‑115 (MQ=255)
                              tAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggt                                                                                                                >  2:877212‑M1/1‑113 (MQ=255)
                              tAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggt                                                                                                                >  2:1046379‑M1/1‑113 (MQ=255)
                                 cAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtt                                                                                                             >  1:65474‑M1/1‑110 (MQ=255)
                                   ataatGAACAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttt                                                                                                           <  1:335391‑M1/151‑44 (MQ=255)
                                   ataatGAACAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttt                                                                                                           >  2:514932‑M1/1‑108 (MQ=255)
                                       tGAACAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgag                                                                                                       <  2:333687‑M1/151‑48 (MQ=255)
                                         aaCAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgagg                                                                                                      <  1:674170‑M1/150‑49 (MQ=255)
                                           cAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaaattactgatttagtgtatgatggtgtttttgaggtg                                                                                                    >  2:905757‑M1/1‑100 (MQ=255)
                                             attattTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgct                                                                                                  >  2:171068‑M1/1‑98 (MQ=255)
                                             attattTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgct                                                                                                  >  2:105377‑M1/1‑98 (MQ=255)
                                                atttatGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccag                                                                                              >  1:772755‑M1/1‑95 (MQ=255)
                                                atttatGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccag                                                                                              >  2:700248‑M1/1‑95 (MQ=255)
                                                            tGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgttt                                                                                  <  2:394832‑M1/151‑69 (MQ=255)
                                                             gACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgttttgaggtgctccagtggcttctgttt                                                                                   >  1:648928‑M1/1‑82 (MQ=255)
                                                               cTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTGTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttcta                                                                               <  1:709499‑M1/151‑72 (MQ=255)
                                                                 tCTATCGAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatc                                                                             <  2:498134‑M1/151‑74 (MQ=255)
                                                                  cTATCGAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatca                                                                            >  2:784521‑M1/1‑77 (MQ=255)
                                                                   tATCGAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcag                                                                           <  2:113556‑M1/151‑76 (MQ=255)
                                                                    aTCGAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagc                                                                          <  2:904894‑M1/151‑77 (MQ=255)
                                                                       gAATAAATCGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcag                                                                           >  1:707757‑M1/1‑72 (MQ=255)
                                                                       gAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcag                                                                           >  2:348149‑M1/1‑72 (MQ=255)
                                                                          taaataGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtcc                                                                     >  2:714523‑M1/1‑69 (MQ=255)
                                                                           aaataGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccct                                                                   <  2:70917‑M1/151‑84 (MQ=255)
                                                                                  tttGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttc                                                            <  1:1034656‑M1/151‑91 (MQ=255)
                                                                                                   gTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtttgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccttcctgttcagctactgacggggtgg                                           >  2:489857‑M1/1‑44 (MQ=255)
                                                                                                   gTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtgg                                           >  2:820046‑M1/1‑44 (MQ=255)
                                                                                                   gTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtgg                                           >  1:86332‑M1/1‑44 (MQ=255)
                                                                                                   gTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtgg                                           >  2:1011971‑M1/1‑44 (MQ=255)
                                                                                                   gTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtgg                                           >  1:183130‑M1/1‑44 (MQ=255)
                                                                                                   gTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtgg                                           >  2:634804‑M1/1‑44 (MQ=255)
                                                                                                   gTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtgg                                           >  1:707857‑M1/1‑44 (MQ=255)
                                                                                                   gTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagccactgacggggtgg                                           >  2:477615‑M1/1‑44 (MQ=255)
                                                                                                        attatATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgt                                      >  1:944480‑M1/1‑39 (MQ=255)
                                                                                                        attatATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgt                                      >  2:1016147‑M1/1‑39 (MQ=255)
                                                                                                        attatATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgt                                      >  2:585421‑M1/1‑39 (MQ=255)
                                                                                                        attatATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtggatgatggtgtttttgaggggctccagtggcttctgtttctattagctgttcctcctgttcagctactgacggggtggtgcgt                                      >  2:39743‑M1/1‑39 (MQ=255)
                                                                                                         ttatATACAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaa                                     <  2:649177‑M1/151‑114 (MQ=255)
                                                                                                           atataCAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaac                                   <  1:534465‑M1/151‑116 (MQ=255)
                                                                                                           atataCAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaa                                    >  1:164834‑M1/1‑36 (MQ=255)
                                                                                                           atataCAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaa                                    >  2:777431‑M1/1‑36 (MQ=255)
                                                                                                           atataCAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgcaaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgagggggtggtgcgtaa                                    >  2:449169‑M1/1‑36 (MQ=255)
                                                                                                             ataCAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggggtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacgg                                 >  2:448908‑M1/1‑34 (MQ=255)
                                                                                                                cAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaa                              >  2:178481‑M1/1‑31 (MQ=255)
                                                                                                                cAACTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagccgtccctcctgttcagctactgacggggtggtgcgtaacggcaa                              >  1:886885‑M1/1‑31 (MQ=255)
                                                                                                                  aCTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaa                            >  2:81709‑M1/1‑29 (MQ=255)
                                                                                                                  aCTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaa                            >  1:532126‑M1/1‑29 (MQ=255)
                                                                                                                  aCTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaa                            >  2:39462‑M1/1‑29 (MQ=255)
                                                                                                                  aCTAAGAGGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaa                            >  1:76044‑M1/1‑29 (MQ=255)
                                                                                                                        agGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcac                        <  2:684069‑M1/149‑127 (MQ=255)
                                                                                                                        agGCGCTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcac                        <  2:177284‑M1/149‑127 (MQ=255)
                                                                                                                             cTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccgga                 >  1:216830‑M1/1‑18 (MQ=255)
                                                                                                                             cTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccgga                 >  1:1055747‑M1/1‑18 (MQ=255)
                                                                                                                             cTTTTTTATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccgg                  <  2:9017‑M1/150‑133 (MQ=255)
                                                                                                                                 tttATCGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatc             <  1:370163‑M1/151‑138 (MQ=255)
                                                                                                                                      cGGCTATTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgc        <  2:630618‑M1/151‑143 (MQ=255)
                                                                                                                                           aTTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctat     >  2:230009‑M1/1‑4 (MQ=255)
                                                                                                                                           aTTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctat     >  1:133400‑M1/1‑4 (MQ=255)
                                                                                                                                            ttCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctc  <  1:881536‑M1/151‑149 (MQ=255)
                                                                                                                                            ttCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctc  <  2:76052‑M1/151‑149 (MQ=255)
                                                                                                                                             tCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatct   >  1:71842‑M1/1‑2 (MQ=255)
                                                                                                                                              cggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatct   >  1:824599‑M1/1‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
TAAAAGCACTCCATTATTTATTAATTCCATTAACAATAATGAACAATTATTTATGGGATTTGACTTCTATCGAATAAATAGATTTGGTAGATTGGAAAAGTATTATTATATACAACTAAGAGGCGCTTTTTTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAATTCAAA  >  NC_000913/3582063‑3582353

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