Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,554,088 +C coding (2/624 nt) yhgKJ ← RNA phosphage cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,554,0881.C100.0% 632.4 / NA 188S133W (TCG→TGG) yhgKJRNA phosphage cyclase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (108/80);  total (108/80)

Only 100 of 219 total aligned reads displayed.

CGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGAGGCGGCGCCGACGGGTTATCGGTCGCCGCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTT  >  NC_000913/3553996‑3554174
                                                                                             |                                                                                       
cGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCg                                                                                          >  1:5124222/1‑93 (MQ=255)
cGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCg                                                                                          >  2:4200908/1‑93 (MQ=255)
 ggCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAAACCgc                                                                                         >  1:5583145/1‑92 (MQ=255)
  gCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCgca                                                                                        <  2:883335/93‑3 (MQ=255)
  gCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCgca                                                                                        >  1:1592153/1‑91 (MQ=255)
    aGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCgcagg                                                                                      >  2:1230993/1‑89 (MQ=255)
    aGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCgcagg                                                                                      >  1:1471493/1‑89 (MQ=255)
    aGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCgcagg                                                                                      >  2:2914527/1‑89 (MQ=255)
    aGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCgcacg                                                                                      >  1:5241663/1‑89 (MQ=255)
    aGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATACgcagg                                                                                      >  2:6104824/1‑89 (MQ=255)
      ggTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAggcg                                                                                    >  2:5217645/1‑93 (MQ=255)
        tAAAAACCGTGGCTTAACAGCGTGTTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCTCCCGGCGGATAAAATCCGCAggcggc                                                                                  <  1:420674/93‑1 (MQ=255)
        tAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAggcggc                                                                                  >  2:593379/1‑93 (MQ=255)
          aaaaCCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGgcgc                                                                                <  1:6352418/93‑1 (MQ=255)
            agccGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCg                                                                              <  2:3999951/91‑1 (MQ=255)
            aaCCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCg                                                                              >  1:47820/1‑93 (MQ=255)
              ccGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAc                                                                            >  2:3860898/1‑93 (MQ=255)
              ccGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAc                                                                            <  2:5354649/93‑1 (MQ=255)
               cGTGGCGTAACAGCGTAGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACg                                                                           <  2:1793650/93‑1 (MQ=255)
                 tGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACggg                                                                         >  1:4864423/1‑93 (MQ=255)
                     gTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTAt                                                                     >  1:6094609/1‑93 (MQ=255)
                     gTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTAt                                                                     <  2:581106/93‑1 (MQ=255)
                     gTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGGTAt                                                                     >  1:828130/1‑93 (MQ=255)
                      tAACAGCGTGGTTTGCTGATGAAGTCCTATTTTCGCCAGCAGCGGCTCCAGCGGCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATc                                                                    <  2:2513232/93‑1 (MQ=255)
                       aaCAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCg                                                                   >  1:3733421/1‑93 (MQ=255)
                       aaCAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCg                                                                   >  2:4908006/1‑93 (MQ=255)
                         cAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGt                                                                 <  2:5164713/93‑1 (MQ=255)
                          aGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGtg                                                                >  1:1991857/1‑92 (MQ=255)
                          aGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGtg                                                                >  2:5761513/1‑92 (MQ=255)
                           gCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGtgc                                                               >  2:4343919/1‑91 (MQ=255)
                           gCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGtgc                                                               >  2:661121/1‑91 (MQ=255)
                           gCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGtgc                                                               <  1:5432398/93‑3 (MQ=255)
                           gCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGtgc                                                               <  1:2860963/93‑3 (MQ=255)
                           gCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGtgc                                                               >  1:3898669/1‑91 (MQ=255)
                            cGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGtgcc                                                              <  1:5749945/93‑4 (MQ=255)
                               ggTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGtgccgcc                                                         >  1:2011230/1‑92 (MQ=255)
                               ggTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTggcgcc                                                           >  2:463756/1‑90 (MQ=255)
                                gTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGtgccgccg                                                        >  2:2438865/1‑91 (MQ=255)
                                  ttGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGtgccgccgct                                                      <  1:1096755/93‑5 (MQ=255)
                                      tgatgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACt                                                   >  2:6020763/1‑93 (MQ=255)
                                      tgatgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACt                                                   >  2:3942883/1‑93 (MQ=255)
                                       gatgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACtt                                                  >  1:2759401/1‑93 (MQ=255)
                                        atgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTc                                                 <  2:4993785/93‑1 (MQ=255)
                                         tgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCa                                                >  1:2547920/1‑93 (MQ=255)
                                         tgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGGTATCGGTGCCGCCGCTCACTTCa                                                >  2:213417/1‑93 (MQ=255)
                                           aaTTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAac                                              >  1:3305822/1‑93 (MQ=255)
                                            aTTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAaca                                             >  2:839653/1‑93 (MQ=255)
                                             ttCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAacac                                            >  1:5463871/1‑93 (MQ=255)
                                             ttCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAacac                                            >  2:2817341/1‑93 (MQ=255)
                                               ccTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACAcgc                                          >  2:4603073/1‑93 (MQ=255)
                                                cTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACAcgcg                                         >  1:2797909/1‑93 (MQ=255)
                                                 tATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGa                                        >  1:5922204/1‑93 (MQ=255)
                                                 tATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGa                                        >  2:3653284/1‑93 (MQ=255)
                                                 tATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGa                                        >  2:2855598/1‑93 (MQ=255)
                                                  aTTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCTCCGCTCACTTCAACACGCGaa                                       <  1:5631147/93‑1 (MQ=255)
                                                  aTTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGaa                                       <  1:3823403/93‑1 (MQ=255)
                                                     ttCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGt                                    <  1:2334048/93‑1 (MQ=255)
                                                       cgcgAGCAGCGGCTCCAGCACCCGGCGGATAAAGTCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTcc                                  <  2:1154584/93‑1 (MQ=255)
                                                       cGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTcc                                  >  1:3394840/1‑93 (MQ=255)
                                                       cGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTcc                                  <  2:5842689/93‑1 (MQ=255)
                                                        gCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCa                                 >  1:1951941/1‑93 (MQ=255)
                                                         ccAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCAt                                >  1:1201592/1‑93 (MQ=255)
                                                          cagcagCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATc                               <  1:3520429/93‑1 (MQ=255)
                                                          cagcagCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATc                               <  1:930756/93‑1 (MQ=255)
                                                           agcagcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCg                              >  2:4802773/1‑93 (MQ=255)
                                                           agcaccGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCg                              <  1:4326327/93‑1 (MQ=255)
                                                              agcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCa                           >  1:5259730/1‑93 (MQ=255)
                                                               gcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCaa                          >  2:2783244/1‑93 (MQ=255)
                                                               gcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCaa                          >  2:1881167/1‑93 (MQ=255)
                                                                cGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCaaa                         <  1:4422891/93‑1 (MQ=255)
                                                                 gccTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAAc                        <  2:6161572/91‑1 (MQ=255)
                                                                  gCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTAACGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAAcc                       >  1:3205282/1‑93 (MQ=255)
                                                                   ctgcAGCCCCCGGCAGATAAAAACCGCCTGCGGGGCCGCCTGGTTTACGGGTCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACca                      <  1:2026094/90‑1 (MQ=255)
                                                                    tCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACcac                     >  2:5504474/1‑93 (MQ=255)
                                                                     ccAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACcaca                    <  2:5653024/93‑1 (MQ=255)
                                                                     ccAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACcaca                    >  2:1522341/1‑93 (MQ=255)
                                                                     ccAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACcaca                    >  1:2148228/1‑93 (MQ=255)
                                                                      cAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAg                   >  1:1318593/1‑93 (MQ=255)
                                                                      cAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGATGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAg                   <  1:5226284/93‑1 (MQ=255)
                                                                       aGCCCCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAgc                  <  2:104702/93‑1 (MQ=255)
                                                                       aGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAgc                  >  1:3513141/1‑93 (MQ=255)
                                                                        gCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAgcg                 >  1:3263347/1‑93 (MQ=255)
                                                                         cACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAgcgc                <  2:6393934/93‑1 (MQ=255)
                                                                         cACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAgcgc                <  2:4778157/93‑1 (MQ=255)
                                                                          ccccGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGcgcg               <  2:457501/92‑1 (MQ=255)
                                                                           cccGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCgg              <  1:5838846/93‑1 (MQ=255)
                                                                            ccGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGGGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGgggg             >  2:1949706/1‑93 (MQ=255)
                                                                             cggcggATAAAACCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgc            <  1:3910415/93‑1 (MQ=255)
                                                                              ggcggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCAATTCAACACGCGAAGGTCCATCGGCAAACCACACCGCGGgca           >  2:273109/1‑93 (MQ=255)
                                                                               gcggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGTCGCTCACTTCAACACGCGACAGTCCATCGGCAAACCACATAGCGGgcag          >  1:1527719/1‑93 (MQ=255)
                                                                               gcggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgcag          >  2:2212730/1‑93 (MQ=255)
                                                                                 ggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgcagca        >  2:318246/1‑93 (MQ=255)
                                                                                  gATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCAc       <  2:1966401/93‑1 (MQ=255)
                                                                                     aaaaTCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGt    >  1:39986/1‑93 (MQ=255)
                                                                                     aaaaTCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGt    >  1:5166347/1‑93 (MQ=255)
                                                                                     aaaaTCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGt    <  1:530792/93‑1 (MQ=255)
                                                                                     aaaaTCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGt    <  2:172815/93‑1 (MQ=255)
                                                                                     aaaaTCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGt    <  2:6344931/93‑1 (MQ=255)
                                                                                      aaaTCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGtt   >  2:3540471/1‑93 (MQ=255)
                                                                                      aaaTCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGtt   <  2:2678416/93‑1 (MQ=255)
                                                                                             |                                                                                       
CGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGAGGCGGCGCCGACGGGTTATCGGTCGCCGCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTT  >  NC_000913/3553996‑3554174

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