Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,135,439 T→G 12.0% T311P (ACC→CCC)  pitB ← phosphate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,135,4390TG12.0% 301.6 / 2.6 126T311P (ACC→CCC) pitBphosphate transporter
Reads supporting (aligned to +/- strand):  ref base T (28/75);  new base G (17/3);  total (48/78)
Fisher's exact test for biased strand distribution p-value = 1.77e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 126 total aligned reads displayed.

GCCATGATGATCCACACCGGAGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATC  >  NC_000913/3135204‑3135671
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                         
cccaggagtatCCCCACCGGGGCGTACTCAATGGTGCTTAACAAATCGCTGGGAAGTTTTGTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGTTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTAtt                                                                                                                                                                                                                                        <  2:548122/230‑1 (MQ=255)
        tttCCACCCCTGGGCGGACCAAAGGGTGCTTACCATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCGTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGGTTTAACGCGCGTAATAGCATCAAAGGTATTTGCCGGAt                                                                                                                                                                                                                                <  2:856186/237‑1 (MQ=255)
        gatCCACCCCGGGGTGTTCTTAATGGTGCTTAACCTATCGCGGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGGTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGAt                                                                                                                                                                                                                                <  2:474394/239‑1 (MQ=255)
               cccGGAGCGTACTCAATGGTGCTTAACATTTCGCTGCGAAGTTTTTTCAGGCGGTTTTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTg                                                                                                                                                                                                                         <  1:700158/238‑1 (MQ=255)
                 cGGTGCGTACTCAATGGTGCTTAACATTTCGCTACGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGaa                                                                                                                                                                                                                       <  1:469186/239‑1 (MQ=255)
                   gAGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGGTACATTCCATATTGCCTGGGCGCATCCGTTTTACGCGCCCAAAAACATACAAGGGATTTTCCGGAGGAAAggaaaa                                                                                                                                                                                                                      >  2:68431/1‑233 (MQ=255)
                    aGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGGAACTTTCCATATTGCCTGGCGGCATCCCTTTAAAGCGCGCAATAGCATCAAAAGTATTTTCCGGCTGCCCGTGaaaa                                                                                                                                                                                                                     >  1:1035969/1‑237 (MQ=255)
                    aGCGTACTCAATGCTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACCTTCCCAATTTCCTGGCAGCATCCGTTTTAATGACGCAATTGCCACCAAAGTACTTGCCCGAGGCCAGTTaaaa                                                                                                                                                                                                                     >  2:883289/1‑237 (MQ=255)
                      cGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGGAAACTTTCCAATTGGCCGGGAGGATCGGTTTAAACGCGCGCAAAACATTCAAAGGAAttttcgggtgtaagtgaaaacct                                                                                                                                                                                                                   >  1:585543/1‑219 (MQ=255)
                            cAATGGTTCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGGCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACt                                                                                                                                                                                                            <  1:1249841/239‑1 (MQ=255)
                               tGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTTGCAGTTTCGCTAGCCTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGc                                                                                                                                                                                                         <  1:99659/239‑1 (MQ=255)
                                ggtacctAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACGCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCg                                                                                                                                                                                                        <  1:1195351/233‑1 (MQ=255)
                                  ggCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGt                                                                                                                                                                                                       <  1:966831/237‑1 (MQ=255)
                                   gCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCCGGCAGCATCGTTTTTACCCGCCCCAAAGGATCAAAGGGATTTTCCGGATGGCAGTGAAACCCTGTTACtttgggg                                                                                                                                                                                                      >  1:534693/1‑233 (MQ=255)
                                    cTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCCAAGGGATTTTCCGGATGACAGTGAAACTCTGTTACTTTCGttc                                                                                                                                                                                                     >  1:752449/1‑236 (MQ=255)
                                     ttAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCa                                                                                                                                                                                                   <  1:1102946/239‑1 (MQ=255)
                                       aaCATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCCTATTTCCCGGCAGCATCGTTTTTAAGCCCGCAATAGCATCAAAGGGATTTGCCGGATGTAAGTGGAACCTCGTTAATTCCGTCCCATc                                                                                                                                                                                                 >  1:101199/1‑239 (MQ=255)
                                       aaCATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAAATTTCCACATTTCCCGGCGACCATCGTTTAAGCGCACCCAAAAGACACAAAGGATTTCCCCGGTGACAAGGAAAcatcttttatttcgggcctac                                                                                                                                                                                                 >  2:656882/1‑219 (MQ=255)
                                            atCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTTCCTGGCCGCATACTTTTTAAGACCGCAATACCCTCAAAAGGATTTGCCGGTTGAACGGGAAAATCTGTTACTTTCGTGCCacacgacc                                                                                                                                                                                            >  2:854218/1‑232 (MQ=255)
                                             taagagaccggTTTTTTCAGGAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGGTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGa                                                                                                                                                                                           <  1:685675/229‑1 (MQ=255)
                                              cGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCCTATTTCCTGGCAGCATCGGTTTAACCCGCGCAAAAGCATCCAAGGGATTTTCCGGATGACCGTGAAACTCTGTTACTTGCCTGCCCCCAGCCGAt                                                                                                                                                                                          >  1:813005/1‑239 (MQ=255)
                                              cGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCCTATTGCCTGGGAGCATCGTTTTAAAGCGCGCACAAGCATCAAAGGGATTTGCCGGATGACAGTGAAAACCTGTTACTTTCGGGCCATCAGTccgt                                                                                                                                                                                          >  1:504197/1‑236 (MQ=255)
                                                 tGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACATCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCCTGACCTTCCCTATTTCCTGGCAGCCCACCTTTTACACGCGCCAAAGCATCTATGCTAGTTGCCGGATGAAAAGGAAACACCGCTACTCCCGTACCAACACTCCAtgca                                                                                                                                                                                       >  2:558780/1‑238 (MQ=255)
                                                  gggAAGTTTTTGCAGCGGGGTCTGGTCTTCTTTACGGACGCCTGGCAGTTTCGGTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGAtgctg                                                                                                                                                                                     <  1:364021/237‑1 (MQ=255)
                                                                             ttttttGCTGGCCCCTGGGAGGTTCGCGCGGTTCGGGGGGGTATCAGAGATGCACTGCATAAGGGGGCCAAGCTGGCTGCGTTGAGTCTGCCTTAACGGCTCGTAACTTGCCATATTGCCTGGCAGCATCGGTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCa                                                                                                                                                          <  2:417822/236‑1 (MQ=255)
                                                                               cttTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGTCCGGCCGCATCGTTTTTAAACGCGCCACAGCATCAAAGGGATTTTCCGGATGACAGTTGAAATCTTTTACTTGCGTGCCCTCAGCCGAGGCTGCCGGCAAAGGGGGGTcccaccacataaa                                                                                                                                                         >  2:517902/1‑227 (MQ=255)
                                                                                       aCGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCCAAGGGATTTGCCGGATGACCGTGAAAATCTGTTACTTGCGTGCCATCAGTCGATTCTGCCGGCCATAGAGGTTCCTACCCAATCAACTTCTGCg                                                                                                                                                 >  1:361870/1‑238 (MQ=255)
                                                                                            tGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTTCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTTCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTCCGGCGGt                                                                                                                                            >  1:54550/1‑238 (MQ=255)
                                                                                              gCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCCTATTTCCTGGCAGCATCGTTTTAACGCGCGCAAAAGCATCAAAGGGATTTTCCGGATGACAGGTAAACACTTTTTATTGCGTGCCATCCAACGAGTCTGCCGGCAAAGGGGGGTCCATACCAATCAACCTCTTCGGCGGTTCa                                                                                                                                         >  2:1251427/1‑239 (MQ=255)
                                                                                                 gTTTCGCTTGCTTCGCGGAGGTATGCGTGATTCCCAGCATAATGCGGGGCAGCTGGCTGCGCTGACCCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGGGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGa                                                                                                                                      <  2:504197/239‑1 (MQ=255)
                                                                                                   ttccctggtttCGCGGTGTTTTCAGAGATGCACAGCTTAAGGCGGCGCGGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATg                                                                                                                                    <  2:24537/230‑1 (MQ=255)
                                                                                                   ttCGCTAGCTTCGCGGAGGTTTCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATg                                                                                                                                    <  1:656882/239‑1 (MQ=255)
                                                                                                   ttCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAAAGCGCGCAATAGCATCAAAGGTATTTTCCGGGTGACAGTGAAACTCTGTTATTTTCGTGCCATCAATCGATGCTGCAGGCAATAGGGGTTCCATCGCAAACAACTTTTGCGGcgagtacgggag                                                                                                                                    >  2:1325071/1‑228 (MQ=255)
                                                                                                   ttCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCAAATTTCCTGGCAGCATCGTTTTAACGCGCGCCATAGCATCAAAAGGATTTTCCGGATGACCGTGAAACTCTGTTAACTTCGGGCCATTCTGCGATGCTGCAGGCAATGGAGGGTCCATCGCAATCaaactttccggcgggtccgggg                                                                                                                                     >  1:423271/1‑218 (MQ=255)
                                                                                                     cGCTGGCTTCGGGGAGGTATCAGAGATGCACAGAATAATGCGGCGCAGCTGGCTGCGCTGACTCACGTTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCt                                                                                                                                  <  2:54550/239‑1 (MQ=255)
                                                                                                        tAGCGTCGCGGAGGTAACAAGGATGCCCAGGATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGtt                                                                                                                               <  1:1122062/239‑1 (MQ=255)
                                                                                                         aGCTTCGGGGGGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTg                                                                                                                              <  1:93061/239‑1 (MQ=255)
                                                                                                              cgcgGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGt                                                                                                                         <  1:287419/239‑1 (MQ=255)
                                                                                                               ggtggagtaacaggagtGCACAGCATAATTGGGCGGAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTa                                                                                                                        <  1:1239886/223‑1 (MQ=255)
                                                                                                                   aGGTATCAGGAAGGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGAATCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGa                                                                                                                                      <  1:780760/221‑1 (MQ=255)
                                                                                                                    ggTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTTGCAGCATCGCTTTAACCCCCGCAATAGCCTCCAAAGTATTTGCCGGGTTGCCGTGAAACTCTGTTACTTTGGTGCCATCAGTCGATGCCGCAGGCAATGGGGGTTTCAACGCAATCAACTTCTTCGGCGGTTCAGGATgtctttggtgggatggtt                                                                                                                   >  2:798099/1‑222 (MQ=255)
                                                                                                                     gTATCAGATATGCAACGCATAATGCGGCCCAGCTTGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCAGCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGtt                                                                                                                   <  2:981128/238‑1 (MQ=255)
                                                                                                                                 cccaGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACGTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTaa                                                                                                       <  1:944500/236‑1 (MQ=255)
                                                                                                                                      cATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTTCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCa                                                                                                                                                                                                   >  1:511355/1‑142 (MQ=255)
                                                                                                                                      cATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTTGCCGCATCGTTTTAAAGCGCGCAATAGCATCCAAAGTATTTGCCGGATGACAGGGAAACTCTGTTACTTGCGTGCCACCAGTCGGTGCATGCGGGAAAAGAGGGTCCATCGCAATCAACATCTTCCGCAGTTTACGAGGCCGTTGCGGGTAGGGTTCGGAGTTGGTAAAggaca                                                                                                 >  2:1066973/1‑236 (MQ=255)
                                                                                                                                      cATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGGATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCCTCGCAATCAACTTCTGCGGCCGTTCAGGATGCTGTTGCAGGTAGTGTTCCCAGTTGGGAAAGGCa                                                                                                  >  1:1080076/1‑238 (MQ=255)
                                                                                                                                      cATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCACCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCa                                                                                                                                                                                                   <  2:511355/142‑1 (MQ=255)
                                                                                                                                       aTAATGGGGCCCAGGTGGCTGGTGGTACTAACGCTTAGCGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATc                                                                                                <  2:361870/239‑1 (MQ=255)
                                                                                                                                        tAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATcg                                                                                               <  1:517902/239‑1 (MQ=255)
                                                                                                                                        tAATGCGGCCCAGGTGGCTGCACTGGCTCACGCTTAGCGGCTCGTACCTGTCCATATTGCCTGGCAGCATCGTTTTAACGCTCGAAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACGTCTGCGGCAGTTCAGGAGGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATc                                                                                                <  2:629998/238‑1 (MQ=255)
                                                                                                                                            gtctcccAGCTTGTTGCGCTTCCTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGGATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGGCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGcggg                                                                                           <  2:375253/233‑1 (MQ=255)
                                                                                                                                             cggcgCAGCCGGCTTCGCTGACTCCCCTTTACCGGCTGCTAAATTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGGGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGcgggt                                                                                          <  2:1261051/239‑1 (MQ=255)
                                                                                                                                              ggcgCAGCTGGCTGCGCTTCCTTCCGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGcgggta                                                                                         <  2:1080076/239‑1 (MQ=255)
                                                                                                                                                cgcAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCCTATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAAATCTGTTACTTGCGTGCCATCCTTCGCTGCTGCAGGCCATGGAGGTTCCATCGCAATCAACTTCTGCGGCCGTTCAGGATGCTGTTGCAAGTGTTGTTCGTAGTTGGTAACGGCATCCcgggtacg                                                                                       >  1:285144/1‑239 (MQ=255)
                                                                                                                                                  cAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAgg                                                                                                                          <  1:597624/202‑1 (MQ=255)
                                                                                                                                                   gggttgggtcggttacgtccGCTTAACGGCTCGTAACTTTCCATAGTGCCTGGCAGCGACGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGcgggtacgggt                                                                                    <  2:101199/220‑1 (MQ=255)
                                                                                                                                                    gctggctgCGCTGGCTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGcgggtacgggta                                                                                   <  1:1325071/239‑1 (MQ=255)
                                                                                                                                                      tggctgCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGcgggtacgggt                                                                                    <  1:215740/236‑1 (MQ=255)
                                                                                                                                                              cTGACTCACGCTTTACGGGTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGc                                                                         <  1:690726/239‑1 (MQ=255)
                                                                                                                                                                  cTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATGGCCGGa                                                                     <  1:1293425/239‑1 (MQ=255)
                                                                                                                                                                   tCACGCTTAACGGTTCTTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGAc                                                                    <  1:1138426/239‑1 (MQ=255)
                                                                                                                                                                              ggCTCGGAACTTTCCATATTGCCTGGCAGCATCGTTTGAACGCGCGAAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACGTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATAtt                                                         <  1:883289/239‑1 (MQ=255)
                                                                                                                                                                                 tCGGAACCTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTTCCGGATGTAAGTGAAACTCTGTTACTTGCGGGCCATCAGTCGATGCTGCAGGCAATGGAGGTTTCATCGGCATCAACTTCTGCGGCAGCTCATGATGCCGTTGCAGGGTGTGGTCGCAGGTGGGGACGGCATACCCGGGACGGGGCATTTTATAGTCGGagggatgtgtaatga                                                       >  2:163221/1‑224 (MQ=255)
                                                                                                                                                                                    gTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGGGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTgacga                                                    <  1:220688/239‑1 (MQ=255)
                                                                                                                                                                                      aaCCTTCCATTTTGCCTGGCCGCCTCGGTTTAACGGCGGCAATAGGATCCAAGGGAATTGCCGGATGGCAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGGTGCAGGCAAGGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGacgacg                                                  <  1:678718/239‑1 (MQ=255)
                                                                                                                                                                                            ccATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGa                                                                     <  1:835550/214‑1 (MQ=255)
                                                                                                                                                                                              atatTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCcagc                                          <  1:991944/239‑1 (MQ=255)
                                                                                                                                                                                               atttGCCTGGCAGGATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCcagca                                         <  1:383010/237‑1 (MQ=255)
                                                                                                                                                                                                 tttccgtgcagcaTCGTTTTAACCCGCGCAATTGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCcagc                                          <  1:404848/229‑1 (MQ=255)
                                                                                                                                                                                                   gCCTGGGCGCATGGGTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGGTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAgggg                                     <  2:1090623/239‑1 (MQ=255)
                                                                                                                                                                                                   gCCTGGCAGTATCGTTTTAACGCGCGAAATAGCCACAAAGGGATTGGCCGGGTGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAgggg                                     <  2:50911/239‑1 (MQ=255)
                                                                                                                                                                                                   gCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGCTGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCCGGTAGGGTTCGAAGTTGGCAACGGCATTCCGGGTACGGGTACTTTCATAGCCGGACGCATTATTATTGAGCACGAACCCCGCAgggg                                     >  1:611019/1‑239 (MQ=255)
                                                                                                                                                                                                     cgggCAGCAACGTTTTAACGCGGGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCa                                   <  1:798099/237‑1 (MQ=255)
                                                                                                                                                                                                        gcagcaTCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGt                                                                                                                         <  1:1152758/150‑1 (MQ=255)
                                                                                                                                                                                                         cagcaTCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATc                                >  1:65154/1‑238 (MQ=255)
                                                                                                                                                                                                             caTCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAAc                                                                                                                                                                                                                     >  2:23399/1‑55 (MQ=255)
                                                                                                                                                                                                             caTCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAAc                                                                                                                                                                                                                     <  1:23399/55‑1 (MQ=255)
                                                                                                                                                                                                               tCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGCGTTCGAAGTTGGTAACGGCATCGCGGGGACGGGTAATTTCATAGCCGGACGCCTTCATATTTCCGACGAAGCCAGCAAGGGCAATCCccacc                          >  1:343141/1‑239 (MQ=255)
                                                                                                                                                                                                                ctgtttttAACGCCCGCAATAGCATCCAAGGTTGTGGCCGGTGTGCAGGGAAACTCTGTTCCTGGCGGCCCTTCCGTCGGTGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCccac                           <  2:269785/235‑1 (MQ=255)
                                                                                                                                                                                                                actgTTTAAACCCGCGAAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGGGCCAGCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTGCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCccac                           <  2:285144/236‑1 (MQ=255)
                                                                                                                                                                                                                    ttttAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACATGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGt                       <  1:473656/239‑1 (MQ=255)
                                                                                                                                                                                                                     tttAACGCGCGCAATTGCATCAAAGGTTTTTGGCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTa                      <  2:534693/239‑1 (MQ=255)
                                                                                                                                                                                                                     tttAACGCGAGCAATAGCCTCAAAGGTATTTGCCGGATGACAGTGAAATGTTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTa                      <  1:311124/239‑1 (MQ=255)
                                                                                                                                                                                                                      ttaaccgcggcaataggtgaaaaGGTATTGCCCGGATGACAGGGAAACTCTGTTACTTCCGTGCCATCAGGCGAGGCTTCAGGTAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTAc                     <  2:982603/219‑1 (MQ=255)
                                                                                                                                                                                                                         ccccgcgcAATAGCATCAAAGGTATTTGCCGGAGGACAGTGAAACTTTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGGAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAg                  <  2:1144226/236‑1 (MQ=255)
                                                                                                                                                                                                                         agggcagaaatcggctcgaaGGTTTTTGGCGGATTGCAGGTAAACCGTGTTACTTGGTTGCCATCGGTCGAGTCTGCAGGGAATGGAGGTTACATTGCAATCAATGTCTGCGGCAGTTCGGGTTGCTGTTGGAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGGAGTCATATTGACGACGAAGACAGCAGGGGCAATCCCCACCAGTTCCAg                  <  2:465170/221‑1 (MQ=255)
                                                                                                                                                                                                                         accccgcaaatggcgtaaaaTGTTTTTGCCGGAGAGCAATGAAACTGTGTTCCGTGCGTGCCATCAGTCGAGGCTGCAGGCAATGGAGGTTCCATGTCAATCAACTTCTGCGGCAGTTCAGTATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCCGGTACGGGTAATTTCATTGCCTGACGTATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCAGCAGGACCAg                  <  2:641761/222‑1 (MQ=255)
                                                                                                                                                                                                                         aCGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCGGCAGGGGCAATCCCCACCAGTACCa                   >  1:1144226/1‑238 (MQ=255)
                                                                                                                                                                                                                              cgcAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCCTACTGACGACGAAGCCAGCAGGGGAAATCCCCACCAGTACCACCAtt              >  1:202767/1‑238 (MQ=255)
                                                                                                                                                                                                                              cgcAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGGGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCCCCAGTCCCCGCCtt              >  1:52309/1‑238 (MQ=255)
                                                                                                                                                                                                                                        tCAAAGGTTTTTTCTGGATGACATTGAAACTCTGTTATCTGCGTGCCATCAGGGGATGCTGAGGGAAAATGAGGTTTCATCGGAATCAACGTACGCGGCAGTTCAGGATGTTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGAt   <  2:370145/239‑1 (MQ=255)
                                                                                                                                                                                                                                        tCAAAGGTATTTGGCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGAt   <  1:1229772/239‑1 (MQ=255)
                                                                                                                                                                                                                                        acAAAGTTGTGTGCCGGATGAAAGTGAATCTCTTTTACTTGCGTGCCATCAGTCGATGCTACAGGCAATGGAGTTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGAt   <  1:558780/238‑1 (MQ=255)
                                                                                                                                                                                                                                         cAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTTACGACGAAGCCAGACGGGGAGATCCCCACCAGGACCAGCATTCCCAggacact   >  1:693344/1‑233 (MQ=255)
                                                                                                                                                                                                                                           aaGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCCATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGCCGCATTCATATTGCCGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCCGCATTACCAGGCCGATc  >  1:502113/1‑237 (MQ=255)
                                                                                                                                                                                                                                           aaGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGa                                                                     <  2:457005/170‑1 (MQ=255)
                                                                                                                                                                                                                                           aaGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATc  >  1:765069/1‑237 (MQ=255)
                                                                                                                                                                                                                                           aaGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGAAATTTCATAGCCGGACGCATTCATATTGACGAAGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCCGGCCGATc  >  1:888180/1‑237 (MQ=255)
                                                                                                                                                                                                                                           aaGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGGACGGGTAATTTCATAGCCGGACGCATTCATATTTACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCCGCATTACCAGGCCGATc  >  1:838691/1‑237 (MQ=255)
                                                                                                                                                                                                                                           aaGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGGGTTCGAAGTTGGTAACGGCACCGCGGGTACGGGTAATTACATAGCCGGACGCATTAATATTGACGCCGAAGACAGCAAGGTCAAACCCCCCCCGAACCGGCATTAcacgccagacc  >  2:965445/1‑227 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                         
GCCATGATGATCCACACCGGAGCGTACTCAATGGTGCTTAACATATCGCTGCGAAGTTTTTTCAGCAGGTTCTGGTCTTCTTTACTGACGCCTGGCAGTTTCGCTAGCTTCGCGGAGGTATCAGAGATGCACAGCATAATGCGGCGCAGCTGGCTGCGCTGACTCACGCTTAACGGCTCGTAACTTTCCATATTGCCTGGCAGCATCGTTTTAACGCGCGCAATAGCATCAAAGGTATTTGCCGGATGACAGTGAAACTCTGTTACTTGCGTGCCATCAGTCGATGCTGCAGGCAATGGAGGTTCCATCGCAATCAACTTCTGCGGCAGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATC  >  NC_000913/3135204‑3135671

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