Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,224,895 C→T L755L (CTG→TTG)  ebgA → evolved beta‑D‑galactosidase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,224,8950CT98.7% 248.6 / NA 74L755L (CTG→TTG) ebgAevolved beta‑D‑galactosidase subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (41/32);  total (42/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

ATTACGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACTTCGGGATGCGCTGCACCTACATCTGGCGCATCGCTGCCGATGGCCAGGTTAACGT  >  NC_000913/3224655‑3225139
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                    
aTTACGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGa                                                                                                                                                                                                                                            <  1:73026/251‑1 (MQ=255)
       tGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCaa                                                                                                                                                                                                                                       >  1:191687/1‑249 (MQ=255)
       tGACCAAAGATTACCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCac                                                                                                                                                                                                                                       >  1:332631/1‑248 (MQ=255)
        gACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGa                                                                                                                                                                                                                                    <  2:151365/251‑1 (MQ=255)
           cAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATc                                                                                                                                                                                                                                  >  1:160244/1‑250 (MQ=255)
            aaaGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCa                                                                                                                                                                                                                                 >  2:1558/1‑250 (MQ=255)
              aGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTGCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTGCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGTCTCGCGAGCCAAAGATCAAct                                                                                                                                                                                                                              <  1:377407/251‑1 (MQ=255)
                        aCCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGcc                                                                                                                                                                                                                    >  1:334942/1‑251 (MQ=255)
                                           gACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACGCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAg                                                                                                                                                                                                 <  1:293203/251‑1 (MQ=255)
                                              ccccTATCGCCACTTATCAGGTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAAGGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTGCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTGCTTCAAGCCGATGATCGACAACCACAAGCAg                                                                                                                                                                                              <  2:168394/250‑1 (MQ=255)
                                                  tATCGCCACTTATCAGTTACCGCTGAAGGAAAACACCGAGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGt                                                                                                                                                                                          >  2:307167/1‑251 (MQ=255)
                                                   aTCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGt                                                                                                                                                                                          >  2:112977/1‑250 (MQ=255)
                                                         aCTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAgg                                                                                                                                                                                   >  1:126965/1‑251 (MQ=255)
                                                         aCTTACCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCTTTTCGCACCAAATAATGCTCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCGGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAgg                                                                                                                                                                                   <  2:253364/251‑1 (MQ=255)
                                                             aTCAGTTCCCGCTGAAGGAAAACCCCGCGCAGCCAGTGCCTGTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCt                                                                                                                                                                                <  2:124046/250‑1 (MQ=255)
                                                             aTCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCtg                                                                                                                                                                               >  1:65679/1‑251 (MQ=255)
                                                                 gTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGc                                                                                                                                                                           <  1:130810/251‑1 (MQ=255)
                                                                  ttCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAACGATCAACTTCTTCAAGCCGATGATCGACACCCACAAGCAGGAGTACGAAGGGCTGTGGc                                                                                                                                                                           >  2:184611/1‑250 (MQ=255)
                                                                  ttCCCGCTGAAGGAAAACACCGCGCAACCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGAACGAAGGGCtgtagc                                                                                                                                                                           >  2:323027/1‑247 (MQ=255)
                                                                    cccGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAAc                                                                                                                                                                        >  1:285576/1‑251 (MQ=255)
                                                                     ccGCTGAAGGAAAACCCCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGTCAAcc                                                                                                                                                                       <  1:126724/251‑1 (MQ=255)
                                                                     ccGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAAcc                                                                                                                                                                       <  2:65679/251‑1 (MQ=255)
                                                                       gCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCg                                                                                                                                                                      >  2:333440/1‑250 (MQ=255)
                                                                          gAAGGAAAGCACCGAGCAGCAAGAGCCTTTCGAACAGCATAATGCGCGTCCGATGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTAGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGAGGAAACCGGACATCCTGGCAGGTGAAAGGCGAATCGCTGTTGACACGCGAGCCAAAGATCAACTTCTTAAAGACGATGATCGACAACACCAAGCCGGAATACGCAGGGCTGTGGCACCCGAATc                                                                                                                                                                  >  2:277749/1‑251 (MQ=255)
                                                                          gAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCCCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAAt                                                                                                                                                                   <  1:137483/250‑1 (MQ=255)
                                                                            aGGAAAACACCGCGCAGCCAGTGCCGTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGGGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTGCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCAt                                                                                                                                                                <  2:369125/251‑1 (MQ=255)
                                                                                aaaCACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGc                                                                                                                                                            <  1:332066/251‑1 (MQ=255)
                                                                                aaaCACCGCGCAGCCAGTGCCTTTCGCACCAAAGAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCACGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGc                                                                                                                                                            <  1:74593/251‑1 (MQ=255)
                                                                                  acacCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGAGCGTTTGAGCTGCACCGTGCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAg                                                                                                                                                          <  2:383893/251‑1 (MQ=255)
                                                                                                gTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGCCTACAACTTCGCGATCACCTTCTCAAACATGAGTGGCAAACCGACATACTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGAACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAAc                                                                                                                                             >  1:214305/1‑250 (MQ=255)
                                                                                                  gCCTTTCGCACCAAATAATGCGCGTCCGCTGAAGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATc                                                                                                                                          >  1:166418/1‑251 (MQ=255)
                                                                                                   ccTTTCGCCGCAAAGAAGGCGCGCCCGCTGACGCTGGAAGGCGAGCGTTTGGGCTGCACCGTGCGCGGCTACAACTTCGCGATCACCTTCTCAAAAAGGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCTAATCGCTGTGGACTCGCGAGTCAAAGATCAACTGCGTCAAGCCGATGATCGACAACCACAAGCAGGATTACGAAGGGCTGGGGCAACCGAATCATGTGCATCTCATGCAGGAACATCt                                                                                                                                         <  2:214305/251‑1 (MQ=255)
                                                                                                    cTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGACTGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGCTGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTg                                                                                                                                        >  1:44677/1‑251 (MQ=255)
                                                                                                     tttCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCt                                                                                                                                         <  1:111044/249‑1 (MQ=255)
                                                                                                       tCGCACCAAATAATGCGCGTCCGCTGACGATGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACAGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGTAGATCATGCAGGAACATCTgcgc                                                                                                                                     >  2:147659/1‑251 (MQ=255)
                                                                                                            ccAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGCTGATCGAAAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCCGGAACATCTGCGCGACtt                                                                                                                                >  2:273023/1‑251 (MQ=255)
                                                                                                              aaaTAATGCGCGTCCGCGGACGCTGGAAGACGATCTTTTGAGCTGCCCCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTg                                                                                                                              <  2:285576/251‑1 (MQ=255)
                                                                                                                 taatGTGCGTCCGCTGACGCTGGAAGCCGATCGTTTGAGCTGAACAGTTAGCGGCTACACCTTCGCGATCCCCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGCATCGCTGTTGACTCGCGAGc                                                                                                                                                                                                                                          >  2:200371/1‑140 (MQ=255)
                                                                                                                 taatGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCGTCTCAAAAATGAGTGTCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGc                                                                                                                                                                                                                                          <  1:200371/140‑1 (MQ=255)
                                                                                                                        cgTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAAc                                                                                                                     >  2:416434/1‑250 (MQ=255)
                                                                                                                             gCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCg                                                                                                               >  2:376216/1‑251 (MQ=255)
                                                                                                                             gCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCg                                                                                                               <  1:112977/251‑1 (MQ=255)
                                                                                                                                         aGACGATCGTTTGAGCTGCACAGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATAAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGACGGGCTGTGGCAACCGACTCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGc                                                                                                   >  1:380191/1‑251 (MQ=255)
                                                                                                                                           aCGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTg                                                                                                 <  1:184611/251‑1 (MQ=255)
                                                                                                                                            cGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGa                                                                                                <  2:99168/251‑1 (MQ=255)
                                                                                                                                             gATCGTTTGAGCTTCACCGTTCGCGGCTACACCTTCTCGATCACCTTATACAACATGCGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCACGCCGCTGATCGCCAACCACCAGCAGGAGCACGACAGGCTGAGGCAACCGAATCCTTTTCAGATCATGCAGGCACATCTGCGCGCCTATACCGCAGACCAGAGCGATGGCGAAGTGCtcat                                                                                               >  2:25321/1‑248 (MQ=255)
                                                                                                                                                cgattGAGATGCACCGTTCGCGGCTACAACTTCGCGAGCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTTAATGGCGAATCGCTGTTGACTCCCGAGCCAAAGATCAACTTCTGCAAGCCGATGAGCGACAACCACAAGCAGGAGGACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTGTGCCGTAGAGCAGATCGAGTGTGAAGTGCGGatcat                                                                                            <  2:210227/248‑1 (MQ=255)
                                                                                                                                                     gAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATACTGGCAGGTGAATGGCGAACCGCTGTTGACGCGCGAGCCAAAGATCAACTTCTTCAAGCCGAGGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGAGCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGc                                                                                        <  2:32801/250‑1 (MQ=255)
                                                                                                                                                             ccGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGACCATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGtt                                                                               >  1:373883/1‑251 (MQ=255)
                                                                                                                                                                 tCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAACCCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTAtt                                                                            <  2:44677/250‑1 (MQ=255)
                                                                                                                                                                                  cgATCAACTTCTCAAAAATGAGTGGCAAACCGATATCCTGGCAGGTGAATGGCGAATCGATGTTGGCTCGCGAGCCAAAGATCAACTTCTTAAAGCCGATGATCGACAACCACAAGCAGGAGTACGCAGGGCTGAGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGCCTTTGCCTTAGACAAGAGCGATGGTGAAGTGCTGATCATCAGCCGAACAGTTATTGCCCCGCCGGTGTTTgcc                                                          >  2:75686/1‑249 (MQ=255)
                                                                                                                                                                                   gATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACt                                                         >  1:430042/1‑251 (MQ=255)
                                                                                                                                                                                       aCCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGCGTTTGACTTCgg                                                     >  2:344604/1‑251 (MQ=255)
                                                                                                                                                                                         cTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACTTCggg                                                    >  2:339398/1‑250 (MQ=255)
                                                                                                                                                                                           tctcAAAAATGAGTGGCAACCCGACATCCTGGCAGGGGAATGGCGAACCGCTGTTGACTCGCGACCCAAAGACCAACTCCTCCAAGCCGATGACCGACACCCACAAGCAGGAGTACGAAGGGCTGTGGCACCCGAATCATTTGCAGATCATGCAGGACCATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACTTCGGGAt                                                  <  2:373883/250‑1 (MQ=255)
                                                                                                                                                                                              cAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACTTCGGGATg                                                 <  1:1558/248‑1 (MQ=255)
                                                                                                                                                                                              aaaaaaTGGGTGGCAAACCTACAGCCTGGCGGGGGAATGGCGAACCGCTGTGGCCTCGGGAGCAAATGATCAACTTCTTCAAGCCGATGAGCGACACCCACAAGCAGGAGTACGAAGGGCTGTGGCATCCGAATCATTTGCAGGTCATGCAGGAACATCTGCGCGACTTGGCCGTAGAACAGAGGGATGGTTAATTGCTGAGCATCAGCCGCACATTTATTGCCCCGCCGGTGTTTGACTTCGGGATgcgc                                              <  2:22343/250‑1 (MQ=255)
                                                                                                                                                                                                    tGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACTTCGGGATGCGCTGCAcc                                        >  1:399781/1‑251 (MQ=255)
                                                                                                                                                                                                      aGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACTTCGGGATGCGCTGCACCt                                       >  1:315491/1‑250 (MQ=255)
                                                                                                                                                                                                      aGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACTTCGGGATGCGCTGCACCt                                       >  1:130591/1‑250 (MQ=255)
                                                                                                                                                                                                        tGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGAGTTTGACTTCGGGATGCGCTGCACCTAc                                     >  2:160502/1‑250 (MQ=255)
                                                                                                                                                                                                         ggCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACTTCGGGATGCGCTGCACCTACAt                                   >  2:152343/1‑251 (MQ=255)
                                                                                                                                                                                                            aacccGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGAATCGCGAGACAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCCCAGTTATTGCCCCGCCGGTGTTTGACTTCGGGATGCGCTGCACCTACATCTg                                >  2:228642/4‑251 (MQ=255)
                                                                                                                                                                                                             aaCCGCCATCCTGGCCAGTGAATGGCGAATCGCTGTTGAATCGCGAGCCAAAGATCAACTTATACAAGCCGATGATCGACAACCACAAGCAGGCGTACGCAGGGCTGTGGCAACCGAATCATTTGAAGATCCTGCAGGAACATCAGCGCGACTTTGCAGTGGACCAGAGCGATGGTGAAGAGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGAGTTTGACTTCGGGATGCGCTGCACCTCCATCTgg                               >  2:464633/1‑251 (MQ=255)
                                                                                                                                                                                                             aaCCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGCCTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTAATGCCCCGCCGGTGTTTGACTTCGGTATTCGCTGCACCGACATCTgg                               >  1:314221/1‑251 (MQ=255)
                                                                                                                                                                                                              aCCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACTTCGTGATGCGCTGCACCTACATCTGgc                              >  1:278958/1‑251 (MQ=255)
                                                                                                                                                                                                               ccGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCACCTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACCGTTATTGCCCCGCCGGTGTTTGACTTCGGGATGCGCTGCACCTACATCTGgcg                             >  2:106924/1‑251 (MQ=255)
                                                                                                                                                                                                                cGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGTGAGCCAAAGACCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGGGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGGGAAGTGCTGATCATCAGCCGCACAGTTATGGCCCCGCCGGTGTTTGACTTCGGGATGCGCTGCCCCTACATCTGgcgc                            <  1:222905/251‑1 (MQ=255)
                                                                                                                                                                                                                 gACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACTTCg                                                      >  1:374347/1‑224 (MQ=255)
                                                                                                                                                                                                                     aTCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACTTCGGGATGCGCTGCACCTACATCTGGCGCATCg                        >  1:114628/1‑251 (MQ=255)
                                                                                                                                                                                                                           gagggggAATGGCGATTCGCTGTTGACTCCGGAGCCAAAGACCAACTTCTCCAAGCCGAGTACCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCGGATCATGCAGGAACATCTGCGCGACTTTGCCGTAAAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACTTCGGGATGCGCTGCACCTACATCTGGCGCATCGCTGCCg                  <  1:88285/245‑1 (MQ=255)
                                                                                                                                                                                                                           gCAGGGGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTCCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACTTCGGGATGCGCTGCACCTACATCTGGCGCATCGCTGCCg                  <  1:438439/251‑1 (MQ=255)
                                                                                                                                                                                                                                    tGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACTTCGGGATGCGCTGCACCTACATCTGGCGCATCGCTGCCGATGGCCAgg         >  1:104169/1‑251 (MQ=255)
                                                                                                                                                                                                                                           acGCTGTTGACTCGCGAGCCAAAGATCACCTTCTTCAAGCCGATGACCAACACCCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGTAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACTTCGGGATGCGCTGCACCTACATCTGGCGCATCGCTGCCGATGGCCAGGTTAACGt  <  1:133770/250‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                    
ATTACGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACTTCGGGATGCGCTGCACCTACATCTGGCGCATCGCTGCCGATGGCCAGGTTAACGT  >  NC_000913/3224655‑3225139

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