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% 410.9 / NA 132T61A (ACC→GCC) selDselenophosphate synthase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (69/63);  total (69/63)

Only 100 of 137 total aligned reads displayed.

CTGATGGCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTA  >  NC_000913/1845764‑1845940
                                                                                        |                                                                                        
cTGATGGCGTTAGTCGCCGCAATGCGGCCAAAATAGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCa                                                                                      >  1:973306/1‑93 (MQ=255)
 tGATGGCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCAt                                                                                     <  2:195653/93‑1 (MQ=255)
   aTGGCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTg                                                                                   >  1:1820845/1‑93 (MQ=255)
    tGGCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGc                                                                                  <  2:1235581/93‑1 (MQ=255)
        gTTATTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAg                                                                              <  2:4066945/93‑1 (MQ=255)
        gTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAg                                                                              >  1:2851707/1‑93 (MQ=255)
         ttAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGa                                                                             <  2:1824270/93‑1 (MQ=255)
           aGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATc                                                                           >  1:2874388/1‑93 (MQ=255)
            gTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCg                                                                          >  1:421473/1‑93 (MQ=255)
            gTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCg                                                                          >  1:3861898/1‑93 (MQ=255)
             tCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGt                                                                         >  1:2717397/1‑93 (MQ=255)
              cgccgcAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTa                                                                        >  1:1474532/1‑93 (MQ=255)
              cgccgcAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTa                                                                        <  2:239984/93‑1 (MQ=255)
               gccgcAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTac                                                                       <  2:2577089/93‑1 (MQ=255)
               gccgcAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTac                                                                       >  1:2432266/1‑93 (MQ=255)
               gccgcAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTac                                                                       <  2:464751/93‑1 (MQ=255)
               gccgcAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTac                                                                       <  1:493984/93‑1 (MQ=255)
                ccgcAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGACATTGCCCAGATCGTaca                                                                      >  2:4207208/1‑93 (MQ=255)
                  gcAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACAcc                                                                    >  1:3254838/1‑93 (MQ=255)
                  gcAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACAcc                                                                    >  1:2946416/1‑93 (MQ=255)
                    aaTGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACAccgc                                                                  <  1:3092391/93‑1 (MQ=255)
                    aaTGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACAccgc                                                                  <  2:1784367/93‑1 (MQ=255)
                     aTGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACAccgcc                                                                 >  1:2234912/1‑93 (MQ=255)
                          gCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCAtc                                                            <  2:1142985/93‑1 (MQ=255)
                          gCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCAtc                                                            <  2:2742140/93‑1 (MQ=255)
                           ccAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCAtcg                                                           >  1:4235919/1‑93 (MQ=255)
                            cAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCAtcgt                                                          >  1:624313/1‑93 (MQ=255)
                            cAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCAtcgt                                                          <  1:1649753/93‑1 (MQ=255)
                            cAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGCGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCAtcgt                                                          >  1:1026934/1‑93 (MQ=255)
                              aaaTCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCAtcgtcg                                                        >  2:1447545/1‑93 (MQ=255)
                                 tCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCgg                                                     <  2:431996/93‑1 (MQ=255)
                                   gAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGt                                                   <  2:833484/93‑1 (MQ=255)
                                     aaGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGttt                                                 >  2:589469/1‑93 (MQ=255)
                                     aaGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGttt                                                 >  2:1939398/1‑93 (MQ=255)
                                     aaGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGttt                                                 <  2:2200306/93‑1 (MQ=255)
                                     aaGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGttt                                                 >  1:1199127/1‑93 (MQ=255)
                                      aGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTc                                                >  2:4254397/1‑93 (MQ=255)
                                       ggATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCa                                               >  1:1470763/1‑93 (MQ=255)
                                       ggATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCa                                               <  2:3240392/93‑1 (MQ=255)
                                        gATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCAt                                              <  1:4408884/93‑1 (MQ=255)
                                          ttATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACGCCGCCGCATCGTCGCGGGTTTCATTa                                            >  2:3794457/1‑93 (MQ=255)
                                           tATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTAc                                           >  1:4278817/1‑93 (MQ=255)
                                           tATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTAc                                           <  1:2639467/93‑1 (MQ=255)
                                           tATCAACAATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTAc                                           <  2:3123113/93‑1 (MQ=255)
                                            aTCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTAcc                                          <  2:692198/93‑1 (MQ=255)
                                             tCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTAccc                                         >  2:1994209/1‑93 (MQ=255)
                                             tCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTAccc                                         >  2:1062362/1‑93 (MQ=255)
                                              cAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCca                                        <  2:1870126/93‑1 (MQ=255)
                                               aaCGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCcac                                       >  1:1858980/1‑93 (MQ=255)
                                                aCGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCcaca                                      >  2:1131027/1‑93 (MQ=255)
                                                 cGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACaa                                     >  2:2872198/1‑93 (MQ=255)
                                                 cGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACaa                                     <  1:2890584/93‑1 (MQ=255)
                                                 cGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACaa                                     <  2:1537946/93‑1 (MQ=255)
                                                 cGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACaa                                     <  1:214117/93‑1 (MQ=255)
                                                 cGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACaa                                     <  1:362203/93‑1 (MQ=255)
                                                 cGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACaa                                     <  1:3726957/93‑1 (MQ=255)
                                                   aTCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGc                                   >  1:2800172/1‑93 (MQ=255)
                                                    tCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCa                                  >  2:2672748/1‑93 (MQ=255)
                                                     cGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCaa                                 <  2:398169/93‑1 (MQ=255)
                                                      ggCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCaaa                                <  2:1891391/93‑1 (MQ=255)
                                                      ggCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCaaa                                <  2:3992050/93‑1 (MQ=255)
                                                       gCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAAt                               >  1:1262343/1‑93 (MQ=255)
                                                       gCATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAAt                               >  2:508768/1‑93 (MQ=255)
                                                        cATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAAtt                              <  1:629099/93‑1 (MQ=255)
                                                        cATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAAtt                              >  1:4486376/1‑93 (MQ=255)
                                                        cATAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACACACAAGCAAAtt                              >  2:2284117/1‑93 (MQ=255)
                                                         aTAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCTCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTc                             >  2:2851151/1‑93 (MQ=255)
                                                          tAAAGAAGTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCg                            >  2:1090154/1‑93 (MQ=255)
                                                            aagaagTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGa                          <  2:1736723/93‑1 (MQ=255)
                                                            aagaagTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGa                          >  2:3592415/1‑93 (MQ=255)
                                                            aagaagTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGa                          >  2:3464448/1‑93 (MQ=255)
                                                             agaagTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGAt                         >  1:953495/1‑93 (MQ=255)
                                                             agaagTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGAt                         >  1:2217088/1‑93 (MQ=255)
                                                              gaagTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAGATTCGGAta                        >  1:4386710/1‑92 (MQ=255)
                                                              gaagTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATc                        <  1:373983/93‑1 (MQ=255)
                                                              aaagTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCAATCCCCACAAGCAAATTCGGATc                        <  2:3225948/92‑1 (MQ=255)
                                                                agTCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATcaa                      >  1:109852/1‑93 (MQ=255)
                                                                  tCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATcaaca                    <  2:2332801/93‑1 (MQ=255)
                                                                  tCGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATcaaca                    >  2:79885/1‑93 (MQ=255)
                                                                   cGGTGGTACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATcaacaa                   <  1:3063919/93‑1 (MQ=255)
                                                                    ggtggtACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACaaa                  >  1:1167575/1‑93 (MQ=255)
                                                                       ggtACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACtt               >  2:1523203/1‑93 (MQ=255)
                                                                       ggtACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACtt               <  2:4605395/93‑1 (MQ=255)
                                                                        gtACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTc              <  1:3419668/93‑1 (MQ=255)
                                                                        gtACTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTc              <  1:2459355/93‑1 (MQ=255)
                                                                        gtACTGAGAACGCTGGCGCCATTGCCCAGATCATACACCGCCCCATCGTCGCAGGTTTCATTCCCCACAAGCAAATTCGGACCAACAAACTTc              <  1:221936/93‑1 (MQ=255)
                                                                           cTGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGcc           <  1:3461025/93‑1 (MQ=255)
                                                                            tGATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCt          <  1:707065/93‑1 (MQ=255)
                                                                             gATAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTg         >  2:3197507/1‑93 (MQ=255)
                                                                              aTAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGc        <  2:2321034/93‑1 (MQ=255)
                                                                              aTAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGc        <  1:1553897/93‑1 (MQ=255)
                                                                              aTAACGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGc        <  2:3518067/93‑1 (MQ=255)
                                                                               aaaCGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCt       <  1:1731339/92‑1 (MQ=255)
                                                                                 aCGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCa     >  1:2734458/1‑93 (MQ=255)
                                                                                  cGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCAc    <  1:4540267/93‑1 (MQ=255)
                                                                                  cGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCAc    >  2:366989/1‑93 (MQ=255)
                                                                                  cGCTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCAc    <  1:3755390/93‑1 (MQ=255)
                                                                                    cTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTa  <  1:1004048/93‑1 (MQ=255)
                                                                                    cTGGCGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTa  <  2:1559368/93‑1 (MQ=255)
                                                                                    cTGGCGCCATTGCACAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTa  <  1:2677390/93‑1 (MQ=255)
                                                                                        |                                                                                        
CTGATGGCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTA  >  NC_000913/1845764‑1845940

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: