Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,960,578 2 bp→CG coding (1266‑1267/1485 nt) gppA ← guanosine pentaphosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,960,5780GC100.0% 758.1 / NA 241L423V (CTC→GTC) gppAguanosine pentaphosphatase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (111/130);  total (111/130)
*NC_0009133,960,5790CG100.0% 760.9 / NA 241R422R (CGG→CGCgppAguanosine pentaphosphatase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (111/130);  total (111/130)

Only 100 of 257 total aligned reads displayed.

GTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGAGCCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCAGCAATTTTTT  >  NC_000913/3960491‑3960665
                                                                                       ||                                                                                      
gTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGAcgcgcg                                                                                    >  1:2361162/1‑93 (MQ=255)
gTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGAcgcgcg                                                                                    >  1:575125/1‑93 (MQ=255)
gTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGAcgcgcg                                                                                    <  2:1345246/93‑1 (MQ=255)
gTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGAcgcgcg                                                                                    >  2:3666326/1‑93 (MQ=255)
 tAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCgcgg                                                                                   >  2:127702/1‑93 (MQ=255)
  aaTGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCgcggc                                                                                  <  1:1890879/93‑1 (MQ=255)
  aaTGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGATGGTCTGGGACGCgcggc                                                                                  >  2:418696/1‑93 (MQ=255)
   aTGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATCGCCAGGCGTAGTAAACGGCAGAGGTGTTCTGCGACGCgcggcg                                                                                 >  1:2950871/1‑93 (MQ=255)
     gTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGt                                                                               <  2:3327146/93‑1 (MQ=255)
      tCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTa                                                                              <  1:326087/93‑1 (MQ=255)
      tCATCTCTGGCACGAGATAGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTa                                                                              <  2:3893831/93‑1 (MQ=255)
       cATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTAc                                                                             <  2:3587249/93‑1 (MQ=255)
       cATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTAc                                                                             <  1:1811681/93‑1 (MQ=255)
         tctcTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACgg                                                                           >  1:4229344/1‑93 (MQ=255)
         tctcTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACgg                                                                           >  1:3111247/1‑93 (MQ=255)
          ctctGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGc                                                                          >  1:2427941/1‑93 (MQ=255)
           tctGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCa                                                                         >  2:1478544/1‑93 (MQ=255)
            ctGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCAt                                                                        >  1:2899546/1‑93 (MQ=255)
             tGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCAtt                                                                       <  1:930662/93‑1 (MQ=255)
             tGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCAtt                                                                       <  2:1591297/93‑1 (MQ=255)
              ggCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCAttt                                                                      <  2:641829/93‑1 (MQ=255)
                 aCGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTgt                                                                   >  2:246424/1‑92 (MQ=255)
                 aCGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGc                                                                   >  2:4122083/1‑93 (MQ=255)
                 aCGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGc                                                                   >  2:699146/1‑93 (MQ=255)
                  cGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCt                                                                  >  2:4496874/1‑93 (MQ=255)
                  cGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCt                                                                  >  1:3825062/1‑93 (MQ=255)
                    agaTCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGa                                                                >  1:3963540/1‑93 (MQ=255)
                    agaTCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGa                                                                <  2:806572/93‑1 (MQ=255)
                    agaTCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGa                                                                <  1:4500682/93‑1 (MQ=255)
                    agaTCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGa                                                                <  2:2781344/93‑1 (MQ=255)
                      aTCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATg                                                              >  1:34813/1‑93 (MQ=255)
                      aTCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATg                                                              >  2:2228794/1‑93 (MQ=255)
                       tCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGc                                                             <  2:1260338/93‑1 (MQ=255)
                        cGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCa                                                            <  2:1131138/93‑1 (MQ=255)
                        cGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCa                                                            <  2:2256268/93‑1 (MQ=255)
                         gTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAg                                                           >  1:4369739/1‑93 (MQ=255)
                          tCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGc                                                          <  1:4049988/93‑1 (MQ=255)
                          tCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGc                                                          >  2:4344374/1‑93 (MQ=255)
                          tCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGc                                                          <  2:2821394/93‑1 (MQ=255)
                            aCGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGa                                                        <  2:3930228/93‑1 (MQ=255)
                            aCGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGa                                                        <  2:4462187/93‑1 (MQ=255)
                               gCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATga                                                     <  2:3805952/93‑1 (MQ=255)
                                 gACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATgaga                                                   >  2:468652/1‑93 (MQ=255)
                                 gACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATgaga                                                   <  2:4099311/93‑1 (MQ=255)
                                   cGGCTGGCAAAAATGATGGCCAGTCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATgagaga                                                 <  2:2215746/93‑1 (MQ=255)
                                   cGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATgagaga                                                 <  1:493481/93‑1 (MQ=255)
                                   cGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATgagaga                                                 <  1:3405214/93‑1 (MQ=255)
                                    ggCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGAt                                                <  2:1806195/93‑1 (MQ=255)
                                    ggCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGAt                                                <  2:1401532/93‑1 (MQ=255)
                                    ggCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGAt                                                <  1:1071105/93‑1 (MQ=255)
                                     gCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATc                                               >  2:4162828/1‑93 (MQ=255)
                                     gCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATc                                               >  1:1499970/1‑93 (MQ=255)
                                       tGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGa                                             <  1:2066467/93‑1 (MQ=255)
                                         gCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGAcc                                           >  1:1978141/1‑93 (MQ=255)
                                          cAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCg                                          <  2:1867579/93‑1 (MQ=255)
                                             aaaTGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGAt                                       <  1:1535846/93‑1 (MQ=255)
                                                tgatgGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAg                                    <  2:203418/93‑1 (MQ=255)
                                                 gatgGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGt                                   <  2:3940851/93‑1 (MQ=255)
                                                   tgGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCGTTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCt                                 <  2:293808/93‑1 (MQ=255)
                                                    ggCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGGAGCGATGAGAGATCGACCGGATTAGTCTg                                >  1:2280346/1‑93 (MQ=255)
                                                    ggCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTg                                >  1:2535986/1‑93 (MQ=255)
                                                     gCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTgg                               <  1:2185169/93‑1 (MQ=255)
                                                      ccAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGt                              <  2:4041446/93‑1 (MQ=255)
                                                      ccAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGt                              <  1:1163301/93‑1 (MQ=255)
                                                       cAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGtt                             >  2:2396549/1‑93 (MQ=255)
                                                       cAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGtt                             >  1:3828195/1‑93 (MQ=255)
                                                        aGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTg                            >  1:3012643/1‑93 (MQ=255)
                                                         ggCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGa                           <  2:1997046/93‑1 (MQ=255)
                                                         ggCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGa                           <  1:2852294/93‑1 (MQ=255)
                                                         ggCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGa                           <  1:1710174/93‑1 (MQ=255)
                                                         ggCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGa                           <  2:648973/93‑1 (MQ=255)
                                                         ggCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGa                           >  1:4590319/1‑93 (MQ=255)
                                                         ggCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGa                           >  2:2980270/1‑93 (MQ=255)
                                                         ggCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGa                           <  1:1123766/93‑1 (MQ=255)
                                                         ggCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGa                           >  2:3494084/1‑93 (MQ=255)
                                                            gtagtaAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCa                        >  1:3669450/1‑93 (MQ=255)
                                                            gtagtaAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCa                        <  2:3297364/93‑1 (MQ=255)
                                                             tagtaAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAg                       >  2:1175279/1‑93 (MQ=255)
                                                             tagtaAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAg                       >  1:4254926/1‑93 (MQ=255)
                                                               gtaAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTa                     <  1:1768615/93‑1 (MQ=255)
                                                               gtaAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTa                     <  2:1650653/93‑1 (MQ=255)
                                                                taAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAg                    >  1:3883052/1‑93 (MQ=255)
                                                                 aaaCGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGc                   <  1:25919/93‑1 (MQ=255)
                                                                  aaCGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCg                  >  2:942989/1‑93 (MQ=255)
                                                                  aaCGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCg                  >  2:2117452/1‑93 (MQ=255)
                                                                      gCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGc              >  2:565314/1‑93 (MQ=255)
                                                                      gCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGc              >  2:3142219/1‑93 (MQ=255)
                                                                       cAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATTGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGcc             >  2:2153252/1‑93 (MQ=255)
                                                                       cAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGcc             <  2:3350733/93‑1 (MQ=255)
                                                                          agTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCAGc          >  2:265646/1‑93 (MQ=255)
                                                                           gttgttCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCAGCa         >  1:3522379/1‑93 (MQ=255)
                                                                            ttgttCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCAGCaa        <  1:811029/93‑1 (MQ=255)
                                                                            ttgttCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCAGCaa        >  1:4460898/1‑93 (MQ=255)
                                                                             tgttCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCAGCAAt       >  2:1098993/1‑93 (MQ=255)
                                                                             tgttCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCAGCAAt       >  1:2134859/1‑93 (MQ=255)
                                                                              gttCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCAGCAAtt      <  1:71158/93‑1 (MQ=255)
                                                                              gttCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCAGCAAtt      >  1:1589450/1‑93 (MQ=255)
                                                                                tCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCAGCAAtttt    >  1:1008032/1‑93 (MQ=255)
                                                                                tCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCAGCAAtttt    <  1:1246718/93‑1 (MQ=255)
                                                                                  tGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCAGCAAtttttt  >  1:1388653/1‑93 (MQ=255)
                                                                                       ||                                                                                      
GTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGAGCCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCAGCAATTTTTT  >  NC_000913/3960491‑3960665

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