Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,407,017 C→T 12.9% I83I (ATC→ATT yjfP → acyl CoA esterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,407,0170CT12.9% 925.6 / 93.8 303I83I (ATC→ATTyjfPacyl CoA esterase
Reads supporting (aligned to +/- strand):  ref base C (116/148);  new base T (13/26);  total (129/174)
Fisher's exact test for biased strand distribution p-value = 2.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 304 total aligned reads displayed.

CACGGCTTTACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACGGCACTGGGGATTACCGCTCGCCACCCCACGGTGAGAT  >  CP009273/4406868‑4407165
                                                                                                                                                     |                                                                                                                                                    
cACGGCTTTACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGCTAAATCAATTCTGGCAAATc                                                                                                                                                      >  2:663289/1‑150 (MQ=255)
 aCGGCTTTACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCt                                                                                                                                                     >  2:4063067/1‑150 (MQ=255)
  cGGCTTTACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCtt                                                                                                                                                    >  1:1672245/1‑150 (MQ=255)
   ggCTTTACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATTTTg                                                                                                                                                   >  2:501083/1‑150 (MQ=255)
    gCTTTACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGTCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGc                                                                                                                                                  <  2:2174342/150‑1 (MQ=255)
    gCTTTACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATTTTGc                                                                                                                                                  <  1:1448363/150‑1 (MQ=255)
      tttACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTa                                                                                                                                                >  1:12671/1‑150 (MQ=255)
      tttACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTa                                                                                                                                                <  1:4402329/150‑1 (MQ=255)
      tttACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTa                                                                                                                                                <  1:4229599/150‑1 (MQ=255)
      tttACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTa                                                                                                                                                >  1:2395756/1‑150 (MQ=255)
      tttACTTCATCCAGTCTGGTGTATAGCTATTGTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTa                                                                                                                                                <  2:1953949/150‑1 (MQ=255)
        tACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACa                                                                                                                                              <  1:603020/150‑1 (MQ=255)
          cTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACaaa                                                                                                                                            <  2:2908798/150‑1 (MQ=255)
           ttCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAg                                                                                                                                           <  1:3451912/150‑1 (MQ=255)
           ttCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAg                                                                                                                                           <  2:1193294/150‑1 (MQ=255)
             cATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTa                                                                                                                                         >  2:3778611/1‑150 (MQ=255)
              aTCCGGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTAt                                                                                                                                        >  2:3188514/1‑150 (MQ=255)
              aTCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTAt                                                                                                                                        <  1:4011483/150‑1 (MQ=255)
               tCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATg                                                                                                                                       <  1:3712796/150‑1 (MQ=255)
               tCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATg                                                                                                                                       <  1:2723447/150‑1 (MQ=255)
                 cAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATTTTGCTACAAAGTATGCa                                                                                                                                     <  2:1933606/150‑1 (MQ=255)
                       ggTGTATAGCTATTTTGCCGGTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAAtt                                                                                                                               <  2:1356155/150‑1 (MQ=255)
                          gtATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCac                                                                                                                            <  2:535925/150‑1 (MQ=255)
                          gtATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCac                                                                                                                            <  2:3896420/150‑1 (MQ=255)
                            ataGCTATTGTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCacta                                                                                                                          <  2:516830/150‑1 (MQ=255)
                             taGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATTTTGCTACAAAGTATGCAGGAATTCactac                                                                                                                         >  1:3558792/1‑150 (MQ=255)
                               gCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATTTTGCTACAAAGTATGCAGGAATTCACTACtt                                                                                                                       <  1:2545486/150‑1 (MQ=255)
                                cTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACttt                                                                                                                      <  1:3304025/150‑1 (MQ=255)
                                 tATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATTTTGCTACAAAGTATGCAGGAATTCACTACTTTa                                                                                                                     >  1:4367546/1‑150 (MQ=255)
                                   ttttGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACg                                                                                                                   >  2:2989896/1‑150 (MQ=255)
                                     ttGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTg                                                                                                                 >  2:3271761/1‑150 (MQ=255)
                                      tGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGc                                                                                                                >  2:1499107/1‑150 (MQ=255)
                                         cGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGc                                                                                                             >  1:3279168/1‑150 (MQ=255)
                                             gcgcTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATTTTGCTACAAAGTATGCAGGAATTCACTACTTTCCGTGCGGCAATa                                                                                                         <  2:425659/150‑1 (MQ=255)
                                               gcTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGc                                                                                                       <  1:2632861/150‑1 (MQ=255)
                                                cTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGcc                                                                                                      <  2:400757/150‑1 (MQ=255)
                                                cTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGcc                                                                                                      <  2:1432224/150‑1 (MQ=255)
                                                 tGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCg                                                                                                     >  2:72042/1‑150 (MQ=255)
                                                    cgcAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCgaag                                                                                                  >  1:3591022/1‑150 (MQ=255)
                                                    cgcAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCgaag                                                                                                  >  1:305436/1‑150 (MQ=255)
                                                     gcAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATTTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCgaaga                                                                                                 <  1:1956621/150‑1 (MQ=255)
                                                         gCTGGTTTGCGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAAc                                                                                             <  1:86310/150‑1 (MQ=255)
                                                          cTGGTTTGCGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACt                                                                                            <  1:480725/150‑1 (MQ=255)
                                                           tGGTTTGCGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTg                                                                                           <  2:3105515/150‑1 (MQ=255)
                                                           tGGTTTGCGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTg                                                                                           <  2:3704383/150‑1 (MQ=255)
                                                             gTTTGCGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGgc                                                                                         >  1:1149592/1‑150 (MQ=255)
                                                              tttGCGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGgct                                                                                        <  1:3955477/150‑1 (MQ=255)
                                                                  cGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCtt                                                                                    <  1:413656/150‑1 (MQ=255)
                                                                  cGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCtt                                                                                    <  1:1750296/150‑1 (MQ=255)
                                                                  cGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCtt                                                                                    <  2:3962708/150‑1 (MQ=255)
                                                                   gggTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTtg                                                                                   <  1:3621904/150‑1 (MQ=255)
                                                                      tGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTtgatg                                                                                >  1:1611223/1‑150 (MQ=255)
                                                                        aTCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGAc                                                                              <  2:222442/150‑1 (MQ=255)
                                                                         tCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGAcc                                                                             <  2:2429617/150‑1 (MQ=255)
                                                                         tCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAAATTTGCTCCAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGACGAAAACTGGCTGCTTGATGAcc                                                                             >  2:4041130/1‑150 (MQ=255)
                                                                             gCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCt                                                                         <  2:229119/150‑1 (MQ=255)
                                                                              ccGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTg                                                                        <  2:2773622/150‑1 (MQ=255)
                                                                              ccGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTg                                                                        <  2:4497632/150‑1 (MQ=255)
                                                                                       cccGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGt                                                               >  2:941137/1‑150 (MQ=255)
                                                                                       cccGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGt                                                               <  1:1882870/150‑1 (MQ=255)
                                                                                        ccGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTg                                                              <  2:1380108/150‑1 (MQ=255)
                                                                                        ccGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTg                                                              <  1:4298416/150‑1 (MQ=255)
                                                                                         cGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATTTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTgg                                                             <  1:1409796/150‑1 (MQ=255)
                                                                                           aTCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATTTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGgcg                                                           <  2:4451526/150‑1 (MQ=255)
                                                                                           aTCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGgcg                                                           >  2:3550420/1‑150 (MQ=255)
                                                                                             cACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGcgcg                                                         <  1:4155286/150‑1 (MQ=255)
                                                                                                   aGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATg                                                   >  2:1113346/1‑150 (MQ=255)
                                                                                                      cGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGgc                                                >  1:2239758/1‑150 (MQ=255)
                                                                                                       gTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGgcg                                               >  1:2243479/1‑150 (MQ=255)
                                                                                                       gTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGgcg                                               >  2:1054320/1‑150 (MQ=255)
                                                                                                         tttAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGcgcg                                             <  1:3189546/150‑1 (MQ=255)
                                                                                                          ttAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGa                                            >  2:171088/1‑150 (MQ=255)
                                                                                                          ttAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGa                                            >  2:1745548/1‑150 (MQ=255)
                                                                                                              tggtgACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATTTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGAc                                        <  1:3245640/150‑1 (MQ=255)
                                                                                                              tggtgACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATTTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGAc                                        <  2:2038517/150‑1 (MQ=255)
                                                                                                                gtgACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACgg                                      >  1:4596286/1‑150 (MQ=255)
                                                                                                                    cGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACGGCACt                                  >  1:1978546/1‑150 (MQ=255)
                                                                                                                    cGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACGGCACt                                  >  2:420528/1‑150 (MQ=255)
                                                                                                                     gCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACGGCACTg                                 >  1:246126/1‑150 (MQ=255)
                                                                                                                      cAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACGGCACTgg                                >  1:181500/1‑150 (MQ=255)
                                                                                                                      cAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACGGCACTgg                                >  1:2321595/1‑150 (MQ=255)
                                                                                                                       aGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACGGCACTggg                               >  1:4519247/1‑150 (MQ=255)
                                                                                                                        gcgcGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACGGCACTgggg                              <  1:3896790/150‑1 (MQ=255)
                                                                                                                           cggcggTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGTAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACGGCACTGGGGAtt                           <  2:1765572/150‑1 (MQ=255)
                                                                                                                             gcggTTAAATCAATTCTGGCAAATTTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACGGCACTGGGGATTAc                         <  2:1468137/150‑1 (MQ=255)
                                                                                                                             gcggTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACGGCACTGGGGATTAc                         <  2:818811/150‑1 (MQ=255)
                                                                                                                                gTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACGGCACTGGGGATTACCGc                      >  1:790178/1‑150 (MQ=255)
                                                                                                                                gTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACGGCACTGGGGATTACCGc                      <  1:194945/150‑1 (MQ=255)
                                                                                                                                 ttAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACGGCACTGGGGATTACCGCt                     >  1:2891614/1‑150 (MQ=255)
                                                                                                                                 ttAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACGGCACTGGGGATTACCGCt                     >  2:1720458/1‑150 (MQ=255)
                                                                                                                                  tAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACGGCACTGGGGATTACCGCTc                    <  1:3977712/150‑1 (MQ=255)
                                                                                                                                  tAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACGGCACTGGGGATTACCGCTc                    <  2:2495625/150‑1 (MQ=255)
                                                                                                                                   aaaTCAATTCTGGCAAATTTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACGGCACTGGGGATTACCGCTCg                   >  2:2700498/1‑150 (MQ=255)
                                                                                                                                    aaTCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACGGCACTGGGGATTACCGCTCGc                  >  2:584686/1‑150 (MQ=255)
                                                                                                                                        aaTTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACGGCACTGGGGATTACCGCTCGCCAcc              >  2:3378388/1‑150 (MQ=255)
                                                                                                                                             tGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACGGCACTGGGGATTACCGCTCGCCACCCCACg         <  1:1759623/150‑1 (MQ=255)
                                                                                                                                               gCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACGGCACTGGGGATTACCGCTCGCCACCCCACGGt       >  1:2153385/1‑150 (MQ=255)
                                                                                                                                                 aaaTCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACGGCACTGGGGATTACCGCTCGCCACCCCACGGTga     <  1:368176/150‑1 (MQ=255)
                                                                                                                                                   aTCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACGGCACTGGGGATTACCGCTCGCCACCCCACGGTgaga   >  1:2822263/1‑150 (MQ=255)
                                                                                                                                                   aTCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACGGCACTGGGGATTACCGCTCGCCACCCCACGGTgaga   >  1:4075047/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CACGGCTTTACTTCATCCAGTCTGGTGTATAGCTATTTTGCCGTTGCGCTGGCGCAGGCTGGTTTGCGGGTGATCATGCCGGATGCGCCCGATCACGGTAGCCGTTTTAGTGGTGACGCAGCGCGGCGGTTAAATCAATTCTGGCAAATCTTGCTACAAAGTATGCAGGAATTCACTACTTTACGTGCGGCAATAGCCGAAGAAAACTGGCTGCTTGATGACCGTCTGGCAGTCGGTGGCGCGTCGATGGGCGCGATGACGGCACTGGGGATTACCGCTCGCCACCCCACGGTGAGAT  >  CP009273/4406868‑4407165

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: