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.6% 243.1 / NA 72L755L (CTG→TTG) ebgAevolved beta‑D‑galactosidase subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (47/24);  total (48/24)
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.98e-01

ACATTACGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACTTCGGGATGCGCTGCACCTACATCTGGCGCATCGCTGCCGATGGCCAGGTTAACGTG  >  NC_000913/3224653‑3225140
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     
aCATTACGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTcgc                                                                                                                                                                                                                                               >  1:330116/1‑251 (MQ=255)
   ttACGGTGACCAAAGATTCCCGACCCCGCTACAGAGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACCCCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGCAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACCTCCTGGCAGGTGAATGGCGAATCGCTGTTGACTcgcg                                                                                                                                                                                                                                              >  1:123293/1‑249 (MQ=255)
    tACGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTTCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAg                                                                                                                                                                                                                                            <  2:349220/250‑1 (MQ=255)
     aCGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGc                                                                                                                                                                                                                                           <  1:105816/250‑1 (MQ=255)
     aCGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGGCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGc                                                                                                                                                                                                                                           <  1:410576/250‑1 (MQ=255)
             cAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGACTGTCGCACCAAATAAGGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAAGCGCTGTTGACTCGCGAGCCAAAGATCa                                                                                                                                                                                                                                  <  1:78239/251‑1 (MQ=255)
                        gCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGACTGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCa                                                                                                                                                                                                                         >  2:140802/1‑249 (MQ=255)
                           cccGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGGCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAAGGAGTTGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTGCTTCAAGCCg                                                                                                                                                                                                                    <  2:287886/251‑1 (MQ=255)
                            ccGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAACCCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCg                                                                                                                                                                                                                    >  1:164513/1‑250 (MQ=255)
                             cGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTAGCGATCACCTTCTCAAAAATGAGTGGCAAACCGCCATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCCACTTCTTCAAGCCgat                                                                                                                                                                                                                  >  1:309065/1‑251 (MQ=255)
                              gCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCgatg                                                                                                                                                                                                                 >  1:401979/1‑251 (MQ=255)
                                 aCAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATc                                                                                                                                                                                                              <  2:221123/251‑1 (MQ=255)
                                 aCAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATc                                                                                                                                                                                                              >  2:398419/1‑251 (MQ=255)
                                      gAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACaa                                                                                                                                                                                                         <  2:330116/251‑1 (MQ=255)
                                       aaGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGCCAAc                                                                                                                                                                                                        >  2:14019/1‑251 (MQ=255)
                                       aaGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAAc                                                                                                                                                                                                        <  2:195282/251‑1 (MQ=255)
                                          ccGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACCCCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCGGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTGCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACcac                                                                                                                                                                                                     <  2:240139/251‑1 (MQ=255)
                                              acacCCTATCGCCCCTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTGTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGGCTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGc                                                                                                                                                                                                 <  1:45338/251‑1 (MQ=255)
                                              acacCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGc                                                                                                                                                                                                 >  1:360467/1‑251 (MQ=255)
                                                   cTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGAGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAAATTCTTCAAGCCGATGATCGTCAACCACAAGCAgg                                                                                                                                                                                              >  1:10896/1‑249 (MQ=255)
                                                           ccTTAGCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAgg                                                                                                                                                                                    <  2:342290/250‑1 (MQ=255)
                                                                          cTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGACAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGGGTACGAAGGGCTGTGGCAACCGa                                                                                                                                                                      >  1:9570/1‑250 (MQ=255)
                                                                                   aaccccGCGCGGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGGGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCa                                                                                                                                                            <  2:360467/247‑1 (MQ=255)
                                                                                   aaCACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAACATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGCCAACCACAAGCAGGGGTACGCAGGGCTGTGGCAACCGAATCATTTGc                                                                                                                                                             >  1:38472/1‑250 (MQ=255)
                                                                                   aaCACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCa                                                                                                                                                            <  1:43974/251‑1 (MQ=255)
                                                                                    acacCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGc                                                                                                                                                             >  1:6692/1‑249 (MQ=255)
                                                                                     cacCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGc                                                                                                                                                             >  2:131983/1‑248 (MQ=255)
                                                                                           gcAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGTAACCGAATCATTTGCAGATCATGc                                                                                                                                                    >  1:308586/1‑251 (MQ=255)
                                                                                                     ccTTGCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCTTTTGAGCTGCACCGTTCGCGTCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGTCAACCGAATCATTTGCAGATCAGGCAGGAACATCt                                                                                                                                          <  2:321614/251‑1 (MQ=255)
                                                                                                        ttCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTg                                                                                                                                         <  2:259650/249‑1 (MQ=255)
                                                                                                             aCCAAATAATGCGCGGCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGCACATCTGCGCGACt                                                                                                                                  >  2:95090/1‑251 (MQ=255)
                                                                                                               cAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTGCGCGGCTAAAACTTCGCGAGCACCGTCTCAAAAATGAGTGGCAAACCGACATCCAGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCATCTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCCGCGCGACttt                                                                                                                                <  2:391533/251‑1 (MQ=255)
                                                                                                                   taatGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCCCCTTCTCAAATATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACAACTGCGCGCCTTTGCCg                                                                                                                            >  2:143204/1‑251 (MQ=255)
                                                                                                                   taatGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGCATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCg                                                                                                                            >  2:149319/1‑251 (MQ=255)
                                                                                                                         gcgTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAAc                                                                                                                      >  2:304353/1‑251 (MQ=255)
                                                                                                                             ccGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACag                                                                                                                    >  2:123096/1‑249 (MQ=255)
                                                                                                                                         gAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGTGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGt                                                                                                      >  2:84740/1‑251 (MQ=255)
                                                                                                                                           aGACGATCGTTTGAGCTGCCCCGTTCGCGGCTACAACTTCGCGATCAACTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCCGGTGAATGGCGAATCGCTGTTGACTCACGAGCCAAAGATCAACTTCTTCAAGCCGCTGATCGACAACCACAAGACGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGc                                                                                                    >  1:284336/1‑251 (MQ=255)
                                                                                                                                             aCGATCTTTTGAAATGCACCGTAAGCGGCTACAACTTCGCGCTCCCCTTCCCAGAAATGAGCGGCACACAGAAATCCTGGAAGGTGAATGGCGCATCGCTGTTAACTCGCGAGCCAACGATCAACTTCTACAAGCCGATGATCGACAACCACAAGCAGGAGTACGACGGGCTGTGTCCACCGAAACGTTTGCAGATCATGCAGGCAAATCTGCGCGACTTTGCCGGAGAACAGAGCGATGGTGAAgaactg                                                                                                  >  2:312221/1‑246 (MQ=255)
                                                                                                                                             aCGATCGTGTGAGCTGCAACGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCACACCGACATCCTGGCAGGTGAATGGCGCACCGCTGTTGACTCGCGAGACAAAGATCAACTTCTTCAAGCCGATGATCCACAACCACAAGCTGGAGCACGATGGGCTGTGGCAACAGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAgggtgg                                                                                                  >  2:145926/1‑246 (MQ=255)
                                                                                                                                                      tGAGCTGCCCCGTTCGCGGCTACAACTTCGCGATCCCCTTCTCAAAAATGAGGGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAg                                                                                          <  1:91720/250‑1 (MQ=255)
                                                                                                                                                           tGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCTTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGcac                                                                                    >  2:216251/1‑251 (MQ=255)
                                                                                                                                                           tGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGcac                                                                                    >  1:120989/1‑251 (MQ=255)
                                                                                                                                                               cctttCGCGGCTACAACTTCGCGAGCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTGCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGtt                                                                                <  2:272420/248‑1 (MQ=255)
                                                                                                                                                                cGTTCGCGGCTACAACGTCGCGATCACCTTCTCAAAAATGAGTGGCAACCCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCTCAGTTa                                                                               <  2:301772/251‑1 (MQ=255)
                                                                                                                                                                 gTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTAt                                                                              >  1:406407/1‑251 (MQ=255)
                                                                                                                                                                      aggCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGGGAAGTGCTGATCATCAGCCGCACAGTTATTGccc                                                                         >  2:246429/2‑251 (MQ=255)
                                                                                                                                                                       ggCTACCACTGCGCGACCACATTCTCAAAAATGAGTGGCAAAACGACATCCTAGCAGGCGAACGGCGAATCGCTGCTGACCCGCGAGCCAAAAATCAACATCTTCAAGACGTTGCTTGACAACCACACGCAGGAGTACGAAGGGCTGTGGCCACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCACCCGCAACGTTATTGcccc                                                                        >  2:50801/1‑251 (MQ=255)
                                                                                                                                                                                ttCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCCGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGAGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCACACAGTTATTGCCCCGCCggattt                                                               >  1:408430/1‑247 (MQ=255)
                                                                                                                                                                                 tCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGttt                                                              >  1:244965/1‑251 (MQ=255)
                                                                                                                                                                                   gcgATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGACGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTg                                                             >  2:406806/1‑250 (MQ=255)
                                                                                                                                                                                         aCCTTCTCAAGAATGAGTAGCAAACCGACATCCTGGCCGGTGCACGGCTAATCGCTGTTGACTCGCGAGCCACACATCAACTTCTTCAAGCCGATGATAGACAACCACAAGCAGGAGTACGAAGGGATGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGAATTTGCCGTAGTCCAGAGCGATGGTGAAGTGCTGATCCTCAGCCGCACCGTTATTGCCCCGCCGGTGGTTGACTTCgg                                                      >  2:74196/1‑251 (MQ=255)
                                                                                                                                                                                                    aaggAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGCCTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACTTCGGGATGCGCTGCa                                           <  1:202752/248‑1 (MQ=255)
                                                                                                                                                                                                    aaTGAGCGGCAAACAGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGAGGCAACCGAATCATTTGCAGATCATGCAGGGACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGACGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACTTCGGGATGCGCTGc                                            >  2:57770/1‑250 (MQ=255)
                                                                                                                                                                                                           ggCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCCCCAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGATATTGCACCGCCGGTGTTTGACTTCGGGATGCGCTGCCCCTACAt                                    >  2:269841/1‑251 (MQ=255)
                                                                                                                                                                                                              aaaCCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCACCTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACTTCGGGATGCGCGGCACCTACATCTg                                 <  1:131983/251‑1 (MQ=255)
                                                                                                                                                                                                              aaaCCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCACGCCGATGATCGACAAACACACGCAGGAGAACGCAGGGCTGTGGCACCCGAATCATTTACAGATCATGCAGGACCATCTGCGCGACTTTGCCGTAGCACAGAGCGATGGGGAAGTGCTGATCATCAGCCGCACAGTTATTCCCCCGCCGGTGTTTGACTTCGGGATGCGCTGCACCTACATCTg                                 >  1:429848/1‑251 (MQ=255)
                                                                                                                                                                                                              aaaCCGAAATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGACGAGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACTTCGGGATGCGCTGCACCTACATCTg                                 >  2:189434/1‑251 (MQ=255)
                                                                                                                                                                                                                  cGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCTATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACTTCGGGATGCGCTGCACCTACATCTGgcgc                             >  1:216006/1‑251 (MQ=255)
                                                                                                                                                                                                                    aCATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGGGTTTGACTCCGGGTTGCGCTCCCCCTACATCTGGCGCac                           >  1:48291/1‑250 (MQ=255)
                                                                                                                                                                                                                    aCATCCGGGCAGGTGAAGGGCGAAGCGCTGTGGCCGCGCGAGCCAAAGATCAACTCCTTCAAGCCGATGATCGACACCCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATGGGCGGGACTTGGCCGTAGAGCAGAGCGATGGTGAAGTGCTGAGCATCAGCCGCACATTTATTGCCCCGCCGGTGTTTGACTTCGGGATGCGGTGCACCTACATCGGGCGCAt                           <  2:433078/251‑1 (MQ=255)
                                                                                                                                                                                                                     cATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGTTGGTGAAGTGATGATCATCAGCCGCACAGTTTTTGCCCCGCCGGTGTTTGACTTCGGGATACGCTGCACCTACATCTGGCGCATc                          >  2:19007/1‑251 (MQ=255)
                                                                                                                                                                                                                        ccTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGcc                                                                          >  1:81753/1‑200 (MQ=255)
                                                                                                                                                                                                                        ccTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGcc                                                                          <  2:81753/200‑1 (MQ=255)
                                                                                                                                                                                                                         cTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGTTCGACAACCACCAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACTTCGGGATGCGCTGCACCTACATCTGGCGCATCGCTg                      >  1:33326/1‑251 (MQ=255)
                                                                                                                                                                                                                            gCAGGGGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACTTCGGGATGCGCTGCCCCTACATCTGGCGCCTCGCTGCCg                   >  2:335862/1‑251 (MQ=255)
                                                                                                                                                                                                                                  gAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACTTCGGGATGCGCTGCACCTACATCTGGCGCATCGCTGCCGATGGcc             >  1:11605/1‑251 (MQ=255)
                                                                                                                                                                                                                                   aaTGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGCCTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACTTCGGGATGCGCTGCACCTACATCTGGCGCATCGCTGCCGATGGcc             >  1:239009/1‑250 (MQ=255)
                                                                                                                                                                                                                                    aTGGCGAATCGCTGTTGACTCGCGCGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGttt                                                              >  2:365037/1‑200 (MQ=255)
                                                                                                                                                                                                                                    aTGGCGAATCGCTGTTGACTCGCGAGCAAAAGATCAACTTCTCCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCGGCGCGACTTTGCCGTAGAACAGAGGGATGGTGAAGTGCTGATCATCAGCCGCACATTTATTGCCCCGCCGGTGttt                                                              <  1:365037/200‑1 (MQ=255)
                                                                                                                                                                                                                                         gAATCGCTGTTGACTCGCGAGCCAAAGCTCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACTTCGGGATGCGCTGCACCTACCTCTGGCGCATCGCTGCCGATGGCCAGGTTa       >  2:428366/1‑250 (MQ=255)
                                                                                                                                                                                                                                             cGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACTGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACTTCGGGATGCGCTGCACCTACATCTGGCGCATCGCTGCCGATGGCCAGGTTAACGTg  >  1:254267/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     
ACATTACGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACTTCGGGATGCGCTGCACCTACATCTGGCGCATCGCTGCCGATGGCCAGGTTAACGTG  >  NC_000913/3224653‑3225140

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