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

ATTACGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCG‑ACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACTTCGGGATGCGCTGCACCTACATCTGGCGCATCGCTGCCGATGGCCAGGTTAACGT  >  NC_000913/3224655‑3225139
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                    
aTTACGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCG‑ACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGa                                                                                                                                                                                                                                            <  1:73026/251‑1 (MQ=255)
       tGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCG‑ACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCaa                                                                                                                                                                                                                                       >  1:191687/1‑249 (MQ=255)
       tGACCAAAGATTACCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCG‑ACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCac                                                                                                                                                                                                                                       >  1:332631/1‑248 (MQ=255)
        gACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCG‑ACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGa                                                                                                                                                                                                                                    <  2:151365/251‑1 (MQ=255)
           cAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCG‑ACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATc                                                                                                                                                                                                                                  >  1:160244/1‑250 (MQ=255)
            aaaGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCG‑ACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCa                                                                                                                                                                                                                                 >  2:1558/1‑250 (MQ=255)
              aGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTGCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTGCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCG‑ACATCCTGGCAGGTGAATGGCGAATCGCTGTTGTCTCGCGAGCCAAAGATCAAct                                                                                                                                                                                                                              <  1:377407/251‑1 (MQ=255)
                        aCCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCG‑ACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGcc                                                                                                                                                                                                                    >  1:334942/1‑251 (MQ=255)
                                           gACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCG‑ACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACGCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAg                                                                                                                                                                                                 <  1:293203/251‑1 (MQ=255)
                                              ccccTATCGCCACTTATCAGGTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAAGGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTGCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCG‑ACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTGCTTCAAGCCGATGATCGACAACCACAAGCAg                                                                                                                                                                                              <  2:168394/250‑1 (MQ=255)
                                                  tATCGCCACTTATCAGTTACCGCTGAAGGAAAACACCGAGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCG‑ACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGt                                                                                                                                                                                          >  2:307167/1‑251 (MQ=255)
                                                   aTCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCG‑ACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGt                                                                                                                                                                                          >  2:112977/1‑250 (MQ=255)
                                                         aCTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCG‑ACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAgg                                                                                                                                                                                   >  1:126965/1‑251 (MQ=255)
                                                         aCTTACCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCTTTTCGCACCAAATAATGCTCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCG‑ACATCCTGGCAGGTGAATGGCGAATCGCGGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAgg                                                                                                                                                                                   <  2:253364/251‑1 (MQ=255)
                                                             aTCAGTTCCCGCTGAAGGAAAACCCCGCGCAGCCAGTGCCTGTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCG‑ACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCt                                                                                                                                                                                <  2:124046/250‑1 (MQ=255)
                                                             aTCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCG‑ACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCtg                                                                                                                                                                               >  1:65679/1‑251 (MQ=255)
                                                                 gTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCG‑ACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGc                                                                                                                                                                           <  1:130810/251‑1 (MQ=255)
                                                                  ttCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCG‑ACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAACGATCAACTTCTTCAAGCCGATGATCGACACCCACAAGCAGGAGTACGAAGGGCTGTGGc                                                                                                                                                                           >  2:184611/1‑250 (MQ=255)
                                                                  ttCCCGCTGAAGGAAAACACCGCGCAACCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCG‑ACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGAACGAAGGGCtgtagc                                                                                                                                                                           >  2:323027/1‑247 (MQ=255)
                                                                    cccGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCG‑ACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAAc                                                                                                                                                                        >  1:285576/1‑251 (MQ=255)
                                                                     ccGCTGAAGGAAAACCCCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCG‑ACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGTCAAcc                                                                                                                                                                       <  1:126724/251‑1 (MQ=255)
                                                                     ccGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCG‑ACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAAcc                                                                                                                                                                       <  2:65679/251‑1 (MQ=255)
                                                                       gCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCG‑ACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCg                                                                                                                                                                      >  2:333440/1‑250 (MQ=255)
                                                                          gAAGGAAAGCACCGAGCAGCAAGAGCCTTTCGAACAGCATAATGCGCGTCCGATGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTAGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGAGG‑AAACCGGACATCCTGGCAGGTGAAAGGCGAATCGCTGTTGACACGCGAGCCAAAGATCAACTTCTTAAAGACGATGATCGACAACACCAAGCCGGAATACGCAGGGCTGTGGCACCCGAATc                                                                                                                                                                  >  2:277749/1‑251 (MQ=255)
                                                                          gAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCCCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCG‑ACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAAt                                                                                                                                                                   <  1:137483/250‑1 (MQ=255)
                                                                            aGGAAAACACCGCGCAGCCAGTGCCGTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGGGGCAAACCG‑ACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTGCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCAt                                                                                                                                                                <  2:369125/251‑1 (MQ=255)
                                                                                aaaCACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCG‑ACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGc                                                                                                                                                            <  1:332066/251‑1 (MQ=255)
                                                                                aaaCACCGCGCAGCCAGTGCCTTTCGCACCAAAGAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCACGATCACCTTCTCAAAAATGAGTGGCAAACCG‑ACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGc                                                                                                                                                            <  1:74593/251‑1 (MQ=255)
                                                                                  acacCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGAGCGTTTGAGCTGCACCGTGCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCG‑ACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAg                                                                                                                                                          <  2:383893/251‑1 (MQ=255)
                                                                                                gTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGCCTACAACTTCGCGATCACCTTCTCAAACATGAGTGGCAAACCG‑ACATACTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGAACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAAc                                                                                                                                             >  1:214305/1‑250 (MQ=255)
                                                                                                  gCCTTTCGCACCAAATAATGCGCGTCCGCTGAAGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCG‑ACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATc                                                                                                                                          >  1:166418/1‑251 (MQ=255)
                                                                                                   ccTTTCGCCGCAAAGAAGGCGCGCCCGCTGACGCTGGAAGGCGAGCGTTTGGGCTGCACCGTGCGCGGCTACAACTTCGCGATCACCTTCTCAAAAAGGAGTGGCAAACCG‑ACATCCTGGCAGGTGAATGGCTAATCGCTGTGGACTCGCGAGTCAAAGATCAACTGCGTCAAGCCGATGATCGACAACCACAAGCAGGATTACGAAGGGCTGGGGCAACCGAATCATGTGCATCTCATGCAGGAACATCt                                                                                                                                         <  2:214305/251‑1 (MQ=255)
                                                                                                    cTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCG‑ACATCCTGGCAGGTGACTGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGCTGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTg                                                                                                                                        >  1:44677/1‑251 (MQ=255)
                                                                                                     tttCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCG‑ACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCt                                                                                                                                         <  1:111044/249‑1 (MQ=255)
                                                                                                       tCGCACCAAATAATGCGCGTCCGCTGACGATGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACAG‑ACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGTAGATCATGCAGGAACATCTgcgc                                                                                                                                     >  2:147659/1‑251 (MQ=255)
                                                                                                            ccAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCG‑ACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGCTGATCGAAAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCCGGAACATCTGCGCGACtt                                                                                                                                >  2:273023/1‑251 (MQ=255)
                                                                                                              aaaTAATGCGCGTCCGCGGACGCTGGAAGACGATCTTTTGAGCTGCCCCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCG‑ACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTg                                                                                                                              <  2:285576/251‑1 (MQ=255)
                                                                                                                 taatGTGCGTCCGCTGACGCTGGAAGCCGATCGTTTGAGCTGAACAGTTAGCGGCTACACCTTCGCGATCCCCTTCTCAAAAATGAGTGGCAAACCG‑ACATCCTGGCAGGTGAATGGCGCATCGCTGTTGACTCGCGAGc                                                                                                                                                                                                                                          >  2:200371/1‑140 (MQ=255)
                                                                                                                 taatGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCGTCTCAAAAATGAGTGTCAAACCG‑ACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGc                                                                                                                                                                                                                                          <  1:200371/140‑1 (MQ=255)
                                                                                                                        cgTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCG‑ACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAAc                                                                                                                     >  2:416434/1‑250 (MQ=255)
                                                                                                                             gCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCG‑ACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCg                                                                                                               >  2:376216/1‑251 (MQ=255)
                                                                                                                             gCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCG‑ACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCg                                                                                                               <  1:112977/251‑1 (MQ=255)
                                                                                                                                         aGACGATCGTTTGAGCTGCACAGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCG‑ACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATAAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGACGGGCTGTGGCAACCGACTCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGc                                                                                                   >  1:380191/1‑251 (MQ=255)
                                                                                                                                           aCGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCG‑ACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTg                                                                                                 <  1:184611/251‑1 (MQ=255)
                                                                                                                                            cGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCG‑ACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGa                                                                                                <  2:99168/251‑1 (MQ=255)
                                                                                                                                             gATCGTTTGAGCTTCACCGTTCGCGGCTACACCTTCTCGATCACCTTATACAACATGCGTGGCAAACCG‑ACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCACGCCGCTGATCGCCAACCACCAGCAGGAGCACGACAGGCTGAGGCAACCGAATCCTTTTCAGATCATGCAGGCACATCTGCGCGCCTATACCGCAGACCAGAGCGATGGCGAAGTGCtcat                                                                                               >  2:25321/1‑248 (MQ=255)
                                                                                                                                                cgattGAGATGCACCGTTCGCGGCTACAACTTCGCGAGCACCTTCTCAAAAATGAGTGGCAAACCG‑ACATCCTGGCAGGTTAATGGCGAATCGCTGTTGACTCCCGAGCCAAAGATCAACTTCTGCAAGCCGATGAGCGACAACCACAAGCAGGAGGACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTGTGCCGTAGAGCAGATCGAGTGTGAAGTGCGGatcat                                                                                            <  2:210227/248‑1 (MQ=255)
                                                                                                                                                     gAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCG‑ACATACTGGCAGGTGAATGGCGAACCGCTGTTGACGCGCGAGCCAAAGATCAACTTCTTCAAGCCGAGGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGAGCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGc                                                                                        <  2:32801/250‑1 (MQ=255)
                                                                                                                                                             ccGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCG‑ACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGACCATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGtt                                                                               >  1:373883/1‑251 (MQ=255)
                                                                                                                                                                 tCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAACCCG‑ACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTAtt                                                                            <  2:44677/250‑1 (MQ=255)
                                                                                                                                                                                  cgATCAACTTCTCAAAAATGAGTGGCAAACCG‑ATATCCTGGCAGGTGAATGGCGAATCGATGTTGGCTCGCGAGCCAAAGATCAACTTCTTAAAGCCGATGATCGACAACCACAAGCAGGAGTACGCAGGGCTGAGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGCCTTTGCCTTAGACAAGAGCGATGGTGAAGTGCTGATCATCAGCCGAACAGTTATTGCCCCGCCGGTGTTTgcc                                                          >  2:75686/1‑249 (MQ=255)
                                                                                                                                                                                   gATCACCTTCTCAAAAATGAGTGGCAAACCG‑ACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACt                                                         >  1:430042/1‑251 (MQ=255)
                                                                                                                                                                                       aCCTTCTCAAAAATGAGTGGCAAACCG‑ACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGCGTTTGACTTCgg                                                     >  2:344604/1‑251 (MQ=255)
                                                                                                                                                                                         cTTCTCAAAAATGAGTGGCAAACCG‑ACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACTTCggg                                                    >  2:339398/1‑250 (MQ=255)
                                                                                                                                                                                           tctcAAAAATGAGTGGCAACCCG‑ACATCCTGGCAGGGGAATGGCGAACCGCTGTTGACTCGCGACCCAAAGACCAACTCCTCCAAGCCGATGACCGACACCCACAAGCAGGAGTACGAAGGGCTGTGGCACCCGAATCATTTGCAGATCATGCAGGACCATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACTTCGGGAt                                                  <  2:373883/250‑1 (MQ=255)
                                                                                                                                                                                              cAAAAATGAGTGGCAAACCG‑ACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACTTCGGGATg                                                 <  1:1558/248‑1 (MQ=255)
                                                                                                                                                                                              aaaaaaTGGGTGGCAAACCT‑ACAGCCTGGCGGGGGAATGGCGAACCGCTGTGGCCTCGGGAGCAAATGATCAACTTCTTCAAGCCGATGAGCGACACCCACAAGCAGGAGTACGAAGGGCTGTGGCATCCGAATCATTTGCAGGTCATGCAGGAACATCTGCGCGACTTGGCCGTAGAACAGAGGGATGGTTAATTGCTGAGCATCAGCCGCACATTTATTGCCCCGCCGGTGTTTGACTTCGGGATgcgc                                              <  2:22343/250‑1 (MQ=255)
                                                                                                                                                                                                    tGAGTGGCAAACCG‑ACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACTTCGGGATGCGCTGCAcc                                        >  1:399781/1‑251 (MQ=255)
                                                                                                                                                                                                      aGTGGCAAACCG‑ACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACTTCGGGATGCGCTGCACCt                                       >  1:315491/1‑250 (MQ=255)
                                                                                                                                                                                                      aGTGGCAAACCG‑ACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACTTCGGGATGCGCTGCACCt                                       >  1:130591/1‑250 (MQ=255)
                                                                                                                                                                                                        tGGCAAACCG‑ACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGAGTTTGACTTCGGGATGCGCTGCACCTAc                                     >  2:160502/1‑250 (MQ=255)
                                                                                                                                                                                                         ggCAAACCG‑ACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACTTCGGGATGCGCTGCACCTACAt                                   >  2:152343/1‑251 (MQ=255)
                                                                                                                                                                                                            aacccG‑ACATCCTGGCAGGTGAATGGCGAATCGCTGTTGAATCGCGAGACAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCCCAGTTATTGCCCCGCCGGTGTTTGACTTCGGGATGCGCTGCACCTACATCTg                                >  2:228642/4‑251 (MQ=255)
                                                                                                                                                                                                             aaCCG‑CCATCCTGGCCAGTGAATGGCGAATCGCTGTTGAATCGCGAGCCAAAGATCAACTTATACAAGCCGATGATCGACAACCACAAGCAGGCGTACGCAGGGCTGTGGCAACCGAATCATTTGAAGATCCTGCAGGAACATCAGCGCGACTTTGCAGTGGACCAGAGCGATGGTGAAGAGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGAGTTTGACTTCGGGATGCGCTGCACCTCCATCTgg                               >  2:464633/1‑251 (MQ=255)
                                                                                                                                                                                                             aaCCG‑ACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGCCTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTAATGCCCCGCCGGTGTTTGACTTCGGTATTCGCTGCACCGACATCTgg                               >  1:314221/1‑251 (MQ=255)
                                                                                                                                                                                                              aCCG‑ACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACTTCGTGATGCGCTGCACCTACATCTGgc                              >  1:278958/1‑251 (MQ=255)
                                                                                                                                                                                                               ccG‑ACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCACCTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACCGTTATTGCCCCGCCGGTGTTTGACTTCGGGATGCGCTGCACCTACATCTGgcg                             >  2:106924/1‑251 (MQ=255)
                                                                                                                                                                                                                cG‑ACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGTGAGCCAAAGACCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGGGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGGGAAGTGCTGATCATCAGCCGCACAGTTATGGCCCCGCCGGTGTTTGACTTCGGGATGCGCTGCCCCTACATCTGgcgc                            <  1:222905/251‑1 (MQ=255)
                                                                                                                                                                                                                 g‑ACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACTTCg                                                      >  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)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                    
ATTACGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCG‑ACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACTTCGGGATGCGCTGCACCTACATCTGGCGCATCGCTGCCGATGGCCAGGTTAACGT  >  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: ‑