Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,260,920 G→T 5.8% R112L (CGT→CTT)  yjbR → DUF419 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,260,9200GT 5.8% 1127.8 / 27.4 299R112L (CGT→CTT) yjbRDUF419 family protein
Reads supporting (aligned to +/- strand):  ref base G (139/142);  new base T (13/5);  total (152/147)
Fisher's exact test for biased strand distribution p-value = 8.71e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 299 total aligned reads displayed.

CAGCACAGCGATGTGCGTCCAAGCCGCCATCTGAATAAAGCGCACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGAGGATCTCGCCGCCAC  >  CP009273/4260772‑4261068
                                                                                                                                                    |                                                                                                                                                    
 agcaCAGCGATGTGCGTCCAAGCCGCCATCTGAATAAAGCGCACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTa                                                                                                                                                    <  2:3076908/150‑1 (MQ=255)
 agcaCAGCGATGTGCGTCCAAGCCGCCATCTGAATAAAGCGCACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTa                                                                                                                                                    >  1:1287256/1‑150 (MQ=255)
 agcaCAGCGATGTGCGTCCAAGCCGCCATCTGAATAAAGCGCACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTa                                                                                                                                                    >  1:4486478/1‑150 (MQ=255)
 agcaCAGCGATGTGCGTCCAAGCCGCCATCTGAATAAAGCGCACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTa                                                                                                                                                    <  2:794973/150‑1 (MQ=255)
 agcaCAGCGATGTGCGTCCAAGCCGCCATCTGAATAAAGCGCACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTa                                                                                                                                                    <  1:3493042/150‑1 (MQ=255)
 agcaCAGCGATGTGCGTCCAAGCCGCCATCTGAATAAAGCGCACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTa                                                                                                                                                    <  1:3296480/150‑1 (MQ=255)
   cacaGCGATGTGCGTCCAAGCCGCCATCTGAATAAAGCGCACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTaaa                                                                                                                                                  >  2:4280141/1‑150 (MQ=255)
    acaGCGATGTGCGTCCAAGCCGCCATCTGAATAAAGCGCACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAAt                                                                                                                                                 <  2:3853914/150‑1 (MQ=255)
     caGCGATGTGCGTCCAAGCCGCCATCTGAATAAAGCGCACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAAtt                                                                                                                                                <  2:2852599/150‑1 (MQ=255)
      aGCGATGTGCGTCCAAGCCGTCATCTGAATAAAGCGCACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTg                                                                                                                                               >  2:2895650/1‑150 (MQ=255)
      aGCGATGTGCGTCCAAGCCGCCATCTGAATAAAGCGCACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTg                                                                                                                                               >  1:161600/1‑150 (MQ=255)
                gTCCAAGCCGCCATCTGAATAAAGCGCACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAAc                                                                                                                                     <  2:3847136/150‑1 (MQ=255)
                 tCCAAGCCGCCATCTGAATAAAGCGCACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAAct                                                                                                                                    <  1:2970831/150‑1 (MQ=255)
                  acAAGCCTCCATCTGAATAAAGCGCACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAActc                                                                                                                                   >  1:3173706/2‑150 (MQ=255)
                   cAAGCCGCCATCTGAATAAAGCGCACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAActct                                                                                                                                  <  2:4132759/150‑1 (MQ=255)
                     aGCCGCCATCTGAATAAAGCGCACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGa                                                                                                                                <  2:4381880/150‑1 (MQ=255)
                     aGCCGCCATCTGAATAAAGCGCACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCCGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGa                                                                                                                                >  2:70000/1‑150 (MQ=255)
                       ccgccATCTGAATAAAGCGCACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGaaa                                                                                                                              <  1:597637/150‑1 (MQ=255)
                        cgccATCTGAATAAAGCGCACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAg                                                                                                                             >  2:3361748/1‑150 (MQ=255)
                           cATCTGAATAAAGCGCACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGaa                                                                                                                          >  2:3821595/1‑150 (MQ=255)
                           cATCTGAATAAAGCGCACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGaa                                                                                                                          >  2:1458876/1‑150 (MQ=255)
                            aTCTGAATAAAGCGCACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGaaa                                                                                                                         >  2:2329179/1‑150 (MQ=255)
                            aTCTGAATAAAGCGCACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGaaa                                                                                                                         >  2:1149660/1‑150 (MQ=255)
                            aTCTGAATAAAGCGCACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGaaa                                                                                                                         >  1:1000590/1‑150 (MQ=255)
                            aTCTGAATAAAGCGCACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGaaa                                                                                                                         <  1:376123/150‑1 (MQ=255)
                                 aaTAAAGCGCACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGcc                                                                                                                    <  1:274139/150‑1 (MQ=255)
                                   tAAAGCGCACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCCGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGc                                                                                                                  >  2:1630013/1‑150 (MQ=255)
                                   tAAAGCGCACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAATCGTACATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGc                                                                                                                  >  2:659463/1‑150 (MQ=255)
                                   tAAAGCGCACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGc                                                                                                                  >  1:3993358/1‑150 (MQ=255)
                                        cgcACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGa                                                                                                             <  2:2966566/150‑1 (MQ=255)
                                        cgcACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGa                                                                                                             <  2:2727560/150‑1 (MQ=255)
                                          cACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGaaa                                                                                                           <  1:2641186/150‑1 (MQ=255)
                                            cTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGc                                                                                                         <  1:3071396/150‑1 (MQ=255)
                                            cTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGc                                                                                                         <  1:3071338/150‑1 (MQ=255)
                                             tGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCg                                                                                                        >  2:1729901/1‑150 (MQ=255)
                                             tGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCg                                                                                                        >  1:1835152/1‑150 (MQ=255)
                                             tGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCg                                                                                                        <  1:810125/150‑1 (MQ=255)
                                             tGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCg                                                                                                        <  1:809989/150‑1 (MQ=255)
                                              ggAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACTTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCgg                                                                                                       <  2:3981798/150‑1 (MQ=255)
                                               gAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGc                                                                                                      <  1:255682/150‑1 (MQ=255)
                                                  cACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCtt                                                                                                   >  1:2070510/1‑150 (MQ=255)
                                                   aCCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTa                                                                                                  >  2:1964791/1‑150 (MQ=255)
                                                     cGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTTTTATCTGGTGGATGCGTCTTATCAGCTGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAAc                                                                                                <  1:3793223/150‑1 (MQ=255)
                                                           tctcGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTAc                                                                                          <  1:131217/150‑1 (MQ=255)
                                                            ctcGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACa                                                                                         <  1:4334734/150‑1 (MQ=255)
                                                            ctcGACGGTTCGCGGCCAGATTCGCAACTCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACa                                                                                         <  2:3620989/150‑1 (MQ=255)
                                                              cGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGc                                                                                       >  2:572704/1‑150 (MQ=255)
                                                               gACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCa                                                                                      >  2:4364479/1‑150 (MQ=255)
                                                                 cGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATc                                                                                    <  1:2229624/150‑1 (MQ=255)
                                                                  ggTTCGCTGCCAGATTCGCAAATCTTTTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCg                                                                                   <  1:474128/150‑1 (MQ=255)
                                                                  ggTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCg                                                                                   >  1:2000413/1‑150 (MQ=255)
                                                                      cGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCtt                                                                               <  1:2027057/150‑1 (MQ=255)
                                                                       gCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTa                                                                              <  2:2490101/150‑1 (MQ=255)
                                                                        cTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTaa                                                                             >  1:228577/1‑150 (MQ=255)
                                                                        cTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTaa                                                                             >  1:3304644/1‑150 (MQ=255)
                                                                         tGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAg                                                                            >  1:439742/1‑150 (MQ=255)
                                                                         tGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAg                                                                            <  2:3135702/150‑1 (MQ=255)
                                                                          gCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACTTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAgg                                                                           >  2:1401612/1‑150 (MQ=255)
                                                                          gCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACTTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAgg                                                                           <  2:492890/150‑1 (MQ=255)
                                                                            cAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGaa                                                                         <  2:4383770/150‑1 (MQ=255)
                                                                             aGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAg                                                                        >  2:390095/1‑150 (MQ=255)
                                                                             aGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAg                                                                        >  2:3400813/1‑150 (MQ=255)
                                                                              gATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGc                                                                       >  1:2899461/1‑150 (MQ=255)
                                                                                 tCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGATAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTg                                                                    <  2:1087988/150‑1 (MQ=255)
                                                                                   gCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGcc                                                                  <  2:1615128/150‑1 (MQ=255)
                                                                                       aTCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCgtgt                                                              >  1:908121/1‑150 (MQ=255)
                                                                                          tattatCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACTTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGtgtgtg                                                           <  2:4390190/150‑1 (MQ=255)
                                                                                          tattatCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGtgtgtg                                                           <  2:896056/150‑1 (MQ=255)
                                                                                          tattatCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGtgtgtg                                                           <  1:4168467/150‑1 (MQ=255)
                                                                                           attatCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGa                                                          >  2:2686523/1‑150 (MQ=255)
                                                                                           attatCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGa                                                          >  2:3743531/1‑150 (MQ=255)
                                                                                               tCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACTTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCt                                                      >  2:4530798/1‑150 (MQ=255)
                                                                                                cTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCtt                                                     >  1:2816502/1‑150 (MQ=255)
                                                                                                  ggtggATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCg                                                   <  2:766001/150‑1 (MQ=255)
                                                                                                   gtggATGCGTCTTATCAGCAGGCGGTGAATTGACTGCCGGAAGAAAAACGTAAAGTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGc                                                  <  2:3149309/150‑1 (MQ=255)
                                                                                                     ggATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCAc                                                >  1:3869637/1‑150 (MQ=255)
                                                                                                           gTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCcgcg                                          >  2:3503200/1‑150 (MQ=255)
                                                                                                              ttATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACTTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACg                                       >  1:1914880/1‑150 (MQ=255)
                                                                                                               tATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACgg                                      <  2:298061/150‑1 (MQ=255)
                                                                                                                aTCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGt                                     <  1:1962902/150‑1 (MQ=255)
                                                                                                                aTCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGt                                     <  1:3972703/150‑1 (MQ=255)
                                                                                                                aTCAGCAGGCGGTGAATTTACTGCAGGAAGAAAAACGTAAAATGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGt                                     >  1:1080161/1‑150 (MQ=255)
                                                                                                                 tCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGtt                                    <  2:960742/150‑1 (MQ=255)
                                                                                                                  cagcagGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGttt                                   >  1:4115373/1‑150 (MQ=255)
                                                                                                                      agGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTgg                               >  2:3635928/1‑150 (MQ=255)
                                                                                                                             gAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCg                        <  2:3305230/150‑1 (MQ=255)
                                                                                                                              aaTTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCgg                       >  1:2625240/1‑150 (MQ=255)
                                                                                                                                 ttACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGgaga                    <  1:1454073/150‑1 (MQ=255)
                                                                                                                                  tACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGAc                   >  1:976532/1‑150 (MQ=255)
                                                                                                                                  tACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGAc                   >  1:498935/1‑150 (MQ=255)
                                                                                                                                   aCTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACg                  >  1:3631741/1‑150 (MQ=255)
                                                                                                                                    cTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTATAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGa                 >  1:2884098/1‑150 (MQ=255)
                                                                                                                                     tGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGAg                >  2:1796851/1‑150 (MQ=255)
                                                                                                                                       ccGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGAGGa              >  1:4171502/1‑150 (MQ=255)
                                                                                                                                        cGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGAGGAt             >  1:501390/1‑150 (MQ=255)
                                                                                                                                           aagaaAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGAGGAtctc          >  2:1236538/1‑150 (MQ=255)
                                                                                                                                             gaaAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGAGGATCTcgc        >  2:1119763/1‑150 (MQ=255)
                                                                                                                                              aaaaaCGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGAGGATCTcgcc       <  2:4091170/150‑1 (MQ=255)
                                                                                                                                                 aaCTTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGAGGATCTCgccgcc    >  1:4083774/1‑150 (MQ=255)
                                                                                                                                                   cGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGAGGATCTCGCCGccac  >  1:1162409/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CAGCACAGCGATGTGCGTCCAAGCCGCCATCTGAATAAAGCGCACTGGAGCACCGTGTATCTCGACGGTTCGCTGCCAGATTCGCAAATCTATTATCTGGTGGATGCGTCTTATCAGCAGGCGGTGAATTTACTGCCGGAAGAAAAACGTAAATTGCTGGTGCAACTCTGAAAGGAAAAGGCCGCTCAGAAAGCGGCCTTAACGATTACAGCATCGGCTTAAGGAAGCGTGCCGTGTGTGATGCTTCGCACTCCGCGACGGTTTCTGGCGTACCGGAGACGAGGATCTCGCCGCCAC  >  CP009273/4260772‑4261068

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