Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_bme_pgi 4,182,934 C→T P564S (CCT→TCT)  rpoB → RNA polymerase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi4,182,9340CT98.5% 668.3 / ‑3.0 198P564S (CCT→TCT) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base T (93/102);  total (95/103)
Fisher's exact test for biased strand distribution p-value = 6.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

Only 100 of 200 total aligned reads displayed.

CGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACC  >  NC_000913_3_bme_pgi/4182837‑4183029
                                                                                                 |                                                                                               
 gCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACctct                                                                                               >  2:664503/1‑99 (MQ=255)
    cTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGATACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCTCTGAAg                                                                                           >  2:3135726/1‑100 (MQ=255)
      cGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCTCTGAAGGt                                                                                         >  2:3708531/1‑100 (MQ=255)
       ggCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCTCTGaa                                                                                            <  2:2145110/96‑1 (MQ=255)
       ggCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCTCTGAAgg                                                                                          <  2:4852618/98‑1 (MQ=255)
       ggCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTc                                                                                        >  1:1001731/1‑100 (MQ=255)
       ggCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTc                                                                                        >  1:2378360/1‑100 (MQ=255)
        gCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCTCTGAAgg                                                                                          >  2:2699735/1‑97 (MQ=255)
        gCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTcc                                                                                       <  1:3013182/100‑1 (MQ=255)
        gCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTcc                                                                                       >  1:4394715/1‑100 (MQ=255)
         cccAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCTCCGAAGGTcc                                                                                       >  2:2498534/1‑99 (MQ=255)
           cAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTCCGaa                                                                                    >  1:194460/1‑100 (MQ=255)
            aGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTCCGaa                                                                                    >  1:3473003/1‑99 (MQ=255)
            aGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTCCGaa                                                                                    >  2:1251839/1‑99 (MQ=255)
               cGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTCCGAAc                                                                                   <  1:2629599/97‑1 (MQ=255)
               cGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTCCGAACATc                                                                                <  2:72922/100‑1 (MQ=255)
                 gTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTCCGAACATCgg                                                                              <  1:1734414/100‑1 (MQ=255)
                 gTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTCCGAACATCgg                                                                              <  2:3758129/100‑1 (MQ=255)
                 gTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACACCgg                                                                              <  1:3797892/100‑1 (MQ=255)
                  tCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTCCGAACATCGGt                                                                             >  1:4529955/1‑100 (MQ=255)
                   cTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCTCTGACGGTCCGAACATCGGt                                                                             >  2:5040814/1‑99 (MQ=255)
                    tGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACGTCTGAAGGTCCGAACATCGGTCt                                                                           >  1:5045279/1‑100 (MQ=255)
                     gACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTCCGAACATCGGTCt                                                                           >  1:356932/1‑99 (MQ=255)
                      aCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTCCGAACATCGGTc                                                                            <  1:3965873/97‑1 (MQ=255)
                      aCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTCCGAACATCGGTCTg                                                                          <  2:77705/99‑1 (MQ=255)
                       cccGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTCCGAACATCGGTCTGAt                                                                        <  1:1653502/100‑1 (MQ=255)
                        ccGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTCCGAACATCGGTCTGATc                                                                       <  1:53641/100‑1 (MQ=255)
                        ccGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTCCGAACATCGGTCTGATc                                                                       <  1:220082/100‑1 (MQ=255)
                           tGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCCATCGAACCCTCTGAAGGTCCGAACATCGGTCTGATCAAc                                                                    <  2:4791885/100‑1 (MQ=255)
                             aaCGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTCCGAACATCGGTCTGATCAActc                                                                  <  1:76852/100‑1 (MQ=255)
                                gTGCAGGCTTCGAAGTTCGAGACGTGCACCCGACTCACTACGGTCGCGTATGTCCAAGCGAAACCTCTGAAGGTCCGAACATCGGTCTGATCAActctc                                                                <  1:1770305/99‑1 (MQ=255)
                                gTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTCCGAACATCGGTCTGATCAActctc                                                                <  2:3395782/99‑1 (MQ=255)
                                   cAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGt                                                             <  1:4285703/99‑1 (MQ=255)
                                   cAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTc                                                            <  1:3569269/100‑1 (MQ=255)
                                    aGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTcc                                                           <  1:2830456/100‑1 (MQ=255)
                                     ggCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTCCGACCATCGGTCTGATCAACTCTCTGTcc                                                           >  2:1114473/1‑99 (MQ=255)
                                     ggCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCg                                                          <  1:1540048/100‑1 (MQ=255)
                                       cTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTCCGAACATCGGTCTGATCACCTCTCTGTCCgtg                                                        >  2:4001829/1‑100 (MQ=255)
                                       cTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCgtg                                                        >  1:375316/1‑100 (MQ=255)
                                       cTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCgt                                                         >  1:2214992/1‑99 (MQ=255)
                                       cTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTCCGAACATCGGTCTGATCAACTCGCTGTCCgtg                                                        >  1:1111324/1‑100 (MQ=255)
                                           gAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACg                                                    >  2:237859/1‑100 (MQ=255)
                                             aGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGca                                                  <  2:4751059/100‑1 (MQ=255)
                                             aGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGca                                                  <  1:2635017/100‑1 (MQ=255)
                                              gTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGcac                                                 <  1:3846565/100‑1 (MQ=255)
                                              gTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGcac                                                 <  2:5023669/100‑1 (MQ=255)
                                               ttCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGca                                                  >  2:774941/1‑98 (MQ=255)
                                                 cGAGTCGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGa                                              >  2:4414780/1‑100 (MQ=255)
                                                  gagaCGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGa                                              >  2:3597327/1‑99 (MQ=255)
                                                   agaCGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACt                                            <  1:510376/100‑1 (MQ=255)
                                                   agaCGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACt                                            <  1:3187599/100‑1 (MQ=255)
                                                   agaCGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACt                                            <  1:486712/100‑1 (MQ=255)
                                                       gTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACg                                        <  1:2578750/100‑1 (MQ=255)
                                                       gTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACg                                        <  1:1422112/100‑1 (MQ=255)
                                                          cacCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGaa                                      >  2:4814718/1‑99 (MQ=255)
                                                          cacCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAAt                                     <  2:529795/100‑1 (MQ=255)
                                                          cacCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACACACAGACTAACGaa                                      <  1:4710913/99‑1 (MQ=255)
                                                           acCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATa                                    >  1:637893/1‑100 (MQ=255)
                                                           acCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATa                                    >  1:4138783/1‑100 (MQ=255)
                                                            cccGACTCACTACGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATAc                                   >  1:1498551/1‑100 (MQ=255)
                                                            cccGACTCACTACGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATAc                                   >  2:1429690/1‑100 (MQ=255)
                                                            cccGACTCACTACGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATAc                                   >  1:254294/1‑100 (MQ=255)
                                                             ccGACTCACTACGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTCGGAACATCGGTCTGATCAACGCTCTGTCCGTGTACGCACAGACTAACGAATAc                                   <  2:498431/99‑1 (MQ=255)
                                                              cGACTCACTACGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACgg                                 >  2:2312699/1‑100 (MQ=255)
                                                              cGACTCACTACGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACgg                                 <  1:3640497/100‑1 (MQ=255)
                                                                   cACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTcc                            >  2:2087129/1‑100 (MQ=255)
                                                                    aCTACGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCt                           >  2:4503705/1‑100 (MQ=255)
                                                                     cTACGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCtt                          <  1:1675395/100‑1 (MQ=255)
                                                                       aCGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCtt                          >  2:3817711/1‑98 (MQ=255)
                                                                       aCGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGGGTACGCACAGACTAACGAATACGGCTTCCt                           >  1:2106541/1‑97 (MQ=255)
                                                                        cGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCtt                          >  1:504966/1‑97 (MQ=255)
                                                                        cGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTgag                       <  1:3600061/100‑1 (MQ=255)
                                                                        cGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTg                         <  2:5176137/98‑1 (MQ=255)
                                                                        cGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTg                         <  1:1714039/98‑1 (MQ=255)
                                                                        cGGTCGCGTATGTCCAATCGAAACCTCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTg                         <  1:5061052/98‑1 (MQ=255)
                                                                          gTCGCGTATGTCCAATCGAAACCTCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGAc                     <  2:254392/100‑1 (MQ=255)
                                                                             gcgTATGTCCAATCGAAACCTCTGAAGGTCCGAACATCCGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTc                             >  2:2170707/1‑89 (MQ=255)
                                                                              cgTATGTCCAATCGAAACCTCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCg                 >  2:464168/1‑100 (MQ=255)
                                                                                tATGTCCAATCGAAACCTTTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGGCCGTGTACGCACTGATTAACGAATACGGCGTCCTTgag                       <  2:4738602/92‑1 (MQ=255)
                                                                                tATGTCCAATCGAAACCTCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTa               <  2:1749446/100‑1 (MQ=255)
                                                                                   gTCCAATCGAAACCTCTGAAGGTCCGACCATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCg            >  1:4052578/1‑100 (MQ=255)
                                                                                   gTCCAATCGAAACCTCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCg            >  1:3366711/1‑100 (MQ=255)
                                                                                   gTCCAATCGAAACCTCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCg            >  1:1693598/1‑100 (MQ=255)
                                                                                    tCCAATCGAAACCTCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCg            <  1:285797/99‑1 (MQ=255)
                                                                                    tCCAATCGAAACCTCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCg            >  1:2495060/1‑99 (MQ=255)
                                                                                     ccAATCGAAACCTCTGAAGGTACGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCtt                          >  2:4118012/1‑84 (MQ=255)
                                                                                      cAATCGAAACCTCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGt           <  2:3451507/98‑1 (MQ=255)
                                                                                      cAATCGAAACCTCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACCCCGTATCGTa          >  1:808389/1‑99 (MQ=255)
                                                                                        aTCGAAACCTCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAg       <  1:3902632/100‑1 (MQ=255)
                                                                                        aTCGAAACCTCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAg       <  2:490857/100‑1 (MQ=255)
                                                                                         tCGAAACCTCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGt      >  2:4146908/1‑100 (MQ=255)
                                                                                          cGAAACCTCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGt      >  2:1721190/1‑99 (MQ=255)
                                                                                          cGAAACCTCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGt      >  2:2902314/1‑99 (MQ=255)
                                                                                           gAAACCTCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTTTGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTg     <  2:5123831/99‑1 (MQ=255)
                                                                                           gAAACCTCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAg       <  2:2966951/97‑1 (MQ=255)
                                                                                           gAAACCTCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGa    >  2:2023605/1‑100 (MQ=255)
                                                                                            aaaCCTCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGAc   >  1:5322871/1‑100 (MQ=255)
                                                                                             aaCCTCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGAcc  <  2:2863927/100‑1 (MQ=255)
                                                                                             aaCCTCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGAcc  <  1:1301354/100‑1 (MQ=255)
                                                                                             aaCCTCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGAcc  <  1:4551935/100‑1 (MQ=255)
                                                                                                 |                                                                                               
CGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACC  >  NC_000913_3_bme_pgi/4182837‑4183029

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: