Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,750,424 G→T 15.4% E303* (GAA→TAA)  intA → CP4‑57 prophage; integrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,750,4240GT15.4% 901.0 / 93.2 277E303* (GAA→TAA) intACP4‑57 prophage; integrase
Reads supporting (aligned to +/- strand):  ref base G (121/113);  new base T (27/16);  total (148/129)
Fisher's exact test for biased strand distribution p-value = 1.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 286 total aligned reads displayed.

TGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCTTGCGTTCTATCGCCAGTACGGCACTCAATGAGGAAGGATTT  >  CP009273/2750276‑2750573
                                                                                                                                                    |                                                                                                                                                     
     cTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTa                                                                                                                                                 <  1:4303775/150‑1 (MQ=255)
       cGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTTATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATc                                                                                                                                               >  1:3476908/1‑150 (MQ=255)
       cGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATc                                                                                                                                               <  1:1668107/150‑1 (MQ=255)
        gATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGATAATTTATCt                                                                                                                                              >  2:705049/1‑150 (MQ=255)
        gATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCt                                                                                                                                              <  2:4304540/150‑1 (MQ=255)
        gATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCt                                                                                                                                              <  1:3613583/150‑1 (MQ=255)
         aTGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGATAATTTATCtt                                                                                                                                             >  1:1536525/1‑150 (MQ=255)
          tGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCttt                                                                                                                                            >  2:3195876/1‑150 (MQ=255)
          tGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCttt                                                                                                                                            >  2:622665/1‑150 (MQ=255)
             gATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAATATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGATATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCt                                                                                                                                         >  1:491235/1‑150 (MQ=255)
             gATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCt                                                                                                                                         <  1:3377441/150‑1 (MQ=255)
              aTGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTa                                                                                                                                        <  1:4359892/150‑1 (MQ=255)
              aTGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTa                                                                                                                                        <  1:4586477/150‑1 (MQ=255)
               tGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAg                                                                                                                                       <  2:217860/150‑1 (MQ=255)
                  gaTCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCg                                                                                                                                    >  2:2345893/1‑150 (MQ=255)
                  gaTCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCg                                                                                                                                    >  1:2304565/1‑150 (MQ=255)
                  gaTCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCg                                                                                                                                    <  2:2971407/150‑1 (MQ=255)
                     cGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTAt                                                                                                                                 >  2:1496945/1‑150 (MQ=255)
                     cGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTAt                                                                                                                                 >  2:3822613/1‑150 (MQ=255)
                      gATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGATAATTTATCTTTCCTAGCCGTATc                                                                                                                                >  1:2382767/1‑150 (MQ=255)
                       aTTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCa                                                                                                                               >  2:4143434/1‑150 (MQ=255)
                       aTTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCa                                                                                                                               <  2:4012691/150‑1 (MQ=255)
                       aTTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGAACATTCGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCa                                                                                                                               <  2:195444/150‑1 (MQ=255)
                         ttCAATGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAg                                                                                                                             <  2:2598472/150‑1 (MQ=255)
                            aaTGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGccc                                                                                                                          >  1:2851445/1‑150 (MQ=255)
                            aaTGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGccc                                                                                                                          <  1:1872833/150‑1 (MQ=255)
                              tGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCccaa                                                                                                                        >  2:2513068/1‑150 (MQ=255)
                                   gCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCccaaccaac                                                                                                                   >  2:1836475/1‑150 (MQ=255)
                                     gAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGATAATTTATCTTTCCTAGCCGTATCAAGCccaaccaacca                                                                                                                 >  1:2061858/1‑150 (MQ=255)
                                     gAATGGAAAAGTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCccaaccaacca                                                                                                                 <  1:2468819/150‑1 (MQ=255)
                                      aaTGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCccaaccaaccaa                                                                                                                >  2:1078984/1‑150 (MQ=255)
                                       aTGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAAt                                                                                                               <  2:1771341/150‑1 (MQ=255)
                                        tGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATg                                                                                                              <  2:2379771/150‑1 (MQ=255)
                                          gAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGaa                                                                                                            <  2:463954/150‑1 (MQ=255)
                                            aaaTTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATa                                                                                                          <  1:543683/150‑1 (MQ=255)
                                             aaTTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAg                                                                                                         >  2:176400/1‑150 (MQ=255)
                                             aaTTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAg                                                                                                         >  2:176496/1‑150 (MQ=255)
                                             aaTTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAg                                                                                                         >  2:2735079/1‑150 (MQ=255)
                                                   tGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGATAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAAc                                                                                                   <  2:592542/150‑1 (MQ=255)
                                                     cAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGATAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAg                                                                                                 >  2:4436987/1‑150 (MQ=255)
                                                       gCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTg                                                                                               <  2:1899890/150‑1 (MQ=255)
                                                       gCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTg                                                                                               >  1:2969003/1‑150 (MQ=255)
                                                           gAATTAAGATGAACAGGGACCATACGGTTCCACTATCTGATGGGGCTCTTTCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGATAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATg                                                                                           <  2:1788470/150‑1 (MQ=255)
                                                            aaTGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATgc                                                                                          <  2:3424241/150‑1 (MQ=255)
                                                                aaGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATgcagca                                                                                      >  1:1311890/1‑150 (MQ=255)
                                                                aaGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATgcagca                                                                                      <  2:820589/150‑1 (MQ=255)
                                                                  gATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACt                                                                                    >  1:707610/1‑150 (MQ=255)
                                                                  gATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCCGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACt                                                                                    >  2:3027715/1‑150 (MQ=255)
                                                                     gAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCaa                                                                                 <  1:228381/150‑1 (MQ=255)
                                                                        ccGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCg                                                                              <  2:3316920/150‑1 (MQ=255)
                                                                         cGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGt                                                                             >  1:4050275/1‑150 (MQ=255)
                                                                            gACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCt                                                                          <  2:1343974/150‑1 (MQ=255)
                                                                             aCCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTg                                                                         <  2:1039748/150‑1 (MQ=255)
                                                                              ccATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTgg                                                                        >  1:1867129/1‑150 (MQ=255)
                                                                                aTACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCt                                                                      >  1:724281/1‑150 (MQ=255)
                                                                                 tACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGATAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCtt                                                                     <  1:1961896/150‑1 (MQ=255)
                                                                                  aCGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTa                                                                    <  1:1079889/150‑1 (MQ=255)
                                                                                    ggTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTagg                                                                  <  1:2016245/150‑1 (MQ=255)
                                                                                    ggTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTagg                                                                  >  1:2537342/1‑150 (MQ=255)
                                                                                    ggTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTagg                                                                  <  1:1439766/150‑1 (MQ=255)
                                                                                    ggTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTagg                                                                  >  2:2601385/1‑150 (MQ=255)
                                                                                       tCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTaggagg                                                               >  1:179252/1‑150 (MQ=255)
                                                                                       tCCACTATCTGAATGTGATATTGCTATTCTGGAAATGATGAATCCTCTCATTTGTGTCCGAGAATTTATCTTTCCTAGACGTATCAAGCCCAACCAACCAATGAATATCAAAACAGTGAATTCAGCACTCAATCTTGCTGGCTTaggagt                                                               >  2:3085364/1‑149 (MQ=255)
                                                                                         cACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAAGAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGgtg                                                             <  2:2834003/150‑1 (MQ=255)
                                                                                           cTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTa                                                           <  2:3864505/150‑1 (MQ=255)
                                                                                             aTCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACt                                                         >  2:64292/1‑150 (MQ=255)
                                                                                               cTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTg                                                       >  1:2138793/1‑150 (MQ=255)
                                                                                               cTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTg                                                       >  1:158418/1‑150 (MQ=255)
                                                                                                tGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGATAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGt                                                      >  2:3991611/1‑150 (MQ=255)
                                                                                                 gATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGtt                                                     <  2:3677191/150‑1 (MQ=255)
                                                                                                 gATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGtt                                                     >  1:1816272/1‑150 (MQ=255)
                                                                                                   tGGGGCTCTTGCCATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTc                                                   >  1:4347718/1‑150 (MQ=255)
                                                                                                     gggCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGATAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTcac                                                 <  1:635245/150‑1 (MQ=255)
                                                                                                        ctctTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGATAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACg                                              >  1:1674432/1‑150 (MQ=255)
                                                                                                         tctTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACgg                                             >  2:1684095/1‑150 (MQ=255)
                                                                                                          ctTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGc                                            >  1:3449008/1‑150 (MQ=255)
                                                                                                           ttGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCt                                           <  1:3656074/150‑1 (MQ=255)
                                                                                                            tGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCtt                                          <  2:119688/150‑1 (MQ=255)
                                                                                                            tGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCtt                                          <  2:1273991/150‑1 (MQ=255)
                                                                                                             gCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCTTg                                         >  1:3815519/1‑150 (MQ=255)
                                                                                                              cTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCTTGc                                        <  2:483815/150‑1 (MQ=255)
                                                                                                                  tCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCTTGCGTTc                                    <  1:1007641/150‑1 (MQ=255)
                                                                                                                   cTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCTTGCGTTCt                                   >  2:1571014/1‑150 (MQ=255)
                                                                                                                   cTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCTTGCGTTCt                                   >  2:4204975/1‑150 (MQ=255)
                                                                                                                      gAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCTTGCGTTCTATc                                <  1:2620453/150‑1 (MQ=255)
                                                                                                                            atgaAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCTTGCGTTCTATCGCCAGt                          >  2:2847932/1‑150 (MQ=255)
                                                                                                                              gaAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCTTGCGTTCTATCGCCAGTAc                        <  2:3758895/150‑1 (MQ=255)
                                                                                                                                   ctctCAGTGGTGGCCGATAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCTTGCGTTCTATCGCCAGTACGGCAc                   >  1:3003078/1‑150 (MQ=255)
                                                                                                                                   ctctCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCTTGCGTTCTATCGCCAGTACGGCAc                   >  2:3656671/1‑150 (MQ=255)
                                                                                                                                    tctcAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCTTGCGTTCTATCGCCAGTACGGCACt                  >  2:4347609/1‑150 (MQ=255)
                                                                                                                                     ctcAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCTTGCGTTCTATCGCCAGTACGGCACTc                 >  2:1131624/1‑150 (MQ=255)
                                                                                                                                     ctcAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCTTGCGTTCTATCGCCAGTACGGCACTc                 >  1:897547/1‑150 (MQ=255)
                                                                                                                                     ctcAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCTTGCGTTCTATCGCCAGTACGGCACTc                 >  2:3173390/1‑150 (MQ=255)
                                                                                                                                      tcAGTGGTGGCCGATAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCTTGCGTTCTATCGCCAGTACGGCACTCa                <  1:116781/150‑1 (MQ=255)
                                                                                                                                        aGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCTTGCGTTCTATCGCCAGTACGGCACTCAAt              <  1:2919901/150‑1 (MQ=255)
                                                                                                                                        aGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCTTGCGTTCTATCGCCAGTACGGCACTCAAt              >  1:2677977/1‑150 (MQ=255)
                                                                                                                                          tggtggCCGATAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCTTGCGTTCTATCGCCAGTACGGCACTCAATGa            >  1:300367/1‑150 (MQ=255)
                                                                                                                                                   agaATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCTTGCGTTCTATCGCCAGTACGGCACTCAATGAGGAAGGAtt   >  2:1474428/1‑150 (MQ=255)
                                                                                                                                                   agaATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCTTGCGTTCTATCGCCAGTACGGCACTCAATGAGGAAGGAtt   >  2:88590/1‑150 (MQ=255)
                                                                                                                                                    gaATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCTTGCGTTCTATCGCCAGTACGGCACTCAATGAGGAAGGAttt  >  1:223114/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGAGGCTCGATGGGATGAGATCGATTTCAATGCTAGCGAATGGAAAATTCCTGCAGCTCGAATGAAGATGAACCGGGACCATACGGTTCCACTATCTGATGGGGCTCTTGCTATTCTGGAAATGATGAAGCCTCTCAGTGGTGGCCGAGAATTTATCTTTCCTAGCCGTATCAAGCCCAACCAACCAATGAATAGCCAAACAGTGAATGCAGCACTCAAGCGTGCTGGCTTAGGAGGTGTACTTGTTTCACACGGCTTGCGTTCTATCGCCAGTACGGCACTCAATGAGGAAGGATTT  >  CP009273/2750276‑2750573

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: