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.7% 950.9 / NA 290E90G (GAA→GGA) ebgCevolved beta‑D‑galactosidase beta subunit
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (155/134);  total (155/135)
Fisher's exact test for biased strand distribution p-value = 4.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 292 total aligned reads displayed.

CCGGGCATCGACGCTATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGAATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATA  >  NC_000913/3223541‑3223720
                                                                                          |                                                                                         
ccGGGCATCGACGCTATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGat                                                                                         >  1:2868448/1‑93 (MQ=255)
 cGGGCATCGACGCTATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAACAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGata                                                                                        <  2:4701383/93‑1 (MQ=255)
      aTCGACGCTATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATc                                                                                   <  2:323585/93‑1 (MQ=255)
       tCGACGCTATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCg                                                                                  >  1:3682799/1‑93 (MQ=255)
       tCGACGCTATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCg                                                                                  >  2:3114778/1‑93 (MQ=255)
        cGACGCTATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGt                                                                                 <  2:4073869/93‑1 (MQ=255)
        cGACGCTATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGt                                                                                 >  2:1972226/1‑93 (MQ=255)
        cGACGCTATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGt                                                                                 <  1:5135897/93‑1 (MQ=255)
         gACGCTATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGtg                                                                                >  2:6467830/1‑93 (MQ=255)
         gACGCTATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGtg                                                                                >  1:3192559/1‑93 (MQ=255)
         gACGCTATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGtg                                                                                <  1:5158401/93‑1 (MQ=255)
          aCGCTATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGtga                                                                               >  2:1617187/1‑93 (MQ=255)
           cGCTATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGtgat                                                                              <  2:4764854/93‑1 (MQ=255)
           cGCTATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGtgat                                                                              >  2:3891629/1‑93 (MQ=255)
           cGCTATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGtgat                                                                              <  2:4118904/93‑1 (MQ=255)
              tATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGaa                                                                           >  2:2480905/1‑93 (MQ=255)
              tATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGaa                                                                           >  1:1947910/1‑93 (MQ=255)
               aTTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGaaa                                                                          >  1:2172531/1‑93 (MQ=255)
                ttttGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAAc                                                                         <  2:2934618/93‑1 (MQ=255)
                ttttGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAAc                                                                         <  2:4246797/93‑1 (MQ=255)
                ttttGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAAc                                                                         >  1:1191573/1‑93 (MQ=255)
                 tttGAAGTGCATTAATACCTGCGAGGGCAGCAATAACTTGAATATGCAAAGAAAGAGACATTACAGGAAGTGGACTATTATCGTGATGAAACt                                                                        <  1:6725025/93‑1 (MQ=255)
                  ttGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTg                                                                       <  2:615317/93‑1 (MQ=255)
                   tGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGa                                                                      <  2:6297138/93‑1 (MQ=255)
                   tGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGa                                                                      <  1:5696677/93‑1 (MQ=255)
                      aGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGGGCCGAAAGAGACATTACAGGTAGGGGGATATTATCTTGATGAAACTGACCg                                                                   >  1:2430465/1‑93 (MQ=255)
                       gTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGCAAGAGACATTACAGGTAGCGGGCTATTATCGTGCTGAAACTGACCgg                                                                  >  2:5393371/1‑92 (MQ=255)
                       gTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGt                                                                  >  1:4324744/1‑93 (MQ=255)
                          cATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGaa                                                               >  2:3040068/1‑93 (MQ=255)
                          cATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGaa                                                               <  2:2143454/93‑1 (MQ=255)
                           aTTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAat                                                              >  2:1324330/1‑93 (MQ=255)
                           aTTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAat                                                              <  2:2376826/93‑1 (MQ=255)
                           aTTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAat                                                              >  1:5995065/1‑93 (MQ=255)
                           aTTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAat                                                              >  2:4842448/1‑93 (MQ=255)
                              actacCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATAtt                                                           <  2:1588902/93‑1 (MQ=255)
                               ctacCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATAttt                                                          <  1:4414180/93‑1 (MQ=255)
                                tacCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTa                                                         >  1:3359608/1‑93 (MQ=255)
                                tacCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTa                                                         <  1:5818419/93‑1 (MQ=255)
                                tacCTGCAAGGGCAGCAAAAAATTGAATATGAGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACTGTGAATATTTa                                                         >  2:6648490/1‑93 (MQ=255)
                                    tGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAg                                                     <  2:5569421/93‑1 (MQ=255)
                                     gCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAgg                                                    <  1:1781457/93‑1 (MQ=255)
                                      cAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGc                                                   <  2:5516415/93‑1 (MQ=255)
                                      cAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGc                                                   <  2:2314444/93‑1 (MQ=255)
                                      cAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGc                                                   <  1:1367697/93‑1 (MQ=255)
                                      cAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGc                                                   <  1:6282075/93‑1 (MQ=255)
                                      cAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGc                                                   <  2:4014915/93‑1 (MQ=255)
                                       aaGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGcg                                                  >  2:4369083/1‑92 (MQ=255)
                                       aaGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCt                                                  <  1:4755447/93‑1 (MQ=255)
                                        aGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTg                                                 >  1:4976420/1‑93 (MQ=255)
                                        aGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTg                                                 >  2:1374381/1‑93 (MQ=255)
                                         gggCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGc                                                <  2:2491925/93‑1 (MQ=255)
                                         gggCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGc                                                <  1:2281764/93‑1 (MQ=255)
                                         gggCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGc                                                >  1:391194/1‑93 (MQ=255)
                                         gggCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGc                                                >  2:4636566/1‑93 (MQ=255)
                                         gggCAGCAAAAAATTGAATATGCGCAGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGTATATTTAAAAGGCTGc                                                <  1:4558041/93‑1 (MQ=255)
                                           gcagcaAAAAATTGAATATGCGCCGAAAGAGACATTACAAGTAGTGGGACATTATCGTGATGAAACTGACCGTGAATATTGCAAAGGCCGCgg                                              >  1:4251834/1‑93 (MQ=255)
                                            cagcaAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGga                                             >  2:4486954/1‑93 (MQ=255)
                                            cagcaAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGga                                             <  2:1863/93‑1 (MQ=255)
                                             agcaAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGgag                                            >  1:4091626/1‑93 (MQ=255)
                                              gcgaaaaaTTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGgaga                                           >  2:5081453/4‑93 (MQ=255)
                                               caAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGaa                                          <  2:5843389/93‑1 (MQ=255)
                                               caAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGaa                                          >  1:5579704/1‑93 (MQ=255)
                                                 aaaaaTTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAAc                                        >  1:602720/1‑93 (MQ=255)
                                                   aaaTTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCg                                      >  2:931595/1‑93 (MQ=255)
                                                   aaaTTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCg                                      >  1:2735955/1‑93 (MQ=255)
                                                     aTTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGtt                                    >  2:2587983/1‑93 (MQ=255)
                                                      ttGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTg                                   >  1:6149923/1‑93 (MQ=255)
                                                      ttGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTg                                   >  2:861573/1‑93 (MQ=255)
                                                       tGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGa                                  <  1:341801/93‑1 (MQ=255)
                                                       tGAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGa                                  <  1:1250484/93‑1 (MQ=255)
                                                        gAATATGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAg                                 >  1:3264884/1‑93 (MQ=255)
                                                        gAATATACGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAg                                 <  2:901768/93‑1 (MQ=255)
                                                          atatGCGCCGAAAGTGACATTACAGGCAGTGGGATATTAACGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGATAAACCGTTGAggc                               >  1:3937727/1‑92 (MQ=255)
                                                           tatGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTACAAGGCTGCGGAGAAACCGTTGAGGTc                              >  1:2961731/1‑93 (MQ=255)
                                                           tatGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTc                              >  2:1231030/1‑93 (MQ=255)
                                                            atGCGCCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTcc                             <  1:5287044/93‑1 (MQ=255)
                                                              gcgcCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCAc                           >  2:3260114/1‑93 (MQ=255)
                                                                gcCGAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGa                         >  1:4692215/1‑93 (MQ=255)
                                                                   gAAAGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAgg                      >  2:374777/1‑93 (MQ=255)
                                                                    aaaGAGACATTACAGGTAGTGGGATATTATTGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAggg                     <  1:3123240/93‑1 (MQ=255)
                                                                    aaaGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAggg                     >  1:1515301/1‑93 (MQ=255)
                                                                    aaaGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAggg                     >  2:1672406/1‑93 (MQ=255)
                                                                     aaGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGc                    >  2:6013788/1‑93 (MQ=255)
                                                                     aaGAGACATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGc                    >  2:350517/1‑93 (MQ=255)
                                                                        agaCATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCaaa                 <  1:1111444/93‑1 (MQ=255)
                                                                        agaCATTACAGGTAGTGAGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCaaa                 <  1:6197596/93‑1 (MQ=255)
                                                                         gaCATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAAt                >  1:6369390/1‑93 (MQ=255)
                                                                         gaCATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAAt                >  1:888332/1‑93 (MQ=255)
                                                                         gaCATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAAt                <  2:114907/93‑1 (MQ=255)
                                                                          aCATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATc               >  1:6583194/1‑93 (MQ=255)
                                                                          aCATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATc               <  2:3322119/93‑1 (MQ=255)
                                                                           cATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCg              >  1:3627899/1‑93 (MQ=255)
                                                                           cATTACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCg              >  1:2587413/1‑93 (MQ=255)
                                                                              tACAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTa           <  1:416851/93‑1 (MQ=255)
                                                                               aCAGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTAt          <  1:1290703/93‑1 (MQ=255)
                                                                                 aGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTAttt        <  2:5423409/93‑1 (MQ=255)
                                                                                 aGGTAGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTAttt        <  1:2036059/93‑1 (MQ=255)
                                                                                     aGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGa    <  2:711395/93‑1 (MQ=255)
                                                                                     aGTGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCCGCGGAGAAACCGTTGCGGTCCACGAAGGGCAAATCGTTATTTGCGa    <  1:3858336/93‑1 (MQ=255)
                                                                                       tGGGATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGata  <  2:4971338/93‑1 (MQ=255)
                                                                                          |                                                                                         
CCGGGCATCGACGCTATTTTGAAGTGCATTACTACCTGCAAGGGCAGCAAAAAATTGAATATGCGCCGAAAGAGACATTACAGGTAGTGGAATATTATCGTGATGAAACTGACCGTGAATATTTAAAAGGCTGCGGAGAAACCGTTGAGGTCCACGAAGGGCAAATCGTTATTTGCGATA  >  NC_000913/3223541‑3223720

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