Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,845,852 T→C T61A (ACC→GCC)  selD ← selenophosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,845,8520TC100.0% 565.2 / NA 176T61A (ACC→GCC) selDselenophosphate synthase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (96/80);  total (96/80)

Only 100 of 180 total aligned reads displayed.

GCTGATGGCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTA  >  NC_000913/1845763‑1845940
                                                                                         |                                                                                        
   gATGGCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCAtt                                                                                    <  1:1143051/93‑1 (MQ=255)
   gATGGCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCAtt                                                                                    >  2:5864355/1‑93 (MQ=255)
    aTGGCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTg                                                                                   >  1:3556861/1‑93 (MQ=255)
     tGGCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGc                                                                                  >  1:2970525/1‑93 (MQ=255)
      ggCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGcc                                                                                 <  1:3441017/93‑1 (MQ=255)
        cGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCa                                                                               <  2:3136680/93‑1 (MQ=255)
         gTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAg                                                                              >  2:777680/1‑93 (MQ=255)
          ttAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGa                                                                             >  2:1534193/1‑93 (MQ=255)
           tAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGAt                                                                            <  1:6034166/93‑1 (MQ=255)
            aGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATc                                                                           >  1:2113614/1‑93 (MQ=255)
            aGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATc                                                                           <  1:2315227/93‑1 (MQ=255)
             gTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCg                                                                          >  1:1386701/1‑93 (MQ=255)
             gTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCg                                                                          <  2:2164913/93‑1 (MQ=255)
              tCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGt                                                                         >  2:3896226/1‑93 (MQ=255)
              tCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGt                                                                         >  1:5615840/1‑93 (MQ=255)
              tCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGt                                                                         >  2:3926373/1‑93 (MQ=255)
                gccgcAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTac                                                                       <  2:4148193/93‑1 (MQ=255)
                gccgcAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTac                                                                       <  2:1430561/93‑1 (MQ=255)
                 ccgcAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTaca                                                                      >  2:426189/1‑93 (MQ=255)
                 ccgcAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTaca                                                                      >  1:1855678/1‑93 (MQ=255)
                   gcAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACAcc                                                                    <  2:4373599/93‑1 (MQ=255)
                     aaTGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACAccgc                                                                  <  1:4045538/93‑1 (MQ=255)
                      aTGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACAccgcc                                                                 >  1:5891442/1‑93 (MQ=255)
                      aTGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACAccgcc                                                                 >  2:1110879/1‑93 (MQ=255)
                       tGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACAccgccg                                                                <  2:4514047/93‑1 (MQ=255)
                        gCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACcgccgc                                                               <  1:5099505/93‑1 (MQ=255)
                          ggCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCAt                                                             <  1:3884321/93‑1 (MQ=255)
                          ggCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCAt                                                             <  1:5569528/93‑1 (MQ=255)
                          ggCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCAt                                                             >  1:4542449/1‑93 (MQ=255)
                           gCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCAtc                                                            >  2:445516/1‑93 (MQ=255)
                             cAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCAtcgt                                                          >  1:491219/1‑93 (MQ=255)
                             cAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCAtcgt                                                          <  1:4526432/93‑1 (MQ=255)
                             cAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCAtcgt                                                          <  2:2968727/93‑1 (MQ=255)
                              aaaaTCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCAtcgtc                                                         <  1:4912492/93‑1 (MQ=255)
                              aaaaTCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCAtcgtc                                                         <  1:3937006/93‑1 (MQ=255)
                                  tCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCgg                                                     <  1:371807/93‑1 (MQ=255)
                                  tCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCgg                                                     >  1:5706664/1‑93 (MQ=255)
                                    gAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGt                                                   <  1:1671254/93‑1 (MQ=255)
                                    gAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGt                                                   <  2:613207/93‑1 (MQ=255)
                                      aaGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGttt                                                 >  1:4800312/1‑93 (MQ=255)
                                      aaGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGttt                                                 >  2:6043489/1‑93 (MQ=255)
                                       aGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTc                                                >  1:2600453/1‑93 (MQ=255)
                                        ggATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCa                                               <  1:81106/93‑1 (MQ=255)
                                         gATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCTTTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCAt                                              >  2:4093261/1‑93 (MQ=255)
                                         gATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCAt                                              <  1:73736/93‑1 (MQ=255)
                                           ttATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGGTTCATTa                                            >  2:368420/1‑93 (MQ=255)
                                           ttATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCAAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTa                                            <  2:5672313/93‑1 (MQ=255)
                                            tATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTAc                                           <  2:2188735/93‑1 (MQ=255)
                                            tATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTAc                                           <  1:3802628/93‑1 (MQ=255)
                                             aTCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTAcc                                          >  1:3469669/1‑93 (MQ=255)
                                              tCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTAccc                                         >  2:127636/1‑93 (MQ=255)
                                               cAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCca                                        <  1:2946945/93‑1 (MQ=255)
                                                aaCGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGGCCAGATCGTACACCGCCGCGTCGTCGCGGGTTTCATTACCcac                                       <  2:1836458/93‑1 (MQ=255)
                                                aaCGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCcac                                       >  1:2062023/1‑93 (MQ=255)
                                                 aCGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCcaca                                      <  1:5694106/93‑1 (MQ=255)
                                                  cGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACaa                                     <  1:3132363/93‑1 (MQ=255)
                                                  cGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACaa                                     <  2:4492126/93‑1 (MQ=255)
                                                  cGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACaa                                     <  2:3260584/93‑1 (MQ=255)
                                                  cGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACaa                                     <  2:681251/93‑1 (MQ=255)
                                                   gATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAg                                    <  2:965277/93‑1 (MQ=255)
                                                   gATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAg                                    >  2:5289497/1‑93 (MQ=255)
                                                   gATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAg                                    >  2:6223520/1‑93 (MQ=255)
                                                     tCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCa                                  >  1:5830150/1‑93 (MQ=255)
                                                     tCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCa                                  >  2:3473279/1‑93 (MQ=255)
                                                       ggCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCaaa                                <  2:5535604/93‑1 (MQ=255)
                                                        gCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAAt                               >  1:5860595/1‑93 (MQ=255)
                                                        gCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAAt                               <  2:1210121/93‑1 (MQ=255)
                                                         cATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAAtt                              <  2:25626/93‑1 (MQ=255)
                                                          aTAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTc                             <  2:2934935/93‑1 (MQ=255)
                                                          aTAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTc                             >  2:3395510/1‑93 (MQ=255)
                                                          aTAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTc                             <  2:36145/93‑1 (MQ=255)
                                                          aTAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTc                             >  1:4726379/1‑93 (MQ=255)
                                                           tAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCg                            >  1:1764773/1‑93 (MQ=255)
                                                           tAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCg                            >  2:5626241/1‑93 (MQ=255)
                                                            aaaGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCgg                           >  2:2333005/1‑93 (MQ=255)
                                                              agaagTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGAt                         <  2:918623/93‑1 (MQ=255)
                                                              agaagTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGAt                         >  1:1206931/1‑93 (MQ=255)
                                                              agaagTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGAt                         >  1:1638884/1‑93 (MQ=255)
                                                              agaagTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGAt                         <  2:5082699/93‑1 (MQ=255)
                                                               gaagTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATc                        <  1:3312201/93‑1 (MQ=255)
                                                               gaagTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATc                        <  1:5980800/93‑1 (MQ=255)
                                                               gaagTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATc                        <  1:6089783/93‑1 (MQ=255)
                                                               gaagTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATc                        >  1:6196130/1‑93 (MQ=255)
                                                               gaagTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATc                        >  2:5159846/1‑93 (MQ=255)
                                                               gaagTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATc                        <  2:4537411/93‑1 (MQ=255)
                                                                     ggtggtACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACaaa                  <  2:4651049/93‑1 (MQ=255)
                                                                     ggtggtACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACaaa                  <  2:1693144/93‑1 (MQ=255)
                                                                      gtggtACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAAc                 >  1:6276307/1‑93 (MQ=255)
                                                                       tggtACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACt                <  1:5354107/93‑1 (MQ=255)
                                                                        ggtACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACtt               <  1:2006296/93‑1 (MQ=255)
                                                                        ggaaCTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACtt               >  1:112605/4‑93 (MQ=255)
                                                                         gtACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTc              <  1:4358627/93‑1 (MQ=255)
                                                                         gtACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTc              >  2:631486/1‑93 (MQ=255)
                                                                           aCTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGc            <  1:4311187/93‑1 (MQ=255)
                                                                            cTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGcc           <  2:4988396/93‑1 (MQ=255)
                                                                              gATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACCAACTTCGCCTg         >  2:2758822/1‑93 (MQ=255)
                                                                                tAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCt       >  1:1628904/1‑93 (MQ=255)
                                                                                tAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCt       >  1:3577450/1‑93 (MQ=255)
                                                                                 aaCGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTc      >  2:2049428/1‑93 (MQ=255)
                                                                                     cTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTa  <  1:1709214/93‑1 (MQ=255)
                                                                                         |                                                                                        
GCTGATGGCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTA  >  NC_000913/1845763‑1845940

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: