Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,370,778 A→G V145A (GTA→GCA)  nanA ← N‑acetylneuraminate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,370,7780AG99.7% 1210.2 / ‑5.3 386V145A (GTA→GCA) nanAN‑acetylneuraminate lyase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (197/188);  total (198/188)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.84e-01

Only 100 of 396 total aligned reads displayed.

TCCATCTGATAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGTTCTTCA  >  NC_000913/3370688‑3370867
                                                                                          |                                                                                         
    tCTGATAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCAc                                                                                     >  2:1190785/1‑93 (MQ=255)
    tCTGATAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCAc                                                                                     <  2:4958379/93‑1 (MQ=255)
    tCTGATAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCAc                                                                                     >  1:2786190/1‑93 (MQ=255)
      tGATAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTc                                                                                   >  2:6495736/1‑93 (MQ=255)
      tGATAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTc                                                                                   >  1:659784/1‑93 (MQ=255)
       gATAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCa                                                                                  <  1:7338722/93‑1 (MQ=255)
       gATAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCa                                                                                  >  1:284845/1‑93 (MQ=255)
        aTAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAg                                                                                 >  1:4089006/1‑93 (MQ=255)
         tAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAgg                                                                                <  2:4877867/93‑1 (MQ=255)
         tAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAgg                                                                                <  2:1071370/93‑1 (MQ=255)
         tAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAgg                                                                                >  1:4414395/1‑93 (MQ=255)
           gagaTCGCCAGAGGTCTGTTTCATCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGc                                                                              <  2:4284258/93‑1 (MQ=255)
           gagaTCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGc                                                                              <  1:5852944/93‑1 (MQ=255)
           gagaTCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGc                                                                              >  2:1334870/1‑93 (MQ=255)
           gagaTCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGc                                                                              >  2:5086889/1‑93 (MQ=255)
             gaTCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTg                                                                            >  1:6808948/1‑93 (MQ=255)
              aTCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTgg                                                                           >  1:3188514/1‑93 (MQ=255)
               tCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAAGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGa                                                                          >  2:2573312/1‑93 (MQ=255)
                cGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGaa                                                                         <  2:5673145/93‑1 (MQ=255)
                cGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGaa                                                                         <  2:4581139/93‑1 (MQ=255)
                 gCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAAt                                                                        <  2:6784399/93‑1 (MQ=255)
                 gCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAAt                                                                        >  2:5507178/1‑93 (MQ=255)
                   cAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAAtgt                                                                      <  1:4656627/93‑1 (MQ=255)
                    agagGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAAtgtt                                                                     <  2:4024797/93‑1 (MQ=255)
                    agagGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAAtgtt                                                                     >  1:7168514/1‑93 (MQ=255)
                    agagGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAAtgtt                                                                     <  2:3057503/93‑1 (MQ=255)
                    agagGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAAtgtt                                                                     >  1:5477855/1‑93 (MQ=255)
                    agagGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAAtgtt                                                                     >  2:5410796/1‑93 (MQ=255)
                     gagGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAAtgttg                                                                    <  2:7359561/93‑1 (MQ=255)
                      agGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAAtgttgt                                                                   >  1:7584366/1‑93 (MQ=255)
                        gTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTac                                                                 >  2:3198112/1‑93 (MQ=255)
                        gTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTac                                                                 <  1:5677067/93‑1 (MQ=255)
                         tCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTaca                                                                >  1:887528/1‑93 (MQ=255)
                         tCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTaca                                                                >  1:3897470/1‑93 (MQ=255)
                         tCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTaca                                                                >  2:166042/1‑93 (MQ=255)
                                cAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCAt                                                         <  2:1953730/93‑1 (MQ=255)
                                 aGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATc                                                        <  2:6007632/93‑1 (MQ=255)
                                 aGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATc                                                        >  2:7184347/1‑93 (MQ=255)
                                    gcACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGc                                                     >  2:5770516/1‑93 (MQ=255)
                                     cACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCa                                                    >  1:4298719/1‑93 (MQ=255)
                                     cACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGCCa                                                    >  1:1360846/1‑93 (MQ=255)
                                       ccTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCACCGGCaaa                                                  <  2:4160821/93‑1 (MQ=255)
                                         tACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAAcc                                                <  1:7714029/93‑1 (MQ=255)
                                         tACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAAcc                                                >  1:3757709/1‑93 (MQ=255)
                                          aCGCCAGGCAATGTAACAAGTGTTTTGATCTGATCCAGGGTCAGGTTTAACCCACTCAGGGCTGGAATTCTGTAAACCACCATCGGCAAAcct                                               >  2:2717041/1‑92 (MQ=255)
                                          aCGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCa                                               <  2:3297063/93‑1 (MQ=255)
                                             ccAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATcc                                            <  1:2150830/93‑1 (MQ=255)
                                              cAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATccg                                           >  1:630488/1‑93 (MQ=255)
                                              cAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATccg                                           >  1:2637686/1‑93 (MQ=255)
                                               aGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATccgc                                          <  1:502836/93‑1 (MQ=255)
                                               aGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACAACCATCGGCAAACCATccgc                                          >  2:6571561/1‑93 (MQ=255)
                                                ggCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATccgcc                                         <  1:1471966/93‑1 (MQ=255)
                                                  cAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGa                                       >  1:483038/1‑93 (MQ=255)
                                                   aaTGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGaa                                      >  2:3904454/1‑93 (MQ=255)
                                                     tGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATc                                    >  1:6058697/1‑93 (MQ=255)
                                                       tAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCaa                                  >  1:7406092/1‑93 (MQ=255)
                                                       tAACAAGTCTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCaa                                  <  1:4344688/93‑1 (MQ=255)
                                                        aaCAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCaat                                 <  2:3610279/93‑1 (MQ=255)
                                                        aaCAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCaat                                 >  1:6915812/1‑93 (MQ=255)
                                                         aCAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCaata                                <  2:4372659/93‑1 (MQ=255)
                                                         aCAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCaata                                <  1:5208059/93‑1 (MQ=255)
                                                         aCAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCaata                                <  1:5989191/93‑1 (MQ=255)
                                                          cAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCaataa                               >  1:2446538/1‑93 (MQ=255)
                                                           aaGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCaataat                              <  2:1276016/93‑1 (MQ=255)
                                                             gtgtgtTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTg                            >  1:1571732/1‑93 (MQ=255)
                                                              tgtgtTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGc                           <  2:2938758/93‑1 (MQ=255)
                                                              tgtgtTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGc                           >  2:7421763/1‑93 (MQ=255)
                                                              tgtgtTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGc                           <  2:6093282/93‑1 (MQ=255)
                                                               gtgtTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGcc                          <  2:4443722/93‑1 (MQ=255)
                                                                tgtTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGccc                         <  1:3715175/93‑1 (MQ=255)
                                                                 gtTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCg                        >  1:1968413/1‑93 (MQ=255)
                                                                 gtTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCg                        >  2:1850132/1‑93 (MQ=255)
                                                                  ttGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGa                       >  1:1163089/1‑93 (MQ=255)
                                                                  ttGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGa                       >  2:7896997/1‑93 (MQ=255)
                                                                  ttGATCTGATCCAGGGGCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCACCCGCCGAATCAATAATTGCCCGa                       <  2:1564904/93‑1 (MQ=255)
                                                                   tgatctgatcCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGAt                      <  2:2838303/93‑1 (MQ=255)
                                                                    gatctgatcCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATa                     <  1:1812308/93‑1 (MQ=255)
                                                                    gatctgatcCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATa                     <  2:2191395/93‑1 (MQ=255)
                                                                    gatctgatcCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATa                     >  1:307459/1‑93 (MQ=255)
                                                                    gatctgatcCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATa                     <  1:6418901/93‑1 (MQ=255)
                                                                      tctgatcCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGt                   >  2:3426785/1‑93 (MQ=255)
                                                                        tgatcCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGa                 >  2:4759852/1‑93 (MQ=255)
                                                                        tgatcCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGa                 <  2:6135340/93‑1 (MQ=255)
                                                                         gatcCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGAt                <  2:6330886/93‑1 (MQ=255)
                                                                          atcCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATc               >  1:1686339/1‑93 (MQ=255)
                                                                          atcCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATc               >  2:6633655/1‑93 (MQ=255)
                                                                           tcCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCg              >  2:866923/1‑93 (MQ=255)
                                                                           tcCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCg              >  2:382348/1‑93 (MQ=255)
                                                                           tcCAGGGTCAGTTTTGCCCCACTCAGGACTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCg              <  1:3557310/93‑1 (MQ=255)
                                                                            ccAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGc             <  1:7864218/93‑1 (MQ=255)
                                                                             cAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCa            >  1:3793082/1‑93 (MQ=255)
                                                                              aGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAg           <  1:5402438/93‑1 (MQ=255)
                                                                               gggTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAgt          <  2:3090384/93‑1 (MQ=255)
                                                                               gggTCAATTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAgt          <  2:7031988/93‑1 (MQ=255)
                                                                                ggTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAgtg         <  1:6668368/93‑1 (MQ=255)
                                                                                ggTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAgtg         >  2:7644513/1‑93 (MQ=255)
                                                                                 gTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAgtgt        >  1:5559721/1‑93 (MQ=255)
                                                                                 gTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAgtgt        <  2:7269743/93‑1 (MQ=255)
                                                                                      ttttGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGttcttc   <  1:7106348/93‑1 (MQ=255)
                                                                                      ttttGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGttcttc   >  2:6851043/1‑93 (MQ=255)
                                                                                          |                                                                                         
TCCATCTGATAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGTTCTTCA  >  NC_000913/3370688‑3370867

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