Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,303,614 T→A 44.1% *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,6140TA44.1% 127.8 / inf 302*150R (TGA→AGA) rpiBribose 5‑phosphate isomerase B/allose 6‑phosphate isomerase
Reads supporting (aligned to +/- strand):  ref base T (81/87);  new base A (49/84);  total (131/171)
Fisher's exact test for biased strand distribution p-value = 6.07e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

Only 100 of 303 total aligned reads displayed.

CGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACG  >  CP009273/4303466‑4303763
                                                                                                                                                    |                                                                                                                                                     
cGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTg                                                                                                                                                      >  2:20717/1‑150 (MQ=255)
 gACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAga                                                                                                                                                     >  2:1410303/1‑150 (MQ=255)
    acCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGAt                                                                                                                                                  <  2:3867860/150‑1 (MQ=255)
     cctaCGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGTCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGAtt                                                                                                                                                 <  2:4755508/147‑1 (MQ=255)
       aaCGTGCTGGCTTTTTGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCa                                                                                                                                               <  2:4154068/150‑1 (MQ=255)
       aaCGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCa                                                                                                                                               >  2:737835/1‑150 (MQ=255)
        aCGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCAt                                                                                                                                              <  2:446280/150‑1 (MQ=255)
          gTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATcc                                                                                                                                            <  1:1051969/150‑1 (MQ=255)
             cTGGCTTTTGGGTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCact                                                                                                                                         <  2:4219927/150‑1 (MQ=255)
              tGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCacta                                                                                                                                        <  1:2612079/150‑1 (MQ=255)
                 cTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACtt                                                                                                                                     <  2:918797/150‑1 (MQ=255)
                 cTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACtt                                                                                                                                     <  2:5118035/150‑1 (MQ=255)
                  ttttGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTg                                                                                                                                    <  1:1413202/150‑1 (MQ=255)
                  ttttGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTg                                                                                                                                    <  2:1593207/150‑1 (MQ=255)
                  ttttGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTg                                                                                                                                    >  1:2188765/1‑150 (MQ=255)
                         tCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGAt                                                                                                                             <  2:4653076/150‑1 (MQ=255)
                          cACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATg                                                                                                                            >  2:3161882/1‑150 (MQ=255)
                          cACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATg                                                                                                                            >  1:1287263/1‑150 (MQ=255)
                              aGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTa                                                                                                                        >  2:5081639/1‑150 (MQ=255)
                              aGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTa                                                                                                                        <  2:530048/150‑1 (MQ=255)
                                tGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAAt                                                                                                                      <  1:1495858/150‑1 (MQ=255)
                                   ttGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGAt                                                                                                                   <  1:163898/150‑1 (MQ=255)
                                      gCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTaa                                                                                                                <  2:4552437/150‑1 (MQ=255)
                                        cTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAAGGATTAAtg                                                                                                              >  1:3781777/1‑150 (MQ=255)
                                         tCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAAtgt                                                                                                             >  2:978094/1‑150 (MQ=255)
                                           gAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTgg                                                                                                           >  2:4826423/1‑150 (MQ=255)
                                           gAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTgg                                                                                                           <  1:4420661/150‑1 (MQ=255)
                                            aaCTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGa                                                                                                          >  2:1521455/1‑150 (MQ=255)
                                            aaCTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGAGTAATGTGGa                                                                                                          >  1:3610350/1‑150 (MQ=255)
                                              cTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATa                                                                                                        >  2:4706345/1‑150 (MQ=255)
                                              cTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATa                                                                                                        >  2:2102205/1‑150 (MQ=255)
                                              cTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATa                                                                                                        <  1:506599/150‑1 (MQ=255)
                                                 gCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATagag                                                                                                     <  2:1345707/150‑1 (MQ=255)
                                                   aaaaaTGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGAGGTCGTCATCAACAACTCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATACATCCACTACTTGCATGGATGAGTAATGATTAATTTGGATAGAGtt                                                                                                   >  2:2590257/1‑150 (MQ=255)
                                                      aaTGATTGTGGATGCGTGGCTGGGCGCACCGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCt                                                                                                >  2:187976/1‑150 (MQ=255)
                                                      aaTGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCt                                                                                                >  1:3275727/1‑150 (MQ=255)
                                                          aTTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCttttt                                                                                            >  1:3151682/1‑150 (MQ=255)
                                                           ttGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTg                                                                                           >  1:4124572/1‑150 (MQ=255)
                                                           ttGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTg                                                                                           <  1:1911739/150‑1 (MQ=255)
                                                             gtgGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAg                                                                                         <  2:681767/150‑1 (MQ=255)
                                                              tgGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAACTTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAgg                                                                                        >  1:3034002/1‑150 (MQ=255)
                                                                gATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGtt                                                                                      >  2:3565967/1‑150 (MQ=255)
                                                                gATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGtt                                                                                      >  2:4292054/1‑150 (MQ=255)
                                                                 aTGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTg                                                                                     >  1:2334071/1‑150 (MQ=255)
                                                                 aTGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTg                                                                                     <  1:630022/150‑1 (MQ=255)
                                                                 aTGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTg                                                                                     <  1:875072/150‑1 (MQ=255)
                                                                    cGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCt                                                                                  >  2:430955/1‑150 (MQ=255)
                                                                    cGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCt                                                                                  >  1:2837416/1‑150 (MQ=255)
                                                                      tGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTaa                                                                                >  1:460861/1‑150 (MQ=255)
                                                                       ggCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAAt                                                                               <  2:3232687/150‑1 (MQ=255)
                                                                        gCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATg                                                                              <  1:4978131/150‑1 (MQ=255)
                                                                          tGGGAGCACATGACGAAGGCGGTCGTCATCAACATCGCGTGGAGGCGATTACGGCAATAGAGCTGCTGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGaa                                                                            <  2:4820175/150‑1 (MQ=255)
                                                                            ggCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTACTGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAAc                                                                          <  1:517380/150‑1 (MQ=255)
                                                                             gcgcACAGTACGAAGGCGGTCGTCATCAACTACGCGTGGAGGCTATTACGGCTATAGATCTTCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTTATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACg                                                                         <  1:3543822/150‑1 (MQ=255)
                                                                               gcACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCt                                                                       <  1:540048/150‑1 (MQ=255)
                                                                               gcACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCt                                                                       >  1:401519/1‑150 (MQ=255)
                                                                                    gTACGAAGGCGGTCGTCATCAACAACGCGTGGATGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAATCGCTTTCCg                                                                  <  1:1082153/150‑1 (MQ=255)
                                                                                    gTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCACTAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCg                                                                  <  2:2184825/150‑1 (MQ=255)
                                                                                     tACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGc                                                                 >  1:4658027/1‑150 (MQ=255)
                                                                                      aCGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCt                                                                <  1:1964513/150‑1 (MQ=255)
                                                                                       cGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTa                                                               >  1:2743861/1‑150 (MQ=255)
                                                                                       cGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTa                                                               <  2:859223/150‑1 (MQ=255)
                                                                                        gAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTAt                                                              >  2:3740227/1‑150 (MQ=255)
                                                                                           ggCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTc                                                           >  1:3967833/1‑150 (MQ=255)
                                                                                             cGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCtt                                                         <  1:476976/150‑1 (MQ=255)
                                                                                              ggTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCttt                                                        <  2:4401783/150‑1 (MQ=255)
                                                                                              ggTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCttt                                                        >  2:17458/1‑150 (MQ=255)
                                                                                              ggTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCttt                                                        >  2:877870/1‑150 (MQ=255)
                                                                                              ggTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCttt                                                        >  2:3015050/1‑150 (MQ=255)
                                                                                                tcgtcATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTa                                                      <  1:1853344/150‑1 (MQ=255)
                                                                                                   tcatcaACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTc                                                   >  1:486338/1‑150 (MQ=255)
                                                                                                   tcatcaACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTc                                                   >  1:4277891/1‑150 (MQ=255)
                                                                                                     atcaACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCAc                                                 >  2:4978460/1‑150 (MQ=255)
                                                                                                       caacaaCGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTACTGATTAAGGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTCTGGCTGTTATTCAccc                                               >  2:2673979/1‑150 (MQ=255)
                                                                                                        aacaacGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCt                                              <  1:2491493/150‑1 (MQ=255)
                                                                                                         acaacGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTg                                             <  1:812657/150‑1 (MQ=255)
                                                                                                          ccacGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGc                                            <  2:3233131/148‑1 (MQ=255)
                                                                                                          caacGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGc                                            >  2:3045553/1‑150 (MQ=255)
                                                                                                           aacGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCt                                           <  2:2869175/150‑1 (MQ=255)
                                                                                                           aacGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGAATAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCt                                           <  2:1050482/150‑1 (MQ=255)
                                                                                                              gcgTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCAc                                        <  2:2408409/150‑1 (MQ=255)
                                                                                                                 tGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCt                                     >  1:3841109/1‑150 (MQ=255)
                                                                                                                  ggAGGCGATTACGGCAATAGCGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGAGGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTg                                    <  2:3983008/150‑1 (MQ=255)
                                                                                                                       cGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCAc                               >  2:3660943/1‑150 (MQ=255)
                                                                                                                        gATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCAcc                              >  1:3335616/1‑150 (MQ=255)
                                                                                                                             cGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTcc                         >  1:3572576/1‑150 (MQ=255)
                                                                                                                                   tAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGTGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCa                   >  2:1914121/1‑150 (MQ=255)
                                                                                                                                     gagCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAgg                 >  2:3489299/1‑150 (MQ=255)
                                                                                                                                      agcagcGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGc                >  1:1998949/1‑150 (MQ=255)
                                                                                                                                        cagcGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCg              <  2:1711383/150‑1 (MQ=255)
                                                                                                                                        cagcGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCg              <  1:1041876/150‑1 (MQ=255)
                                                                                                                                         agcGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGa             <  1:234247/150‑1 (MQ=255)
                                                                                                                                           cGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGata           >  2:3807012/1‑150 (MQ=255)
                                                                                                                                           cGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGata           >  1:1721775/1‑150 (MQ=255)
                                                                                                                                           cGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGata           >  1:943483/1‑150 (MQ=255)
                                                                                                                                            ggAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGatat          >  1:922699/1‑150 (MQ=255)
                                                                                                                                                aaaTAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTATTGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATc      <  2:3524822/150‑1 (MQ=255)
                                                                                                                                                 aaTTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATcg     <  1:2134749/150‑1 (MQ=255)
                                                                                                                                                 aaTTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATcg     >  1:5197739/1‑150 (MQ=255)
                                                                                                                                                  aTTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATcga    >  2:1243777/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGACG  >  CP009273/4303466‑4303763

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: