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% 1939.4 / NA 559L423V (CTC→GTC) gppAguanosine pentaphosphatase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (257/302);  total (257/302)
*NC_0009133,960,5790CG100.0% 1937.4 / NA 560R422R (CGG→CGC) gppAguanosine pentaphosphatase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (257/303);  total (257/303)

Only 100 of 579 total aligned reads displayed.

GTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGAGCCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCAGCAATTTTTT  >  NC_000913/3960491‑3960665
                                                                                       ||                                                                                      
gTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGAcgcgcg                                                                                    >  1:7666294/1‑93 (MQ=255)
   aTGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCgcggcg                                                                                 >  1:7147090/1‑93 (MQ=255)
    tGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTATTAAACGGCAGAGTTGTTCTGCGACGCGcggcgg                                                                                >  2:3483270/1‑93 (MQ=255)
    tGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGcggcgg                                                                                >  1:4902147/1‑93 (MQ=255)
    agTCATCTCTGGCACGAGATCGGCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGcggcgg                                                                                <  1:6647280/92‑1 (MQ=255)
       cATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTAc                                                                             >  1:5862499/1‑93 (MQ=255)
         tctcTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACgg                                                                           >  1:6372483/1‑93 (MQ=255)
         tctcTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACgg                                                                           >  2:3937756/1‑93 (MQ=255)
         tctcTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTAAgg                                                                           <  1:7484076/93‑1 (MQ=255)
         tctcTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGGACgg                                                                           >  1:6176341/1‑93 (MQ=255)
          ctctGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGc                                                                          >  1:6754540/1‑93 (MQ=255)
          ctctGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGc                                                                          >  2:4696966/1‑93 (MQ=255)
           tctGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCa                                                                         <  1:5787245/93‑1 (MQ=255)
            ctGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTAAGGCAt                                                                        >  2:2050200/1‑93 (MQ=255)
             tGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCAtt                                                                       <  1:5062612/93‑1 (MQ=255)
             tGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCAtt                                                                       >  2:6486916/1‑93 (MQ=255)
              ggCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCAttt                                                                      <  2:2932110/93‑1 (MQ=255)
              ggCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCAttt                                                                      >  1:8458135/1‑93 (MQ=255)
                 aCGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGc                                                                   <  2:3358187/93‑1 (MQ=255)
                 aCGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGGACGGCATTTTgt                                                                   >  2:5159778/1‑92 (MQ=255)
                  cGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCt                                                                  >  1:3748494/1‑93 (MQ=255)
                   gagaTCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTg                                                                 <  1:5976612/93‑1 (MQ=255)
                    agaTCGTCACGGCGACGGCTGGCAAAAGTGATGGCCAGGCGTAGTAAACGGAAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGa                                                                <  1:6986931/93‑1 (MQ=255)
                    agaTCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGa                                                                >  1:946999/1‑93 (MQ=255)
                     gaTCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGAt                                                               >  2:5963763/1‑93 (MQ=255)
                       tCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGc                                                             <  2:7577447/93‑1 (MQ=255)
                        cGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCa                                                            <  1:5364890/93‑1 (MQ=255)
                        cGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCa                                                            <  2:2675549/93‑1 (MQ=255)
                         gTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAg                                                           <  1:3864662/93‑1 (MQ=255)
                          tCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGc                                                          >  1:1263072/1‑93 (MQ=255)
                          tCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGc                                                          >  1:3330545/1‑93 (MQ=255)
                          tCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGc                                                          >  2:5612001/1‑93 (MQ=255)
                          tCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGc                                                          >  2:3641313/1‑93 (MQ=255)
                            aCGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGa                                                        <  1:5688204/93‑1 (MQ=255)
                            aCGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGa                                                        <  2:7371451/93‑1 (MQ=255)
                            aCGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGa                                                        <  2:1303328/93‑1 (MQ=255)
                              ggCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATg                                                      <  2:1046165/93‑1 (MQ=255)
                               gCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATga                                                     <  1:2352650/93‑1 (MQ=255)
                               gCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATga                                                     >  2:7976998/1‑93 (MQ=255)
                                cGACGGCCGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGATGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATgag                                                    >  1:7339996/1‑93 (MQ=255)
                                 gACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATgaga                                                   <  2:3763831/93‑1 (MQ=255)
                                 gACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATgaga                                                   >  1:8650467/1‑93 (MQ=255)
                                 gACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATgaga                                                   >  1:7942122/1‑93 (MQ=255)
                                  aCGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATgagag                                                  >  2:6955583/1‑93 (MQ=255)
                                  aCGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATgagag                                                  <  2:545150/93‑1 (MQ=255)
                                   cGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATgagaga                                                 <  2:1198942/93‑1 (MQ=255)
                                    ggCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGAt                                                >  1:1423801/1‑93 (MQ=255)
                                    ggCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGAt                                                <  1:2158839/93‑1 (MQ=255)
                                    ggCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGAt                                                <  2:1699024/93‑1 (MQ=255)
                                     gCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATc                                               <  2:6714935/93‑1 (MQ=255)
                                     gCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATc                                               <  1:2510135/93‑1 (MQ=255)
                                     gCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATc                                               <  1:6283121/93‑1 (MQ=255)
                                       tGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCgg                                             >  2:4566949/1‑92 (MQ=255)
                                       tGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGa                                             >  2:7900191/1‑93 (MQ=255)
                                       tGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGa                                             >  1:5552073/1‑93 (MQ=255)
                                        ggCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGAc                                            >  1:2980163/1‑93 (MQ=255)
                                         gCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGAcc                                           <  1:4350013/93‑1 (MQ=255)
                                          cAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCg                                          <  2:4081728/93‑1 (MQ=255)
                                          cAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCg                                          <  1:269236/93‑1 (MQ=255)
                                           aaaaaTGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCgg                                         <  2:8159334/93‑1 (MQ=255)
                                             aaaTGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGAt                                       >  2:2449006/1‑93 (MQ=255)
                                             aaaTGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGAt                                       >  2:7118780/1‑93 (MQ=255)
                                              aaTGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGAtt                                      <  2:4187232/93‑1 (MQ=255)
                                              aaTGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGAtt                                      <  1:7238965/93‑1 (MQ=255)
                                               atgatgGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTa                                     >  1:4047661/1‑93 (MQ=255)
                                                tgatgGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAg                                    >  1:2414251/1‑93 (MQ=255)
                                                 gatgGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGt                                   >  1:4248319/1‑93 (MQ=255)
                                                   tgGCCAGGCGTAGTAAACGGCAGCGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCt                                 <  2:6271759/93‑1 (MQ=255)
                                                   tgGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCt                                 <  2:2223538/93‑1 (MQ=255)
                                                     gCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTgg                               <  2:5086802/93‑1 (MQ=255)
                                                      ccAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGt                              >  2:4968658/1‑93 (MQ=255)
                                                      ccAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGt                              >  2:8647882/1‑93 (MQ=255)
                                                      ccAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGt                              <  2:2596114/93‑1 (MQ=255)
                                                       cAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGtt                             >  1:4534807/1‑93 (MQ=255)
                                                       cAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGtt                             >  2:6390818/1‑93 (MQ=255)
                                                        aGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTg                            <  1:8363779/93‑1 (MQ=255)
                                                        aGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTg                            >  1:1045217/1‑93 (MQ=255)
                                                        aGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTg                            <  2:8569070/93‑1 (MQ=255)
                                                         ggCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGa                           >  1:3223244/1‑93 (MQ=255)
                                                           cGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGc                         <  1:8747968/93‑1 (MQ=255)
                                                            gtagtaAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCa                        >  1:1599927/1‑93 (MQ=255)
                                                             tagtaAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAg                       >  1:1880661/1‑93 (MQ=255)
                                                             tagtaAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAg                       >  1:4697996/1‑93 (MQ=255)
                                                              agtaAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGt                      <  1:8045831/93‑1 (MQ=255)
                                                               gtaAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTa                     <  2:4270919/93‑1 (MQ=255)
                                                                taAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAg                    >  1:2269323/1‑93 (MQ=255)
                                                                taAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAg                    <  2:7754518/93‑1 (MQ=255)
                                                                taAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAg                    <  2:4788276/93‑1 (MQ=255)
                                                                taAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAg                    <  1:3449276/93‑1 (MQ=255)
                                                                  aaCGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCg                  >  1:7821156/1‑93 (MQ=255)
                                                                  aaCGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCg                  >  2:6562080/1‑93 (MQ=255)
                                                                   aCGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCgc                 >  1:7887550/1‑92 (MQ=255)
                                                                   aCGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGt                 <  1:2042648/93‑1 (MQ=255)
                                                                   aCGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGt                 <  2:4397163/93‑1 (MQ=255)
                                                                    cGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTc                >  2:8351970/1‑93 (MQ=255)
                                                                        agagTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCa            <  2:2757535/93‑1 (MQ=255)
                                                                          agTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCAGc          >  1:6551305/1‑93 (MQ=255)
                                                                          agTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCAGc          >  2:186709/1‑93 (MQ=255)
                                                                                tCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTCGCGTCGCCAGCAAtttt    >  2:3109764/1‑93 (MQ=255)
                                                                                  tGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCAGCAAtttttt  <  1:8189254/93‑1 (MQ=255)
                                                                                       ||                                                                                      
GTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGAGCCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCAGCAATTTTTT  >  NC_000913/3960491‑3960665

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