Predicted mutation
evidence seq id position mutation annotation gene description
RA AC_000091 4,585,134 (A)7→8 coding (1009/1395 nt) hsdS ← specificity determinant for hsdM and hsdR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AC_0000914,585,1271.A100.0% 476.3 / NA 116S339L (TCA→TTA) hsdSspecificity determinant for hsdM and hsdR
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base A (48/68);  total (48/68)

Only 100 of 146 total aligned reads displayed.

CTTGGGATTTGATATCTTTTCCTGAAATACCTTTTTGACCAGAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATG‑AAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGAT  >  AC_000091/4585038‑4585211
                                                                                          |                                                                                    
cTTGGGATTTGATATCTTTTCCTGAAATACCTTTTTGACCAGAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATg                                                                                       <  1:144486/90‑1 (MQ=255)
cTTGGGATTTGATATCTTTTCCTGAAATACCTTTTTGACCAGAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATg                                                                                       <  2:3158608/90‑1 (MQ=255)
 ttGGGATTTGATATCTTTTCCTGAAATACCTTTTTGACCAGAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATG‑a                                                                                     >  1:2640248/1‑90 (MQ=255)
 ttGGGATTTGATATCTTTTCCTGAAATACCTTTTTGACCAGAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATG‑a                                                                                     >  1:627540/1‑90 (MQ=255)
   gggATTTGATATCTTTTCCTGAAATACCTTTTTGACCAGAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATG‑aaa                                                                                   <  1:2012001/90‑1 (MQ=255)
    ggATTTGATATCTTTTCCTGAAATACCTTTTTGACCAGAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATG‑aaaa                                                                                  <  1:2968421/90‑1 (MQ=255)
    ggATTTGATATCTTTTCCTGAAATACCTTTTTGACCAGAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATG‑aaaa                                                                                  <  1:814028/90‑1 (MQ=255)
    ggATTTGATATCTTTTCCTGAAATACCTTTTTGACCAGAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATG‑aaaa                                                                                  <  1:2686074/90‑1 (MQ=255)
    ggATTTGATATCTTTTCCTGAAATACCTTTTTGACCAGAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATG‑aaaa                                                                                  >  2:1704950/1‑90 (MQ=255)
     gATTTGATATCTTTTCCTGAAATACCTTTTTGACCAGAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATG‑aaaaa                                                                                 <  2:2175386/90‑1 (MQ=255)
     gATTTGATATCTTTTCCTGAAATACCTTTTTGACCAGAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATG‑aaaaa                                                                                 <  1:2305614/90‑1 (MQ=255)
      aTTTGATATCTTTTCCTGAAATACCTTTTTGACCAGAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATG‑aaaaaa                                                                                <  1:314076/90‑1 (MQ=255)
      aTTTGATATCTTTTCCTGAAATACCTTTTTGACCAGAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATG‑aaaaaa                                                                                >  1:3256327/1‑90 (MQ=255)
       tttGATATCTTTTCCTGAAATACCTTTTTGACCAGAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATG‑aaaaaaa                                                                               >  2:1118205/1‑90 (MQ=255)
        ttGATATCTTTTCCTGAAATACCTTTTTGACCAGAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATG‑aaaaaaaa                                                                              >  1:1115091/1‑89 (MQ=255)
        ttGATATCTTTTCCTGAAATACCTTTTTGACCAGAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATG‑aaaaaaaa                                                                              <  2:2462255/90‑2 (MQ=255)
        ttGATATCTTTTCCTGAAATACCTTTTTGACCAGAAGTTGTCTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATG‑aaaaaaaa                                                                              >  1:464791/1‑89 (MQ=255)
         tGATATCTTTTCCTGAAATACCTTTTTGACCAGAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATG‑aaaaaaaat                                                                             <  1:1214449/90‑3 (MQ=255)
         tGATATCTTTTCCTGAAATACCTTTTTGACCAGAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATG‑aaaaaaaat                                                                             >  2:1661130/1‑88 (MQ=255)
         tGATATCTTTTCCTGAAATACCTTTTTGACCAGAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATG‑aaaaaaaat                                                                             >  2:3242774/1‑88 (MQ=255)
          gATATCTTTTCCTGAAATACCTTTTTGACCAGAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATG‑aaaaaaaata                                                                            <  1:471989/90‑4 (MQ=255)
             atCTTTTCCTGAAATACCTTTTTGACCAGAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAAATAttt                                                                          <  1:599714/90‑1 (MQ=255)
             atCTTTTCCTGAAATACCTTTTTGACCAGAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAAATAttt                                                                          >  2:18780/1‑90 (MQ=255)
              tCTTTTCCTGAAATACCTTTTTGACCAGAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAAATATTTc                                                                         >  2:2483426/1‑90 (MQ=255)
                 tttCCTGAAATACCTTTTTGACCAGAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAAATATTTCGat                                                                      <  2:2230864/90‑1 (MQ=255)
                 tttCCTGAAATACCTTTTTGACCAGAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAAATATTTCGat                                                                      >  2:2297104/1‑90 (MQ=255)
                    ccTGAAATACCTTTTTGACCAGAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAAATATTTCGatata                                                                   <  2:1176316/90‑1 (MQ=255)
                     cTGAAATACCTTTTTGACCAGAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAAATATTTCGatatat                                                                  <  2:451013/90‑1 (MQ=255)
                      tGAAATACCTTTTTGACCAGAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAAATATTTCGATATAtt                                                                 <  1:1823073/90‑1 (MQ=255)
                      tGAAATACCTTTTTGACCAGAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAAATATTTCGATATAtt                                                                 >  2:3057504/1‑90 (MQ=255)
                         aaTACCTTTTTGACCAGAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAAATATTTCGATATATTCTg                                                              <  1:2990173/90‑1 (MQ=255)
                         aaTACCTTTTTGACCAGAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAAATATTTCGATATATTCTg                                                              <  1:3084790/90‑1 (MQ=255)
                          aTACCTTTTTGACCAGAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAAATATTTCGATATATTCTgg                                                             >  1:2906393/1‑90 (MQ=255)
                           tACCTTTTTGACCAGAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAAATATTTCGATATATTCTGGt                                                            >  2:1928489/1‑90 (MQ=255)
                           tACCTTTTTGACCAGAAGTTGGTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAAATATTTCGATATATTCTGGt                                                            <  1:2050328/90‑1 (MQ=255)
                             ccTTTTTGACCAGAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAAATATTTCGATATATTCTGGTaa                                                          <  2:2630563/90‑1 (MQ=255)
                             ccTTTTTGACCAGAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAAATATTTCGATATATTCTGGTaa                                                          <  2:1243445/90‑1 (MQ=255)
                              cTTTTTGACCAGAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAAATATTTCGATATATTCTGGTaaa                                                         >  2:1291709/1‑90 (MQ=255)
                                ttttGGCCAGAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAAATATTTCGATATATTCTGGTAAAGc                                                       <  2:549544/90‑1 (MQ=255)
                                  ttGACCAGAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAAATATTTCGATATATTCTGGTAAAGCAt                                                     <  1:677476/90‑1 (MQ=255)
                                  ttGACCAGAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAAATATTTCGATATATTCTGGTAAAGCAt                                                     <  1:2533055/90‑1 (MQ=255)
                                    gACCAGAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAAATATTTCGATATATTCTGGTAAAGCATCt                                                   >  1:2910644/1‑90 (MQ=255)
                                    gACCAGAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAAATATTTCGATATATTCTGGTAAAGCATCt                                                   >  1:1175247/1‑90 (MQ=255)
                                    gACCAGAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAAATATTTCGATATATTCTGGTAAAGCATCt                                                   >  2:1338341/1‑90 (MQ=255)
                                     aCCAGAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAAATATTTCGATATATTCTGGTAAAGCATCtt                                                  <  1:1529687/90‑1 (MQ=255)
                                     aCCAGAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAAATATTTCGATATATTCTGGTAAAGCATCtt                                                  >  1:1869524/1‑90 (MQ=255)
                                      ccAGAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAAATATTTCGATATATTCTGGTAAAGCATCttt                                                 <  2:254506/90‑1 (MQ=255)
                                      ccAGAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAAATATTTCGATATATTCTGGTAAAGCATCttt                                                 <  1:1984986/90‑1 (MQ=255)
                                      ccAGAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAAATATTTCGATATATTCTGGTAAAGCATCttt                                                 <  1:3471033/90‑1 (MQ=255)
                                        aGAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTgg                                               <  1:2783812/90‑1 (MQ=255)
                                         gAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGt                                              >  2:2581365/1‑90 (MQ=255)
                                         gAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGt                                              <  2:2970410/90‑1 (MQ=255)
                                         gAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGt                                              <  2:1461306/90‑1 (MQ=255)
                                          aaGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGtt                                             <  1:3232399/90‑1 (MQ=255)
                                           aGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTa                                            <  2:532171/90‑1 (MQ=255)
                                            gttgttTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTaa                                           <  2:1215638/90‑1 (MQ=255)
                                             ttgttTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAAt                                          <  2:1327342/90‑1 (MQ=255)
                                             ttgttTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAAt                                          >  1:2545967/1‑90 (MQ=255)
                                              tgttTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATc                                         <  1:107580/90‑1 (MQ=255)
                                              tgttTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATc                                         >  1:1038684/1‑90 (MQ=255)
                                              tgttTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATc                                         <  1:1954151/90‑1 (MQ=255)
                                                ttttCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGa                                       >  2:3271017/1‑90 (MQ=255)
                                                ttttCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGa                                       <  2:1439695/90‑1 (MQ=255)
                                                ttttCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGa                                       >  1:3365602/1‑90 (MQ=255)
                                                 tttCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAg                                      <  2:647768/90‑1 (MQ=255)
                                                  ttCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGc                                     >  2:2340910/1‑90 (MQ=255)
                                                    cACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTc                                   <  1:281123/90‑1 (MQ=255)
                                                      cGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGa                                 <  1:579535/90‑1 (MQ=255)
                                                      cGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGa                                 >  1:1719911/1‑90 (MQ=255)
                                                       gCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGaa                                <  1:1412573/90‑1 (MQ=255)
                                                        cAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAAt                               <  1:3436483/90‑1 (MQ=255)
                                                         aGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATa                              >  2:3367897/1‑90 (MQ=255)
                                                           ttCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAg                            >  2:2247433/1‑90 (MQ=255)
                                                            tcatcaTTGCATTTCGTGCTGAGGGGGATGAAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGt                           <  1:2372187/90‑1 (MQ=255)
                                                             catcatTGCATTTCGTGCTGAGGGGGATGAAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGtt                          >  2:851176/1‑90 (MQ=255)
                                                               tcatTGCATTTCGTGCTGAGGGGGATGAAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTa                        <  1:3380020/90‑1 (MQ=255)
                                                                  ttGCATTTCGTGCTGAGGGGGATGAAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCa                     >  1:354420/1‑90 (MQ=255)
                                                                    gCATTTCGTGCTGAGGGGGATGAAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAgg                   <  2:3296049/90‑1 (MQ=255)
                                                                    gCATTTCGTGCTGAGGGGGATGAAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAgg                   >  1:1338829/1‑90 (MQ=255)
                                                                     cATTTCGTGCTGAGGGGGATGAAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGa                  <  1:2927031/90‑1 (MQ=255)
                                                                        ttCGTGCTGAGGGGGATGAAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGatat               >  1:1360509/1‑90 (MQ=255)
                                                                        ttCGTGCTGAGGGGGATGAAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGatat               <  2:3466271/90‑1 (MQ=255)
                                                                         tCGTGCTGAGGGGGATGAAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGAtata              <  1:3263729/90‑1 (MQ=255)
                                                                         tCGTGCTGAGGGGGATGAAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGAtata              <  1:1618714/90‑1 (MQ=255)
                                                                         tCGTGCTGAGGGGGATGAAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGAtata              >  1:2905654/1‑90 (MQ=255)
                                                                          cGTGCTGAGGGGGATGAAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAg             <  2:2162015/90‑1 (MQ=255)
                                                                           gTGCTGAGGGGGATGAAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGc            >  1:1900724/1‑90 (MQ=255)
                                                                           gTGCTGAGGGGGATGAAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGc            <  2:2184443/90‑1 (MQ=255)
                                                                            tGCTGAGGGGGATGAAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCa           >  1:1812576/1‑90 (MQ=255)
                                                                             gCTGAGGGGGATGAAAAAAAATATGTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCaa          <  1:3430467/90‑1 (MQ=255)
                                                                              cTGAGGGGGATGAAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCaaa         >  2:3285749/1‑90 (MQ=255)
                                                                               tGAGGGGGATGAAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAAt        >  1:1537158/1‑90 (MQ=255)
                                                                               tGAGGGGGATGAAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAAt        <  1:681867/90‑1 (MQ=255)
                                                                               tGAGGGGGATGAAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAAt        >  2:3379115/1‑90 (MQ=255)
                                                                                gAGGGGGATGAAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAAtt       <  1:1802001/90‑1 (MQ=255)
                                                                                 aGGGGGATGAAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAAttt      <  1:3465170/90‑1 (MQ=255)
                                                                                 aGGGGGATGAAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAAttt      >  2:2054458/1‑90 (MQ=255)
                                                                                  gggggATGAAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAAtttt     <  1:234803/90‑1 (MQ=255)
                                                                                   ggggATGAAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTg    <  1:987433/90‑1 (MQ=255)
                                                                                     ggATGAAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGAt  <  1:38734/90‑1 (MQ=255)
                                                                                          |                                                                                    
CTTGGGATTTGATATCTTTTCCTGAAATACCTTTTTGACCAGAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATG‑AAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTTTATCAGGATATAGCAAATTTTGAT  >  AC_000091/4585038‑4585211

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