Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,106,773 A→G V120V (GTA→GTG sbp → periplasmic sulfate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,106,7730AG100.0% 755.5 / NA 233V120V (GTA→GTGsbpperiplasmic sulfate‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (117/116);  total (117/116)

Only 100 of 240 total aligned reads displayed.

ATTGCGGAACGCGGGCGGATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTACGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGAAAAGCTCTGG  >  NC_000913/4106684‑4106862
                                                                                         |                                                                                         
aTTGCGGAACGCGGGCGGATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGt                                                                                        >  2:4902567/1‑93 (MQ=255)
aTTGCGGAACGCGGGCGGATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGt                                                                                        >  1:1319514/1‑93 (MQ=255)
 ttGCGGAACGCGGGCGGATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTa                                                                                       >  2:2169949/1‑93 (MQ=255)
   gCGTAACGCGGGCGGATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAg                                                                                     <  2:5693677/93‑1 (MQ=255)
   gCGGACCGCGGGCGGATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAg                                                                                     <  1:5668994/93‑1 (MQ=255)
    cGGAACGCGGGCGGATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAgg                                                                                    >  2:6010112/1‑93 (MQ=255)
    cGGAACGCGGGCGGATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAgg                                                                                    <  1:2242168/93‑1 (MQ=255)
     ggAACGCGGGCGGATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAggg                                                                                   <  2:33433/93‑1 (MQ=255)
     ggAACGCGGGCGGATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAggg                                                                                   <  1:4709666/93‑1 (MQ=255)
     ggAACGCGGGCGGATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAggg                                                                                   <  1:5085143/93‑1 (MQ=255)
       aaCGCGGGCGGATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGaa                                                                                 >  2:5145611/1‑93 (MQ=255)
       aaCGCGGGCGGATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGaa                                                                                 <  2:6247139/93‑1 (MQ=255)
        aCGCGGGCGGATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGaaa                                                                                >  1:1793633/1‑93 (MQ=255)
          gcgggcggATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATc                                                                              <  2:5453407/93‑1 (MQ=255)
          gcgggcggATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATc                                                                              <  1:2461687/93‑1 (MQ=255)
            gggcggATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCg                                                                            >  1:2730236/1‑93 (MQ=255)
             ggcggATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGa                                                                           <  1:5412657/93‑1 (MQ=255)
              gcggATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCCCCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGaa                                                                          <  1:555511/93‑1 (MQ=255)
                ggATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGc                                                                        <  1:4862438/93‑1 (MQ=255)
                ggATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGc                                                                        >  2:3964684/1‑93 (MQ=255)
                 gATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCa                                                                       >  1:1945083/1‑93 (MQ=255)
                  aTTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAg                                                                      <  1:6111617/93‑1 (MQ=255)
                   ttGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGa                                                                     >  2:1202810/1‑93 (MQ=255)
                      aTAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATcc                                                                  <  1:4126547/93‑1 (MQ=255)
                       tAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCa                                                                 >  1:2999848/1‑93 (MQ=255)
                       tAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCa                                                                 <  2:3535391/93‑1 (MQ=255)
                         aaGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATg                                                               <  2:2364215/93‑1 (MQ=255)
                         aaGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATg                                                               <  1:4503761/93‑1 (MQ=255)
                          agagTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGa                                                              >  1:2534409/1‑93 (MQ=255)
                           gagTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGAc                                                             <  2:5195950/93‑1 (MQ=255)
                           gagTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGAc                                                             <  1:313103/93‑1 (MQ=255)
                             gTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTg                                                           <  2:4979364/93‑1 (MQ=255)
                              tGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTgg                                                          >  1:2137181/1‑93 (MQ=255)
                              tGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTgg                                                          <  2:2691970/93‑1 (MQ=255)
                               ggATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGa                                                         >  2:5215547/1‑93 (MQ=255)
                               ggATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGa                                                         <  1:2778994/93‑1 (MQ=255)
                                  tCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACg                                                      >  1:3638962/1‑93 (MQ=255)
                                     aaCGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATc                                                   <  2:4588038/93‑1 (MQ=255)
                                      aCGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGCTTTCCTGTTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCt                                                  >  1:4440144/1‑93 (MQ=255)
                                       cGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTg                                                 <  1:2603415/93‑1 (MQ=255)
                                        gTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGa                                                <  1:4278487/93‑1 (MQ=255)
                                          cTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGAtt                                              <  1:1251858/93‑1 (MQ=255)
                                           tGCCTGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTa                                             <  2:2568288/93‑1 (MQ=255)
                                           tGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTa                                             <  1:1859857/93‑1 (MQ=255)
                                           tGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTa                                             >  2:5865057/1‑93 (MQ=255)
                                           tGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTa                                             <  2:5523356/93‑1 (MQ=255)
                                           tGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTa                                             <  1:4370518/93‑1 (MQ=255)
                                            gCCGGATAATTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGCAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTaa                                            >  2:4500636/1‑93 (MQ=255)
                                            gCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTaa                                            >  2:4604708/1‑93 (MQ=255)
                                            gCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTaa                                            >  2:5722493/1‑93 (MQ=255)
                                              cGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAAc                                          <  2:3446105/93‑1 (MQ=255)
                                               ggATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAAcc                                         >  1:6182181/1‑93 (MQ=255)
                                               ggATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAAcc                                         >  1:5333480/1‑93 (MQ=255)
                                                gATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCg                                        >  2:3665105/1‑93 (MQ=255)
                                                gATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCg                                        >  2:468528/1‑93 (MQ=255)
                                                gATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCg                                        >  2:6144374/1‑93 (MQ=255)
                                                gATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCg                                        <  1:3278424/93‑1 (MQ=255)
                                                 aTAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCgg                                       >  2:2804134/1‑93 (MQ=255)
                                                 aTAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCgg                                       <  2:5903501/93‑1 (MQ=255)
                                                  tagctcctcACCGTCCACTTCCTCCATTGTTTCCCTGGTGCGTAAGGTAAAGCCAAAGCATACCCATGACTGTAACGATCTGATTAATCCggg                                      <  2:1898336/90‑1 (MQ=255)
                                                  tAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCggg                                      >  1:4602748/1‑93 (MQ=255)
                                                   aaCTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGgt                                     <  1:2834408/93‑1 (MQ=255)
                                                        cGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTc                                >  1:3721542/1‑93 (MQ=255)
                                                        cGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTc                                >  2:1551260/1‑93 (MQ=255)
                                                             cGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGa                           <  2:2209699/93‑1 (MQ=255)
                                                             cGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGa                           >  2:611739/1‑93 (MQ=255)
                                                              gTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGAt                          >  2:116528/1‑93 (MQ=255)
                                                               tACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATc                         <  2:4302422/93‑1 (MQ=255)
                                                               tACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATc                         >  2:1299386/1‑93 (MQ=255)
                                                               tACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGCTTCGGTGATc                         >  2:282772/1‑93 (MQ=255)
                                                                acacTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCa                        <  2:4221336/93‑1 (MQ=255)
                                                                acacTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCa                        >  2:3149047/1‑93 (MQ=255)
                                                                acacTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCa                        <  1:5733071/93‑1 (MQ=255)
                                                                acacTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCa                        <  1:733648/93‑1 (MQ=255)
                                                                 cacTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCAc                       <  2:1467846/93‑1 (MQ=255)
                                                                   cTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGc                     >  1:4966227/1‑93 (MQ=255)
                                                                    ttCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGcc                    <  1:5167460/93‑1 (MQ=255)
                                                                      ccaccaTTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTa                  <  1:1004741/93‑1 (MQ=255)
                                                                       caccaTTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTaa                 >  2:1961067/1‑93 (MQ=255)
                                                                       caccaTTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTaa                 <  1:4661007/93‑1 (MQ=255)
                                                                        accaTTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAAt                <  2:3015585/93‑1 (MQ=255)
                                                                        accaTTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAAt                >  2:212260/1‑93 (MQ=255)
                                                                        accaTTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCAATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAAt                >  1:3380879/1‑93 (MQ=255)
                                                                             tGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGaa           >  2:1779134/1‑93 (MQ=255)
                                                                             tGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGaa           >  1:2006067/1‑93 (MQ=255)
                                                                              gTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGaaa          >  1:6272135/1‑93 (MQ=255)
                                                                              gTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGaaa          <  2:5616950/93‑1 (MQ=255)
                                                                               ttttCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGaaaa         <  1:545565/93‑1 (MQ=255)
                                                                               ttttCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGaaaa         <  1:1202727/93‑1 (MQ=255)
                                                                                 ttCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGAAAAGc       <  2:4057588/93‑1 (MQ=255)
                                                                                 ttCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGAAAAGc       <  2:3761473/93‑1 (MQ=255)
                                                                                 ttCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGAAAAGc       >  1:1633211/1‑93 (MQ=255)
                                                                                  tCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGAAAAGct      >  1:1172710/1‑93 (MQ=255)
                                                                                   ccTGTTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGAAAAGctc     <  1:974834/93‑1 (MQ=255)
                                                                                   ccTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGAAAAGctc     <  2:633922/93‑1 (MQ=255)
                                                                                    cTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGAAAAGctct    <  1:1955476/93‑1 (MQ=255)
                                                                                      ggTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGAAAAGCTCtgg  >  2:3881816/1‑93 (MQ=255)
                                                                                      ggTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGAAAAGCTCtgg  <  2:250460/93‑1 (MQ=255)
                                                                                      ggTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGAAAAGCTCtgg  >  2:1039640/1‑93 (MQ=255)
                                                                                      ggTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGAAAAGCTCtgg  <  1:58389/93‑1 (MQ=255)
                                                                                         |                                                                                         
ATTGCGGAACGCGGGCGGATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTACGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGAAAAGCTCTGG  >  NC_000913/4106684‑4106862

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