Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,806,842 C→T R226C (CGC→TGC)  kdtA → 3‑deoxy‑D‑manno‑octulosonic‑acid transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,806,8420CT100.0% 938.4 / NA 256R226C (CGC→TGC) kdtA3‑deoxy‑D‑manno‑octulosonic‑acid transferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (133/123);  total (133/123)

Only 100 of 275 total aligned reads displayed.

GCAAAAAATAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTA  >  NC_000913/3806752‑3806934
                                                                                          |                                                                                            
 cAAAAAATAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCgctgcc                                                                                           >  2:4024901/1‑93 (MQ=255)
  aaaaaaTAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCa                                                                                          <  2:1407443/93‑1 (MQ=255)
   aaaaaTAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAg                                                                                         >  2:570047/1‑93 (MQ=255)
    aaaaTAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGt                                                                                        >  1:3954942/1‑93 (MQ=255)
    aaaaTAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGt                                                                                        <  1:1766387/93‑1 (MQ=255)
     aaaTAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTg                                                                                       <  1:3724874/93‑1 (MQ=255)
     aaaTAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTg                                                                                       <  1:1876928/93‑1 (MQ=255)
      aataatCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTgg                                                                                      <  2:3049534/93‑1 (MQ=255)
      aataatCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTgg                                                                                      <  2:1024678/93‑1 (MQ=255)
        taatCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGc                                                                                    >  1:3069142/1‑93 (MQ=255)
          atCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGcac                                                                                  <  2:704171/93‑1 (MQ=255)
           tCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGcacc                                                                                 >  1:2533610/1‑93 (MQ=255)
              ggTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACcaca                                                                              >  1:3194488/1‑93 (MQ=255)
               gTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCacac                                                                             <  2:2163485/93‑1 (MQ=255)
                  aCCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGc                                                                          <  2:3154223/93‑1 (MQ=255)
                  aCCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGc                                                                          >  1:4657151/1‑93 (MQ=255)
                  aCCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGc                                                                          >  2:208754/1‑93 (MQ=255)
                  aCCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGc                                                                          >  1:2781801/1‑93 (MQ=255)
                   ccGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGcc                                                                         <  2:3710561/93‑1 (MQ=255)
                   ccGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGcc                                                                         <  1:223608/93‑1 (MQ=255)
                    cGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGccc                                                                        <  1:2408082/93‑1 (MQ=255)
                    cGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGccc                                                                        <  2:1461818/93‑1 (MQ=255)
                     gTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCg                                                                       >  1:1130587/1‑93 (MQ=255)
                      ttACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCgg                                                                      >  2:1677335/1‑93 (MQ=255)
                      ttACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCgg                                                                      >  1:3338487/1‑93 (MQ=255)
                      ttACCGGTAGCCTGAAATTCGATATTTCTGCAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACcccgcccgt                                                                      >  2:152412/1‑92 (MQ=255)
                           ggTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATgg                                                                 <  2:3419165/93‑1 (MQ=255)
                           ggTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATgg                                                                 >  2:3639913/1‑93 (MQ=255)
                             tAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGAt                                                               >  1:2264671/1‑93 (MQ=255)
                             tAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGAt                                                               >  2:112526/1‑93 (MQ=255)
                             tAGCCTGAAATTCGACATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACTGCCCGGTATGGAt                                                               <  2:4605166/93‑1 (MQ=255)
                              aGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGAtt                                                              >  2:2251237/1‑93 (MQ=255)
                               gCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTg                                                             <  2:3365368/93‑1 (MQ=255)
                                ccTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGc                                                            <  1:9271/93‑1 (MQ=255)
                                 cTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGcc                                                           <  2:1182012/93‑1 (MQ=255)
                                  tGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGcca                                                          >  1:1638942/1‑93 (MQ=255)
                                  tGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGcca                                                          >  1:2708971/1‑93 (MQ=255)
                                   gAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGccac                                                         >  2:3387312/1‑93 (MQ=255)
                                     aaTTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGccacca                                                       >  1:367179/1‑93 (MQ=255)
                                     aaTTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGccacca                                                       >  2:1877232/1‑93 (MQ=255)
                                     aaTTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGccacca                                                       >  1:260169/1‑93 (MQ=255)
                                     aaTTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGccacca                                                       >  1:315254/1‑93 (MQ=255)
                                       ttCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGGACCACACCGCCCGGTATGGATTGCCCCCAGc                                                     >  2:423682/1‑93 (MQ=255)
                                       ttCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTTCcagcagg                                                     >  2:4297093/1‑92 (MQ=255)
                                        tCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCa                                                    <  2:622544/93‑1 (MQ=255)
                                         cGATATTTCTGTAACGCCGCAGTTGCCTTCTAAAGCCTGTACGCTGCGCTGCCAGTGGGCACCCCGCCGCCCGGTATGGATTGCCACCAGCAc                                                   <  1:2896912/93‑1 (MQ=255)
                                           atatTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTc                                                 >  2:746821/1‑93 (MQ=255)
                                           atatTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTc                                                 >  1:762311/1‑93 (MQ=255)
                                                tCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGaa                                            <  1:1320164/93‑1 (MQ=255)
                                                tCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGaa                                            <  1:3240089/93‑1 (MQ=255)
                                                tCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGaa                                            <  2:68699/93‑1 (MQ=255)
                                                tCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGaa                                            <  1:4083157/93‑1 (MQ=255)
                                                  tGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAgg                                          >  1:2077958/1‑93 (MQ=255)
                                                  tGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAgg                                          <  2:2061083/93‑1 (MQ=255)
                                                   gTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGc                                         >  1:2549089/1‑93 (MQ=255)
                                                    tAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCg                                        <  2:2568221/93‑1 (MQ=255)
                                                     aaCGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGa                                       >  1:427535/1‑93 (MQ=255)
                                                        gccgcAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAaga                                    <  2:2358119/93‑1 (MQ=255)
                                                         ccgcAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAagag                                   >  2:1213386/1‑93 (MQ=255)
                                                         ccgcAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAagag                                   >  1:1223832/1‑93 (MQ=255)
                                                          cgcAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAagaga                                  <  1:1252781/93‑1 (MQ=255)
                                                          cgcACTAGGCTGTTGAAGCTGGGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCAATCACGAAGGCGAagaga                                  <  2:802573/93‑1 (MQ=255)
                                                           gcAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAagagag                                 >  1:4000767/1‑93 (MQ=255)
                                                            cAGTTGGCTGCTAAAGCCGTGACGCTGCTCTGCCAGTGGGCACCACCCCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAgt                                <  2:3779079/93‑1 (MQ=255)
                                                            cAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAgt                                <  2:3099078/93‑1 (MQ=255)
                                                            cAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAgt                                <  1:1352986/93‑1 (MQ=255)
                                                            cAGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAgt                                >  2:3456920/1‑93 (MQ=255)
                                                             aGTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAgtg                               <  2:4654939/93‑1 (MQ=255)
                                                              gTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAgtgt                              <  1:4408516/93‑1 (MQ=255)
                                                              gTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAgtgt                              >  1:728855/1‑93 (MQ=255)
                                                              gTTGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAgtgg                              >  1:503512/1‑92 (MQ=255)
                                                               ttGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGtgtg                             >  1:559408/1‑93 (MQ=255)
                                                                tGGCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTgtgg                            <  2:2025041/93‑1 (MQ=255)
                                                                 ggCTGCTAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTgtggg                           >  1:4477642/1‑92 (MQ=255)
                                                                   ctgctAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGa                         <  1:151550/93‑1 (MQ=255)
                                                                    tgctAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGGACTCACGAAGGCGAAGAGAGTGGGGTGAt                        >  2:12803/1‑93 (MQ=255)
                                                                    tgctAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGAt                        <  1:3419691/93‑1 (MQ=255)
                                                                    tgctAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGAt                        >  2:3231339/1‑93 (MQ=255)
                                                                    tgctAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGAt                        >  1:2947740/1‑93 (MQ=255)
                                                                     gctAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGCGTGTGGTGATc                       >  2:300000/1‑93 (MQ=255)
                                                                       tAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATcgc                     >  1:3780142/1‑93 (MQ=255)
                                                                       tAAAGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATcgc                     <  1:82327/93‑1 (MQ=255)
                                                                        aaaGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATcgcc                    <  1:641493/93‑1 (MQ=255)
                                                                         aaGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATcgccg                   >  1:1397851/1‑93 (MQ=255)
                                                                          aGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATcgccgc                  >  1:1036519/1‑93 (MQ=255)
                                                                          aGCCGTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATcgccgc                  >  1:1443531/1‑93 (MQ=255)
                                                                              gTGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACAt              <  1:283215/93‑1 (MQ=255)
                                                                               tGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATc             <  1:1581687/93‑1 (MQ=255)
                                                                               tGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATc             <  1:948620/93‑1 (MQ=255)
                                                                               tGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATc             <  2:4683571/93‑1 (MQ=255)
                                                                               tGACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATc             >  2:3965672/1‑93 (MQ=255)
                                                                               tGACGCTCCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATc             <  1:2021469/93‑1 (MQ=255)
                                                                                gACGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCa            <  1:2294569/93‑1 (MQ=255)
                                                                                 aCGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAg           <  1:1282183/93‑1 (MQ=255)
                                                                                 aCGCTGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAg           <  2:2911151/93‑1 (MQ=255)
                                                                                   gctgcgctgcCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGc         <  2:2458021/93‑1 (MQ=255)
                                                                                     tGCGCTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCAt       <  2:268798/93‑1 (MQ=255)
                                                                                      gcgcTGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCAtt      <  2:4435800/93‑1 (MQ=255)
                                                                                       cgctgcGAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTg     <  2:1664338/93‑1 (MQ=255)
                                                                                       cgcggcCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTg     <  2:2197319/93‑1 (MQ=255)
                                                                                          |                                                                                            
GCAAAAAATAATCAGGTGACCGTTACCGGTAGCCTGAAATTCGATATTTCTGTAACGCCGCAGTTGGCTGCTAAAGCCGTGACGCTGCGCCGCCAGTGGGCACCACACCGCCCGGTATGGATTGCCACCAGCACTCACGAAGGCGAAGAGAGTGTGGTGATCGCCGCACATCAGGCATTGTTA  >  NC_000913/3806752‑3806934

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