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

CTCCCCGCTCCGCCAATGACTACCGAGGCAATAACCGACACAATAAGCGCGATGGCGGGCATCATCTGTTCTGCTTTAAAGCGCGGCAGCACTCCAGGAATAGCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCCGCGCATGGGGGTTTGCTTACCGT  >  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)
                           gCAATAACCGACACAATAAGCGCGATGGCGGGCATCATCTGTTCTGCTTTAAAGCGCGGCAGCACTCCAGGAATAGCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTAttt                                                                                                                                                                                                                >  1:1120598/1‑239 (MQ=255)
                                       acaATAAGCGCGATGGCGGGCATCATCTGTTCTGCTTTAAAGCGCGGCAGCACTCCAGGAATAGCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCg                                                                                                                                                                                                    <  2:942782/239‑1 (MQ=255)
                                         aaTAAGCGCGATGGCGGGCATCATCTGTTCTGCTTTAAAGCGCGGCAGCACTCCAGGAATAGCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCa                                                                                                                                                                                                  <  1:1024978/239‑1 (MQ=255)
                                         aaTAAGCGCGATGGCGGGCATCATCTGTTCTGCTTTAAAGCGCGGCAGCACTCCAGGAATAGCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCa                                                                                                                                                                                                  <  1:1250702/239‑1 (MQ=255)
                                                gcgATGGCGGGCATCATCTGTTCTGCTTTAAAGCGCGGCAGCACTCCAGGAATAGCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATtct                                                                                                                                                                                             <  2:832308/237‑1 (MQ=255)
                                                   aTGGCGGGCATCATCTGTTCTGCTTTAAAGCGCGGCAGCACTCCAGGAATAGCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATg                                                                                                                                                                                         >  1:522235/1‑238 (MQ=255)
                                                   aTGGCGGGCATCATCTGTTCTGCTTTAAAGCGCGGCAGCACTCCAGGAATAGCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATg                                                                                                                                                                                         >  2:339651/1‑238 (MQ=255)
                                                        gggCATCATCTGTTCTGCTTTAAAGCGCGGCAGCACTCCAGGAATAGCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGa                                                                                                                                                                                    <  1:920036/238‑1 (MQ=255)
                                                        gggCATCATCTGTTCTGCTTTAAAGCGCGGCAGCACTCCAGGAATAGCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGa                                                                                                                                                                                    <  1:756712/238‑1 (MQ=255)
                                                        gggCATCATCTGTTCTGCTTTAAAGCGCGGCAGCACTCCAGGAATAGCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGa                                                                                                                                                                                    <  2:1088687/238‑1 (MQ=255)
                                                              catcTGTTCTGCTTTAAAGCGCGGCAGCACTCCAGGAATAGCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCa                                                                                                                                                                                >  1:353096/1‑236 (MQ=255)
                                                                   gTTCTGCTTTAAAGCGCGGCAGCACTCCAGGAATAGCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACa                                                                                                                                                                        <  2:627322/239‑1 (MQ=255)
                                                                    ttCTGCTTTAAAGCGCGGCAGCACTCCAGGAATAGCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACAAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGAACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACa                                                                                                                                                                        >  1:361669/1‑238 (MQ=255)
                                                                     tCTGCTTTAAAGCGCGGCAGCACTCCAGGAATAGCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAAc                                                                                                                                                                      <  1:382561/239‑1 (MQ=255)
                                                                     tCTGCTTTAAAGCGCGGCAGCACTCCAGGAATAGCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAAc                                                                                                                                                                      <  1:488948/239‑1 (MQ=255)
                                                                     tCTGCTTTAAAGCGCGGCAGCACTCCAGGAATAGCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAAc                                                                                                                                                                      <  2:1011219/239‑1 (MQ=255)
                                                                     tCTGCTTTAAAGCGCGGCAGCACTCCAGGAATAGCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATGCTCATGGCTGATTCATCTGTACAAc                                                                                                                                                                      <  1:563813/239‑1 (MQ=255)
                                                                        gCTTTAAAGCGCGGCAGCACTCCAGGAATAGCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAAccag                                                                                                                                                                   <  2:292347/239‑1 (MQ=255)
                                                                         cTTTAAAGCGCGGCAGCACTCCAGGAATAGCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAAccag                                                                                                                                                                   >  1:893053/1‑238 (MQ=255)
                                                                         cTTTAAAGCGCGGCAGCACTCCAGGAATAGCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTCGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAAccag                                                                                                                                                                   >  2:1248654/1‑238 (MQ=255)
                                                                          tttAAAGCGCGGCAGCACTCCAGGAATAGCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAAccagc                                                                                                                                                                  <  1:519951/238‑1 (MQ=255)
                                                                                  gcgGCAGCACTCCAGGAATAGCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAg                                                                                                                                                         <  2:896651/239‑1 (MQ=255)
                                                                                   cgGCAACACTCCAGGAATAGCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGa                                                                                                                                                        >  1:2034/1‑239 (MQ=255)
                                                                                          aCTCCAGAAATAGCAACCTGCACCAGACAAGGCCCGATCAACACGCCCGTTGAAAATTGCTCACTGAACCAGTCCGCAAGAAAAGGCCAGAAACCCAGACAATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGGGGCGGAAAAAAAAGGTACCCTGGAAATGCTATTTGGTCGTTTCCCGAAATTCTAATGGCTGATTCATCTGTACACCCAGCCATTGATAGATTACCGc                                                                                                                                                 <  2:354448/239‑1 (MQ=255)
                                                                                            tCCAGGAATAGCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATAAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTCCCCTGGAAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTAAAACCAGCCAGTGATAGATTACCGCCa                                                                                                                                               <  1:880670/239‑1 (MQ=255)
                                                                                            cccAGAAATAGCAAAGGTCAGCATACAAGGCACGATCAACACGCCCGTGGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTCCCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGTTGATTCATCTGTACACCCAGCCAGTGATAGATTACCGCCa                                                                                                                                               <  2:493151/238‑1 (MQ=255)
                                                                                            ccaaGGAATAGCCAACGTCAGCAAACAAGCCACGACCAACACGCCCGTGGAGAATTGCTCACTAACCCAGCCCGAAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGAAAATAAAGGTACCCGGGGAATGCTATTTTGTCGTTTTCCGAAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCa                                                                                                                                               <  2:736508/236‑1 (MQ=255)
                                                                                              cAGGAATAGCCAACGACAGCATACAAGGCCCGATCAACACGCCCGTTGAGAATTGCTCACTGACCCAGTCCGCAAGAAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACg                                                                                                                                             <  1:566159/239‑1 (MQ=255)
                                                                                                    tAGCCACCGTCAGCATACAAGGCACAATCAACCCGCCCGTTGAGAATTGCTCACGAAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATaa                                                                                                                                       <  1:718979/239‑1 (MQ=255)
                                                                                                    tAGCCAACGCCAGCAAAAAAGGCACGATCAACACGCCCTTTGAAAATTGCTCACTGAACCAGTCCGCAAGAAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGAAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATaa                                                                                                                                       <  2:1120598/239‑1 (MQ=255)
                                                                                                     aGCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAAc                                                                                                                                      >  2:221013/1‑239 (MQ=255)
                                                                                                      gacAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCAGTGAACCAGTCCGCAAGCAAAGGCCAGAAATCAAGACTATCAATACTGACCGAATCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAACGGTAACCTGGGAACGCTATTTTGTCGTTTTCCGCAATCCCCTAGACAGATTCATCTGTACAGCCAGCCAGTGAAAGATTACCGCCACGATATAACg                                                                                                                                     >  2:281734/3‑239 (MQ=255)
                                                                                                      gCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAGGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACg                                                                                                                                     >  1:976750/1‑239 (MQ=255)
                                                                                                      gCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATACCg                                                                                                                                     >  2:1051881/1‑239 (MQ=255)
                                                                                                      gCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACg                                                                                                                                     >  1:395588/1‑239 (MQ=255)
                                                                                                      gCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACg                                                                                                                                     >  1:947413/1‑239 (MQ=255)
                                                                                                         aaCGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACgttg                                                                                                                                  <  1:279611/239‑1 (MQ=255)
                                                                                                                 cATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGAAAGCAAAGGCCAGAAATCCAGCCTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGttctt                                                                                                                          <  1:66415/239‑1 (MQ=255)
                                                                                                                     caagacaagaTCAACCCGCCCTTTGAGAATTGTCCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAATCAAGATAGCATGTTGGCTAACGGCGGCGGAAAAAAAAGGTCCCCTGGGAATGCTATTTTTTCGTTTTCCGAAATTCTCAGGGCTGATTCATCTGTACACCCAGCAAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCt                                                                                                                      <  1:281734/239‑1 (MQ=255)
                                                                                                                     cAAGGCACAACCAACACGCCCGTGGAGAATGGCTCACTGAACCAGCCCGCAAGCAAAGGCAAGAAACCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGAAAATAAAGGTACCTTGGGAAGCCTATTTTGTCGTTTCCCGCAATTCTCATGGCTGATTCATCTGTACACCCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCt                                                                                                                      <  1:70498/239‑1 (MQ=255)
                                                                                                                       aGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTa                                                                                                                    >  1:687477/1‑239 (MQ=255)
                                                                                                                       aGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGAAAGCAAAGGCCAAAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTa                                                                                                                    <  1:221013/239‑1 (MQ=255)
                                                                                                                             gaacAACACCCCTGTTAAGATTTGCTCACTGAACCATTCAGAAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCAACCAAAGATACCATTTTCCCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATTCTATTTTGTCGTTTTCCGCAATGCTCATGGCTGATTCATCTGTACAACCAGCCAGTAATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGaaa                                                                                                              <  2:1157388/236‑1 (MQ=255)
                                                                                                                             gATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGaa                                                                                                               >  1:639553/1‑238 (MQ=255)
                                                                                                                              aTCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAac                                                                                                             >  2:892982/1‑239 (MQ=255)
                                                                                                                                caccacGCCGGTTGAGAATGCCCCACCGACCCAGTCCGAAAGCAAAGGCCAGAAAACCAGACTACCAATACTGACCCAACCAAGATAGCATGTTCGCTGCCAGCGGCGGAAAATAAAGGTCCCTGGGGAATGTTTTTGTGTCGTTTTCCGCAATCCGCATGGCTTATTCTTCTGTACAACCAGCCAGTGATAGATGACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAacac                                                                                                           <  2:361669/236‑1 (MQ=255)
                                                                                                                                 aaCACGCCCGTTGAAAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACAcc                                                                                                          <  1:1063773/239‑1 (MQ=255)
                                                                                                                                   cacGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGAATGTTCGCTGAAAGCGGCGGAAAACAAAGGTACCCTGGGAATGATATTTTGACGTTTTCCGCAATACTCATGGCTGATTCATCTGTACAACAAGCCAGTGATAGATTACAGCACCGATATAACGATGTAGTTCTTCGCTTAACGCAacaa                                                                                                           >  1:671690/1‑235 (MQ=255)
                                                                                                                                      gcccttgaagaATTGCTCACTAACCCGGTCCGCAAACAAAGCCCAAAAATCCAGCCTACCAATACGGACCGAACCAAGATAGAATGTGCGCTGACAGCGCGGGCAAATAAAGGTACCCTGGAAATGCTATTTTGACGTTTTCCGCAATTCCCATGGCTGATTCATCTGTACAACCAGCCAGTAATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCg                                                                                                     <  2:105692/231‑1 (MQ=255)
                                                                                                                                      ccccTTTGAGAATTGCTCACTAAACCAGCCCGAAAGAAAAGGCAAAAAATCAAGCCTATAAATACTGACCAAACCAAGATAGCATGTCCGGTAACAGCGGCGGAAAATAAAGGTCCCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACACCCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCg                                                                                                     <  2:177706/238‑1 (MQ=255)
                                                                                                                                       cccttcgaaaaTTGCTCACTGAACCAGTCCGCAAGCAAAGGCAAGAAAGCCAGACTAACAATACTGACCGAACCAAGATAGCATGTTCGCTGCAAGCGGCGGCAAATAAAGGTCCCTTGGGAATGCTATTTTGTCGTTTGCCGCAATTCTCATGGCTGATTCATCTGTACAACCGGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCgg                                                                                                    <  2:715278/230‑1 (MQ=255)
                                                                                                                                           ttGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGAACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATACATTACCGCCACGTTATACCGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGCAATg                                                                                                >  2:227492/1‑239 (MQ=255)
                                                                                                                                               aaaTTGCTCACTGAACCAGTCCGAAACCAAAGGCCAGAAACCAAGTCTATAAATACTGACCAAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAAAAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATCCTCATGGCTGATTCATCTGAACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTAt                                                                                            <  2:522235/238‑1 (MQ=255)
                                                                                                                                                   tGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCGGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCa                                                                                        <  1:339651/239‑1 (MQ=255)
                                                                                                                                                      tCACTGAACCAGTCCGCAAGCAAGGCCCAGAAATCAAGTCTATAAATCCTGACCGAACCAAGAATCCATGTTCGTTGACAGCGGGGTCAAATAAAGGTACCCTGGGAATCTTTTTTCGTCGTTTTCCACAATCCTCGTGCCAAATTCCTCTGTACACCCAGCAAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACtt                                                                                     <  2:137931/239‑1 (MQ=255)
                                                                                                                                                              cctgcccGAAAACAAAGGCCAGAAATCAAGATTACAAATACCGCCCGAACCAAGATAGCATTTTCGCTGACAGCGGAGAAAAATAAGGGTCCCCTGGAAATGCTTTTTTGTCGTTTTCCGAATTTCTCATGGCTGATTCATCTGTACAGCCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGAc                                                                             <  2:960967/234‑1 (MQ=255)
                                                                                                                                                               cAGTCGCCAAGCAAAGACAAAAATACCGGCTTATCAATCCTGACCGAACAAAGCTACCATGTCCGATGCCGGCGGCGCAAAATAAAGTTCCCCGGGGAATGTTTTTTTGTCTTTTCCCGCATTTCTCATGTTTTTTTCTTCTGTACAACCACCCAGTGATATATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACa                                                                            <  2:1077866/239‑1 (MQ=255)
                                                                                                                                                                      caagcaaAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGAAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCGATTATCATGGCTGATTCATCTGTACAACCAGCCAGCGATACCTTACCGCCACGATATACCGTTGTTGTTGTTCGCTTAACGAAACACCGTGCGGACTGTTATGCCACTTAGCCAGaacaccccg                                                                      >  2:509436/1‑230 (MQ=255)
                                                                                                                                                                                           gACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAAccagc                                                                                                                                                                  >  2:487428/1‑125 (MQ=255)
                                                                                                                                                                                           gACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAAccagc                                                                                                                                                                  <  1:487428/125‑1 (MQ=255)
                                                                                                                                                                                              tATCAATACTGACCGAACCAAGATAGCATGTTCGCGACCAGCGGCGGAAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTGCCGCAATTCTCATGGCTGATTCACCTGTAAACCCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTGCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGa                                             <  1:509436/239‑1 (MQ=255)
                                                                                                                                                                                               ataaaTACTGACGGACCCAAGAAAGTGTGTTCGTTGACAGCGGGGGCAAATAAAGTTTCCTGGGAAATGTTTTTTTGTCGTTTCCCGCAATCCTCATGGTTGATTCATCTTTCCATCCACCCGGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTGCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGGCAGCAGGTCGCGGTGGCGTGCTGAg                                            <  2:893053/236‑1 (MQ=255)
                                                                                                                                                                                                      cTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATaa                                      >  1:695911/1‑238 (MQ=255)
                                                                                                                                                                                                        gACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAac                                    >  1:704133/1‑238 (MQ=255)
                                                                                                                                                                                                        gACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAac                                    >  1:35160/1‑238 (MQ=255)
                                                                                                                                                                                                        gACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGAGCTGAGCGATAAac                                    >  2:865474/1‑238 (MQ=255)
                                                                                                                                                                                                          ccGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCa                                                                                                                                                                                                                                            <  2:59551/37‑1 (MQ=255)
                                                                                                                                                                                                          ccGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCa                                                                                                                                                                                                                                            >  1:59551/1‑37 (MQ=255)
                                                                                                                                                                                                             aaCCAAGATAGCATGTTCCCGGACGGCGGCGGAAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTCCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGa                              <  2:687477/239‑1 (MQ=255)
                                                                                                                                                                                                                   gATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGccc                         >  1:56854/1‑238 (MQ=255)
                                                                                                                                                                                                                        cATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCCGCGCAt                   >  1:1134214/1‑239 (MQ=255)
                                                                                                                                                                                                                         aTGTTCGCTGACAGGGGCGGAAAAAAAAGGTCCCCTGAGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCAGGTAAAAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCCGCGCa                    <  2:138941/237‑1 (MQ=255)
                                                                                                                                                                                                                            ttCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCCGCGCATgggg               >  1:826586/1‑239 (MQ=255)
                                                                                                                                                                                                                               acttccggcggcgCAAAAAAAAGGTACCCGGTGATTCCTTTTTTGTCTTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAAATTCCCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACGCCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCCGCGCATGGGGGtt            <  2:2034/232‑1 (MQ=255)
                                                                                                                                                                                                                                       cggcggCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGCTAAACACCGGATGCCCGCGCATGGGGGTTTGCTTAcc    >  1:404812/1‑239 (MQ=255)
                                                                                                                                                                                                                                         gcggcAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCCGCGCATGGGGGTTTGCTTACCg   >  1:166004/1‑238 (MQ=255)
                                                                                                                                                                                                                                          cggcAAATAAAGTTCCCCGGGGATGGTATTTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAAATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCCGCGCATGGGGGTTTGCTTACCGt  <  1:5840/238‑1 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                           
CTCCCCGCTCCGCCAATGACTACCGAGGCAATAACCGACACAATAAGCGCGATGGCGGGCATCATCTGTTCTGCTTTAAAGCGCGGCAGCACTCCAGGAATAGCCAACGTCAGCATACAAGGCACGATCAACACGCCCGTTGAGAATTGCTCACTGAACCAGTCCGCAAGCAAAGGCCAGAAATCCAGACTATCAATACTGACCGAACCAAGATAGCATGTTCGCTGACAGCGGCGGCAAATAAAGGTACCCTGGGAATGCTATTTTGTCGTTTTCCGCAATTCTCATGGCTGATTCATCTGTACAACCAGCCAGTGATAGATTACCGCCACGATATAACGTTGTTGTTCTTCGCTTAACGAAACACCCTGCGGAATGTTATGCCACTTAGCCAGACAACCGCGACAGCAGGTCGCGGTGGCGTGCTGAGCGATAAACACCGGATGCCCGCGCATGGGGGTTTGCTTACCGT  >  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: