Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,473,900 C→T G52R (GGA→AGA)  bdcA ← c‑di‑GMP‑binding biofilm dispersal mediator protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,473,9000CT90.3% 295.5 / 15.1 103G52R (GGA→AGA) bdcAc‑di‑GMP‑binding biofilm dispersal mediator protein
Reads supporting (aligned to +/- strand):  ref base C (4/6);  new base T (41/52);  total (45/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.62e-01

Only 100 of 104 total aligned reads displayed.

ATTCACGGAGCCGATGATTAAGATGCGCCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCT  >  NC_000913/4473652‑4474144
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     
aTTCACGGAGCCGATGATTAAGATGCGCCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTcc                                                                                                                                                                                                                                                       >  1:689843/1‑249 (MQ=255)
   cacggagcgggctatTAAGGTGCGCCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCAGGATAAGGAGCATTAATATTGAGTTTGAAAAGGCGATCAATAGCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATTCCTGCATTTACCACCAGGATAGCCAATGCACCGCTCTTACGAACGACATCAATGACAGCTTCTCTGGCAGCACTATCTGTGAATACTGCTGTCGCGCTA‑Gtct                                                                                                                                                                                                                                                 <  2:390523/238‑1 (MQ=255)
         gCCGATGATTAAGATGCGCCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTtgtg                                                                                                                                                                                                                                           >  1:382921/1‑251 (MQ=255)
          ccGATGATTAAGATGCGCCCGGCGTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATTCCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTAGCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGc                                                                                                                                                                                                                                          <  2:392077/251‑1 (MQ=255)
           cGATGATTAAGATGCGCCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGcc                                                                                                                                                                                                                                         >  2:866612/1‑251 (MQ=255)
           cGATGAGTAAGATGCGCCTGCCTTTGGGCATATGCCGGGCAGCTTCAACAGAGGCATGATAAGGATAATTAATATTGATTGTGAAAAGGCGATCAATATTGTCGACATTTAATGCCAGGGCCGCGCCAGAGACGCAAATACCTGCATTTACCCCCAGGATATCCAATGCGCCGTTCCTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCCCTA‑GTCTCTTTTGcc                                                                                                                                                                                                                                         <  2:620599/251‑1 (MQ=255)
            gatgatTAAGATGCGCCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCa                                                                                                                                                                                                                                        <  1:876216/251‑1 (MQ=255)
               gatTAAGATGCGCCCGCCGTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCGCTA‑GTCTCTTGTGCCAGGc                                                                                                                                                                                                                                     <  2:963710/251‑1 (MQ=255)
                   aaGATGCGCCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGttt                                                                                                                                                                                                                                 <  1:641835/251‑1 (MQ=255)
                    aGATGCGCCCACCTTCAGGCATCTGCCGGGCGGCTTCAACCGCGGCATGATAAGAAGCATGAATCTTGACTTTGAAAATGAGATCAATATCGTCGGCATTTAATTCCAGGGCATCGCCAAAAACGCCAATACCTGCATTTTCCACCAGGATATCCAATGCCCCGCCCTTACGAACGACATCAATGACAGCGTCTCTGTCAACAATCTCTGTGACTACTGCTTACGATCTACGT‑TCTTGTGCCAGGCGtttg                                                                                                                                                                                                                                >  2:922444/1‑250 (MQ=255)
                          gcccgccTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTTCCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCa                                                                                                                                                                                                                          <  1:350313/251‑1 (MQ=255)
                           cctgccTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAATGAGCATGAAGATTGATTTGGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATTCCTGCATTTTCCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATTCTGCTGTCGCGCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCAt                                                                                                                                                                                                                         <  1:736808/248‑1 (MQ=255)
                            ccgccTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATc                                                                                                                                                                                                                        <  1:430812/251‑1 (MQ=255)
                              gccTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAATGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCtt                                                                                                                                                                                                                      >  2:772496/1‑251 (MQ=255)
                              gccTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCtt                                                                                                                                                                                                                      >  2:688416/1‑251 (MQ=255)
                               ccTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCCGGCGTTTAGCGGAATCttt                                                                                                                                                                                                                     >  2:175931/1‑251 (MQ=255)
                                cttagggCATCTGCCGGACGGCTTCAACAGAGGCATGATAAGGAGCCTGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTACTTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCAACAGGATATCCAATGCGCAGCTCTTACGAACGACCTCCATGACAGCGTCTCTGTCAGCACTATCTGTGAATCATGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGGTTAGCGGCATCTTTc                                                                                                                                                                                                                    >  2:181005/5‑251 (MQ=255)
                                           tGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCAt                                                                                                                                                                                                          >  2:462773/1‑250 (MQ=255)
                                                ggcggcTTCAACAGAGGCATGATAAGGAGCATGAATATTGAGTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGa                                                                                                                                                                                                    <  2:218784/251‑1 (MQ=255)
                                                   ggcTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATc                                                                                                                                                                                                 <  1:505778/251‑1 (MQ=255)
                                                    gcTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGACTCg                                                                                                                                                                                                >  2:136128/1‑251 (MQ=255)
                                                            agagGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATTAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATTCCTGCATTTACCACCAGGATATCCAATGCGCCGCTCGTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTgg                                                                                                                                                                                        <  2:723245/251‑1 (MQ=255)
                                                            agagGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCTCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCCCTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTgg                                                                                                                                                                                        >  1:375857/1‑251 (MQ=255)
                                                                    gATAAGGAGCATGAATATTGATTTTGAAAAGGCAATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCA‑GTCTCTTGTGCCAGGCGt                                                                                                                                                                                                                                   <  1:399182/200‑1 (MQ=255)
                                                                        aGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTg                                                                                                                                                                             >  1:131623/1‑250 (MQ=255)
                                                                        aGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGa                                                                                                                                                                            <  1:417939/251‑1 (MQ=255)
                                                                             cATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGaa                                                                                                                                                                        >  1:592263/1‑250 (MQ=255)
                                                                              aTGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAac                                                                                                                                                                      <  1:461528/251‑1 (MQ=255)
                                                                               tGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAacg                                                                                                                                                                     >  1:156281/1‑251 (MQ=255)
                                                                               tGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGACTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAacg                                                                                                                                                                     >  1:484406/1‑251 (MQ=255)
                                                                                gAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAacg                                                                                                                                                                     >  1:206771/1‑250 (MQ=255)
                                                                                   tatTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAacgaagt                                                                                                                                                                 >  2:803424/1‑248 (MQ=255)
                                                                                     ttGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGGCGCTCCA‑GTCTCTTGTGCCAGTCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTa                                                                                                                                                                <  2:689843/250‑1 (MQ=255)
                                                                                         ttttGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGAt                                                                                                                                                            >  1:16891/1‑250 (MQ=255)
                                                                                          tttGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGAt                                                                                                                                                            >  1:276052/1‑249 (MQ=255)
                                                                                           ttGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAg                                                                                                                                                          <  1:140484/250‑1 (MQ=255)
                                                                                            tGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATTCCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGc                                                                                                                                                         <  1:669022/250‑1 (MQ=255)
                                                                                                aaGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCAc                                                                                                                                                    <  1:402673/251‑1 (MQ=255)
                                                                                                 aGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCACTGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTTGCCCCATCGGTGACGAAACGACGTACGATAGCGGcccc                                                                                                                                                   >  2:66597/1‑248 (MQ=255)
                                                                                                    cGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGAt                                                                                                                                                >  2:8638/1‑251 (MQ=255)
                                                                                                     gATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGAt                                                                                                                                                >  1:265036/1‑250 (MQ=255)
                                                                                                            atCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGAGCCGCTCTTACGAACGACATCCATGACAGCGTCTCTGTCAGCCCTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGACTCGTACATTGGCCCCATCGGTGACGGAACGACGTACGATAGCGGACCCGATAcccagac                                                                                                                                        >  2:407873/1‑246 (MQ=255)
                                                                                                             gcggcgGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGGCAGCACTAGCTGTGAATACTGCTGGCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTGCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACt                                                                                                                                       <  2:319379/247‑1 (MQ=255)
                                                                                                              cgtcgGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACt                                                                                                                                       <  1:862949/250‑1 (MQ=255)
                                                                                                                          aTTCCAGGGCCCCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTAGCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATg                                                                                                                          <  2:75899/251‑1 (MQ=255)
                                                                                                                                      cGCCAAAGCCGCCAATGCCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCTTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTa                                                                                                              <  2:587306/251‑1 (MQ=255)
                                                                                                                                        ccAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCAACCGCGGATGAGAACTGTCTTAcc                                                                                                            >  2:726707/1‑251 (MQ=255)
                                                                                                                                         cAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCt                                                                                                           >  1:477383/1‑251 (MQ=255)
                                                                                                                                         cAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCt                                                                                                           >  1:608374/1‑251 (MQ=255)
                                                                                                                                                gCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTaaa                                                                                                      >  1:6165/1‑249 (MQ=255)
                                                                                                                                                   aaTACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGGCACATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAgcgc                                                                                                 <  1:411012/251‑1 (MQ=255)
                                                                                                                                                       ccTGCCTTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCAt                                                                                             <  2:365245/251‑1 (MQ=255)
                                                                                                                                                         tGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATa                                                                                            >  2:320792/1‑250 (MQ=255)
                                                                                                                                                          gCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATa                                                                                            <  1:439183/249‑1 (MQ=255)
                                                                                                                                                          gCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAg                                                                                          <  1:749526/251‑1 (MQ=255)
                                                                                                                                                          gCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAg                                                                                          <  1:585523/251‑1 (MQ=255)
                                                                                                                                                           cATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAgg                                                                                         <  1:799981/251‑1 (MQ=255)
                                                                                                                                                           cATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACAGATACCACGACTGCCACCGAGGAAGAGAACTGTCTTACCTGTACAAGCGCCCATAAgg                                                                                         >  2:266092/1‑251 (MQ=255)
                                                                                                                                                              ttACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGCCAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCTTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACCGCCACCGAGGATGAGATCCGTCTTACCTGAAATAGCGCCCATAAgaccc                                                                                      >  1:801278/1‑247 (MQ=255)
                                                                                                                                                              ttACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGAc                                                                                       <  1:989837/250‑1 (MQ=255)
                                                                                                                                                                  cacccGGATATCCAATTCGCCGCTCCTACGAACGACATCAATGCCAGCGTCTCTGTCAGCACTATCTGTGAATGCTGCTGTCGCTCTA‑GGCTCTTGTGCCAGGCGTTTAGCGGCAGTTTTCGAGCTCGCATAAGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACTTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCtt                                                                                  <  2:960467/251‑1 (MQ=255)
                                                                                                                                                                    ccaGGATTTCCAATGCGCCGATCTTACGAACTACATCAATGACAGCGTCTCTGTCAGCACTATCTGGGAATACTGCGGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTGAGCGGCATCTTTCGACCCCGCATAGGTTAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCTATAAGGACTCCtt                                                                                  <  2:326884/249‑1 (MQ=255)
                                                                                                                                                                     cagagattCCAATGCGCCGCGCTTACGAACGACATCAATGACAGCGTCGCTGGTAGCACTATCTGTGAATACTGCTGACGCTCTA‑GTCTCTTGTGCTAGGCGTTTAGCTGCATCTTTCGCCCCCGCATAGGTGAATCGTACATTGGCCACATCGGTGACGAAACGACTTACGATAGCGGCACCGATACCACGACTGCCACCTAGGATGATAACTGTCTTACCTGTAAAAGCGCCCATAAGGCCTCCTTGAt                                                                               <  2:115198/244‑1 (MQ=255)
                                                                                                                                                                      aGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGAtt                                                                              <  1:739538/251‑1 (MQ=255)
                                                                                                                                                                      aGGATATCCAATGCGCCGCTCGTACGAACGACAACAATGACAGCTGCTCTGTCAGCACTATCTGTGAATACTGCTGGCGCTCTA‑GTCGCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTAACCTGTAAAAGCGCCCATAAGGACTCCTTGAtt                                                                              <  2:492579/251‑1 (MQ=255)
                                                                                                                                                                           atCCAATGCGACGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCCTCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGTTGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGCGGATGAGAACAGTCTTACCTGTAAAAGCGCCCATAAGAACTCCTTGAtttttt                                                                          >  2:623839/1‑247 (MQ=255)
                                                                                                                                                                            tCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTtattat                                                                        >  2:119093/1‑251 (MQ=255)
                                                                                                                                                                             ccAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATg                                                                       >  1:638418/1‑251 (MQ=255)
                                                                                                                                                                              cAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGt                                                                      <  1:246918/251‑1 (MQ=255)
                                                                                                                                                                              cAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACC‑AGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTa                                                                     <  1:11447/251‑1 (MQ=255)
                                                                                                                                                                               aaTGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTAGATTTATTATGTa                                                                     <  1:462773/251‑1 (MQ=255)
                                                                                                                                                                                  gcgcCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAAc                                                                   <  1:591474/250‑1 (MQ=255)
                                                                                                                                                                                   ctccgcTCTTACGAACGACATCAATGCCAGCGTCTCTGTCAGCACTCTCTGTGAATACTGCTGTCTCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCTATACCACGACTGCCACCGAGGCTGAGAACTGTCTTACCTGTACACGCGCCCATAAGGACTCCTTGATTTATTATGTAACAt                                                                 >  2:611662/3‑251 (MQ=255)
                                                                                                                                                                                    gccgcTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATg                                                                >  1:33874/1‑251 (MQ=255)
                                                                                                                                                                                     ccgcTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGc                                                               <  1:304744/251‑1 (MQ=255)
                                                                                                                                                                                        ctctTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCAtt                                                            >  1:349373/1‑251 (MQ=255)
                                                                                                                                                                                            tACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACa                                                         >  1:233834/1‑250 (MQ=255)
                                                                                                                                                                                            tACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGGCGCTCCA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTa                                                           <  2:427146/248‑1 (MQ=255)
                                                                                                                                                                                             acgaacgaCATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCCTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACaa                                                        >  1:856901/1‑250 (MQ=255)
                                                                                                                                                                                              cgaacgaCATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACaaaa                                                      <  1:866612/251‑1 (MQ=255)
                                                                                                                                                                                              cgaacgaCATCAAGGTCAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTTAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACaaaa                                                      <  2:405101/251‑1 (MQ=255)
                                                                                                                                                                                              cgaacgaCAGCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTGTCGACCCCGCATAGGTGAAGCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGGAACATGCATTGCaaaa                                                      <  2:461590/251‑1 (MQ=255)
                                                                                                                                                                                                aacgaCATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACt                                                    >  1:424345/1‑251 (MQ=255)
                                                                                                                                                                                                      aTCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTa                                              <  1:771717/251‑1 (MQ=255)
                                                                                                                                                                                                         aaTGACAGTGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGTCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAACCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACGGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACt                                           <  2:866879/251‑1 (MQ=255)
                                                                                                                                                                                                            gACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGCAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTc                                        >  1:493285/1‑251 (MQ=255)
                                                                                                                                                                                                             aCAGCGTTTCTTGCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGTCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCTCCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCt                                       <  2:709671/251‑1 (MQ=255)
                                                                                                                                                                                                              cAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGccc                                                                                               <  2:202771/194‑1 (MQ=255)
                                                                                                                                                                                                              cAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGccc                                                                                               >  1:202771/1‑194 (MQ=255)
                                                                                                                                                                                                               aGCGTCTCTGTCAGCACTATCTGTGACTACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATATTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGCTAGCGGCACCGCTACCACGACTGCCACCGAGGATGAGCACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGTTTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGt                                     >  2:900878/1‑251 (MQ=255)
                                                                                                                                                                                                               aGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGt                                     <  1:987100/251‑1 (MQ=255)
                                                                                                                                                                                                               aGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGt                                     >  1:956940/1‑251 (MQ=255)
                                                                                                                                                                                                                           aGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAATCTGTTTTAACTTTCTGTCAACAGGTTTTg                         <  1:210133/251‑1 (MQ=255)
                                                                                                                                                                                                                           aGCACTATCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTg                         <  1:688416/251‑1 (MQ=255)
                                                                                                                                                                                                                                 aTCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAAt                     <  1:136128/249‑1 (MQ=255)
                                                                                                                                                                                                                                 aTCTGTGAATACTGCTGTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATg                    <  1:607269/250‑1 (MQ=255)
                                                                                                                                                                                                                                     gtgAATTCTGCTGGCGCTCCA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTGCGACCCCGCATAGGTGAATCCTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCAc               <  2:856901/251‑1 (MQ=255)
                                                                                                                                                                                                                                              gctgTCGCTCTA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACa      <  1:129977/251‑1 (MQ=255)
                                                                                                                                                                                                                                                tgTCGCTCCA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGCCGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACCAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAAt    >  2:480407/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   cGCTCTA‑GTCTCTTGCGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAAtct  <  1:932034/250‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     
ATTCACGGAGCCGATGATTAAGATGCGCCCGCCTTCGGGCATCTGCCGGGCGGCTTCAACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCA‑GTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTATGTAACATGCATTACAAAACTGTTTTAACTTTCTGTCAACAGGTTTTGTAATGGTCACTAAAAAACAATCT  >  NC_000913/4473652‑4474144

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑