Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,876,347 A→G intergenic (‑18/‑234) cysD ← / → iap sulfate adenylyltransferase subunit 2/alkaline phosphatase isozyme conversion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,876,3470AG100.0% 1117.0 / NA 268intergenic (‑18/‑234)cysD/iapsulfate adenylyltransferase subunit 2/alkaline phosphatase isozyme conversion protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (128/140);  total (128/140)

Only 100 of 268 total aligned reads displayed.

GAATTCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACG  >  NC_000913/2876249‑2876444
                                                                                                  |                                                                                                 
gaATTCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATGcc                                                                                                 >  2:3718834/1‑101 (MQ=255)
gaATTCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATGcc                                                                                                 >  1:3718835/1‑101 (MQ=255)
    tCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTa                                                                                             <  2:3718321/101‑1 (MQ=255)
        ccgccACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTAtttt                                                                                         <  1:3718512/101‑1 (MQ=255)
         cgccACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTc                                                                                        >  2:3718839/1‑101 (MQ=255)
         cgccACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTc                                                                                        <  2:3718574/101‑1 (MQ=255)
           ccACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCtt                                                                                      >  2:3718841/1‑101 (MQ=255)
           ccACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCtt                                                                                      <  2:3718260/101‑1 (MQ=255)
           ccACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCtt                                                                                      >  1:3718844/1‑101 (MQ=255)
            cACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTg                                                                                     >  2:3718847/1‑101 (MQ=255)
               cTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCa                                                                                  >  1:3718849/1‑101 (MQ=255)
                 cgcgAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATc                                                                                >  1:3718851/1‑101 (MQ=255)
                  gcgAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCa                                                                               >  2:3718857/1‑101 (MQ=255)
                  gcgAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCa                                                                               >  1:3718854/1‑101 (MQ=255)
                     aataatGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGAt                                                                            <  1:3718155/101‑1 (MQ=255)
                     aataatGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGAt                                                                            >  2:3718860/1‑101 (MQ=255)
                     aataatGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGAt                                                                            <  1:3718390/101‑1 (MQ=255)
                      ataatGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATg                                                                           >  1:3718861/1‑101 (MQ=255)
                      ataatGTAGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATg                                                                           <  2:3718420/101‑1 (MQ=255)
                       taatGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGTCATCAGATGt                                                                          >  1:3718863/1‑101 (MQ=255)
                       taatGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGt                                                                          <  2:3718311/101‑1 (MQ=255)
                       taatGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGt                                                                          <  1:3718325/101‑1 (MQ=255)
                        aatGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGtt                                                                         <  2:3718588/101‑1 (MQ=255)
                            tgGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGa                                                                     <  2:3718403/101‑1 (MQ=255)
                              gATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCAACt                                                                   <  1:3718348/101‑1 (MQ=255)
                               aTGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACta                                                                  <  2:3718183/101‑1 (MQ=255)
                                tGCTTTCCGCCTCCAGTTGCCGCAGGTAAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACtat                                                                 <  1:3718494/101‑1 (MQ=255)
                                 gCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACtata                                                                >  2:3718868/1‑101 (MQ=255)
                                 gCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACtata                                                                <  2:3718525/101‑1 (MQ=255)
                                  cTTTCCGCCTCCAGTTGTCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAg                                                               <  2:3718754/101‑1 (MQ=255)
                                  cTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAg                                                               >  1:3718869/1‑101 (MQ=255)
                                    ttCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAggg                                                             <  2:3718474/101‑1 (MQ=255)
                                     tCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGa                                                            <  2:3718481/101‑1 (MQ=255)
                                      ccGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAg                                                           >  1:3718873/1‑101 (MQ=255)
                                       cGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGc                                                          >  1:3718877/1‑101 (MQ=255)
                                       cGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGc                                                          >  2:3718876/1‑101 (MQ=255)
                                        gCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGa                                                            >  2:3718880/1‑98 (MQ=255)
                                         ccTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGt                                                        <  2:3718500/101‑1 (MQ=255)
                                          cTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTa                                                       <  1:3718470/101‑1 (MQ=255)
                                          cTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTa                                                       <  1:3718463/101‑1 (MQ=255)
                                           tCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTaa                                                      <  1:3718590/101‑1 (MQ=255)
                                           tCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTaa                                                      <  2:3718444/101‑1 (MQ=255)
                                            ccAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAag                                                     <  2:3718554/101‑1 (MQ=255)
                                            ccAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAag                                                     >  2:3718885/1‑101 (MQ=255)
                                             cAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAaga                                                    >  1:3718886/1‑101 (MQ=255)
                                                ttGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGaa                                                 >  1:3718890/1‑101 (MQ=255)
                                                 tGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAAc                                                <  2:3718634/101‑1 (MQ=255)
                                                 tGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAAc                                                <  1:3718601/101‑1 (MQ=255)
                                                   ccGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGa                                              <  1:3718519/101‑1 (MQ=255)
                                                    ctcAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGaa                                             >  1:3718893/3‑101 (MQ=255)
                                                    cGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGaa                                             >  2:3718894/1‑101 (MQ=255)
                                                     gCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAAt                                            <  1:3718270/101‑1 (MQ=255)
                                                       aGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGa                                          >  1:3718897/1‑101 (MQ=255)
                                                         gTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGaaa                                        <  1:3718486/101‑1 (MQ=255)
                                                          tGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAAt                                       >  1:3718900/1‑101 (MQ=255)
                                                           gAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAAtt                                      >  1:3718903/1‑101 (MQ=255)
                                                           gAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAAtt                                      <  1:3718376/101‑1 (MQ=255)
                                                               aaGTCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCa                                  >  2:3718906/1‑101 (MQ=255)
                                                               aaGTCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCa                                  <  1:3718413/101‑1 (MQ=255)
                                                               aaGTCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCa                                  >  1:3718907/1‑101 (MQ=255)
                                                                aGTCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCaa                                 <  1:3718668/101‑1 (MQ=255)
                                                                 gTCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAAt                                <  1:3718543/101‑1 (MQ=255)
                                                                  tCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCa                                  >  1:3718911/1‑98 (MQ=255)
                                                                  tCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAAtt                               <  1:3718498/101‑1 (MQ=255)
                                                                   cGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTa                              <  1:3718482/101‑1 (MQ=255)
                                                                   cGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTATTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTa                                     <  2:3718379/94‑1 (MQ=255)
                                                                    gTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAAtt                               <  2:3718598/99‑1 (MQ=255)
                                                                     tATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGa                            >  1:3718918/1‑101 (MQ=255)
                                                                     tATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGa                            >  1:3718915/1‑101 (MQ=255)
                                                                     tATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGa                            >  2:3718916/1‑101 (MQ=255)
                                                                        ttGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATg                         <  1:3718507/101‑1 (MQ=255)
                                                                        ttGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATg                         >  1:3718922/1‑101 (MQ=255)
                                                                          gATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGag                       <  1:3718452/101‑1 (MQ=255)
                                                                           aTCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGagt                      <  1:3718529/101‑1 (MQ=255)
                                                                           aTCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGagt                      >  2:3718925/1‑101 (MQ=255)
                                                                            tCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGagta                     >  2:3718929/1‑101 (MQ=255)
                                                                            tCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGagta                     >  1:3718928/1‑101 (MQ=255)
                                                                            tCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGacta                     <  1:3718660/101‑1 (MQ=255)
                                                                             ccATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGagtag                    <  1:3718556/101‑1 (MQ=255)
                                                                             ccATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGagtag                    >  2:3718932/1‑101 (MQ=255)
                                                                             ccATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGagtag                    >  1:3718934/1‑101 (MQ=255)
                                                                               aTAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGtt                  <  1:3718182/101‑1 (MQ=255)
                                                                               aTAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGtt                  >  1:3718937/1‑101 (MQ=255)
                                                                                 aaCCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTcc                >  1:3718942/1‑101 (MQ=255)
                                                                                  aCCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCt               >  1:3718944/1‑101 (MQ=255)
                                                                                   ccGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCtt              >  2:3718945/1‑101 (MQ=255)
                                                                                    cGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTa             <  1:3718573/101‑1 (MQ=255)
                                                                                     gTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTaa            >  2:3718948/1‑101 (MQ=255)
                                                                                      ttCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAAc           <  2:3718468/101‑1 (MQ=255)
                                                                                      ttCCTTCGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATg                         <  1:3718314/87‑1 (MQ=255)
                                                                                        ccTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACgg         >  2:3718952/1‑101 (MQ=255)
                                                                                         cTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGa        <  1:3718595/101‑1 (MQ=255)
                                                                                         cTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGa        <  1:3718719/101‑1 (MQ=255)
                                                                                         cTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGa        >  1:3718953/1‑101 (MQ=255)
                                                                                            tGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATa     <  2:3718692/101‑1 (MQ=255)
                                                                                            tGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATa     <  1:3718364/101‑1 (MQ=255)
                                                                                            tGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATa     <  1:3718767/101‑1 (MQ=255)
                                                                                             gCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATaa    <  1:3718417/101‑1 (MQ=255)
                                                                                             gCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATaa    >  1:3718958/1‑101 (MQ=255)
                                                                                             gCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATaa    >  2:3718959/1‑101 (MQ=255)
                                                                                                  |                                                                                                 
GAATTCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACG  >  NC_000913/2876249‑2876444

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