Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,611,919 C→A 16.7% intergenic (‑65/+47) yneF ← / ← yneG putative membrane‑bound diguanylate cyclase/DUF4186 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,611,9190CA16.7% 232.7 / 4.1 82intergenic (‑65/+47)yneF/yneGputative membrane‑bound diguanylate cyclase/DUF4186 family protein
Reads supporting (aligned to +/- strand):  ref base C (32/35);  new base A (2/13);  total (34/48)
Fisher's exact test for biased strand distribution p-value = 1.93e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCCCCGCTCCGCCAATGACTACCGAGGCAATAACCGACACAATAAGCGCGATGGCGGGCATCATCTGTTCTGCTTTAAAGCGCGGCAGCACTCCAGGAATAGCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTA‑TTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCCGCGCATGGGGGTTTGCTTACCGT  >  NC_000913/1611682‑1612153
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                           
cTCCCCGCTCCGCCAATGACTACCGAGGCAATAACCGACACAATAAGCGCGATGGCGGGCATCATCTGTTCTGCTTTAAAGCGCGGCAGCACTCCAGGAATAGCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCa                                                                                                                                                                                                                                            <  2:312189/239‑1 (MQ=255)
  ccccGCTCCGCCAATGACTACCGAGGCAATAACCGACACAATAAGCGCGATGGCGGGCATCATCTGTTCTGCTTTAAAGCGCGGCAGCACTCCAGGAATAGCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCaaa                                                                                                                                                                                                                                          <  1:1120459/239‑1 (MQ=255)
        tCCGCCAATGACTACCGAGGCAATAACCGACACAATAAGCGCGATGGCGGGCATCATCTGTTCTGCTTTAAAGCGCGGCAGCACTCCAGGAATAGCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAgg                                                                                                                                                                                                                                    <  1:1051191/239‑1 (MQ=255)
         ccGCCAATGACTACCGAGGCAATAACCGACACAATAAGCGCGATGGCGGGCATCATCTGTTCTGCTTTAAAGCGCGGCAGCACTCCAGGAATAGCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGt                                                                                                                                                                                                                                   >  2:488948/1‑239 (MQ=255)
            ccAATGACTACCGAGGCAATAACCGACACAATAAGCGCGATGGCGGGCATCATCTGTTCTGCTTTAAAGCGCGGCAGCACTCCAGGAATAGCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTAc                                                                                                                                                                                                                                 >  1:505373/1‑238 (MQ=255)
                    tACCGAGGCAATAACCGACACAATAAGCGCGATGGCGGGCATCATCTGTTCTGCTTTAAAGCGCGGCAGCACTCCAGGAATAGCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTAACCTGGGAAt                                                                                                                                                                                                                        >  1:1157388/1‑239 (MQ=255)
                    tACCGAGGCAATAACCGACACAATAAGCGCGATGGCGGGCATCATCTGTTCTGCTTTAAAGCGCGGCAGCACTCCAGGAATAGCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTAACCTGGGAAt                                                                                                                                                                                                                        >  1:225864/1‑239 (MQ=255)
                           gCAATAACCGACACAATAAGCGCGATGGCGGGCATCATCTGTTCTGCTTTAAAGCGCGGCAGCACTCCAGGAATAGCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTA‑ttt                                                                                                                                                                                                                >  1:1120598/1‑239 (MQ=255)
                                       acaATAAGCGCGATGGCGGGCATCATCTGTTCTGCTTTAAAGCGCGGCAGCACTCCAGGAATAGCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTA‑TTTTGTCGTTTTCCg                                                                                                                                                                                                    <  2:942782/239‑1 (MQ=255)
                                         aaTAAGCGCGATGGCGGGCATCATCTGTTCTGCTTTAAAGCGCGGCAGCACTCCAGGAATAGCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTA‑TTTTGTCGTTTTCCGCa                                                                                                                                                                                                  <  1:1024978/239‑1 (MQ=255)
                                         aaTAAGCGCGATGGCGGGCATCATCTGTTCTGCTTTAAAGCGCGGCAGCACTCCAGGAATAGCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTA‑TTTTGTCGTTTTCCGCa                                                                                                                                                                                                  <  1:1250702/239‑1 (MQ=255)
                                                gcgATGGCGGGCATCATCTGTTCTGCTTTAAAGCGCGGCAGCACTCCAGGAATAGCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTA‑TTTTGTCGTTTTCCGCAATtct                                                                                                                                                                                             <  2:832308/237‑1 (MQ=255)
                                                   aTGGCGGGCATCATCTGTTCTGCTTTAAAGCGCGGCAGCACTCCAGGAATAGCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTA‑TTTTGTCGTTTTCCGCAATTCTCATg                                                                                                                                                                                         >  1:522235/1‑238 (MQ=255)
                                                   aTGGCGGGCATCATCTGTTCTGCTTTAAAGCGCGGCAGCACTCCAGGAATAGCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTA‑TTTTGTCGTTTTCCGCAATTCTCATg                                                                                                                                                                                         >  2:339651/1‑238 (MQ=255)
                                                        gggCATCATCTGTTCTGCTTTAAAGCGCGGCAGCACTCCAGGAATAGCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTA‑TTTTGTCGTTTTCCGCAATTCTCATGGCTGa                                                                                                                                                                                    <  1:920036/238‑1 (MQ=255)
                                                        gggCATCATCTGTTCTGCTTTAAAGCGCGGCAGCACTCCAGGAATAGCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTA‑TTTTGTCGTTTTCCGCAATTCTCATGGCTGa                                                                                                                                                                                    <  1:756712/238‑1 (MQ=255)
                                                        gggCATCATCTGTTCTGCTTTAAAGCGCGGCAGCACTCCAGGAATAGCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTA‑TTTTGTCGTTTTCCGCAATTCTCATGGCTGa                                                                                                                                                                                    <  2:1088687/238‑1 (MQ=255)
                                                              catcTGTTCTGCTTTAAAGCGCGGCAGCACTCCAGGAATAGCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTA‑TTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCa                                                                                                                                                                                >  1:353096/1‑236 (MQ=255)
                                                                   gTTCTGCTTTAAAGCGCGGCAGCACTCCAGGAATAGCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTA‑TTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACa                                                                                                                                                                        <  2:627322/239‑1 (MQ=255)
                                                                    ttCTGCTTTAAAGCGCGGCAGCACTCCAGGAATAGCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACAAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGAACCCTGGGAATGCTA‑TTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACa                                                                                                                                                                        >  1:361669/1‑238 (MQ=255)
                                                                     tCTGCTTTAAAGCGCGGCAGCACTCCAGGAATAGCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTA‑TTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAAc                                                                                                                                                                      <  1:382561/239‑1 (MQ=255)
                                                                     tCTGCTTTAAAGCGCGGCAGCACTCCAGGAATAGCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTA‑TTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAAc                                                                                                                                                                      <  1:488948/239‑1 (MQ=255)
                                                                     tCTGCTTTAAAGCGCGGCAGCACTCCAGGAATAGCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTA‑TTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAAc                                                                                                                                                                      <  2:1011219/239‑1 (MQ=255)
                                                                     tCTGCTTTAAAGCGCGGCAGCACTCCAGGAATAGCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTA‑TTTTGTCGTTTTCCGCAATGCTCATGGCTGATTCATCTGTACAAc                                                                                                                                                                      <  1:563813/239‑1 (MQ=255)
                                                                        gCTTTAAAGCGCGGCAGCACTCCAGGAATAGCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTA‑TTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAAccag                                                                                                                                                                   <  2:292347/239‑1 (MQ=255)
                                                                         cTTTAAAGCGCGGCAGCACTCCAGGAATAGCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTA‑TTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAAccag                                                                                                                                                                   >  1:893053/1‑238 (MQ=255)
                                                                         cTTTAAAGCGCGGCAGCACTCCAGGAATAGCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTA‑TTTCGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAAccag                                                                                                                                                                   >  2:1248654/1‑238 (MQ=255)
                                                                          tttAAAGCGCGGCAGCACTCCAGGAATAGCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTA‑TTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAAccagc                                                                                                                                                                  <  1:519951/238‑1 (MQ=255)
                                                                                  gcgGCAGCACTCCAGGAATAGCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTA‑TTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAg                                                                                                                                                         <  2:896651/239‑1 (MQ=255)
                                                                                   cgGCAACACTCCAGGAATAGCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTA‑TTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGa                                                                                                                                                        >  1:2034/1‑239 (MQ=255)
                                                                                          aCTCCAGAAATAGCAACCTGCACCAGACAAGGCCCGATCAACACGCCCGTTGAAAATTGCTCACTGAACCAGTCCGCAAGAAAAGGCCAGAAACCCAGACAATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGGGGCGGAAAAAAAAGGTACCCTGGAAATGCTA‑TTTGGTCGTTTCCCGAAATTCTAATGGCTGATTCATCTGTACACCCAGCCATTGATAGATTACCGc                                                                                                                                                 <  2:354448/239‑1 (MQ=255)
                                                                                            tCCAGGAATAGCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATAAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTCCCCTGGAAATGCTA‑TTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTAAAACCAGCCAGTGATAGATTACCGCCa                                                                                                                                               <  1:880670/239‑1 (MQ=255)
                                                                                            cccAGAAATAGCAAAGGTCAGCATACAAGGCACGATCAACACGCCCGTGGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTCCCCTGGGAATGCTA‑TTTTGTCGTTTTCCGCAATTCTCATGGTTGATTCATCTGTACACCCAGCCAGTGATAGATTACCGCCa                                                                                                                                               <  2:493151/238‑1 (MQ=255)
                                                                                            ccaaGGAATAGCCAACGTCAGCAAACAAGCCACGACCAACACGCCCGTGGAGAATTGCTCACTAACCCAGCCCGAAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGAAAATAAAGGTACCCGGGGAATGCTA‑TTTTGTCGTTTTCCGAAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCa                                                                                                                                               <  2:736508/236‑1 (MQ=255)
                                                                                              cAGGAATAGCCAACGACAGCATACAAGGCCCGATCAACACGCCCGTTGAGAATTGCTCACTGACCCAGTCCGCAAGAAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTA‑TTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACg                                                                                                                                             <  1:566159/239‑1 (MQ=255)
                                                                                                    tAGCCACCGTCAGCATACAAGGCACAATCAACCCGCCCGTTGAGAATTGCTCACGAAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTA‑TTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATaa                                                                                                                                       <  1:718979/239‑1 (MQ=255)
                                                                                                    tAGCCAACGCCAGCAAAAAAGGCACGATCAACACGCCCTTTGAAAATTGCTCACTGAACCAGTCCGCAAGAAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGAAAATAAAGGTACCCTGGGAATGCTA‑TTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATaa                                                                                                                                       <  2:1120598/239‑1 (MQ=255)
                                                                                                     aGCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTA‑TTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAAc                                                                                                                                      >  2:221013/1‑239 (MQ=255)
                                                                                                      gacAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCAGTGAACCAGTCCGCAAGCAAAGGCCAGAAATCAAGACTATCAATACTGACCGAATCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAACGGTAACCTGGGAACGCTA‑TTTTGTCGTTTTCCGCAATCCCCTAGACAGATTCATCTGTACAGCCAGCCAGTGAAAGATTACCGCCACGATATAACg                                                                                                                                     >  2:281734/3‑239 (MQ=255)
                                                                                                      gCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAGGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTA‑TTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACg                                                                                                                                     >  1:976750/1‑239 (MQ=255)
                                                                                                      gCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTA‑TTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATACCg                                                                                                                                     >  2:1051881/1‑239 (MQ=255)
                                                                                                      gCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTA‑TTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACg                                                                                                                                     >  1:395588/1‑239 (MQ=255)
                                                                                                      gCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTA‑TTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACg                                                                                                                                     >  1:947413/1‑239 (MQ=255)
                                                                                                         aaCGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTA‑TTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACgttg                                                                                                                                  <  1:279611/239‑1 (MQ=255)
                                                                                                                 cATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGAAAGCAAAGGCCAGAAATCCAGCCTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTA‑TTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGttctt                                                                                                                          <  1:66415/239‑1 (MQ=255)
                                                                                                                     caagacaagaTCAACCCGCCCTTTGAGAATTGTCCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAATCAAGATAGCATGTTGGCTAACGGCGGCGGAAAAAAAAGGTCCCCTGGGAATGCTA‑TTTTTTCGTTTTCCGAAATTCTCAGGGCTGATTCATCTGTACACCCAGCAAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCt                                                                                                                      <  1:281734/239‑1 (MQ=255)
                                                                                                                     cAAGGCACAACCAACACGCCCGTGGAGAATGGCTCACTGAACCAGCCCGCAAGCAAAGGCAAGAAACCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGAAAATAAAGGTACCTTGGGAAGCCTA‑TTTTGTCGTTTCCCGCAATTCTCATGGCTGATTCATCTGTACACCCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCt                                                                                                                      <  1:70498/239‑1 (MQ=255)
                                                                                                                       aGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTA‑TTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTa                                                                                                                    >  1:687477/1‑239 (MQ=255)
                                                                                                                       aGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGAAAGCAAAGGCCAAAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTA‑TTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTa                                                                                                                    <  1:221013/239‑1 (MQ=255)
                                                                                                                             gaacAACACCCCTGTTAAGATTTGCTCACTGAACCATTCAGAAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCAACCAAAGATACCATTTTCCCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATTCTA‑TTTTGTCGTTTTCCGCAATGCTCATGGCTGATTCATCTGTACAACCAGCCAGTAATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGaaa                                                                                                              <  2:1157388/236‑1 (MQ=255)
                                                                                                                             gATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTA‑TTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGaa                                                                                                               >  1:639553/1‑238 (MQ=255)
                                                                                                                              aTCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTA‑TTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAac                                                                                                             >  2:892982/1‑239 (MQ=255)
                                                                                                                                caccacGCCGGTTGAGAATGCCCCACCGACCCAGTCCGAAAGCAAAGGCCAGAAAACCAGACTACCAATACTGACCCAACCAAGATAGCATGTTCGCTGCCAGCGGCGGAAAATAAAGGTCCCTGGGGAATGTTT‑TTGTGTCGTTTTCCGCAATCCGCATGGCTTATTCTTCTGTACAACCAGCCAGTGATAGATGACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAacac                                                                                                           <  2:361669/236‑1 (MQ=255)
                                                                                                                                 aaCACGCCCGTTGAAAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTA‑TTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACAcc                                                                                                          <  1:1063773/239‑1 (MQ=255)
                                                                                                                                   cacGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGAATGTTCGCTGAAAGCGGCGGAAAACAAAGGTACCCTGGGAATGATA‑TTTTGACGTTTTCCGCAATACTCATGGCTGATTCATCTGTACAACAAGCCAGTGATAGATTACAGCACCGATATAACGATGTAGTTCTTCGCTTAACGCAacaa                                                                                                           >  1:671690/1‑235 (MQ=255)
                                                                                                                                      gcccttgaagaATTGCTCACTAACCCGGTCCGCAAACAAAGCCCAAAAATCCAGCCTACCAATACGGACCGAACCAAGATAGAATGTGCGCTGACAGCGCGGGCAAATAAAGGTACCCTGGAAATGCTA‑TTTTGACGTTTTCCGCAATTCCCATGGCTGATTCATCTGTACAACCAGCCAGTAATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCg                                                                                                     <  2:105692/231‑1 (MQ=255)
                                                                                                                                      ccccTTTGAGAATTGCTCACTAAACCAGCCCGAAAGAAAAGGCAAAAAATCAAGCCTATAAATACTGACCAAACCAAGATAGCATGTCCGGTAACAGCGGCGGAAAATAAAGGTCCCCTGGGAATGCTA‑TTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACACCCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCg                                                                                                     <  2:177706/238‑1 (MQ=255)
                                                                                                                                       cccttcgaaaaTTGCTCACTGAACCAGTCCGCAAGCAAAGGCAAGAAAGCCAGACTAACAATACTGACCGAACCAAGATAGCATGTTCGCTGCAAGCGGCGGCAAATAAAGGTCCCTTGGGAATGCTA‑TTTTGTCGTTTGCCGCAATTCTCATGGCTGATTCATCTGTACAACCGGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCgg                                                                                                    <  2:715278/230‑1 (MQ=255)
                                                                                                                                           ttGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGAACCCTGGGAATGCTA‑TTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATACATTACCGCCACGTTATACCGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGCAATg                                                                                                >  2:227492/1‑239 (MQ=255)
                                                                                                                                               aaaTTGCTCACTGAACCAGTCCGAAACCAAAGGCCAGAAACCAAGTCTATAAATACTGACCAAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAAAAAAGGTACCCTGGGAATGCTA‑TTTTGTCGTTTTCCGCAATCCTCATGGCTGATTCATCTGAACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTAt                                                                                            <  2:522235/238‑1 (MQ=255)
                                                                                                                                                   tGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCGGGGAATGCTA‑TTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCa                                                                                        <  1:339651/239‑1 (MQ=255)
                                                                                                                                                      tCACTGAACCAGTCCGCAAGCAAGGCCCAGAAATCAAGTCTATAAATCCTGACCGAACCAAGAATCCATGTTCGTTGACAGCGGGGTCAAATAAAGGTACCCTGGGAATCTTT‑TTTCGTCGTTTTCCACAATCCTCGTGCCAAATTCCTCTGTACACCCAGCAAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACtt                                                                                     <  2:137931/239‑1 (MQ=255)
                                                                                                                                                              cctgcccGAAAACAAAGGCCAGAAATCAAGATTACAAATACCGCCCGAACCAAGATAGCATTTTCGCTGACAGCGGAGAAAAATAAGGGTCCCCTGGAAATGCTT‑TTTTGTCGTTTTCCGAATTTCTCATGGCTGATTCATCTGTACAGCCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGAc                                                                             <  2:960967/234‑1 (MQ=255)
                                                                                                                                                               cAGTCGCCAAGCAAAGACAAAAATACCGGCTTATCAATCCTGACCGAACAAAGCTACCATGTCCGATGCCGGCGGCGCAAAATAAAGTTCCCCGGGGAATGTTT‑TTTTGTCTTTTCCCGCATTTCTCATGTTTTTTTCTTCTGTACAACCACCCAGTGATATATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACa                                                                            <  2:1077866/239‑1 (MQ=255)
                                                                                                                                                                      caagcaaAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGAAAATAAAGGTACCCTGGGAATGCTA‑TTTTGTCGTTTTCCGCGATTATCATGGCTGATTCATCTGTACAACCAGCCAGCGATACCTTACCGCCACGATATACCGTTGTTGTTGTTCGCTTAACGAAACACCGTGCGGACTGTTATGCCACTTAGCCAGaacaccccg                                                                      >  2:509436/1‑230 (MQ=255)
                                                                                                                                                                                           gACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTA‑TTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAAccagc                                                                                                                                                                  >  2:487428/1‑125 (MQ=255)
                                                                                                                                                                                           gACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTA‑TTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAAccagc                                                                                                                                                                  <  1:487428/125‑1 (MQ=255)
                                                                                                                                                                                              tATCAATACTGACCGAACCAAGATAGCATGTTCGCGACCAGCGGCGGAAAATAAAGGTACCCTGGGAATGCTA‑TTTTGTCGTTTGCCGCAATTCTCATGGCTGATTCACCTGTAAACCCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTGCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGa                                             <  1:509436/239‑1 (MQ=255)
                                                                                                                                                                                               ataaaTACTGACGGACCCAAGAAAGTGTGTTCGTTGACAGCGGGGGCAAATAAAGTTTCCTGGGAAATGTTT‑TTTTGTCGTTTCCCGCAATCCTCATGGTTGATTCATCTTTCCATCCACCCGGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTGCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGGCAGCAGGTCGCGGTGGCGTGCTGAg                                            <  2:893053/236‑1 (MQ=255)
                                                                                                                                                                                                      cTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTA‑TTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATaa                                      >  1:695911/1‑238 (MQ=255)
                                                                                                                                                                                                        gACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTA‑TTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAac                                    >  1:704133/1‑238 (MQ=255)
                                                                                                                                                                                                        gACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTA‑TTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAac                                    >  1:35160/1‑238 (MQ=255)
                                                                                                                                                                                                        gACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTA‑TTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGAGCTGAGCGATAAac                                    >  2:865474/1‑238 (MQ=255)
                                                                                                                                                                                                          ccGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCa                                                                                                                                                                                                                                            <  2:59551/37‑1 (MQ=255)
                                                                                                                                                                                                          ccGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCa                                                                                                                                                                                                                                            >  1:59551/1‑37 (MQ=255)
                                                                                                                                                                                                             aaCCAAGATAGCATGTTCCCGGACGGCGGCGGAAAATAAAGGTACCCTGGGAATGCTA‑TTTTGTCGTTTCCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGa                              <  2:687477/239‑1 (MQ=255)
                                                                                                                                                                                                                   gATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTA‑TTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGccc                         >  1:56854/1‑238 (MQ=255)
                                                                                                                                                                                                                        cATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTA‑TTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCCGCGCAt                   >  1:1134214/1‑239 (MQ=255)
                                                                                                                                                                                                                         aTGTTCGCTGACAGGGGCGGAAAAAAAAGGTCCCCTGAGAATGCTA‑TTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCAGGTAAAAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCCGCGCa                    <  2:138941/237‑1 (MQ=255)
                                                                                                                                                                                                                            ttCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTA‑TTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCCGCGCATgggg               >  1:826586/1‑239 (MQ=255)
                                                                                                                                                                                                                               acttccggcggcgCAAAAAAAAGGTACCCGGTGATTCCTT‑TTTTGTCTTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAAATTCCCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACGCCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCCGCGCATGGGGGtt            <  2:2034/232‑1 (MQ=255)
                                                                                                                                                                                                                                       cggcggCAAATAAAGGTACCCTGGGAATGCTA‑TTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGCTAAACACCGGATGCCCGCGCATGGGGGTTTGCTTAcc    >  1:404812/1‑239 (MQ=255)
                                                                                                                                                                                                                                         gcggcAAATAAAGGTACCCTGGGAATGCTA‑TTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCCGCGCATGGGGGTTTGCTTACCg   >  1:166004/1‑238 (MQ=255)
                                                                                                                                                                                                                                          cggcAAATAAA‑GTTCCCCGGGGATGGTATTTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAAATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCCGCGCATGGGGGTTTGCTTACCGt  <  1:5840/238‑1 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                           
CTCCCCGCTCCGCCAATGACTACCGAGGCAATAACCGACACAATAAGCGCGATGGCGGGCATCATCTGTTCTGCTTTAAAGCGCGGCAGCACTCCAGGAATAGCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTA‑TTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCCGCGCATGGGGGTTTGCTTACCGT  >  NC_000913/1611682‑1612153

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