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,6140TA97.6% 1003.5 / 9.8 291*150R (TGA→AGA) rpiBribose 5‑phosphate isomerase B/allose 6‑phosphate isomerase
Reads supporting (aligned to +/- strand):  ref base T (5/2);  new base A (148/136);  total (153/138)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 301 total aligned reads displayed.

GACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGAC  >  CP009273/4303467‑4303762
                                                                                                                                                   |                                                                                                                                                    
gACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAga                                                                                                                                                    >  1:2108780/1‑150 (MQ=255)
 acacCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATagag                                                                                                                                                   <  1:2909283/150‑1 (MQ=255)
 acacCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATagag                                                                                                                                                   <  2:1383634/150‑1 (MQ=255)
 acacCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATagag                                                                                                                                                   <  2:1547506/150‑1 (MQ=255)
 acacCAACGCGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATagag                                                                                                                                                   <  1:579220/150‑1 (MQ=255)
  cacCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAgaga                                                                                                                                                  >  1:1748239/1‑150 (MQ=255)
  cacCAACGTGCTGGCTGTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAgaga                                                                                                                                                  >  2:3482730/1‑150 (MQ=255)
       aCGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCAt                                                                                                                                             <  2:125467/150‑1 (MQ=255)
         gTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATcc                                                                                                                                           >  1:3130433/1‑150 (MQ=255)
          tGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCa                                                                                                                                          <  1:709855/150‑1 (MQ=255)
              ggCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCactac                                                                                                                                      <  1:3725149/150‑1 (MQ=255)
               gCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTTGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGAGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCactact                                                                                                                                     >  2:4475560/1‑150 (MQ=255)
                cTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACtt                                                                                                                                    >  1:629390/1‑150 (MQ=255)
                 ttttGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCACACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTg                                                                                                                                   <  2:1535814/150‑1 (MQ=255)
                   ttGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCa                                                                                                                                 >  1:1018073/1‑150 (MQ=255)
                    tGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCAt                                                                                                                                <  1:4585996/150‑1 (MQ=255)
                        tCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGAt                                                                                                                            >  1:4052926/1‑150 (MQ=255)
                          aCGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGa                                                                                                                          >  2:2955152/1‑150 (MQ=255)
                           cGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGAGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAg                                                                                                                         <  1:1392594/150‑1 (MQ=255)
                            gAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGt                                                                                                                        >  2:652556/1‑150 (MQ=255)
                            gAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGt                                                                                                                        <  1:2343346/150‑1 (MQ=255)
                            gAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGt                                                                                                                        >  2:1038231/1‑150 (MQ=255)
                             aGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTa                                                                                                                       >  1:2890961/1‑150 (MQ=255)
                               tGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAAt                                                                                                                     >  2:3779314/1‑150 (MQ=255)
                                 gTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGa                                                                                                                   <  1:2584548/150‑1 (MQ=255)
                                  ttGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGAt                                                                                                                  <  2:1096740/150‑1 (MQ=255)
                                  ttGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGAt                                                                                                                  >  1:987249/1‑150 (MQ=255)
                                  ttGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGAt                                                                                                                  >  1:3620980/1‑150 (MQ=255)
                                    ggCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTa                                                                                                                <  2:2371786/150‑1 (MQ=255)
                                          gAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTgg                                                                                                          <  1:858815/150‑1 (MQ=255)
                                           aaCTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGa                                                                                                         <  2:1081051/150‑1 (MQ=255)
                                           aaCTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGa                                                                                                         >  2:2475712/1‑150 (MQ=255)
                                            aCTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGAt                                                                                                        >  1:755207/1‑150 (MQ=255)
                                            aCTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGAt                                                                                                        >  2:1901754/1‑150 (MQ=255)
                                                 cAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGt                                                                                                   >  1:3305040/1‑150 (MQ=255)
                                                   aaaaTGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGttt                                                                                                 <  1:3106784/150‑1 (MQ=255)
                                                      aTGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCtt                                                                                              >  2:2424183/1‑150 (MQ=255)
                                                       tGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCttt                                                                                             <  2:2146546/150‑1 (MQ=255)
                                                        gATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCtttt                                                                                            <  1:904274/150‑1 (MQ=255)
                                                        gATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCtttt                                                                                            >  1:2139557/1‑150 (MQ=255)
                                                           tgtgGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGa                                                                                         <  2:178413/150‑1 (MQ=255)
                                                            gtgGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAg                                                                                        >  1:325658/1‑150 (MQ=255)
                                                              ggATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGt                                                                                      >  1:1422105/1‑150 (MQ=255)
                                                              ggATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGt                                                                                      >  2:3741331/1‑150 (MQ=255)
                                                               gATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGtt                                                                                     <  1:4112187/150‑1 (MQ=255)
                                                               gATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGtt                                                                                     <  2:4281806/150‑1 (MQ=255)
                                                                 tGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTgg                                                                                   <  2:2661782/150‑1 (MQ=255)
                                                                 tGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTgg                                                                                   <  2:83826/150‑1 (MQ=255)
                                                                    gTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTa                                                                                >  2:3910885/1‑150 (MQ=255)
                                                                    gTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTa                                                                                >  1:2476064/1‑150 (MQ=255)
                                                                    gTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTa                                                                                >  1:1211617/1‑150 (MQ=255)
                                                                    gTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTa                                                                                >  1:847034/1‑150 (MQ=255)
                                                                      ggCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAAt                                                                              >  1:4510112/1‑150 (MQ=255)
                                                                        cTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGa                                                                            >  1:2930776/1‑150 (MQ=255)
                                                                        cTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGa                                                                            >  1:3560975/1‑150 (MQ=255)
                                                                        cTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGa                                                                            >  1:4266237/1‑150 (MQ=255)
                                                                         tGGGCGCACAGTACGAAGGCTGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGTAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGaa                                                                           <  2:3339695/150‑1 (MQ=255)
                                                                         tGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTTATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGaa                                                                           >  2:3956864/1‑150 (MQ=255)
                                                                         tGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGaa                                                                           >  1:1770475/1‑150 (MQ=255)
                                                                         tGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGaa                                                                           >  2:4171040/1‑150 (MQ=255)
                                                                           ggCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAAc                                                                         <  2:1638558/150‑1 (MQ=255)
                                                                              gcACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCt                                                                      >  2:4064230/1‑150 (MQ=255)
                                                                               cacaGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCtt                                                                     <  1:1690900/150‑1 (MQ=255)
                                                                                     aCGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCt                                                               <  1:4288259/150‑1 (MQ=255)
                                                                                      cGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTa                                                              >  1:3993479/1‑150 (MQ=255)
                                                                                         aGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTAttt                                                           <  1:332104/150‑1 (MQ=255)
                                                                                         aGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTAttt                                                           >  2:598771/1‑150 (MQ=255)
                                                                                             ggTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCttt                                                       <  1:1328468/150‑1 (MQ=255)
                                                                                             ggTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCttt                                                       <  1:1533441/150‑1 (MQ=255)
                                                                                             ggTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCttt                                                       <  2:2133160/150‑1 (MQ=255)
                                                                                               tcgtcATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTa                                                     <  1:1577400/150‑1 (MQ=255)
                                                                                                cgtcATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTAt                                                    >  2:4640402/1‑150 (MQ=255)
                                                                                                 gtcATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCGGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCAGCTATTTCTTTTAtt                                                   >  1:3683186/1‑150 (MQ=255)
                                                                                                 gtcATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTAtt                                                   >  1:786610/1‑150 (MQ=255)
                                                                                                    atcaACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCAc                                                >  1:613395/1‑150 (MQ=255)
                                                                                                    atcaACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCAc                                                >  1:1450949/1‑150 (MQ=255)
                                                                                                    atcaACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCAc                                                <  2:1705461/150‑1 (MQ=255)
                                                                                                       aacaacGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCt                                             >  2:3268320/1‑150 (MQ=255)
                                                                                                           acGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTc                                         <  2:2108611/150‑1 (MQ=255)
                                                                                                                  gAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGt                                  <  2:1967783/150‑1 (MQ=255)
                                                                                                                     gCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGGTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCa                               >  2:3136778/1‑150 (MQ=255)
                                                                                                                      cGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCAc                              >  1:1831332/1‑150 (MQ=255)
                                                                                                                        aTTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCg                            >  1:3415473/1‑150 (MQ=255)
                                                                                                                           aCGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTc                         <  2:1001274/150‑1 (MQ=255)
                                                                                                                             ggCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCg                       >  1:2210645/1‑150 (MQ=255)
                                                                                                                             ggCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCg                       <  1:2295805/150‑1 (MQ=255)
                                                                                                                             ggCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCg                       >  2:1467079/1‑150 (MQ=255)
                                                                                                                               cAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGc                     <  1:2723694/150‑1 (MQ=255)
                                                                                                                                aaTAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGcc                    <  2:4599226/150‑1 (MQ=255)
                                                                                                                                  tAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCa                  <  2:755207/150‑1 (MQ=255)
                                                                                                                                   agagCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAg                 <  2:3155376/150‑1 (MQ=255)
                                                                                                                                     agcagcGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGc               >  1:3537078/1‑150 (MQ=255)
                                                                                                                                      gcagcGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGcc              >  1:2277643/1‑150 (MQ=255)
                                                                                                                                          cGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGata          >  1:1345858/1‑150 (MQ=255)
                                                                                                                                          cGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGata          >  1:3750600/1‑150 (MQ=255)
                                                                                                                                          cGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTGCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGata          >  1:1126448/1‑150 (MQ=255)
                                                                                                                                              gaAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGAtattat      >  1:1398316/1‑150 (MQ=255)
                                                                                                                                              gaAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGAtattat      >  2:2029101/1‑150 (MQ=255)
                                                                                                                                                aaTAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATcg    >  1:198388/1‑150 (MQ=255)
                                                                                                                                                  ttGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATcgac  <  1:1625223/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGAC  >  CP009273/4303467‑4303762

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: