Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,580,094 C→T 100% D233N (GAT→AAT)  ydeM ← putative YdeN‑specific sulfatase‑maturating enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,580,0940CT100.0% 426.8 / NA 124D233N (GAT→AAT) ydeMputative YdeN‑specific sulfatase‑maturating enzyme
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (59/65);  total (59/65)

Only 100 of 124 total aligned reads displayed.

TCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTAT  >  NC_000913/1579997‑1580189
                                                                                                 |                                                                                               
tCCTTGCAGGACTCCTGGAAAATACTCATGGTATGCCCATTCCGCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGt                                                                                               <  2:3535390/100‑1 (MQ=255)
  cTTGCAGGACTCCTGGAAAATACAACTGGTATGCCAATTCCCCAAAAAACGGCTGACAACGGTTTCAAACTTACGGATGAAAAGTTCACCCACATTGtt                                                                                              >  2:1622518/1‑99 (MQ=255)
     gCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTa                                                                                          >  2:618149/1‑100 (MQ=255)
      cAGGACTCCTGGAAAATACAACTGGTATGCAAAAAAAAAAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTaa                                                                                         <  2:3938344/100‑1 (MQ=255)
       aGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAAc                                                                                        <  2:3682961/100‑1 (MQ=255)
       aGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAAc                                                                                        <  2:29400/100‑1 (MQ=255)
       aGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAAc                                                                                        <  1:3187952/100‑1 (MQ=255)
       aGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAAc                                                                                        <  1:256572/100‑1 (MQ=255)
        ggACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAAcc                                                                                       <  1:3746656/100‑1 (MQ=255)
          aCTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCa                                                                                     >  2:1373669/1‑100 (MQ=255)
           cTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCAt                                                                                    <  1:2766180/100‑1 (MQ=255)
            tCCTGGAAAATACAACTGGTATGCCCATTCCCGAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCAtt                                                                                   >  1:3956170/1‑100 (MQ=255)
             ccTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTg                                                                                  <  1:1256260/100‑1 (MQ=255)
             ccTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTg                                                                                  >  1:3635162/1‑100 (MQ=255)
             ccTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTg                                                                                  >  1:1731137/1‑100 (MQ=255)
             ccTGGAAAATACAACTGGTATGCCCATTCCACAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTg                                                                                  <  2:1284912/100‑1 (MQ=255)
              cTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGc                                                                                 <  1:3591079/100‑1 (MQ=255)
              cTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGc                                                                                 >  1:1684592/1‑100 (MQ=255)
              cTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGc                                                                                 >  2:3030668/1‑100 (MQ=255)
              cTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGc                                                                                 >  1:868577/1‑100 (MQ=255)
               tGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCa                                                                                >  1:3505946/1‑100 (MQ=255)
               tGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCa                                                                                <  2:1833824/100‑1 (MQ=255)
                   aaaTACAACTGGTATGCAAATAAAAAAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATaaa                                                                            <  2:2937451/100‑1 (MQ=255)
                     aTACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATaaaaa                                                                          <  1:2671071/100‑1 (MQ=255)
                     aTACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATaaaaa                                                                          <  1:1291507/100‑1 (MQ=255)
                        cAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATgg                                                                       >  1:907931/1‑100 (MQ=255)
                        cAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACACACATTGTTTTTAACCCATTGCATAAAAATgg                                                                       >  1:1875680/1‑100 (MQ=255)
                         aaCTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGt                                                                      >  1:270938/1‑100 (MQ=255)
                         aaCTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGt                                                                      <  1:3758539/100‑1 (MQ=255)
                          aCTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGtt                                                                     <  1:390582/100‑1 (MQ=255)
                          aCTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGtt                                                                     <  1:1858976/100‑1 (MQ=255)
                          aCTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGtt                                                                     >  1:3851392/1‑100 (MQ=255)
                          aCTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGtt                                                                     >  1:2147595/1‑100 (MQ=255)
                          aCTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGtt                                                                     >  1:2589058/1‑100 (MQ=255)
                           cTGGTATGCACATGCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTg                                                                    <  2:244412/100‑1 (MQ=255)
                            tGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATTAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGa                                                                   <  1:709397/100‑1 (MQ=255)
                            tGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGt                                                                      <  2:3635135/97‑1 (MQ=255)
                            tGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGa                                                                   >  1:2738117/1‑100 (MQ=255)
                             ggTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGAc                                                                  >  2:3972238/1‑100 (MQ=255)
                             ggTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGAc                                                                  <  2:2696289/100‑1 (MQ=255)
                             ggTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGAc                                                                  <  2:2896633/100‑1 (MQ=255)
                             ggTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGAc                                                                  >  1:1498066/1‑100 (MQ=255)
                              gTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACa                                                                 <  2:2139997/100‑1 (MQ=255)
                                 tGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACAtga                                                              >  1:3395221/1‑98 (MQ=255)
                                  gCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACAtgaa                                                             >  1:3472702/1‑97 (MQ=255)
                                   cccATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAAc                                                            <  1:3492012/100‑1 (MQ=255)
                                   cccATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAAc                                                            >  1:2729861/1‑100 (MQ=255)
                                    ccATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACt                                                           <  2:2283407/100‑1 (MQ=255)
                                    ccATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACt                                                           >  1:715835/1‑100 (MQ=255)
                                    ccATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTCACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACt                                                           <  2:2805150/100‑1 (MQ=255)
                                     cATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACtt                                                          >  1:3565130/1‑100 (MQ=255)
                                     cATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACtt                                                          <  1:2368642/100‑1 (MQ=255)
                                      aTTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTg                                                         <  1:1242278/100‑1 (MQ=255)
                                      aTTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTg                                                         <  2:3934157/100‑1 (MQ=255)
                                        tCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGcc                                                       <  2:1668141/100‑1 (MQ=255)
                                        tCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGcc                                                       <  1:3694681/100‑1 (MQ=255)
                                         tctcAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCa                                                      <  2:4097829/97‑1 (MQ=255)
                                         ccccAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCa                                                      >  2:188342/1‑100 (MQ=255)
                                         ccccAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCa                                                      >  2:2509626/1‑100 (MQ=255)
                                          cccAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATACAAATGGTTGACATAAAATTGCCAt                                                     >  1:3078173/1‑100 (MQ=255)
                                          cccAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCAt                                                     >  1:3339039/1‑100 (MQ=255)
                                           ccAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATa                                                    >  1:2465664/1‑100 (MQ=255)
                                           ccAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATa                                                    >  1:471864/1‑100 (MQ=255)
                                            cAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCaaga                                                   >  1:4097861/1‑97 (MQ=255)
                                            cAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATaa                                                   >  2:1000452/1‑100 (MQ=255)
                                              aaaaaCGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGCTGACATAAACTTGCCATAAGc                                                 >  2:2847688/1‑100 (MQ=255)
                                                   cGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGt                                              >  2:185625/1‑98 (MQ=255)
                                                   cGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTgg                                            <  1:128478/100‑1 (MQ=255)
                                                    ggCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTggg                                           >  1:1289417/1‑100 (MQ=255)
                                                       tGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAgg                                        <  1:3591353/100‑1 (MQ=255)
                                                       tGACAAAGCTTTCAAACTGACGGAGGAAAATTTCACCCACATTGGTCTTAACCCATTGCGTAAAAATGGTTTACATAAACTTGCCATAGGCCGTGGGAgg                                        <  2:1216360/100‑1 (MQ=255)
                                                        gACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGc                                       <  1:210167/100‑1 (MQ=255)
                                                        gACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGc                                       <  1:2509641/100‑1 (MQ=255)
                                                         aCAAAGCTTTCAAACTGACGGATGAAAATTTTATCCACATTGTTTTTAAAAAAAAAAAAAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGc                                       <  2:2319353/99‑1 (MQ=255)
                                                         aCAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGca                                      <  2:1022756/100‑1 (MQ=255)
                                                          cAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGcac                                     >  1:419443/1‑100 (MQ=255)
                                                          cAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGcac                                     >  1:962810/1‑100 (MQ=255)
                                                            aaGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAg                                   >  1:2107611/1‑100 (MQ=255)
                                                            aaGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAg                                   <  2:2041895/100‑1 (MQ=255)
                                                            aaGCTTTCAAACTAACGGATGAAAATTTCACGCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAg                                   >  1:2813667/1‑100 (MQ=255)
                                                               cTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGaaa                                >  1:2255556/1‑100 (MQ=255)
                                                               cTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGaaa                                >  2:457966/1‑100 (MQ=255)
                                                                    aaaCTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCa                           <  2:1599438/100‑1 (MQ=255)
                                                                           cGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATAc                    >  1:1567458/1‑100 (MQ=255)
                                                                            aaataaaaaTTTCTCGTACATTTTGTTTAACCCATTGCATAAAAATGGCTAACATAAACTTGCCATAAGCCGTGAGTTTCACAGAAAAATCAATGacac                    <  2:4085912/94‑1 (MQ=255)
                                                                               tGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGaa                >  1:2379546/1‑100 (MQ=255)
                                                                                gAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACGTAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAAt               >  1:2116452/1‑100 (MQ=255)
                                                                                 aaaaTTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGAGAGGCACAGAAAAATCAATGATACGGAAtg              >  2:3108142/1‑100 (MQ=255)
                                                                                    aTTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATgtgt           >  1:880237/1‑100 (MQ=255)
                                                                                      atCACCAACATTCTTTTTAACAAAATGAATAAAAACGGATGACATAAACATGCAATCAGACGTGAAAAGCACAGAAAAATCAATGATACGGaaagaaaac         <  2:3749303/99‑8 (MQ=255)
                                                                                         aCCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCa      <  2:1045760/100‑1 (MQ=255)
                                                                                         aCCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCa      <  2:1045542/100‑1 (MQ=255)
                                                                                         aCCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCa      <  1:1957417/100‑1 (MQ=255)
                                                                                          cccACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCAc     >  1:783226/1‑100 (MQ=255)
                                                                                          cccACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGACGCACAGAAAAATCAATGATACGGAATGTGTTCTCAc     <  1:2943823/100‑1 (MQ=255)
                                                                                           ccACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACt    <  2:12877/100‑1 (MQ=255)
                                                                                            cacaTTGTTTTTAACCCATTGCATACAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATCCGGAATGTGTTCTCACTa   >  2:3869361/1‑100 (MQ=255)
                                                                                            cacaTTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTa   <  2:155528/100‑1 (MQ=255)
                                                                                            cacaTTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTa   <  1:1979905/100‑1 (MQ=255)
                                                                                            cacaTTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTa   >  1:473578/1‑100 (MQ=255)
                                                                                                 |                                                                                               
TCCTTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTAT  >  NC_000913/1579997‑1580189

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