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,7780AG97.7% 533.8 / 1.1 173V145A (GTA→GCA) nanAN‑acetylneuraminate lyase
Reads supporting (aligned to +/- strand):  ref base A (4/0);  new base G (81/88);  total (85/88)
Fisher's exact test for biased strand distribution p-value = 5.62e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.78e-01

Only 100 of 175 total aligned reads displayed.

TCCATCTGATAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGTTCTTC  >  NC_000913/3370688‑3370866
                                                                                          |                                                                                        
tCCATCTGATAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGcc                                                                                        >  1:3598320/1‑93 (MQ=255)
 ccATCTGATAGAGATCGCCAGAGGGCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGccc                                                                                       >  1:3020344/1‑93 (MQ=255)
  cATCTGATAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGcccc                                                                                      <  1:690698/93‑1 (MQ=255)
  cATCTGATAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGcccc                                                                                      >  1:2709656/1‑93 (MQ=255)
   aTCTGATAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCa                                                                                     >  1:88773/1‑93 (MQ=255)
   aTCTGATAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCa                                                                                     <  1:325860/93‑1 (MQ=255)
   aTCTGATAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCa                                                                                     <  1:2100393/93‑1 (MQ=255)
    tCTGATAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCAc                                                                                    >  2:2081354/1‑93 (MQ=255)
      tGATAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTc                                                                                  <  2:102377/93‑1 (MQ=255)
       gATAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCa                                                                                 >  2:4091059/1‑93 (MQ=255)
       gATAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCa                                                                                 <  2:3509704/93‑1 (MQ=255)
       gATAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCa                                                                                 <  1:4669798/93‑1 (MQ=255)
        aTAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAg                                                                                <  1:810154/93‑1 (MQ=255)
        aTAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAg                                                                                <  1:3798906/93‑1 (MQ=255)
          agagATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAggg                                                                              >  1:4650011/1‑93 (MQ=255)
          agagATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAggg                                                                              <  2:1617855/93‑1 (MQ=255)
           gagaTCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGc                                                                             <  1:2825897/93‑1 (MQ=255)
           gagaTCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGc                                                                             <  2:511149/93‑1 (MQ=255)
            agaTCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCt                                                                            <  1:4016752/93‑1 (MQ=255)
               tCGCCAGAGGCCTGTTTCAGCGCACCTACGCCACGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGa                                                                         >  2:1761109/1‑93 (MQ=255)
                     gagGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAAtgttg                                                                   <  1:4602529/93‑1 (MQ=255)
                      agGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAAtgttgt                                                                  >  1:3311872/1‑93 (MQ=255)
                       ggTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTa                                                                 >  2:141319/1‑93 (MQ=255)
                            gTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTAcacca                                                            <  2:2912544/93‑1 (MQ=255)
                               tCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACAccacca                                                         <  1:3907568/93‑1 (MQ=255)
                               tCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACAccacca                                                         >  1:3772337/1‑93 (MQ=255)
                                 aGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCATGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATc                                                       >  1:1355435/1‑93 (MQ=255)
                                 aGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATc                                                       <  2:1986818/93‑1 (MQ=255)
                                  gcgcACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCg                                                      >  1:1489130/1‑93 (MQ=255)
                                   cgcACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCgg                                                     <  1:2533385/93‑1 (MQ=255)
                                    gcACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGc                                                    >  2:1574227/1‑93 (MQ=255)
                                    gcACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGc                                                    >  1:851812/1‑93 (MQ=255)
                                      aCCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCaa                                                  >  1:4379761/1‑93 (MQ=255)
                                       ccTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCaaa                                                 <  1:3508250/93‑1 (MQ=255)
                                       ccTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCaaa                                                 >  2:3591610/1‑93 (MQ=255)
                                        cTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCaagc                                                <  2:2300321/93‑3 (MQ=255)
                                        cTACGCCAGGAAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAAc                                                <  1:1302313/93‑1 (MQ=255)
                                         tACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAAcc                                               >  2:3400780/1‑93 (MQ=255)
                                         tACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAAcc                                               >  2:4225852/1‑93 (MQ=255)
                                            gCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATc                                            <  1:3186447/93‑1 (MQ=255)
                                            gCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATc                                            <  1:581210/93‑1 (MQ=255)
                                             ccAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATcc                                           >  2:4459434/1‑93 (MQ=255)
                                             ccAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATcc                                           >  2:1843636/1‑93 (MQ=255)
                                             ccAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATcc                                           <  2:2440171/93‑1 (MQ=255)
                                             ccAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATcc                                           >  2:728556/1‑93 (MQ=255)
                                             ccAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATcc                                           <  1:1388606/93‑1 (MQ=255)
                                             ccAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATcc                                           <  1:3711967/93‑1 (MQ=255)
                                             ccAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTCGAATGTTGTACACCACCATCGGCAAACCATcc                                           <  2:3941817/93‑1 (MQ=255)
                                                ggCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATccgcc                                        >  1:1406143/1‑93 (MQ=255)
                                                  cAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGa                                      <  2:3613717/93‑1 (MQ=255)
                                                   aaTGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGaa                                     >  2:1723159/1‑93 (MQ=255)
                                                   aaTGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGaa                                     >  2:1581788/1‑93 (MQ=255)
                                                     tGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATc                                   <  2:4329393/93‑1 (MQ=255)
                                                     tGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATc                                   <  2:1738887/93‑1 (MQ=255)
                                                      gTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCa                                  >  2:1919621/1‑93 (MQ=255)
                                                       tAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCaa                                 >  2:2128828/1‑93 (MQ=255)
                                                        aaCAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCaat                                <  1:3436385/93‑1 (MQ=255)
                                                         aCAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCaata                               >  1:1211021/1‑93 (MQ=255)
                                                         aCAAGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCaata                               >  1:3084049/1‑93 (MQ=255)
                                                         aCAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCaata                               >  2:4044637/1‑93 (MQ=255)
                                                           aaGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCaataat                             <  1:2989110/93‑1 (MQ=255)
                                                            aGTGTGTTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCGAACCATCCGCCGAATCAATAAtt                            >  2:2583695/1‑93 (MQ=255)
                                                              tgtgtTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCACACCATCCGCCGCATCAATAATTGc                          >  1:2571431/1‑93 (MQ=255)
                                                              tgtgtTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGc                          <  2:3992253/93‑1 (MQ=255)
                                                              tgtgtTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGc                          <  2:2391883/93‑1 (MQ=255)
                                                              tgtgtTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGc                          >  2:4279823/1‑93 (MQ=255)
                                                              tgtgtTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGc                          <  1:1244645/93‑1 (MQ=255)
                                                               gtgtTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGcc                         <  2:3071871/93‑1 (MQ=255)
                                                               gtgtTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGcc                         <  2:4697213/93‑1 (MQ=255)
                                                               gtgtTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGcc                         >  1:4191790/1‑93 (MQ=255)
                                                                tgtTGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGccc                        <  1:4262280/93‑1 (MQ=255)
                                                                  ttGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGa                      <  2:3642908/93‑1 (MQ=255)
                                                                  ttGATCTGATCCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGa                      >  2:3224119/1‑93 (MQ=255)
                                                                   tgatctgatcCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGAt                     <  1:650951/93‑1 (MQ=255)
                                                                    gatctgatcCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATa                    >  2:191342/1‑93 (MQ=255)
                                                                    gatctgatcCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATa                    <  2:1147579/93‑1 (MQ=255)
                                                                     atctgatcCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAg                   <  1:203022/93‑1 (MQ=255)
                                                                     atctgatcCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAg                   >  2:1308004/1‑93 (MQ=255)
                                                                        tgatcCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGa                <  1:959019/93‑1 (MQ=255)
                                                                        tgatcCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGa                >  1:2953818/1‑93 (MQ=255)
                                                                        tgatcCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGa                <  2:1690381/93‑1 (MQ=255)
                                                                         gatcCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGAt               >  2:2018853/1‑93 (MQ=255)
                                                                         gatcCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGGACACCACCATCGGCAAAGCATCCGCCGAATCAATAATTGCCCGATAGTGAt               >  2:4388044/1‑93 (MQ=255)
                                                                          atcCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATc              >  1:764375/1‑93 (MQ=255)
                                                                          atcCAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATc              <  1:3791296/93‑1 (MQ=255)
                                                                             cAGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCa           <  2:3286685/93‑1 (MQ=255)
                                                                              aGGGTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAg          <  2:4471267/93‑1 (MQ=255)
                                                                               gggTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAgt         <  1:2917262/93‑1 (MQ=255)
                                                                               gggTCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAgt         >  2:3111346/1‑93 (MQ=255)
                                                                                ggTCAGTTTTGCCCCACTCAGGACTGGAATGTTGTACACCACCATCGGCAATCCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAgtg        <  2:2529003/93‑1 (MQ=255)
                                                                                ggTCAGTGTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAgtg        <  1:4179333/93‑1 (MQ=255)
                                                                                  tCAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGtt      >  2:3871156/1‑93 (MQ=255)
                                                                                   cAGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGttc     <  1:1910526/93‑1 (MQ=255)
                                                                                    aGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGttct    <  1:1516019/93‑1 (MQ=255)
                                                                                    aGTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGAAGTGttct    <  1:2155104/93‑1 (MQ=255)
                                                                                     gTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCACACCATCCGCCGACTCAATAATTGCCCGATAGTGATCGCAGTGttctt   >  1:167959/1‑93 (MQ=255)
                                                                                     gTTTTGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGttctt   >  1:167827/1‑93 (MQ=255)
                                                                                      ttttGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGttcttc  >  2:3375142/1‑93 (MQ=255)
                                                                                      ttttGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGttcttc  <  1:2240008/93‑1 (MQ=255)
                                                                                      ttttGCCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGttcttc  >  1:4738142/1‑93 (MQ=255)
                                                                                          |                                                                                        
TCCATCTGATAGAGATCGCCAGAGGTCTGTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGTTCTTC  >  NC_000913/3370688‑3370866

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: