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% 884.8 / NA 271L423V (CTC→GTC) gppAguanosine pentaphosphatase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (120/151);  total (120/151)
*NC_0009133,960,5790CG100.0% 881.7 / NA 271R422R (CGG→CGCgppAguanosine pentaphosphatase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (120/151);  total (120/151)

Only 100 of 278 total aligned reads displayed.

GTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGAGCCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCAGCAATTTTTT  >  NC_000913/3960491‑3960665
                                                                                       ||                                                                                      
gTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGAcgcgcg                                                                                    <  1:1814440/93‑1 (MQ=255)
 tAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCgcgg                                                                                   <  2:5703800/93‑1 (MQ=255)
 tAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCgcgg                                                                                   >  1:2434072/1‑93 (MQ=255)
  aaTGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCgcggc                                                                                  >  1:6292250/1‑93 (MQ=255)
   aTGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCgcggcg                                                                                 >  2:4062385/1‑93 (MQ=255)
   aTGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACAGCAGAGCTGTTCTGCGACGCggggcg                                                                                 >  1:494131/1‑93 (MQ=255)
     gTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGt                                                                               <  1:5745992/93‑1 (MQ=255)
      tCATCTCTGGCACGAGATCGTCACGGCGACGACTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTa                                                                              <  1:2876168/93‑1 (MQ=255)
       cATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTAc                                                                             <  1:4469939/93‑1 (MQ=255)
         tctcTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACgg                                                                           >  1:2765747/1‑93 (MQ=255)
          ctctGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGc                                                                          <  1:52800/93‑1 (MQ=255)
           tctGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCa                                                                         >  2:4188761/1‑93 (MQ=255)
           tctGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCa                                                                         <  2:5230863/93‑1 (MQ=255)
              ggCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCAttt                                                                      >  2:2968384/1‑93 (MQ=255)
               gCACGAGATCGTGACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCAtttt                                                                     <  2:390731/93‑1 (MQ=255)
                 aCGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGc                                                                   <  2:4872603/93‑1 (MQ=255)
                  cGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCt                                                                  <  1:3085673/93‑1 (MQ=255)
                   gagaTCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTg                                                                 <  1:824155/93‑1 (MQ=255)
                   gagaTCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTg                                                                 <  1:1041871/93‑1 (MQ=255)
                    agaTCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGa                                                                >  2:1481949/1‑93 (MQ=255)
                    agaTCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGa                                                                >  2:4023617/1‑93 (MQ=255)
                    agaTCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGa                                                                >  2:546395/1‑93 (MQ=255)
                    agaTCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGa                                                                >  2:5838012/1‑93 (MQ=255)
                    agaTCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGa                                                                <  2:5626464/93‑1 (MQ=255)
                     gaTCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGAt                                                               <  1:5623585/93‑1 (MQ=255)
                      aTCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATg                                                              >  2:6423582/1‑93 (MQ=255)
                      aTCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATg                                                              <  1:6712331/93‑1 (MQ=255)
                          tCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGc                                                          <  1:743836/93‑1 (MQ=255)
                          tCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGc                                                          <  1:6054279/93‑1 (MQ=255)
                          tCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGc                                                          >  2:5895804/1‑93 (MQ=255)
                           cACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCg                                                         >  1:2197429/1‑93 (MQ=255)
                            aCGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGa                                                        >  2:331581/1‑93 (MQ=255)
                            aCGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGa                                                        >  2:2799491/1‑93 (MQ=255)
                              ggCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATg                                                      <  2:172783/93‑1 (MQ=255)
                              ggCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATg                                                      >  2:3178798/1‑93 (MQ=255)
                              ggCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATg                                                      <  1:3596650/93‑1 (MQ=255)
                                 gACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATgaga                                                   <  2:3097452/93‑1 (MQ=255)
                                 gACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATgaga                                                   >  1:5474325/1‑93 (MQ=255)
                                 gACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATgaga                                                   <  2:6045865/93‑1 (MQ=255)
                                 gACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATgaga                                                   <  2:2553258/93‑1 (MQ=255)
                                  aCGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATgagag                                                  >  2:5613469/1‑93 (MQ=255)
                                  aCGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATgagag                                                  <  2:699297/93‑1 (MQ=255)
                                   cGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATgagaga                                                 <  1:51077/93‑1 (MQ=255)
                                    ggCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGAt                                                >  1:1078583/1‑93 (MQ=255)
                                    ggCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGAt                                                >  2:2692304/1‑93 (MQ=255)
                                     gCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGAtg                                               >  1:3350139/1‑92 (MQ=255)
                                     gCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATc                                               >  2:1978927/1‑93 (MQ=255)
                                       tGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGa                                             <  1:2419690/93‑1 (MQ=255)
                                              aaTGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGAtt                                      <  2:4310667/93‑1 (MQ=255)
                                               atgatgGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTa                                     >  2:3816731/1‑93 (MQ=255)
                                               atgatgGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTa                                     >  2:6266941/1‑93 (MQ=255)
                                               atgatgGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCAATGAGAGATCGTCCGGATTa                                     <  1:6572645/93‑1 (MQ=255)
                                                tgatgGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAg                                    >  1:3694517/1‑93 (MQ=255)
                                                tgatgGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAg                                    <  2:2057172/93‑1 (MQ=255)
                                                  atgGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTc                                  <  1:2521335/93‑1 (MQ=255)
                                                   tgGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCt                                 >  1:2452098/1‑93 (MQ=255)
                                                      ccAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGt                              <  1:2267450/93‑1 (MQ=255)
                                                      ccAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGt                              >  2:5495035/1‑93 (MQ=255)
                                                       cAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGtt                             <  2:5919680/93‑1 (MQ=255)
                                                       cAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGtt                             <  1:616392/93‑1 (MQ=255)
                                                        aGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTg                            >  2:2526223/1‑93 (MQ=255)
                                                         ggCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGa                           <  2:1154123/93‑1 (MQ=255)
                                                         ggCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGa                           <  2:5993322/93‑1 (MQ=255)
                                                         ggCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGa                           >  2:4990486/1‑93 (MQ=255)
                                                          gCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAg                          <  2:3483046/93‑1 (MQ=255)
                                                            gtagtaAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCa                        <  2:4750155/93‑1 (MQ=255)
                                                             tagtaAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAg                       >  1:14468/1‑93 (MQ=255)
                                                             tagtaAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAg                       <  2:4656026/93‑1 (MQ=255)
                                                             tagtaAACGGCAGAGTTGTTCTACGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAg                       >  2:3390201/1‑93 (MQ=255)
                                                              agtaAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGt                      <  1:3257192/93‑1 (MQ=255)
                                                              agtaAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGt                      <  1:6330710/93‑1 (MQ=255)
                                                                taAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAg                    <  2:3672013/93‑1 (MQ=255)
                                                                taAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAg                    <  1:5912363/93‑1 (MQ=255)
                                                                 aaaCGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGc                   <  1:3410666/93‑1 (MQ=255)
                                                                  aaCGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCg                  >  2:353760/1‑93 (MQ=255)
                                                                  aaCGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCg                  >  1:5570480/1‑93 (MQ=255)
                                                                  aaCGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCg                  >  2:2331194/1‑93 (MQ=255)
                                                                  aaCGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCg                  >  2:2207350/1‑93 (MQ=255)
                                                                  aaCGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCg                  >  2:19262/1‑93 (MQ=255)
                                                                   aCGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGt                 <  2:828189/93‑1 (MQ=255)
                                                                    cGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTc                <  1:1480387/93‑1 (MQ=255)
                                                                    cGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTc                >  2:4222378/1‑93 (MQ=255)
                                                                    cGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTc                >  2:4584682/1‑93 (MQ=255)
                                                                    cGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTc                <  2:4429714/93‑1 (MQ=255)
                                                                      gCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGc              <  1:2336093/93‑1 (MQ=255)
                                                                       cAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGcc             <  1:415594/93‑1 (MQ=255)
                                                                        agagTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCa            >  2:1378100/1‑93 (MQ=255)
                                                                          agTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCAGc          >  1:99316/1‑93 (MQ=255)
                                                                           gttgttCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCAGCa         <  1:4092318/93‑1 (MQ=255)
                                                                            tggttCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCAGCaa        <  2:3701295/91‑1 (MQ=255)
                                                                             tgttCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCAGCAAt       <  2:952323/93‑1 (MQ=255)
                                                                             tgttCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCAGCAAt       >  1:3952362/1‑93 (MQ=255)
                                                                              gttCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCAGCAAtt      <  1:6238507/93‑1 (MQ=255)
                                                                              gttCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCAGCAAtt      >  2:6515715/1‑93 (MQ=255)
                                                                               ttCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCAGCAAttt     <  2:5401432/93‑1 (MQ=255)
                                                                                tCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCAGCAAtttt    >  2:3949807/1‑93 (MQ=255)
                                                                                 cTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCAGCAAttttt   <  1:2076738/93‑1 (MQ=255)
                                                                                 cTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCAGCAAttttt   <  1:4807517/93‑1 (MQ=255)
                                                                                  tGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCAGCAAtttttt  <  1:5145058/93‑1 (MQ=255)
                                                                                  tGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCAGCAAtttttt  <  2:5260899/93‑1 (MQ=255)
                                                                                       ||                                                                                      
GTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGAGCCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCAGCAATTTTTT  >  NC_000913/3960491‑3960665

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