Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,223,631 A→G E90G (GAA→GGA)  ebgC → evolved beta‑D‑galactosidase beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,223,6310AG99.5% 1277.7 / ‑3.9 404E90G (GAA→GGA) ebgCevolved beta‑D‑galactosidase beta subunit
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (190/212);  total (191/213)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.74e-01

Only 100 of 411 total aligned reads displayed.

TACCGGGCATCGACGCTATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGAATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATAT  >  NC_000913/3223539‑3223721
                                                                                            |                                                                                          
tACCGGGCATCGACGCTATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGa                                                                                            <  1:2625059/93‑1 (MQ=255)
  ccGGGCATCGACGCTATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGat                                                                                          <  2:282699/93‑1 (MQ=255)
   cGGGCATCGACGCTATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGata                                                                                         <  1:4581858/93‑1 (MQ=255)
    gggCATCGACGCTATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGatat                                                                                        >  2:2890796/1‑93 (MQ=255)
    gggCATCGACGCTATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGatat                                                                                        >  1:4146049/1‑93 (MQ=255)
    gggCATCGACGCTATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGatat                                                                                        <  1:2009623/93‑1 (MQ=255)
     ggCATCGACGCTATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGAtatt                                                                                       >  2:2178341/1‑93 (MQ=255)
     ggCATCGACGCTATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGAtatt                                                                                       <  2:4511527/93‑1 (MQ=255)
     ggCATCGACGCTATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGAtatt                                                                                       <  1:4903294/93‑1 (MQ=255)
      gCATCGACGCTATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGAtatta                                                                                      >  2:1592092/1‑93 (MQ=255)
      gCATCGACGCTATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGAtatta                                                                                      <  1:7669300/93‑1 (MQ=255)
       cATCGACGCTATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGAtattat                                                                                     >  1:6731727/1‑93 (MQ=255)
        aTCGACGCTATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATc                                                                                    <  2:7734731/93‑1 (MQ=255)
        aTCGACGCTATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATc                                                                                    >  1:5535303/1‑93 (MQ=255)
         tCGACGCTATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCg                                                                                   >  1:1622694/1‑93 (MQ=255)
          cGACGCTATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGt                                                                                  <  2:1747647/93‑1 (MQ=255)
             cGCTATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGtgat                                                                               >  1:5841197/1‑93 (MQ=255)
              gCTATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGtgatg                                                                              >  2:6501655/1‑93 (MQ=255)
              gCTATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGtgatg                                                                              >  2:7661225/1‑93 (MQ=255)
               cTATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGtgatga                                                                             >  1:6915103/1‑93 (MQ=255)
               cTATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGtgatga                                                                             >  2:6041903/1‑93 (MQ=255)
                tATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGaa                                                                            <  1:4278728/93‑1 (MQ=255)
                tATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGaa                                                                            >  2:7135047/1‑93 (MQ=255)
                 aTTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGaaa                                                                           <  2:2485323/93‑1 (MQ=255)
                 aTTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGaaa                                                                           >  2:6186917/1‑93 (MQ=255)
                 aTTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGaaa                                                                           <  1:5656414/93‑1 (MQ=255)
                  ttttGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAAc                                                                          <  1:2177821/93‑1 (MQ=255)
                   tttGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACt                                                                         >  2:4213973/1‑93 (MQ=255)
                       aaGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGAcc                                                                     >  1:3560646/1‑93 (MQ=255)
                        aGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCg                                                                    >  2:5712447/1‑93 (MQ=255)
                         gTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGt                                                                   >  2:357018/1‑93 (MQ=255)
                         gTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGt                                                                   >  2:3655848/1‑93 (MQ=255)
                         gTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGt                                                                   >  1:730326/1‑93 (MQ=255)
                         gTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGt                                                                   <  1:2973623/93‑1 (MQ=255)
                             aTTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAat                                                               >  2:126366/1‑93 (MQ=255)
                             aTTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAat                                                               >  1:3908892/1‑93 (MQ=255)
                              ttACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAata                                                              <  2:7933153/93‑1 (MQ=255)
                              ttACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCATGAata                                                              >  2:4800903/1‑93 (MQ=255)
                                  tacCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTa                                                          >  2:6876798/1‑93 (MQ=255)
                                  tacCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTa                                                          >  1:70230/1‑93 (MQ=255)
                                     cTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTaaaa                                                       <  1:5904563/93‑1 (MQ=255)
                                     cTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTaaaa                                                       >  2:2461225/1‑93 (MQ=255)
                                     cTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTaaaa                                                       <  1:6300488/93‑1 (MQ=255)
                                      tGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAg                                                      <  2:7433064/93‑1 (MQ=255)
                                        cAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGc                                                    >  2:5532718/1‑93 (MQ=255)
                                         aaGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCt                                                   <  2:5504864/93‑1 (MQ=255)
                                          aGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTg                                                  >  1:5018804/1‑93 (MQ=255)
                                          aGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTg                                                  <  2:7466733/93‑1 (MQ=255)
                                           gggCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGc                                                 <  2:3902823/93‑1 (MQ=255)
                                            ggCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCg                                                >  1:102751/1‑93 (MQ=255)
                                             gcagcaAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCgg                                               >  2:6992136/1‑93 (MQ=255)
                                             gcagcaAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCgg                                               <  2:7271240/93‑1 (MQ=255)
                                             gcagcaAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCgg                                               <  1:3790782/93‑1 (MQ=255)
                                              cagcaAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGga                                              <  2:2020282/93‑1 (MQ=255)
                                              cagcaAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGga                                              >  2:5320358/1‑93 (MQ=255)
                                              cagcaAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGga                                              <  2:238982/93‑1 (MQ=255)
                                              cagcaAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGga                                              <  2:5214539/93‑1 (MQ=255)
                                               agcaAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGgag                                             <  2:4853786/93‑1 (MQ=255)
                                                gcaAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGgaga                                            >  2:5100216/1‑93 (MQ=255)
                                                    aaaaTTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAAcc                                        >  2:1632489/1‑93 (MQ=255)
                                                       aTTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGtt                                     >  2:6303640/1‑93 (MQ=255)
                                                       aTTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGtt                                     <  2:6691202/93‑1 (MQ=255)
                                                        ttGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTg                                    >  1:1862968/1‑93 (MQ=255)
                                                        atGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTg                                    <  1:3641781/92‑1 (MQ=255)
                                                         tGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGa                                   <  2:2681654/93‑1 (MQ=255)
                                                         tGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGa                                   <  2:1866672/93‑1 (MQ=255)
                                                         tGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGa                                   <  1:7407210/93‑1 (MQ=255)
                                                          gAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAg                                  <  2:3736730/93‑1 (MQ=255)
                                                           aaTATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAgg                                 <  1:1686407/93‑1 (MQ=255)
                                                           aaTATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAgg                                 <  1:6560143/93‑1 (MQ=255)
                                                           aaTATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAgg                                 >  2:1510354/1‑93 (MQ=255)
                                                            atatGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGt                                <  1:817946/93‑1 (MQ=255)
                                                             tatGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTc                               >  1:5855233/1‑93 (MQ=255)
                                                             tatGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTc                               >  1:4446134/1‑93 (MQ=255)
                                                              atGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTcc                              <  1:279169/93‑1 (MQ=255)
                                                               tGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCa                             >  1:3035676/1‑93 (MQ=255)
                                                                   ccGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGaa                         <  1:523933/93‑1 (MQ=255)
                                                                   ccGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGaa                         <  1:4701661/93‑1 (MQ=255)
                                                                    cGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAg                        >  2:25351/1‑93 (MQ=255)
                                                                       aaGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAggga                     >  1:3181307/1‑92 (MQ=255)
                                                                         gagaCATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCaa                   >  1:4794168/1‑93 (MQ=255)
                                                                          agaCATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCaaa                  <  1:5489744/93‑1 (MQ=255)
                                                                           gaCATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAAt                 <  2:323465/93‑1 (MQ=255)
                                                                             cATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCg               >  1:1454958/1‑93 (MQ=255)
                                                                              aTTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGt              <  1:1148960/93‑1 (MQ=255)
                                                                              aTTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGCTGAGGTCCACGAAGGGCAAATCGt              >  2:4308360/1‑93 (MQ=255)
                                                                               ttACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGtt             >  2:115458/1‑93 (MQ=255)
                                                                               ttACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGtt             >  1:7189977/1‑93 (MQ=255)
                                                                               ttACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGtt             >  1:2362503/1‑93 (MQ=255)
                                                                               ttACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGtt             >  1:1359496/1‑93 (MQ=255)
                                                                               ttACAGGGAGAGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAATGCTGCGGAGAAACCGTTGAGGTCCACGAAGTGCAAATCGtt             >  2:4463036/1‑93 (MQ=255)
                                                                                tACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTa            <  2:5912807/93‑1 (MQ=255)
                                                                                  cAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTAtt          <  1:5436459/93‑1 (MQ=255)
                                                                                   aGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTAttt         >  2:33387/1‑93 (MQ=255)
                                                                                   aGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTAttt         <  1:7737785/93‑1 (MQ=255)
                                                                                     gTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGc       >  1:961837/1‑93 (MQ=255)
                                                                                     gTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGc       <  1:6169098/93‑1 (MQ=255)
                                                                                       aGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGa     <  2:4664365/93‑1 (MQ=255)
                                                                                        gTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGat    <  2:1362357/93‑1 (MQ=255)
                                                                                         tGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGata   >  2:2144484/1‑93 (MQ=255)
                                                                                            |                                                                                          
TACCGGGCATCGACGCTATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGAATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATAT  >  NC_000913/3223539‑3223721

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