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,6310AG97.6% 783.5 / 2.2 251E90G (GAA→GGA) ebgCevolved beta‑D‑galactosidase beta subunit
Reads supporting (aligned to +/- strand):  ref base A (2/4);  new base G (121/124);  total (123/128)
Fisher's exact test for biased strand distribution p-value = 6.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 253 total aligned reads displayed.

CCGGGCATCGACGCTATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGAATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATA  >  NC_000913/3223541‑3223720
                                                                                          |                                                                                         
 cGGGCATCGACGCTATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGatt                                                                                        >  1:1414989/1‑92 (MQ=255)
  gggTATCGACGCTATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGatat                                                                                       >  1:3586247/1‑93 (MQ=255)
  gggCATCGACGCTATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGatat                                                                                       >  1:4468459/1‑93 (MQ=255)
  gggCATCGACGCTATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGatat                                                                                       >  1:1561192/1‑93 (MQ=255)
    gCATCGACGCTATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGAtatta                                                                                     <  2:3449596/93‑1 (MQ=255)
     cATCGACGCTATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGAtattat                                                                                    <  1:2615403/93‑1 (MQ=255)
      aTCGACGCTATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATc                                                                                   >  1:4701977/1‑93 (MQ=255)
      aTCGACGCTATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATc                                                                                   >  2:2991821/1‑93 (MQ=255)
      aTCGACGCTATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATc                                                                                   <  2:3145019/93‑1 (MQ=255)
      aTCGACGCTATTTTGAAGTGCATAACTACCTGCAAGGGCAGCAAAGAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGAtattaac                                                                                   <  1:4209417/93‑3 (MQ=255)
       tCGACGCTATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCg                                                                                  >  1:4566905/1‑93 (MQ=255)
       tCGACGCTATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCg                                                                                  >  2:1293360/1‑93 (MQ=255)
          aCGCTATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGtga                                                                               >  1:1358392/1‑93 (MQ=255)
          aCGCTATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGtga                                                                               >  1:4165746/1‑93 (MQ=255)
           cGCTATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGtgat                                                                              <  2:2952653/93‑1 (MQ=255)
            gCTATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGtgatg                                                                             >  2:794773/1‑93 (MQ=255)
              tATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGaa                                                                           >  2:1448677/1‑93 (MQ=255)
               aTTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGaaa                                                                          >  1:198351/1‑93 (MQ=255)
                ttttGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAAc                                                                         <  2:3269862/93‑1 (MQ=255)
                ttttGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAAc                                                                         <  2:307360/93‑1 (MQ=255)
                 tttGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACt                                                                        >  2:3354086/1‑93 (MQ=255)
                 tttGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACt                                                                        <  2:2563621/93‑1 (MQ=255)
                   tGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGa                                                                      <  2:1813908/93‑1 (MQ=255)
                      aGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCg                                                                   <  2:1053313/93‑1 (MQ=255)
                       gTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGt                                                                  >  2:486014/1‑93 (MQ=255)
                         gCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGa                                                                <  2:1232810/93‑1 (MQ=255)
                         gCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGa                                                                <  1:1296126/93‑1 (MQ=255)
                         gCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGa                                                                >  1:4357307/1‑93 (MQ=255)
                            ttACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAata                                                             >  2:4605572/1‑93 (MQ=255)
                              actacCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTCATGAAACTGACCGTGAATAtt                                                           >  1:3315006/1‑93 (MQ=255)
                               ctacCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATAttt                                                          >  1:265503/1‑93 (MQ=255)
                               ctacCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATAttt                                                          >  2:2370273/1‑93 (MQ=255)
                                tacCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTa                                                         >  1:3889298/1‑93 (MQ=255)
                                 acCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTaa                                                        <  1:1672623/93‑1 (MQ=255)
                                 acCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTaa                                                        >  1:658507/1‑93 (MQ=255)
                                   cTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCATGAATATTTaaaa                                                      <  2:4076348/93‑1 (MQ=255)
                                    tGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAg                                                     >  2:734756/1‑93 (MQ=255)
                                     gCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAgg                                                    <  1:123920/93‑1 (MQ=255)
                                     gCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAgg                                                    <  2:3688875/93‑1 (MQ=255)
                                      cAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGc                                                   <  1:1774705/93‑1 (MQ=255)
                                      cAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGc                                                   >  1:3812153/1‑93 (MQ=255)
                                      cAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGc                                                   <  2:2207075/93‑1 (MQ=255)
                                       aaGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCt                                                  >  2:4315738/1‑93 (MQ=255)
                                        aGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTg                                                 >  1:1624667/1‑93 (MQ=255)
                                        aGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTg                                                 >  1:897525/1‑93 (MQ=255)
                                        aGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTg                                                 >  1:2748417/1‑93 (MQ=255)
                                         gggCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGc                                                >  2:260518/1‑93 (MQ=255)
                                            cagcaAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGga                                             <  2:56487/93‑1 (MQ=255)
                                            cagcaAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGga                                             >  2:3195418/1‑93 (MQ=255)
                                             agcaAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGgag                                            >  1:2695115/1‑93 (MQ=255)
                                              gcaAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGgaga                                           <  1:2275993/93‑1 (MQ=255)
                                               gaaaaaaTTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATGATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGaa                                          <  1:3104498/92‑1 (MQ=255)
                                               caAAAAATTGAATATGCGCAGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGaa                                          <  1:3288236/93‑1 (MQ=255)
                                                 aaaaaTTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAAc                                        <  2:2835591/93‑1 (MQ=255)
                                                 aaaaaTTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAAc                                        >  2:976097/1‑93 (MQ=255)
                                                  aaaaTTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAAcc                                       >  1:560178/1‑93 (MQ=255)
                                                     aTTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGtt                                    >  1:2567933/1‑93 (MQ=255)
                                                     aTTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGtt                                    >  2:4353337/1‑93 (MQ=255)
                                                          atatGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGt                               >  2:2455616/1‑93 (MQ=255)
                                                          atatGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGt                               <  2:4243123/93‑1 (MQ=255)
                                                          atatGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGt                               >  1:1784816/1‑93 (MQ=255)
                                                          atatGCGCCGAAAGAGACATTACAGGTAGTGGAATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGt                               <  2:1584552/93‑1 (MQ=255)
                                                              gcgcCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCAc                           <  1:2380170/93‑1 (MQ=255)
                                                               cgcCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACg                          >  1:3730095/1‑93 (MQ=255)
                                                                 ccGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGaa                        <  1:2112944/93‑1 (MQ=255)
                                                                  cGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAg                       >  1:1197160/1‑93 (MQ=255)
                                                                  cGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAg                       <  1:4508404/93‑1 (MQ=255)
                                                                  cGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAg                       <  1:2930246/93‑1 (MQ=255)
                                                                   gAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAgg                      <  2:2055328/93‑1 (MQ=255)
                                                                   gAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAgg                      >  2:1728451/1‑93 (MQ=255)
                                                                     aaGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAggga                    >  1:2468911/1‑92 (MQ=255)
                                                                     aaGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGc                    >  2:879470/1‑93 (MQ=255)
                                                                     aaGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGc                    <  1:1034527/93‑1 (MQ=255)
                                                                     aaGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGc                    >  1:3191170/1‑93 (MQ=255)
                                                                     aaGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGc                    <  1:3780063/93‑1 (MQ=255)
                                                                     aaGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGc                    >  1:508050/1‑93 (MQ=255)
                                                                      agagACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCa                   <  1:235703/93‑1 (MQ=255)
                                                                      agagACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCa                   >  2:3562074/1‑93 (MQ=255)
                                                                       gagaCATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCaa                  >  1:2940496/1‑93 (MQ=255)
                                                                       gagaCATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCaa                  >  1:3251273/1‑93 (MQ=255)
                                                                       gagaCATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGGTGAGGTCCACGAAGGGCaa                  >  2:4207498/1‑93 (MQ=255)
                                                                        agaCATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCaaa                 >  1:394721/1‑93 (MQ=255)
                                                                        agaCATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCaaa                 <  2:3822003/93‑1 (MQ=255)
                                                                        agaCATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCaaa                 <  2:2121751/93‑1 (MQ=255)
                                                                         gaCATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAAt                <  2:3274295/93‑1 (MQ=255)
                                                                         gaCATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAAt                <  1:818078/93‑1 (MQ=255)
                                                                          aCATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATc               <  2:1218747/93‑1 (MQ=255)
                                                                          aCATTACAGGTAGTGGGATATTATCGTGAGGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATc               <  2:2211523/93‑1 (MQ=255)
                                                                           cATTACAGGTAGTAGAATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCg              <  1:1853299/93‑1 (MQ=255)
                                                                            aTTGCAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGt             <  2:1907648/93‑1 (MQ=255)
                                                                              tACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTa           <  2:3886179/93‑1 (MQ=255)
                                                                              tACAGGTAGTGCGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAACCCGTTGAGGTCCACGAAGGGCAAATCGTTa           <  2:2913795/93‑1 (MQ=255)
                                                                               aCAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTAt          >  2:1755592/1‑93 (MQ=255)
                                                                               aCAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTAt          >  1:2056760/1‑93 (MQ=255)
                                                                                 aGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTAttt        <  2:1380677/93‑1 (MQ=255)
                                                                                  ggTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTATTTg       <  2:4476781/93‑1 (MQ=255)
                                                                                     aGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGa    >  2:629940/1‑93 (MQ=255)
                                                                                     aGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGa    <  1:2192911/93‑1 (MQ=255)
                                                                                      gTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGat   <  2:4460198/93‑1 (MQ=255)
                                                                                       tGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGata  <  2:701977/93‑1 (MQ=255)
                                                                                          |                                                                                         
CCGGGCATCGACGCTATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGAATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATA  >  NC_000913/3223541‑3223720

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