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.5% 687.0 / NA 219V145A (GTA→GCA) nanAN‑acetylneuraminate lyase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (113/105);  minor base C (0/1);  total (113/106)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 226 total aligned reads displayed.

CCATCTGATAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGTTCTTCA  >  NC_000913/3370689‑3370867
                                                                                         |                                                                                         
 cATCTGATAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGcccc                                                                                       >  1:1965891/1‑93 (MQ=255)
 cATCTGATAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGcccc                                                                                       >  1:4431356/1‑93 (MQ=255)
  aTCTGATAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCa                                                                                      >  2:1066142/1‑93 (MQ=255)
     tGATAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTc                                                                                   >  1:1107641/1‑93 (MQ=255)
     tGATAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTc                                                                                   >  2:2922922/1‑93 (MQ=255)
      gATAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCa                                                                                  >  1:4494802/1‑93 (MQ=255)
       aTAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAg                                                                                 <  2:69746/93‑1 (MQ=255)
       aTAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAg                                                                                 >  2:4638649/1‑93 (MQ=255)
        tAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAgg                                                                                >  2:3033761/1‑93 (MQ=255)
        tAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAgg                                                                                >  2:2457319/1‑93 (MQ=255)
        tAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAgg                                                                                >  2:1760013/1‑93 (MQ=255)
          gagaTCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGAAATGTAACAAGTGGGTTGTTCTGATCCAGGGCCAGTTTTGCCCCACTCAgggt                                                                              >  1:4410253/1‑92 (MQ=255)
            gaTCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTg                                                                            >  1:1165288/1‑93 (MQ=255)
            gaTCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTg                                                                            >  1:2241465/1‑93 (MQ=255)
             aTCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTgg                                                                           >  2:2230212/1‑93 (MQ=255)
               cGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGaa                                                                         <  1:815990/93‑1 (MQ=255)
               cGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGaa                                                                         <  1:169029/93‑1 (MQ=255)
                  cAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAAtgt                                                                      <  2:1217196/93‑1 (MQ=255)
                   agagGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAAtgtt                                                                     >  2:1980710/1‑93 (MQ=255)
                   agagGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAAtgtt                                                                     <  1:2102407/93‑1 (MQ=255)
                   agagGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAAtgtt                                                                     <  2:3594341/93‑1 (MQ=255)
                    gagGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAAtgttg                                                                    >  1:3197295/1‑93 (MQ=255)
                    gagGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGGCAGTTTTCTTCCACTCAGGGCTGGAAtgttg                                                                    >  2:4302338/1‑93 (MQ=255)
                     agGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAAtgttgt                                                                   <  1:3857545/93‑1 (MQ=255)
                      ggTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTa                                                                  <  1:977078/93‑1 (MQ=255)
                        tCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTaca                                                                >  1:1029213/1‑93 (MQ=255)
                          tGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTAcacc                                                              >  1:3814246/1‑93 (MQ=255)
                            tttCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTAcaccac                                                            <  2:599960/93‑1 (MQ=255)
                            tttCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTAcaccac                                                            >  2:2624876/1‑93 (MQ=255)
                              tCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACAccacca                                                          >  1:4126341/1‑93 (MQ=255)
                               cAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCAt                                                         <  1:2706124/93‑1 (MQ=255)
                               cAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCAt                                                         <  1:695207/93‑1 (MQ=255)
                                aGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATc                                                        <  1:2494030/93‑1 (MQ=255)
                                 gcgcACCTCCGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCg                                                       <  1:1794500/93‑1 (MQ=255)
                                  cgcACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCgg                                                      >  1:2456663/1‑93 (MQ=255)
                                      ccTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCTGCaaa                                                  >  1:779040/1‑93 (MQ=255)
                                      ccTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCaaa                                                  >  1:880008/1‑93 (MQ=255)
                                      ccTACGCCAGGCAATGTAACAAGAGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCaaa                                                  >  2:4681578/1‑93 (MQ=255)
                                       cTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAAc                                                 <  2:3253770/93‑1 (MQ=255)
                                          cGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCAt                                              <  2:2211361/93‑1 (MQ=255)
                                           gCCAGGCANTGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATc                                             <  2:4396709/93‑1 (MQ=255)
                                            ccAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATcc                                            >  2:3566007/1‑93 (MQ=255)
                                             cAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATccg                                           >  1:39032/1‑93 (MQ=255)
                                             cAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATccg                                           >  2:572787/1‑93 (MQ=255)
                                              aGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATccgc                                          >  1:2550139/1‑93 (MQ=255)
                                                gCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATccgccg                                        >  1:259059/1‑93 (MQ=255)
                                                  aaTGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGaa                                      <  1:437452/93‑1 (MQ=255)
                                                  aaTGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGaa                                      >  1:564374/1‑93 (MQ=255)
                                                   aTGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAAt                                     <  1:2662851/93‑1 (MQ=255)
                                                    tGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATc                                    <  2:3121265/93‑1 (MQ=255)
                                                    tGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATc                                    <  1:1300971/93‑1 (MQ=255)
                                                    tGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATc                                    >  1:144790/1‑93 (MQ=255)
                                                     gTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCa                                   >  1:2738294/1‑93 (MQ=255)
                                                     gTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCa                                   >  1:3360638/1‑93 (MQ=255)
                                                      tAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCaa                                  >  2:4239456/1‑93 (MQ=255)
                                                      tAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCaa                                  <  2:1562407/93‑1 (MQ=255)
                                                        aCAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCaata                                <  2:36804/93‑1 (MQ=255)
                                                        aCAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCaata                                >  2:3705563/1‑93 (MQ=255)
                                                        aCAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCaata                                <  1:2192771/93‑1 (MQ=255)
                                                        aCAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCaata                                >  1:2001325/1‑93 (MQ=255)
                                                         cAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCaataa                               >  1:888724/1‑93 (MQ=255)
                                                         cAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCaataa                               >  2:62594/1‑93 (MQ=255)
                                                          aaGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCaataat                              <  2:385455/93‑1 (MQ=255)
                                                          aaGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCaataat                              >  1:3072141/1‑93 (MQ=255)
                                                           aGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAAtt                             <  2:4039797/93‑1 (MQ=255)
                                                           aGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAAtt                             >  2:2295117/1‑93 (MQ=255)
                                                           aGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAAtt                             <  2:2324625/93‑1 (MQ=255)
                                                           aGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAAtt                             <  1:1627914/93‑1 (MQ=255)
                                                           aGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAAtt                             <  2:4541650/93‑1 (MQ=255)
                                                            gtgtgtTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTg                            >  1:1848990/1‑93 (MQ=255)
                                                             tgtgtTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGc                           >  2:2847/1‑93 (MQ=255)
                                                              gtgtTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGcc                          <  1:1557569/93‑1 (MQ=255)
                                                               tgtTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGccc                         <  1:415942/93‑1 (MQ=255)
                                                               tgtTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGccc                         <  2:3451466/93‑1 (MQ=255)
                                                                gtTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCg                        >  2:2783055/1‑93 (MQ=255)
                                                                gtTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCg                        >  1:4196129/1‑93 (MQ=255)
                                                                gtTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCg                        >  1:3299506/1‑93 (MQ=255)
                                                                   gatctgatcCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATTTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATa                     >  1:2516589/1‑93 (MQ=255)
                                                                   gatctgatcCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATa                     <  2:1459865/93‑1 (MQ=255)
                                                                     tctgatcCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGt                   <  2:2047311/93‑1 (MQ=255)
                                                                      ctaatcCAGGTTCAGTTTTCCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTg                  <  2:2575961/90‑1 (MQ=255)
                                                                       tgatcCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGa                 <  1:3663988/93‑1 (MQ=255)
                                                                          tcCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCg              >  2:81872/1‑93 (MQ=255)
                                                                          tcCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCg              <  1:2416109/93‑1 (MQ=255)
                                                                           ccAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACGATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGc             >  2:3782437/1‑93 (MQ=255)
                                                                           ccAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGc             >  1:1209263/1‑93 (MQ=255)
                                                                           ccAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGc             >  1:1902474/1‑93 (MQ=255)
                                                                            cAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCa            <  1:2115401/93‑1 (MQ=255)
                                                                              gggTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAgt          >  2:2370341/1‑93 (MQ=255)
                                                                              gggTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAgt          >  2:1899525/1‑93 (MQ=255)
                                                                              gggTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAgt          >  1:4285144/1‑93 (MQ=255)
                                                                              gggTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAgt          >  1:2909136/1‑93 (MQ=255)
                                                                              gggTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAgt          >  2:882591/1‑93 (MQ=255)
                                                                               ggTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAACCCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAgtg         <  1:3436916/93‑1 (MQ=255)
                                                                                gTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAgtgt        >  2:3529959/1‑93 (MQ=255)
                                                                                    gTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGttctt    >  2:4165509/1‑93 (MQ=255)
                                                                                    gTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGttctt    >  1:3573334/1‑93 (MQ=255)
                                                                                    gTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGttctt    <  2:800643/93‑1 (MQ=255)
                                                                                     ttttGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGttcttc   <  2:2505847/93‑1 (MQ=255)
                                                                                      tttGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGTTCTTCa  >  1:3786874/1‑93 (MQ=255)
                                                                                         |                                                                                         
CCATCTGATAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGTTCTTCA  >  NC_000913/3370689‑3370867

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