Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,303,614 T→A 100% *150R (TGA→AGA)  rpiB → ribose 5‑phosphate isomerase B/allose 6‑phosphate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,303,6140TA99.7% 976.6 / ‑5.0 294*150R (TGA→AGA) rpiBribose 5‑phosphate isomerase B/allose 6‑phosphate isomerase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base A (124/169);  total (124/170)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 296 total aligned reads displayed.

GACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCG  >  CP009273/4303467‑4303760
                                                                                                                                                   |                                                                                                                                                  
 acacCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATagag                                                                                                                                                 <  2:1088075/150‑1 (MQ=255)
    ccAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGAtt                                                                                                                                              <  2:4331374/150‑1 (MQ=255)
     cAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTc                                                                                                                                             <  1:4306988/150‑1 (MQ=255)
     cAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTc                                                                                                                                             <  1:1142999/150‑1 (MQ=255)
        cGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATc                                                                                                                                          <  2:4671887/150‑1 (MQ=255)
          tGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCa                                                                                                                                        <  2:131625/150‑1 (MQ=255)
            cTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCact                                                                                                                                      <  1:2738978/150‑1 (MQ=255)
            cTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCact                                                                                                                                      <  1:5057216/150‑1 (MQ=255)
            cTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCact                                                                                                                                      <  2:2477735/150‑1 (MQ=255)
             tGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCacta                                                                                                                                     <  2:3874997/150‑1 (MQ=255)
              ggCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCactac                                                                                                                                    >  1:951411/1‑150 (MQ=255)
                cTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACtt                                                                                                                                  <  1:1762180/150‑1 (MQ=255)
                cTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACtt                                                                                                                                  >  1:3474844/1‑150 (MQ=255)
                    tGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCAt                                                                                                                              <  2:2384470/150‑1 (MQ=255)
                    tGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCAt                                                                                                                              <  2:231749/150‑1 (MQ=255)
                     ggTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATg                                                                                                                             <  2:3653768/150‑1 (MQ=255)
                     ggTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATg                                                                                                                             >  1:762530/1‑150 (MQ=255)
                      gatCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCTATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATgg                                                                                                                            >  1:3272435/3‑150 (MQ=255)
                           cGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAg                                                                                                                       <  1:3912787/150‑1 (MQ=255)
                             aGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTa                                                                                                                     <  1:2228377/150‑1 (MQ=255)
                             aGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTAATTGCATGGATGAGTa                                                                                                                     >  1:159645/1‑150 (MQ=255)
                              gTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTAATCCACTACTTGCATGGATGAGTaa                                                                                                                    <  1:3996493/150‑1 (MQ=255)
                                ggTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATg                                                                                                                  <  1:3054978/150‑1 (MQ=255)
                                ggTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATg                                                                                                                  >  1:4882629/1‑150 (MQ=255)
                                      ccTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAAt                                                                                                            <  2:4582846/150‑1 (MQ=255)
                                         cGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAAtgtg                                                                                                         <  2:377034/150‑1 (MQ=255)
                                             cTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATa                                                                                                     >  1:1996675/1‑150 (MQ=255)
                                             cTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATa                                                                                                     >  1:2482803/1‑150 (MQ=255)
                                              tGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATag                                                                                                    >  2:1557817/1‑150 (MQ=255)
                                              tGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATag                                                                                                    <  2:340231/150‑1 (MQ=255)
                                               ggCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATaga                                                                                                   >  2:3006680/1‑150 (MQ=255)
                                                 cAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGt                                                                                                 <  2:2211765/150‑1 (MQ=255)
                                                 cAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGt                                                                                                 <  1:1371432/150‑1 (MQ=255)
                                                  aaaaaTGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGtt                                                                                                >  1:1027010/1‑150 (MQ=255)
                                                  aaaaaTGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGtt                                                                                                <  2:198085/150‑1 (MQ=255)
                                                  aaaaaTGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGtt                                                                                                >  2:264882/1‑150 (MQ=255)
                                                  aaaaaTGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGtt                                                                                                <  1:3645342/150‑1 (MQ=255)
                                                        gATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCtttt                                                                                          >  1:4026950/1‑150 (MQ=255)
                                                         aTTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCttttt                                                                                         <  2:119465/150‑1 (MQ=255)
                                                         aTTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCttttt                                                                                         >  1:821047/1‑150 (MQ=255)
                                                          ttGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTg                                                                                        <  2:4193328/150‑1 (MQ=255)
                                                           tgtgGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGa                                                                                       <  1:3621661/150‑1 (MQ=255)
                                                            gtgGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAg                                                                                      <  2:1877387/150‑1 (MQ=255)
                                                             tgGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAgg                                                                                     <  2:2447171/150‑1 (MQ=255)
                                                              ggATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGt                                                                                    <  1:3291187/150‑1 (MQ=255)
                                                               gATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGtt                                                                                   <  1:3735028/150‑1 (MQ=255)
                                                                aTGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTg                                                                                  <  2:839704/150‑1 (MQ=255)
                                                                       gCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATg                                                                           <  1:2407763/150‑1 (MQ=255)
                                                                       gCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATg                                                                           <  2:3867930/150‑1 (MQ=255)
                                                                         tGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGaa                                                                         >  2:754497/1‑150 (MQ=255)
                                                                          gggCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGaaa                                                                        <  1:2827544/150‑1 (MQ=255)
                                                                                acaGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCttt                                                                  <  2:81277/150‑1 (MQ=255)
                                                                                 caGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTc                                                                 <  2:3028092/150‑1 (MQ=255)
                                                                                  aGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTcc                                                                <  2:1755743/150‑1 (MQ=255)
                                                                                        aaGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTAtt                                                          <  2:2042154/150‑1 (MQ=255)
                                                                                        aaGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTAtt                                                          >  2:4028554/1‑150 (MQ=255)
                                                                                          ggCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTc                                                        <  1:4629584/150‑1 (MQ=255)
                                                                                          ggCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTc                                                        <  1:424368/150‑1 (MQ=255)
                                                                                            cGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGTGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCATTCCGCGATTTCtt                                                      >  1:1370147/1‑150 (MQ=255)
                                                                                             ggTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCttt                                                     <  1:4345325/150‑1 (MQ=255)
                                                                                             ggTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCttt                                                     <  1:3830186/150‑1 (MQ=255)
                                                                                               tcgtcATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTa                                                   <  2:3133744/150‑1 (MQ=255)
                                                                                                cgtcATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCGATTTCTTTTAt                                                  >  1:4935553/1‑150 (MQ=255)
                                                                                                    atcaACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCAc                                              >  2:4363135/1‑150 (MQ=255)
                                                                                                    atcaACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCAc                                              >  1:2163473/1‑150 (MQ=255)
                                                                                                    atcaACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCAc                                              <  1:1750525/150‑1 (MQ=255)
                                                                                                     tcaACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCAcc                                             <  1:2513771/150‑1 (MQ=255)
                                                                                                      caacaaCGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCAccc                                            >  2:1367234/1‑150 (MQ=255)
                                                                                                       aacaacGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCt                                           <  1:5199706/150‑1 (MQ=255)
                                                                                                       aacaacGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCt                                           >  1:264530/1‑150 (MQ=255)
                                                                                                       aacaacGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCt                                           >  1:4796220/1‑150 (MQ=255)
                                                                                                       aacaacGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCt                                           >  2:1014380/1‑150 (MQ=255)
                                                                                                       aacaacGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCt                                           >  1:2949074/1‑150 (MQ=255)
                                                                                                        acaacGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAAGAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTg                                          <  2:5056454/150‑1 (MQ=255)
                                                                                                         caacGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGTTTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGc                                         <  2:3478963/150‑1 (MQ=255)
                                                                                                           acGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTc                                       <  2:951411/150‑1 (MQ=255)
                                                                                                           acGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTc                                       <  1:4131443/150‑1 (MQ=255)
                                                                                                           acGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGAGAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTc                                       <  1:1467909/150‑1 (MQ=255)
                                                                                                             gcgTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCAc                                     <  1:64056/150‑1 (MQ=255)
                                                                                                                 ggAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTg                                 >  2:877526/1‑150 (MQ=255)
                                                                                                                    ggCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTcc                              >  2:5290190/1‑150 (MQ=255)
                                                                                                                     gctaTTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCa                             <  2:3272435/147‑1 (MQ=255)
                                                                                                                       gATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCAcc                           <  1:5297695/150‑1 (MQ=255)
                                                                                                                        aTTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCg                          >  1:858315/1‑150 (MQ=255)
                                                                                                                        aTTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCg                          >  2:1967456/1‑150 (MQ=255)
                                                                                                                               cAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGc                   >  2:1605774/1‑150 (MQ=255)
                                                                                                                               cAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGc                   >  2:577186/1‑150 (MQ=255)
                                                                                                                               cAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGc                   >  2:5153244/1‑150 (MQ=255)
                                                                                                                                aaTAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGcc                  >  1:1624757/1‑150 (MQ=255)
                                                                                                                                 aTAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGccc                 >  1:2675209/1‑150 (MQ=255)
                                                                                                                                  tAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCa                <  2:4927547/150‑1 (MQ=255)
                                                                                                                                   agagCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAg               >  2:653999/1‑150 (MQ=255)
                                                                                                                                   agagCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAg               >  2:4873285/1‑150 (MQ=255)
                                                                                                                                    gagCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAgg              <  2:2805446/150‑1 (MQ=255)
                                                                                                                                     agcagcGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGc             >  2:3322657/1‑150 (MQ=255)
                                                                                                                                         gcGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGat         <  1:1107500/150‑1 (MQ=255)
                                                                                                                                          cGGAGAAAGTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCAGata        <  2:3049161/150‑1 (MQ=255)
                                                                                                                                            gagaAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGAtatt      >  2:1808569/1‑150 (MQ=255)
                                                                                                                                               aaaTAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATc   >  1:3237390/1‑150 (MQ=255)
                                                                                                                                               aaaTAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATc   <  2:2879658/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCG  >  CP009273/4303467‑4303760

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: