Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,948,349 C→A 6.1% L504L (CTC→CTA ilvD → dihydroxyacid dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,948,3490CA 6.1% 1481.6 / 42.0 375L504L (CTC→CTAilvDdihydroxyacid dehydratase
Reads supporting (aligned to +/- strand):  ref base C (171/180);  new base A (11/13);  total (182/193)
Fisher's exact test for biased strand distribution p-value = 8.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 379 total aligned reads displayed.

TACGAAAGCCAGGACGATGCGGTAGAAGCGATTCTCGGCGGTAAAGTTGTCGCCGGAGATGTGGTAGTAATTCGCTATGAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACCTCTGGTCTTTCCATCGGCCACGTCTCACCGGAAGCGGCAAGCGGCGGCAGCATTGGCCTGATTGAAGATGGTGACCTGATCGCTATCGACATCCCGAACCG  >  CP009273/3948200‑3948498
                                                                                                                                                     |                                                                                                                                                     
  cGAAAGCCAGGACGATGCGGTAGAAGCGATTCTCGGCGGTAAAGTTGTCGCCGGAGATGTGGTAGTAATTCGCTATGAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCgg                                                                                                                                                     <  2:175515/150‑1 (MQ=255)
   gAAAGCCAGGACGATGCGGTAGAAGCGATTCTCGGCGGTAAAGTTGTCGCCGGAGATGTGGTAGTAATTCGCTATGAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGc                                                                                                                                                    >  1:4894131/1‑150 (MQ=255)
    aaaGCCAGGACGATGCGGTAGAAGCGATTCTCGGCGGTAAAGTTGTCGCCGGAGATGTGGTAGTAATTCGCTATGAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCa                                                                                                                                                   >  2:4621390/1‑150 (MQ=255)
    aaaGCCAGGACGATGCGGTAGAAGCGATTCTCGGCGGTAAAGTTGTCGCCGGAGATGTGGTAGTAATTCGCTATGAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCa                                                                                                                                                   <  1:4495008/150‑1 (MQ=255)
     aaGCCAGGACGATGCGGTAGAAGCGATTCTCGGCGGTAAAGTTGTCGCCGGAGATGTGGTAGTAATTCGCTATGAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCaa                                                                                                                                                  <  2:4972253/150‑1 (MQ=255)
     aaGCCAGGACGATGCGGTAGAAGCGATTCTCGGCGGTAAAGTTGTCGCCGGAGATGTGGTAGTAATTCGCTATGAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTAGGCaa                                                                                                                                                  >  2:530539/1‑150 (MQ=255)
       gCCAGGACGATGCGGTAGAAGCGATTCTCGGCGGTAAAGTTGTCGCCGGAGATGTGGTAGTAATTCGCTATGAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAg                                                                                                                                                >  1:91990/1‑150 (MQ=255)
        ccAGGACGATGCGGTAGAAGCGATTCTCGGCGGTAAAGTTGTCGCCGGAGATGTGGTAGTAATTCGCTATGAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGc                                                                                                                                               >  1:4279911/1‑150 (MQ=255)
          aGGACGATGCGGTAGAAGCGATTCTCGGCGGTAAAGTTGTCGCCGGAGATGTGGTAGTAATTCGCTATGAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCt                                                                                                                                             <  1:1920186/150‑1 (MQ=255)
              cGATGCGGTAGAAGCGATTCTCGGCGGTAAAGTTGTCGCCGGAGATGTGGTAGTAATTCGCTATGAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTAGGCAAAGCCTGTgc                                                                                                                                         <  2:2320895/150‑1 (MQ=255)
               gATGCGGTAGAAGCGATTCTCGGCGGTAAAGTTTTCGCCGGAGATGTGGTAGTAATTCGCTATGAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTgcg                                                                                                                                        >  2:864287/1‑150 (MQ=255)
               gATGCGGTAGAAGCGATTCTCGGCGGTAAAGTTGTCGCCGGAGATGTGGTAGTAATTCGCTATGAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTgcg                                                                                                                                        >  2:2905517/1‑150 (MQ=255)
                aTGCGGTAGAAGCGATTCTCGGCGGTAAAGTTGTCGCCGGAGATGTGGTAGTAATTCGCTATGAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTgcgc                                                                                                                                       >  1:3590490/1‑150 (MQ=255)
                 tGCGGTAGAAGCGATTCTCGGCGGTAAAGTTGTCGCCGGAGATGTGGTAGTAATTCGCTATGAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCt                                                                                                                                      <  1:1048428/150‑1 (MQ=255)
                       aGAAGCGATTCTCGGCGGTAAAGTTGTCGCCGGAGATGTGGTAGTAATTCGCTATGAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCAc                                                                                                                                <  2:2941804/150‑1 (MQ=255)
                        gAAGCGATTCTCGGCGGTAAAGTTGTCGCCGGAGATGTGGTAGTAATTCGCTATGAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCAcc                                                                                                                               <  2:3860866/150‑1 (MQ=255)
                          aGCGATTCTCGGCGGTAAAGTTGTCGCCGGAGATGTGGTAGTAATTCGCTATGAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGa                                                                                                                             <  1:991324/150‑1 (MQ=255)
                          aGCGATTCTCGGCGGTAAAGTTGTCGCCGGAGATGTGGTAGTAATTCGCTATGAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGa                                                                                                                             >  1:3154636/1‑150 (MQ=255)
                            cGATTCTCGGCGGTAAAGTTGTCGCCGGAGATGTGGTAGTAATTCGCTATGAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACg                                                                                                                           <  2:1057974/150‑1 (MQ=255)
                            cGATTCTCGGCGGTAAAGTTGTCGCCGGAGATGTGGTAGTAATTCGCTATGAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACg                                                                                                                           >  2:1524796/1‑150 (MQ=255)
                            cGATTCTCGGCGGTAAAGTTGTCGCCGGAGATGTGGTAGTAATTCGCTATGAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACg                                                                                                                           <  2:2752142/150‑1 (MQ=255)
                             gATTCTCGGCGGTAAAGTTGTCGCCGGAGATGTGGTAGTAATTCGCTATGAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACgg                                                                                                                          >  2:473489/1‑150 (MQ=255)
                               ttCTCGGCGGTAAAGTTGTCGCCGGAGATGTGGTAGTAATTCGCTATGAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATTCGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTc                                                                                                                        <  1:467768/150‑1 (MQ=255)
                               ttCTCGGCGGTAAAGTTGTCGCCGGAGATGTGGTAGTAATTCGCTATGAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTc                                                                                                                        <  1:573804/150‑1 (MQ=255)
                               ttCTCGGCGGTAAAGTTGTCGCCGGAGATGTGGTAGTAATTCGCTATGAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTc                                                                                                                        <  2:1757594/150‑1 (MQ=255)
                               ttCTCGGCGGTAAAGTTGTCGCCGGAGATGTGGTAGTAATTCGCTATGAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTAGGCAAAGCCTGTGCGCTGATCACCGACGGTc                                                                                                                        <  1:3750269/150‑1 (MQ=255)
                                tctcGGCGGTAAAGTTGTCGCCGGAGATGTGGTAGTAATTCGCTATGAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCg                                                                                                                       <  1:1003861/150‑1 (MQ=255)
                                 ctcGGCGGTAAAGTTGTCGCCGGAGATGTGGTAGTAATTCGCTATGAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGt                                                                                                                      >  1:1208978/1‑150 (MQ=255)
                                  tcGGCGGTAAAGTTGTCGCCGGAGATGTGGTAGTAATTCGCTATGAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGtt                                                                                                                     >  2:3014392/1‑150 (MQ=255)
                                   cggcggTAAAGTTGTCGCCGGAGATGTGGTAGTAATTCGCTATGAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGttt                                                                                                                    <  2:2038181/150‑1 (MQ=255)
                                   cggcggTAAAGTTGTCGCCGGAGATGTGGTAGTAATTCGCTATGAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGttt                                                                                                                    <  2:488014/150‑1 (MQ=255)
                                    ggcggTAAAGTTGTCGCCGGAGATGTGGTAGTAATTCGCTATGAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTtc                                                                                                                   >  1:3926398/1‑150 (MQ=255)
                                    ggcggTAAAGTTGTCGCCGGAGATGTGGTAGTAATTCGCTATGAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTtc                                                                                                                   >  2:4348911/1‑150 (MQ=255)
                                       ggTAAAGTTGTCGCCGGAGATGTGGTAGTAATTCGCTATGAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTctct                                                                                                                <  2:3905860/150‑1 (MQ=255)
                                       ggTAAAGTTGTCGCCGGAGATGTGGTAGTAATTCGCTATGAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTctct                                                                                                                <  1:2172724/150‑1 (MQ=255)
                                         tAAAGTTGTCGCCGGAGATGTGGTAGTAATTCGCTATGAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCtgg                                                                                                              >  1:2717359/1‑150 (MQ=255)
                                          aaaGTTGTCGCCGGAGATGTGGTAGTAATTCGCTATGAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCtggt                                                                                                             <  1:1533182/150‑1 (MQ=255)
                                            aGTTGTCGCCGGAGATGTGGTAGTAATTCGCTATGAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCtggtgg                                                                                                           >  2:5048582/1‑150 (MQ=255)
                                                 tCGCCGGAGATGTGGTAGTAATTCGCTATGAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACct                                                                                                      >  1:403801/1‑150 (MQ=255)
                                                 tCGCCGGAGATGTGGTAGTAATTCGCTATGAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTAGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACct                                                                                                      >  2:3232310/1‑150 (MQ=255)
                                                    ccGGAGATGTGGTAGTAATTCGCTATGAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACCTCTg                                                                                                   <  1:1581726/150‑1 (MQ=255)
                                                    ccGGAGATGTGGTAGTAATTCGCTATGAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACCTCTg                                                                                                   <  2:91990/150‑1 (MQ=255)
                                                     cGGAGATGTGGTAGTAATTCGCTATGAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACCTCTgg                                                                                                  <  1:2366065/150‑1 (MQ=255)
                                                             tgGTAGTAATTCGCTATGAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACCTCTGGTCTTTCCa                                                                                          >  1:1367813/1‑150 (MQ=255)
                                                             tgGTAGTAATTCGCTATGAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGCTTCTCTGGTGGCACCTCTGGTCTTTCCa                                                                                          <  1:3809905/150‑1 (MQ=255)
                                                              ggTAGTAATTCGCTATGAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACCTCTGGTCTTTCCAt                                                                                         <  1:2732051/150‑1 (MQ=255)
                                                              ggTAGTAATTCGCTATGAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACCTCTGGTCTTTCCAt                                                                                         <  2:2511084/150‑1 (MQ=255)
                                                              ggTAGTAATTCGCTATGAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTAGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACCTCTGGTCTTTCCAt                                                                                         >  1:2328678/1‑150 (MQ=255)
                                                                tagtaATTCGCTATGAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACCTCTGGTCTTTCCATCg                                                                                       >  2:437836/1‑150 (MQ=255)
                                                                tagtaATTCGCTATGAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACCTCTGGTCTTTCCATCg                                                                                       <  2:1597141/150‑1 (MQ=255)
                                                                  gtaATTCGCTATGAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACCTCTGGTCTTTCCATCGGc                                                                                     <  1:4703546/150‑1 (MQ=255)
                                                                  gtaATTCGCTATGAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACCTCTGGTCTTTCCATCGGc                                                                                     >  1:2946216/1‑150 (MQ=255)
                                                                   taATTCGCTATGAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACCTCTGGTCTTTCCATCGGcc                                                                                    <  1:2298436/150‑1 (MQ=255)
                                                                   taATTCGCTATGAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACCTCTGGTCTTTCCATCGGcc                                                                                    >  1:1281353/1‑150 (MQ=255)
                                                                    aaTTCGCTATGAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACCTCTGGTCTTTCCATCGGCCa                                                                                   <  2:4538293/150‑1 (MQ=255)
                                                                    aaTTCGCTATGAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACCTCTGGTCTTTCCATCGGCCa                                                                                   <  2:1931734/150‑1 (MQ=255)
                                                                      ttCGCTATGAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACCTCTGGTCTTTCCATCGGCCACg                                                                                 >  2:3447158/1‑150 (MQ=255)
                                                                          cTATGAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACCTCTGGTCTTTCCATCGGCCACGtctc                                                                             <  2:2496064/150‑1 (MQ=255)
                                                                           tATGAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACCTCTGGTCTTTCCATCGGCCACGTCTCa                                                                            <  2:39157/150‑1 (MQ=255)
                                                                              gAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTTCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACCTCTGGTCTTTCCATCGGCCACGTCTCACCg                                                                         <  1:4407427/150‑1 (MQ=255)
                                                                              gAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACCTCTGGTCTTTCCATCGGCCACGTCTCACCg                                                                         >  1:69106/1‑150 (MQ=255)
                                                                                 ggCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACCTCTGGTCTTTCCATCGGCCACGTCTCACCGGaa                                                                      >  2:4392906/1‑150 (MQ=255)
                                                                                  gCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACCTCTGGTCTTTCCATCGGCCACGTCTCACCGGaag                                                                     <  1:3518984/150‑1 (MQ=255)
                                                                                      gAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACCTCTGGTCTTTCCATCGGCCACGTCTCACCGGaagcggc                                                                 >  1:2987826/1‑150 (MQ=255)
                                                                                        aaGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACCTCTGGTCTTTCCATCGGCCACGTCTCACCGGaagcggcaa                                                               <  1:1136377/150‑1 (MQ=255)
                                                                                        aaGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACCTCTGGTCTTTCCATCGGCCACGTCTCACCGGaagcggcaa                                                               >  1:2852262/1‑150 (MQ=255)
                                                                                          ggCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACCTCTGGTCTTTCCATCGGCCACGTCTCACCGGaagcggcaagc                                                             >  1:848540/1‑150 (MQ=255)
                                                                                          ggCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACCTCTGGTCTTTCCATCGGCCACGTCTCACCGGaagcggcaagc                                                             <  1:4859936/150‑1 (MQ=255)
                                                                                            cGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACCTCTGGTCTTTCCATCGGCCACGTCTCACCGGaagcggcaagcgg                                                           <  2:2658313/150‑1 (MQ=255)
                                                                                            cGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACCTCTGGTCTTTCCATCGGCCACGTCTCACCGGaagcggcaagcgg                                                           >  1:1476482/1‑150 (MQ=255)
                                                                                            cGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACCTCTGGTCTTTCCATCGGCCACGTCTCACCGGaagcggcaagcgg                                                           >  1:743401/1‑150 (MQ=255)
                                                                                             ggTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACCTCTGGTCTTTCCATCGGCCACGTCTCACCGGaagcggcaagcggc                                                          >  2:726658/1‑150 (MQ=255)
                                                                                              gTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACCTCTGGTCTTTCCATCGGCCACGTCTCACCGGAAGCGGCAAgcggcg                                                         >  2:1197827/1‑150 (MQ=255)
                                                                                               tCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACCTCTGGTCTTTCCATCGGCCACGTCTCACCGGAAGCGGCAAGcggcgg                                                        <  1:5059538/150‑1 (MQ=255)
                                                                                               tCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACCTCTGGTCTTTCCATCGGCCACGTCTCACCGGAAGCGGCAAGcggcgg                                                        >  2:4200252/1‑150 (MQ=255)
                                                                                                ccGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACCTCTGGTCTTTCCATCGGCCACGTCTCACCGGAAGCGGCAAGCggcggc                                                       <  2:1828460/150‑1 (MQ=255)
                                                                                                ccGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACCTCTGGTCTTTCCATCGGCCACGTCTCACCGGAAGCGGCAAGCggcggc                                                       <  2:5030087/150‑1 (MQ=255)
                                                                                                  ggggATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACCTCTGGTCTTTCCATCGGCCACGTCTCACCGGAAGCGGCAAGCGGCGgcag                                                     >  1:3342300/1‑150 (MQ=255)
                                                                                                       tGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACCTCTGGTCTTTCCATCGGCCACGTCTCACCGGAAGCGGCAAGCGGCGGCAGCATTg                                                >  1:1811975/1‑150 (MQ=255)
                                                                                                       tGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACCTCTGGTCTTTCCATCGGCCACGTCTCACCGGAAGCGGCAAGCGGCGGCAGCATTg                                                <  1:1618760/150‑1 (MQ=255)
                                                                                                        gCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACCTCTGGTCTTTCCATCGGCCACGTCTCACCGGAAGCGGCAAGCGGCGGCAGCATTgg                                               >  2:683976/1‑150 (MQ=255)
                                                                                                         cAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACCTCTGGTCTTTCCATCGGCCACGTCTCACCGGAAGCGGCAAGCGGCGGCAGCATTGGc                                              >  1:4524289/1‑150 (MQ=255)
                                                                                                          aGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACCTCTGGTCTTTCCATCGGCCACGTCTCACCGGAAGCGGCAAGCGGCGGCAGCATTGGcc                                             >  1:3744481/1‑150 (MQ=255)
                                                                                                                tGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACCTCTGGTCTTTCCATCGGCCACGTCTCACCGGAAGCGGCAAGCGGCGGCAGCATTGGCCTGATTg                                       <  1:3981760/150‑1 (MQ=255)
                                                                                                                     tACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACCTCTGGTCTTTCCATCGGCCACGTCACACCGGAAGCGGCAAGCGGCGGCAGCATTGGCCTGATTGAAGAt                                  >  1:54494/1‑150 (MQ=255)
                                                                                                                      aCCCAACCAGCTTCCTGAAATCAAGGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTGTCGCGGGTGGCACCTCTGGTCTTTCCATCGGCCACGGCTCACCGGAAGCGGCCCGCGGCGGCCGCCGTGGCCGGAGTGCAgcgg                                 >  2:4331077/1‑147 (MQ=255)
                                                                                                                        ccAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACCTCTGGTCTTTCCATCGGCCACGTCTCACCGGAAGCGGCAAGCGGCGGCAGCATTGGCCTGATTGAAGATGGt                               >  2:3508050/1‑150 (MQ=255)
                                                                                                                           aCCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACCTCTGGTCTTTCCATCGGCCACGTCTCACCGGAAGCGGCAAGCGGCGGCAGCATTGGCCTGATTGAAGATGGTGAc                            <  1:423552/150‑1 (MQ=255)
                                                                                                                            ccAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACCTCTGGTCTTTCCATCGGCCACGTCTCACCGGAAGCGGCAAGCGGCGGCAGCATTGGCCTGATTGAAGATGGTGAcc                           >  2:1485822/1‑150 (MQ=255)
                                                                                                                             cAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACCTCTGGTCTTTCCATCGGCCACGTCTCACCGGAAGCGGCAAGCGGCGGCAGCATTGGCCTGATTGAAGATGGTGACCt                          >  2:2209979/1‑150 (MQ=255)
                                                                                                                              aGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACCTCTGGTCTTTCCATCGGCCACGTCTCACCGGAAGCGGCACGCGGCGGCAGCATTGGCCTGAGTGACGATGGTGACCTg                         >  2:3120354/1‑150 (MQ=255)
                                                                                                                              aGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACCTCTGGTCTTTCCATCGGCCACGTCTCACCGGAAGCGGCAAGCGGCGGCAGCATTGGCCTGATTGAAGATGGTGACCTg                         >  2:3771078/1‑150 (MQ=255)
                                                                                                                              aGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACCTCTGGTCTTTCCATCGGCCACGTCTCACCGGAAGCGGCAAGCGGCGGCAGCATTGGCCTGATTGAAGATGGTGACCTg                         >  2:2229804/1‑150 (MQ=255)
                                                                                                                                cTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACCTCTGGTCTTTCCATCGGCCACGTCTCACCGGAAGCGGCAAGCGGCGGCAGCATTGGCCTGATTGAAGATGGTGACCTGAt                       <  2:63219/150‑1 (MQ=255)
                                                                                                                                  tCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACCTCTGGTCTTTCCATCGGCCACGTCTCACCGGAAGCGGCAAGCGGCGGCAGCAGTGGCCTGATTGAAGATGGTGACCTGATCg                     >  2:3696038/1‑150 (MQ=255)
                                                                                                                                    cTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACCTCTGGTCTTTCCATCGGCCACGTCTCACCGGAAGCGGCAAGCGGCGGCAGCATTGGCCTGATTGAAGATGGTGACCTGATCGCt                   <  2:1367813/150‑1 (MQ=255)
                                                                                                                                     tGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACCTCTGGTCTTTCCATCGGCCACGTCTCACCGGAAGCGGCAAGCGGCGGCAGCAGTGGCCTGATTGAAGATGGTGACCTGATCGCTa                  >  2:1251738/1‑150 (MQ=255)
                                                                                                                                            aaTGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACCTCTGGTCTTTCCATCGGCCACGTCTCACCGGAAGCGGCAAGCGGCGGCAGCATTGGCCTGATTGAAGATGGTGACCTGATCGCTATCGACAt           <  1:2594534/150‑1 (MQ=255)
                                                                                                                                              tGGGTCTAGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACCTCTGGTCTTTCCATCGGCCACGTCTCACCGGAAGCGGCAAGCGGCGGCAGCATTGGCCTGATTGAAGATGGTGACCTGATCGCTATCGACATcc         <  1:3162685/150‑1 (MQ=255)
                                                                                                                                                ggTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACCTCTGGTCTTTCCATCGGCCACGTCTCACCGGAAGCGGCAAGCGGCGGCAGCATTGGCCTGATTGAAGATGGTGACCTGATCGCTATCGACATCCCg       <  1:5048334/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TACGAAAGCCAGGACGATGCGGTAGAAGCGATTCTCGGCGGTAAAGTTGTCGCCGGAGATGTGGTAGTAATTCGCTATGAAGGCCCGAAAGGCGGTCCGGGGATGCAGGAAATGCTCTACCCAACCAGCTTCCTGAAATCAATGGGTCTCGGCAAAGCCTGTGCGCTGATCACCGACGGTCGTTTCTCTGGTGGCACCTCTGGTCTTTCCATCGGCCACGTCTCACCGGAAGCGGCAAGCGGCGGCAGCATTGGCCTGATTGAAGATGGTGACCTGATCGCTATCGACATCCCGAACCG  >  CP009273/3948200‑3948498

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: