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% 1212.9 / ‑5.1 380V145A (GTA→GCA) nanAN‑acetylneuraminate lyase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (226/152);  total (228/152)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

Only 100 of 386 total aligned reads displayed.

TCCATCTGATAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGTTCTTCA  >  NC_000913/3370688‑3370867
                                                                                          |                                                                                         
tCCATCTGATAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGcc                                                                                         >  1:1225531/1‑93 (MQ=255)
tCCATCTGATAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGcc                                                                                         >  2:325660/1‑93 (MQ=255)
tCCATCTGATAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGcc                                                                                         >  1:6607447/1‑93 (MQ=255)
 ccATCTGATAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGccc                                                                                        >  2:3691724/1‑93 (MQ=255)
  cATCTGATAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGcccc                                                                                       <  1:6415094/93‑1 (MQ=255)
  cATCTGATAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGAGTCAGTTTTGcccc                                                                                       >  1:187/1‑93 (MQ=255)
    tCTGATAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCAc                                                                                     >  1:7253092/1‑93 (MQ=255)
     cTGATAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACt                                                                                    <  1:105037/93‑1 (MQ=255)
      tGATAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTc                                                                                   >  2:2777910/1‑93 (MQ=255)
      tGATAGAGATCGCCAGAGGCCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTc                                                                                   >  1:7779986/1‑93 (MQ=255)
       gATAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCa                                                                                  >  2:2239959/1‑93 (MQ=255)
       gATAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCa                                                                                  <  1:6364437/93‑1 (MQ=255)
         tAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAgg                                                                                >  1:5002827/1‑93 (MQ=255)
           gagaTCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGc                                                                              >  1:2036574/1‑93 (MQ=255)
            agaTCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCt                                                                             >  1:5496588/1‑93 (MQ=255)
             gaTCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTg                                                                            <  2:5873819/93‑1 (MQ=255)
             gaTCGCCAGAGGCCTGTTTCAACGCACTTACTCCATGCAATGTAAGAAGTGTGTTGATCTGATCCAGGGTCAGGTGTGCTCTACTCAGTGCTg                                                                            >  2:6459158/1‑93 (MQ=255)
               tCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGa                                                                          >  1:8774475/1‑93 (MQ=255)
                cGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGaa                                                                         <  1:240198/93‑1 (MQ=255)
                cGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGaa                                                                         <  1:8045066/93‑1 (MQ=255)
                 gCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAAt                                                                        >  1:9758501/1‑93 (MQ=255)
                 gCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAAt                                                                        >  1:940529/1‑93 (MQ=255)
                 gCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGCGCTGGAAt                                                                        >  2:7684843/1‑93 (MQ=255)
                  ccAGAGGTCTGTTTCAGCGCACCTACGCCAGGCACTGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCTACTCAGGGCTGGAAtg                                                                       >  1:5118849/1‑93 (MQ=255)
                  ccAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAAtg                                                                       >  1:6224372/1‑93 (MQ=255)
                     gagGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAAtgttg                                                                    >  1:3061107/1‑93 (MQ=255)
                        gTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTac                                                                 >  2:531985/1‑93 (MQ=255)
                         tCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTaca                                                                <  2:6650951/93‑1 (MQ=255)
                         tCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTaca                                                                >  1:8392966/1‑93 (MQ=255)
                          cTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTacac                                                               >  2:6322787/1‑93 (MQ=255)
                           tGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTAcacc                                                              >  2:10046557/1‑93 (MQ=255)
                            gTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTAcaccc                                                             >  2:2467472/1‑92 (MQ=255)
                             tttCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTAcaccac                                                            >  1:5662009/1‑93 (MQ=255)
                             tttCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTAcaccac                                                            >  1:4100653/1‑93 (MQ=255)
                             tttCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTAcaccac                                                            >  2:2922569/1‑93 (MQ=255)
                               tCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACAccacca                                                          <  2:1266057/93‑1 (MQ=255)
                               tCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACAccacca                                                          <  2:4231149/93‑1 (MQ=255)
                                cAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCAt                                                         >  1:337481/1‑93 (MQ=255)
                                 aGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATc                                                        >  1:3677412/1‑93 (MQ=255)
                                  gcgcACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCg                                                       >  2:1364191/1‑93 (MQ=255)
                                     cACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCa                                                    >  2:7214655/1‑93 (MQ=255)
                                     cACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCa                                                    >  2:5600878/1‑93 (MQ=255)
                                     cACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCa                                                    <  2:4088744/93‑1 (MQ=255)
                                       ccTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCaaa                                                  <  2:7368705/93‑1 (MQ=255)
                                       ccTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCaaa                                                  <  1:9534603/93‑1 (MQ=255)
                                         tACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAAcc                                                >  2:9693054/1‑93 (MQ=255)
                                          aCGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCa                                               >  1:5295336/1‑93 (MQ=255)
                                            gCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATc                                             <  1:3736872/93‑1 (MQ=255)
                                              cAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATccg                                           >  1:2473837/1‑93 (MQ=255)
                                                 gCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATGGCAAACCAAccgccca                                       >  1:1114668/1‑91 (MQ=255)
                                                   aaTGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGaa                                      >  1:594507/1‑93 (MQ=255)
                                                   aaTGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGaa                                      >  2:1523326/1‑93 (MQ=255)
                                                    aTGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAAt                                     >  2:9247665/1‑93 (MQ=255)
                                                     tGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATc                                    >  2:9504628/1‑93 (MQ=255)
                                                     tGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATc                                    >  1:273182/1‑93 (MQ=255)
                                                      gTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCa                                   >  2:2046336/1‑93 (MQ=255)
                                                       tAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCaa                                  <  1:4774522/93‑1 (MQ=255)
                                                        aaCAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCaat                                 <  1:3146402/93‑1 (MQ=255)
                                                         aCAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCaata                                <  1:7948906/93‑1 (MQ=255)
                                                         aCAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCaata                                <  2:1108771/93‑1 (MQ=255)
                                                          cAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCATTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCaataa                               <  2:6951664/93‑1 (MQ=255)
                                                          cAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCaataa                               <  2:3473807/93‑1 (MQ=255)
                                                           aaGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCaataat                              >  2:3316262/1‑93 (MQ=255)
                                                            aGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAAtt                             >  1:5810257/1‑93 (MQ=255)
                                                            aGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAAtt                             >  1:3966400/1‑93 (MQ=255)
                                                            aGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAAtt                             >  2:4426674/1‑93 (MQ=255)
                                                             gtgtgtTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTg                            >  1:4330980/1‑93 (MQ=255)
                                                             gtgtgtTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTg                            >  2:168361/1‑93 (MQ=255)
                                                              tgtgtTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGc                           >  1:9314858/1‑93 (MQ=255)
                                                              tgtgtTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGc                           >  1:897517/1‑93 (MQ=255)
                                                               gtgtTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGcc                          <  1:4237721/93‑1 (MQ=255)
                                                                   tgatctgatcCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGAt                      <  2:2310872/93‑1 (MQ=255)
                                                                   tgatctgatcCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGAt                      >  1:890372/1‑93 (MQ=255)
                                                                   tgatctgatcCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGAt                      >  1:2836375/1‑93 (MQ=255)
                                                                    gatctgatcCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATa                     <  1:755127/93‑1 (MQ=255)
                                                                    gatctgatcCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATa                     >  2:2081009/1‑93 (MQ=255)
                                                                     atcttatcCAGGGTCAGTTTTGCCCCACGCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAg                    <  1:6821730/93‑1 (MQ=255)
                                                                      tctgatcCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAgg                   >  1:3430923/1‑92 (MQ=255)
                                                                      tctgatcCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGt                   >  1:1532414/1‑93 (MQ=255)
                                                                       ctgatcCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTg                  <  2:9839678/93‑1 (MQ=255)
                                                                       ctgatcCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTg                  <  1:2655856/93‑1 (MQ=255)
                                                                        tgatcCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCACTAATTGCCCGATAGTGa                 <  2:7962294/93‑1 (MQ=255)
                                                                          atcCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATc               <  1:7045771/93‑1 (MQ=255)
                                                                           tcCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCg              >  2:1850097/1‑93 (MQ=255)
                                                                           tcCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCg              >  1:8742779/1‑93 (MQ=255)
                                                                           tcCAGGGGCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCAAATCAATAATTGCCCGATAGAGATCg              >  2:2188621/1‑93 (MQ=255)
                                                                            ccAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGc             >  1:325951/1‑93 (MQ=255)
                                                                              aGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAg           <  2:4656488/93‑1 (MQ=255)
                                                                              aGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAg           <  2:2587284/93‑1 (MQ=255)
                                                                              aGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAg           >  2:8508706/1‑93 (MQ=255)
                                                                               gggTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAgt          <  2:9085124/93‑1 (MQ=255)
                                                                                 gTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAgtgt        >  1:1732723/1‑93 (MQ=255)
                                                                                 gTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAgtgt        <  1:9232426/93‑1 (MQ=255)
                                                                                  tCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAACCCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGtt       <  2:8832153/93‑1 (MQ=255)
                                                                                    aGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGttct     <  1:5478382/93‑1 (MQ=255)
                                                                                    aGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGttct     >  2:3162479/1‑93 (MQ=255)
                                                                                      ttttGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGttcttc   <  1:4438381/93‑1 (MQ=255)
                                                                                      ttttGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGttcttc   >  2:6167365/1‑93 (MQ=255)
                                                                                       tttGCCCCACTCAGGGCTGGAATGTTGTACCCCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGTTCTTCa  <  1:5424194/93‑1 (MQ=255)
                                                                                       tttGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGTTCTTCa  <  2:4954114/93‑1 (MQ=255)
                                                                                          |                                                                                         
TCCATCTGATAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGTTCTTCA  >  NC_000913/3370688‑3370867

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: