Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,992,894 A→G T175A (ACG→GCG)  dapF → diaminopimelate epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,992,8940AG100.0% 684.0 / NA 217T175A (ACG→GCG) dapFdiaminopimelate epimerase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (75/142);  total (75/142)

Only 100 of 231 total aligned reads displayed.

GCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGT  >  NC_000913/3992804‑3992983
                                                                                          |                                                                                         
gccgAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCg                                                                                         >  2:4466815/1‑93 (MQ=255)
 ccgAGCAGACAGTCTTATGCGGCGTGGTGTCGATAGAAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATCCCGCGGCGGTAGAAgcgc                                                                                        <  2:3319067/93‑1 (MQ=255)
  cgAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCt                                                                                       >  1:1235243/1‑93 (MQ=255)
  cgAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCt                                                                                       >  2:1549590/1‑93 (MQ=255)
   gAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCtt                                                                                      <  1:1033577/93‑1 (MQ=255)
   gAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCtt                                                                                      >  1:3650699/1‑93 (MQ=255)
   gAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCtt                                                                                      <  2:679841/93‑1 (MQ=255)
      cAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGt                                                                                   <  1:296783/93‑1 (MQ=255)
          cAATCTTATGCGGCGTGGTGTCGATGGTAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTg                                                                               <  1:655578/93‑1 (MQ=255)
           aaTCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGt                                                                              >  2:1291008/1‑93 (MQ=255)
           aaTCTTAGGCGGCGTTGTGTCGATGGGAAATCCGCCTTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGt                                                                              <  1:55036/93‑1 (MQ=255)
             tCTTATGCGGCGTGGTGTCGATGGTAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTc                                                                            >  1:169841/1‑93 (MQ=255)
             tCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTc                                                                            <  2:2406717/93‑1 (MQ=255)
              cTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCt                                                                           <  2:2931941/93‑1 (MQ=255)
               ttATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTg                                                                          <  1:4050459/93‑1 (MQ=255)
                tATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTgg                                                                         >  1:1130058/1‑93 (MQ=255)
                 aTGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGa                                                                        >  1:692960/1‑93 (MQ=255)
                  tGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGaa                                                                       <  2:1745152/93‑1 (MQ=255)
                  tGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGaa                                                                       <  1:3885695/93‑1 (MQ=255)
                  tGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGaa                                                                       <  2:3568564/93‑1 (MQ=255)
                   gcggcgTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGaaa                                                                      <  2:3036826/93‑1 (MQ=255)
                   gcggcgTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGaaa                                                                      <  1:4282405/93‑1 (MQ=255)
                     ggcgTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGc                                                                    >  2:2206543/1‑93 (MQ=255)
                     ggcgTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGc                                                                    <  2:2532985/93‑1 (MQ=255)
                     cgcgTGGTGTCGATGGGAAATCCGCATTGCGTGATGCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGc                                                                    <  2:2026495/92‑1 (MQ=255)
                      gcgTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGcc                                                                   <  1:3043509/93‑1 (MQ=255)
                       cgTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCa                                                                  >  1:954656/1‑93 (MQ=255)
                       cgTGGTGTCGACGTGAAACCGGCATTGCGTGATTCGGGTCGGTGATGTCGATACCGCGGCGGAAGAAGCGCTTGGTCCAGTCCTGGAAAGCCa                                                                  <  2:2324704/93‑1 (MQ=255)
                          tgtgtCGATGGGAAATCCGCATTGCGTGATTCGGGGCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGa                                                               <  1:67660/92‑1 (MQ=255)
                           gtgtCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAg                                                              >  1:3499459/1‑93 (MQ=255)
                              tCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGt                                                           >  2:1356589/1‑93 (MQ=255)
                              tCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGt                                                           <  1:1016340/93‑1 (MQ=255)
                              gcGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGt                                                           <  2:1455295/92‑1 (MQ=255)
                               cGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGtt                                                          >  2:437491/1‑93 (MQ=255)
                                gATGGGAAATCCGCATTGCGTGATTCGTGTCGATGATCTCGATACCGCGGCGGTAGAAGCGCTTGAACCTGTGCTGGAAAGGCACGAGCGttt                                                         <  2:1517458/93‑1 (MQ=255)
                                gATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTTGTCCTGTTCTGGAAAGCCACGAGCGttt                                                         <  1:2681964/93‑1 (MQ=255)
                                  tGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGTAAAGCCACGAGCGTTTTc                                                       <  1:1459610/93‑1 (MQ=255)
                                  tGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTc                                                       <  2:860719/93‑1 (MQ=255)
                                   gggAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTcc                                                      >  1:1586164/1‑93 (MQ=255)
                                    ggAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCg                                                     >  1:2058718/1‑93 (MQ=255)
                                    ggAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCg                                                     <  1:2569903/93‑1 (MQ=255)
                                      aagtCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGa                                                   <  1:1412767/90‑1 (MQ=255)
                                        acccGCATTGGTTGATTCGGGGCGATGATGCCGATACCGCGGCGGTAGAAGCAGTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAgc                                                 <  2:3828605/91‑1 (MQ=255)
                                        aTCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAgc                                                 <  1:4006534/93‑1 (MQ=255)
                                        aTCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAgc                                                 <  2:1725919/93‑1 (MQ=255)
                                         tCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAgtg                                                >  2:123181/1‑91 (MQ=255)
                                         tCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAgcg                                                >  2:2087631/1‑93 (MQ=255)
                                          ccGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAgcgc                                               <  1:4343613/93‑1 (MQ=255)
                                          ccGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAgcgc                                               <  2:715247/93‑1 (MQ=255)
                                           cGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAgcgcg                                              >  1:1626751/1‑93 (MQ=255)
                                            gCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAgcgcgc                                             >  1:4662031/1‑93 (MQ=255)
                                            gCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGACACGAGCGTTTTCCGGAgcgcgc                                             >  1:2838872/1‑93 (MQ=255)
                                              aTTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCa                                           >  1:723582/1‑93 (MQ=255)
                                               ttGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCaa                                          >  1:2209133/1‑93 (MQ=255)
                                                tGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAat                                         <  1:249611/93‑1 (MQ=255)
                                                 gCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAata                                        <  1:3331304/93‑1 (MQ=255)
                                                  gtggATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAatat                                       <  1:496471/90‑1 (MQ=255)
                                                  cGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAatat                                       <  1:2597833/93‑1 (MQ=255)
                                                    tGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCg                                     <  2:648229/93‑1 (MQ=255)
                                                    tGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCg                                     <  2:3776252/93‑1 (MQ=255)
                                                    tGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCg                                     >  2:3264468/1‑93 (MQ=255)
                                                    tGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCg                                     >  2:537892/1‑93 (MQ=255)
                                                      aTTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGt                                   <  1:1983960/93‑1 (MQ=255)
                                                          gggTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTa                               <  2:2554774/92‑1 (MQ=255)
                                                          gggTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTa                               <  2:3955629/92‑1 (MQ=255)
                                                             tCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGc                            <  2:3551017/93‑1 (MQ=255)
                                                              cGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCa                           <  2:2459844/93‑1 (MQ=255)
                                                               gatgatGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCaa                          <  1:1918202/93‑1 (MQ=255)
                                                               gatgatGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCaa                          <  2:4034306/93‑1 (MQ=255)
                                                               gatgatGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCaa                          <  2:3487892/93‑1 (MQ=255)
                                                                atgatgTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGTGCCAATATCGGTTTTATGCAAg                         >  1:3356164/1‑93 (MQ=255)
                                                                atgatgTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAg                         >  2:210932/1‑93 (MQ=255)
                                                                atgatgTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAg                         >  2:4607625/1‑93 (MQ=255)
                                                                atgatgTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAg                         <  2:4216459/93‑1 (MQ=255)
                                                                atgatgTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAg                         <  2:1528214/93‑1 (MQ=255)
                                                                 tgatgTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGt                        <  1:4411322/93‑1 (MQ=255)
                                                                 tgatgTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGt                        <  1:1518904/93‑1 (MQ=255)
                                                                 tgatgTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGt                        <  2:2797243/93‑1 (MQ=255)
                                                                 tgatgTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGt                        <  1:281680/93‑1 (MQ=255)
                                                                  gatgTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTg                       >  2:3896505/1‑93 (MQ=255)
                                                                   atgTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTgg                      <  2:2240430/93‑1 (MQ=255)
                                                                   atgTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTgg                      >  1:1174163/1‑93 (MQ=255)
                                                                    tgTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGt                     >  1:1810353/1‑93 (MQ=255)
                                                                        gATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAg                 <  2:842801/93‑1 (MQ=255)
                                                                        gATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAg                 <  1:4120287/93‑1 (MQ=255)
                                                                           aCCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAgcgc              >  1:2096813/1‑93 (MQ=255)
                                                                            ccGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGcgcg             <  2:1663443/93‑1 (MQ=255)
                                                                            ccGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGcgcg             <  1:85865/93‑1 (MQ=255)
                                                                             cgcgGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGa            <  1:2337856/93‑1 (MQ=255)
                                                                              gcggcgGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAg           >  1:1075135/1‑93 (MQ=255)
                                                                              gcggcgGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAg           >  1:3590263/1‑93 (MQ=255)
                                                                                 gcggTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCat        <  1:2227896/93‑1 (MQ=255)
                                                                                 gcggTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCat        >  2:27829/1‑93 (MQ=255)
                                                                                   ggTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCatat      >  1:4513028/1‑93 (MQ=255)
                                                                                   ggTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCatat      <  2:295650/93‑1 (MQ=255)
                                                                                    gTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATAtt     <  1:3763467/93‑1 (MQ=255)
                                                                                     tAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTc    >  2:807724/1‑93 (MQ=255)
                                                                                      aGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCg   <  2:2728983/93‑1 (MQ=255)
                                                                                      aGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCg   >  2:3249349/1‑93 (MQ=255)
                                                                                      aGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGACAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCg   <  1:3306762/93‑1 (MQ=255)
                                                                                          |                                                                                         
GCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGT  >  NC_000913/3992804‑3992983

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