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,8950CT97.5% 248.7 / NA 83L755L (CTG→TTG) ebgAevolved beta‑D‑galactosidase subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (42/37);  minor base G (2/0);  total (46/37)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAACATTACGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACTTCGGGATGCGCTGCACCTACATCTGGCGCATCGCTGCCGATGGCCAGGTTA  >  NC_000913/3224651‑3225135
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                
caacacttcGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTAACAGTTCCCGCGGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCTCTGTTGAc                                                                                                                                                                                                                                              <  2:498040/241‑1 (MQ=255)
cAACATTACGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTc                                                                                                                                                                                                                                            <  1:475578/251‑1 (MQ=255)
    aTTACGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACTCCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGa                                                                                                                                                                                                                                        <  1:141572/251‑1 (MQ=255)
         ggTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCa                                                                                                                                                                                                                                    >  1:209261/1‑250 (MQ=255)
            gACCAAAGAGTCCCGCACCCGCTACAGCGAATCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGa                                                                                                                                                                                                                                <  1:517874/251‑1 (MQ=255)
                aaaGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCa                                                                                                                                                                                                                             >  1:340359/1‑250 (MQ=255)
                  aGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAAct                                                                                                                                                                                                                          <  1:265022/251‑1 (MQ=255)
                       cccGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACttcttc                                                                                                                                                                                                                     >  1:254754/1‑251 (MQ=255)
                          ccACCCGCTACAGCGAATCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGGTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGTTGAATGGCGAAGCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAg                                                                                                                                                                                                                  <  2:44056/250‑1 (MQ=255)
                            aCCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGcc                                                                                                                                                                                                                >  1:235196/1‑251 (MQ=255)
                              ccGCTACAGCGAAGCCGGACCCCCTATCGCCCCTTATCAGTTCCCGCTGAAGGAAAACCCCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCga                                                                                                                                                                                                              <  2:544897/251‑1 (MQ=255)
                               cGCTACAGCGACGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGAGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGCATCGCTGGTGACTCGCGAGCCAAAGATAAAATTCTTACAGCCgat                                                                                                                                                                                                             >  2:410899/1‑251 (MQ=255)
                               cGCTACAGCGAAGCCGGCCCCCCTCCCGCCACTTATCAGTCCCCGCTAAAGGAAACCCCCGCGCAGCCAGTGCCTTTCGCACAAAATAAGGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGGGCTGCCCCGTTCGCGGCTACACCTTCGCGATCCCCTTCTCAAAAATGAGTGGCAAACCGACATCCGGGCAGGGGAATGGCGAACCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCgat                                                                                                                                                                                                             <  1:88838/251‑1 (MQ=255)
                                 cTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCgatg                                                                                                                                                                                                            >  1:407810/1‑250 (MQ=255)
                                   aCAGCGAAGCCGGACCCCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCgatgat                                                                                                                                                                                                          <  1:173623/250‑1 (MQ=255)
                                    cAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCg                                                                                                                                                                                                        <  1:216046/251‑1 (MQ=255)
                                          aGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGAACCAAATCATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCAgatgagctaacaac                                                                                                                                                                                                  >  2:258169/1‑242 (MQ=255)
                                               gACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAg                                                                                                                                                                                             <  1:110843/251‑1 (MQ=255)
                                                 cacCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCCCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGc                                                                                                                                                                                            >  1:255458/1‑250 (MQ=255)
                                                  acCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTAGCGGCTACAACTTAGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTATTCAAGCCGATGATCGACAACCACCAGc                                                                                                                                                                                            >  2:462253/1‑249 (MQ=255)
                                                       aTCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGAGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTAGCAAACCGACATCCTGGCAGGTGAATGGCGAGTCGCTGTTGACTCGCGAGCCAAAGATAAACTTCTTCAAGCCGATGATCGTCAAACACAAGCAGGAGt                                                                                                                                                                                      >  1:453399/1‑250 (MQ=255)
                                                         cGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCCCAAGCAGGAGTACg                                                                                                                                                                                   <  1:229810/251‑1 (MQ=255)
                                                          gCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCCAAGATCAACTTCTTCAGGCCGATGATCGAAAACCCCAAGCAGGGGTACg                                                                                                                                                                                   >  2:549782/1‑250 (MQ=255)
                                                            cACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCCCCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCCAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGa                                                                                                                                                                                  >  1:350312/1‑249 (MQ=255)
                                                                 aTCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCtg                                                                                                                                                                           >  1:269899/1‑251 (MQ=255)
                                                                    aGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCCCCAAATAATGCGCGTACGCTGACGCTGGAAGACGTTCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGATTGGCGAATCGCTGTTGACTCGCGATCCAAAGATCCACTTCTTAAAGCCGATGATCGACAACCAAAAGCAGGAGTAAGAAGGGCTGTgg                                                                                                                                                                        >  2:521202/1‑251 (MQ=255)
                                                                      tccccGCGGAAGGAAAACACCGGGCAGCAAGTTCCTTTCGCACCAATGAATGCGCGTCCGCTGACGCTGGAAGCGTATCGTTTGAGCTGCACCGTTCGCGGTTACAACGTCGCGATAACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCAAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGc                                                                                                                                                                       <  2:137458/248‑1 (MQ=255)
                                                                        cccGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCCCAAAATAATGCGCGTCCGCCGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGAACAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGTAGTACGAAGGGCTGTGGCAAc                                                                                                                                                                    <  1:237428/251‑1 (MQ=255)
                                                                         ccgcccttccgatcccccGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCCGACGCTGGAAGACGATCGGTTGACCTGCACCGTTCGCGGCTCCAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGAGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTGCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAAcc                                                                                                                                                                   <  2:435775/235‑1 (MQ=255)
                                                                            cTGACGGAAAACACCGCGCAGACAGTGCCATTCGCACCAAATAATCCGAGACCGCTGACCCGGGCAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTATCAAAAATGCGTGGCAAACCGACATCCTGGCAGGTGACGGGCGAATCGCTGTTGACTCGCGCACCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACTAGCAGGAGTACGAAGGGCtgtaggcaccga                                                                                                                                                                 >  2:320230/1‑241 (MQ=255)
                                                                            cTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGa                                                                                                                                                                 >  1:111960/1‑250 (MQ=255)
                                                                             ggAAGGAAATCCCCGCGTAGCACGTGCCTTTCGCACCAAATATGTCGCGCCCGCTTACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAAGGAGTGGCAAACCGACATCCTGGCAGGTGAAGGGCGAATCGCTGTTGACCCGCGAGCAAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAATCAGGAGTACGAAGGGCTGTGTCATCCGAAt                                                                                                                                                               <  2:368770/250‑1 (MQ=255)
                                                                              gAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATc                                                                                                                                                              >  1:404496/1‑251 (MQ=255)
                                                                               aaGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATc                                                                                                                                                              >  1:117740/1‑250 (MQ=255)
                                                                                agaaaaaCACCGCGCAGCCAGAGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTATCAGCTGCACCGTTCGCGGCTACAACCTCGCGATCACCTTCTCAAAACTGAGTGGCAAACCGACATCCTGGCAGGTGAATTGCGAAGCGCTGTTTACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGCGGCTCGACAACCACAAGAAGGGGTACGAAGGGCcgtgcgaaccggatcct                                                                                                                                                            >  2:16749/4‑238 (MQ=255)
                                                                                        acCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGAGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAAcc                                                                                                                                                                   >  1:435775/1‑236 (MQ=255)
                                                                                           gcccAGCCAGTGTCTTTCGAACCAAACTAGGCGCGTCCGCGTCCGCTGGATGCCAGTCGTCTGACCTGCCCCTTTCGCGGCTACAACTCCGCGATCGCCTTCGCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTAGACCCGCTAGCCAAAGACCAACTTCTTCAAGCCGATGATCGACAACCTCAAGCAGGAGTACTAAGGGCTGTGGCAACCGACTCATTTGCAGATCAt                                                                                                                                                 <  2:396603/248‑1 (MQ=255)
                                                                                            cgccgCCAGGGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACTGACATCCGGGAAGGTGAATGGCGAAACGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATg                                                                                                                                                <  2:497380/247‑1 (MQ=255)
                                                                                            cgcAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATg                                                                                                                                                >  1:422363/1‑251 (MQ=255)
                                                                                            cgcAGCCAGGCCCTTTCGCACCAAATAAGGCGCGCCCGCTGACGCTGGAAGACGATCTTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATg                                                                                                                                                <  2:117740/251‑1 (MQ=255)
                                                                                                gCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAgg                                                                                                                                            >  1:47105/1‑251 (MQ=255)
                                                                                                gCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCAACGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGCGAATGGCGAATCGCTGGTGACTCGCGAGCCAAAGATCAACTTCTTCAATCCGTTGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGCATCATTTGCAGATCAAGCAgg                                                                                                                                            >  2:53363/1‑251 (MQ=255)
                                                                                                 ccggggCCTTTCCCACCAAATAATGCGCGGCCGCTGACGTCTGATGACGATTGGTTGAGTGTCACCGTTCGCGGCTAAAACTTCTCGATCACCTTCTCAAAAATGATGGGCAAACCGACAGCCGTGCAGGGGAAGGGCGAATCGCTTTTGACTCGCGAGCCAAAGACCAACTTCGTCAATCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAAGCATTTGCAGATCATGCAGGa                                                                                                                                           <  2:340359/246‑1 (MQ=255)
                                                                                                   aGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAAc                                                                                                                                         <  1:126924/251‑1 (MQ=255)
                                                                                                       ccTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGACGACGATCGTTTGAGCTGCACCGTTAGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATACTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAACTCAACTTCTTCAAGCCGATGAACGACAACCACAAGCCGGAGAACGAAGGGCTGTGGCAACAGAATCCTTTGCAGATCATGCAGGAACATCt                                                                                                                                     >  2:227973/1‑251 (MQ=255)
                                                                                                                  aaaTAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTg                                                                                                                          <  1:227973/251‑1 (MQ=255)
                                                                                                                            ccggccGCTGCCGCGGGAAGACGCTTGTTATAGCGGCTCCTTTCGCGGCGACATCTTCGCGATCACCGTCCCAAAAATGAGTGGCAAACCCACAGTCTGGAAGGTAATTGGCGAATCGCTGTTGCCCCGCGAGCCAAAGACCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTTCCGTAGAAc                                                                                                                 <  2:533275/247‑1 (MQ=255)
                                                                                                                              gTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGGGCATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCACCTTCTTCAAGCCGATGTTCGCCAACCACAAGCAGGAGTACGAAAGGCTGTGGCAACAGCATACTTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAAc                                                                                                                 >  2:465573/1‑249 (MQ=255)
                                                                                                                               cccGCTGACGCTGGAAGACGACCGTTGGGGCTGCACCGTTCGCGGCTACAACTGCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTAAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACaga                                                                                                              <  2:196941/250‑1 (MQ=255)
                                                                                                                                ccGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCGCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACaga                                                                                                              <  1:522179/250‑1 (MQ=255)
                                                                                                                                 cgctggcgctggAAGACGATCGGTGGAGCGGCACCGTCCGCGGCAACAACGACGCGATCCCCCTCTCAAAAATGGTTGTCAAACCGACATATGGGCAGGTGAATGGCGAATCGCATTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCTCGATTTTGCCTTAGAACagaga                                                                                                            <  2:350312/251‑2 (MQ=255)
                                                                                                                                     gACGCTGGAAGACTAGCGTTTGAGCTGCCCCGTGCGCGGCTACAACTGCGCGATCCCCTTCTCAAAAAGGAGTGGCAAACCGACATCCTGGCAGGGGAATTGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATg                                                                                                        <  2:399319/251‑1 (MQ=255)
                                                                                                                                           tGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGa                                                                                                    >  2:154120/1‑250 (MQ=255)
                                                                                                                                           tGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGa                                                                                                    >  1:402512/1‑250 (MQ=255)
                                                                                                                                                gccGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCGCAAAAATGAGTGGCAAACCGACACCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCt                                                                                              <  1:449072/249‑1 (MQ=255)
                                                                                                                                                 aCGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACCGAGCGATGGGGAAGTGCTg                                                                                             >  2:197661/1‑251 (MQ=255)
                                                                                                                                                  acgacctttgtagcggCACCGTTCGCGGCTACAACTGCGCGATACCTGTCTCAAAAATGAGGGGCAAACCGCAATCCTCGCGGGTGAACGGCTATTCGCTGTTGACTCGCGAGCCAAAGAGCAACTGCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCGGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACATAGCGATGGTGAAGTGCTg                                                                                             <  2:541524/236‑1 (MQ=255)
                                                                                                                                                          ttGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGAtcatca                                                                                      <  1:405552/250‑1 (MQ=255)
                                                                                                                                                               cTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGc                                                                                 >  1:411427/1‑250 (MQ=255)
                                                                                                                                                                 gCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGcac                                                                               >  1:283504/1‑250 (MQ=255)
                                                                                                                                                                      cGTTCGCTGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCTACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTa                                                                          <  1:519496/251‑1 (MQ=255)
                                                                                                                                                                       gttctgggcaccaacttatcaatgaCCTTCACAAAAATGAGTGGCAAACCGACAACCTGGCATGTGAATGGCGAATCGCGGTGGACTCGCGAGCCCAAGATCAACTTCTTCAAGCCGAAGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCAGATCATCAGCCGCACATTTAt                                                                         <  2:471376/227‑1 (MQ=255)
                                                                                                                                                                         tCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGAAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTg                                                                       >  1:383764/1‑251 (MQ=255)
                                                                                                                                                                            cgGCTACAACTTCGCGAAAACATTCTCAAAAATGAGTGGCAAACCGCCATCCTGGCAGGTGAATGGCGAAACGCTGTTGACCCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGCCAACTACAAGCAGGGGTAAGAAGGGGTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTACCGTAGAAAAGAGCGATGGTGAAGTGCTGATCATCACCCCCACAGTTATTGctc                                                                    >  2:261355/1‑249 (MQ=255)
                                                                                                                                                                              gctaaatcttctctgtCACCTTCTAAAAAATGAGTGGCAAACCGACACCCTGGCAGGTGAATGGCGAATCGCTGTGGACTCGCGAGCAAAAGAACAACTCCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACTAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCTTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGccgcg                                                                  <  2:408124/236‑4 (MQ=255)
                                                                                                                                                                                    ccTTCGCGATCACCTACTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAACGGCGCATCGCTGTTGACCCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGAACGAAGGGCTGTGGCAACCGGATCATCTGCAGATCCTGCAGGAACATCTGCGCGACTTTGCCGAAGAAAAGAGCGATGGTGAAGAGCTGCTCATCAACCGCACAGTTATCTCCCCGCCggcg                                                            >  2:259261/2‑249 (MQ=255)
                                                                                                                                                                                              cACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGACCAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACTTCg                                                  >  1:418509/1‑251 (MQ=255)
                                                                                                                                                                                               aCCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGACGAATCGCTGTTGACTCGCGAGCCCACGATCAACTTCTTCAAACCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGACCCGCCGGTGTTTGACTTCgg                                                 >  1:529379/1‑251 (MQ=255)
                                                                                                                                                                                                          aaaTGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCTGGAGTACGAAGGGCTGAGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACa                                                                                                                >  2:18169/1‑178 (MQ=255)
                                                                                                                                                                                                          aaaTGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACa                                                                                                                <  1:18169/178‑1 (MQ=255)
                                                                                                                                                                                                           attGAGTGGCAAACCGCCACCCTGGAAGTGGAATGGCGAACCGCTGTTGCCTCGCGAGCAAAAGATCAATTTCTTCAAGCCGATGATCGACAACCACAAGAAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCCCAGTTATTGCCCCGCCGGTGTTTGACTTCGGGATGCGCTGCa                                      <  1:106180/249‑1 (MQ=255)
                                                                                                                                                                                                                     aaaCCGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGCTCGACAACCACCAGCAGGAGTACGGAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGACCAGAGCGATGGTGAAGTGCTGCTCATCCGCCGCACAGTTATTGCCACGCCGGTGTTTGACTTCGGGATGCGCAGCACCTCCATCTg                            >  2:227015/1‑251 (MQ=255)
                                                                                                                                                                                                                        ccGACATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGGGCTGATCATCAACCGCACAGTTATTGCCCCGCCGGTGTTTGACTTCGGGATGCGCTGCACCTCCTTCTGgcg                         >  2:222374/1‑251 (MQ=255)
                                                                                                                                                                                                                         cGACATCCTGGCAGGTGACTGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACTTCGGGATGCGCTGCACCTACATCTGgcgc                        >  1:307889/1‑251 (MQ=255)
                                                                                                                                                                                                                           aCATCCTGGCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCACGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCAGGTGTTGGACTTCGGGATGCGCTGCACCTCCATCTGGCGCAt                      >  2:176767/1‑251 (MQ=255)
                                                                                                                                                                                                                              atcgggtgggTGAATGGGGAATCGCGGTTAACCCGCGAGCAAAAAACCAATTCCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGTCTGGGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGAAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCCTAGTTATTGCCCCGCCGGTGTTTGACTTCGGGATGCGCTGCACCTACATCTGGCGCATc                     <  1:337139/241‑1 (MQ=255)
                                                                                                                                                                                                                                  ggCAGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTAAAGCCGATGCTCGACAACCAAAAGCAGGAGAACGGCGGGCTGTGGCAACCGAATCATTTGCAGACCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGCGCGATGGAGAAGTGCTGATCATCAACCGCACAGTTATTGCCACGCCGGTGATTGACTTTGGGATGCGCGGCACCTACATCTGGCGCATc                     >  2:337139/1‑245 (MQ=255)
                                                                                                                                                                                                                                     aGGTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACTTCGGGATGCGCTGCACCTACATCTGGCGCATCGCTGCCGAt            >  1:542275/1‑251 (MQ=255)
                                                                                                                                                                                                                                      ggTGAATGGCGAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACTTCGGGATGCGCTGCTCCTACATCTGGCGCATCGCTGCCgctg           >  2:388947/1‑248 (MQ=255)
                                                                                                                                                                                                                                         gAATGGCGAATCGCTGATGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGCTGATCACCAACCACAAGCCGGAGTACGAAGGGCTGCGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGAGCGACTTTGCCGCAGAACACAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACTTCGGTATGCGCTAAACCTACATCtcgcacctccctgacgagggcc        >  2:223851/1‑230 (MQ=255)
                                                                                                                                                                                                                                            gggCGAATCGCTGTTGACTCGCGAGCAAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATAATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGGGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACTTCGGGATGCGCTGCACCTACATCTGGCGCATCGCTGCCGATGGCCAgg     <  1:287788/250‑1 (MQ=255)
                                                                                                                                                                                                                                                gAATCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGGGTTTGACTTCGGGATGCGCTGCACCTACATCTGGCGCATCGCTGCCGATGCCCAGGTTa  >  2:221163/1‑250 (MQ=255)
                                                                                                                                                                                                                                                 aaTCGCTGTTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACTTCGGGATGCGCTGCACCTACATCTGGCGCATCGCTGCCGATGGCCa       >  1:58618/1‑244 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                
CAACATTACGGTGACCAAAGATTCCCGCACCCGCTACAGCGAAGCCGGACACCCTATCGCCACTTATCAGTTCCCGCTGAAGGAAAACACCGCGCAGCCAGTGCCTTTCGCACCAAATAATGCGCGTCCGCTGACGCTGGAAGACGATCGTTTGAGCTGCACCGTTCGCGGCTACAACTTCGCGATCACCTTCTCAAAAATGAGTGGCAAACCGACATCCTGGCAGGTGAATGGCGAATCGCTGCTGACTCGCGAGCCAAAGATCAACTTCTTCAAGCCGATGATCGACAACCACAAGCAGGAGTACGAAGGGCTGTGGCAACCGAATCATTTGCAGATCATGCAGGAACATCTGCGCGACTTTGCCGTAGAACAGAGCGATGGTGAAGTGCTGATCATCAGCCGCACAGTTATTGCCCCGCCGGTGTTTGACTTCGGGATGCGCTGCACCTACATCTGGCGCATCGCTGCCGATGGCCAGGTTA  >  NC_000913/3224651‑3225135

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