Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,303,614 T→A 86.2% *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,6140TA86.2% 734.0 / 90.8 239*150R (TGA→AGA) rpiBribose 5‑phosphate isomerase B/allose 6‑phosphate isomerase
Reads supporting (aligned to +/- strand):  ref base T (12/21);  new base A (87/119);  total (99/140)
Fisher's exact test for biased strand distribution p-value = 5.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 241 total aligned reads displayed.

CGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGAC  >  CP009273/4303466‑4303762
                                                                                                                                                    |                                                                                                                                                    
cGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTg                                                                                                                                                     >  2:2277091/1‑150 (MQ=255)
   cacCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTgaga                                                                                                                                                  >  2:2805656/1‑150 (MQ=255)
    acCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGAt                                                                                                                                                 <  1:1147577/150‑1 (MQ=255)
     ccAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGAtt                                                                                                                                                <  2:421972/150‑1 (MQ=255)
        aCGTGCTGGCTTTTGGTTCACGAGTGGTTGACCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCAt                                                                                                                                             <  2:1779534/150‑1 (MQ=255)
            gCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCac                                                                                                                                         <  1:1331287/150‑1 (MQ=255)
             cTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCact                                                                                                                                        >  2:5081871/1‑150 (MQ=255)
                 cTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACtt                                                                                                                                    >  1:4681407/1‑150 (MQ=255)
                 cTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACtt                                                                                                                                    <  2:4647906/150‑1 (MQ=255)
                  ttttGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTg                                                                                                                                   <  2:3786802/150‑1 (MQ=255)
                   tttGGTTCACGAGTGGTGGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGc                                                                                                                                  <  1:4438033/150‑1 (MQ=255)
                            cGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAg                                                                                                                         <  2:2616250/150‑1 (MQ=255)
                            cGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAg                                                                                                                         >  1:3883877/1‑150 (MQ=255)
                             gAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGt                                                                                                                        >  1:3634455/1‑150 (MQ=255)
                              aGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTa                                                                                                                       >  1:2914298/1‑150 (MQ=255)
                                   ttGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGAt                                                                                                                  <  1:2351391/150‑1 (MQ=255)
                                      gCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTaa                                                                                                               >  1:1578074/1‑150 (MQ=255)
                                        cTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAAtg                                                                                                             >  2:2721605/1‑150 (MQ=255)
                                        cTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAAtg                                                                                                             >  2:4424547/1‑150 (MQ=255)
                                        cTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAAtg                                                                                                             >  2:104619/1‑150 (MQ=255)
                                         tCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAAtgt                                                                                                            <  1:4916821/150‑1 (MQ=255)
                                         tCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAAtgt                                                                                                            <  1:2611173/150‑1 (MQ=255)
                                           gAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTgg                                                                                                          <  2:2781364/150‑1 (MQ=255)
                                           gAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTgg                                                                                                          <  2:1894256/150‑1 (MQ=255)
                                           gAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTgg                                                                                                          <  2:1893968/150‑1 (MQ=255)
                                               tGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATag                                                                                                      >  1:2243274/1‑150 (MQ=255)
                                                ggCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATaga                                                                                                     >  1:4454033/1‑150 (MQ=255)
                                                   aaaaaTGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGtt                                                                                                  <  1:2403606/150‑1 (MQ=255)
                                                    aaaaTGATTGTGGATGCGTGGCTGGGAGAACAGTACGAAGGCGGTCGACATCAACAACGCGTGGAGGCGATTACGGCAAAAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAAAGTGTATAGAGttt                                                                                                 >  1:4257945/1‑150 (MQ=255)
                                                       aTGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCtt                                                                                              <  1:4809152/150‑1 (MQ=255)
                                                         gagtgtgGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGAGTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCtttt                                                                                            <  2:3603470/147‑1 (MQ=255)
                                                          aTTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCttttt                                                                                           <  1:3846540/150‑1 (MQ=255)
                                                                gATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGtt                                                                                     <  2:1582496/150‑1 (MQ=255)
                                                                gATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGtt                                                                                     <  2:175958/150‑1 (MQ=255)
                                                                 aTGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTg                                                                                    <  2:1362313/150‑1 (MQ=255)
                                                                 aTGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTg                                                                                    >  1:405504/1‑150 (MQ=255)
                                                                  tGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATAACGGCAATAGAGCAGCGGAGAACGTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTgg                                                                                   >  2:1733425/1‑150 (MQ=255)
                                                                     gTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTa                                                                                >  2:3317877/1‑150 (MQ=255)
                                                                      tGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCTACAACGCGTGGTGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTaa                                                                               <  2:4929314/150‑1 (MQ=255)
                                                                      tGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTaa                                                                               <  2:1379640/150‑1 (MQ=255)
                                                                         cTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGa                                                                            >  1:1548933/1‑150 (MQ=255)
                                                                         cTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGa                                                                            <  1:2436109/150‑1 (MQ=255)
                                                                         cTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGa                                                                            <  2:1826145/150‑1 (MQ=255)
                                                                         cTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATACATACACTACTTGCATGGATGAGTAATTATTAATGTGGATAGAGTTTATTTTTGAGGTTGGCTAATGa                                                                            >  1:1179411/1‑150 (MQ=255)
                                                                          tGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGaa                                                                           >  1:4196456/1‑150 (MQ=255)
                                                                           gggCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGaaa                                                                          <  2:3213086/150‑1 (MQ=255)
                                                                            ggCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAAc                                                                         >  2:2216098/1‑150 (MQ=255)
                                                                             gcgcACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACg                                                                        <  2:1492481/150‑1 (MQ=255)
                                                                             gcgcACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACg                                                                        <  1:4903643/150‑1 (MQ=255)
                                                                             gcgcACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACg                                                                        <  2:112921/150‑1 (MQ=255)
                                                                              cgcACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGc                                                                       >  2:2066801/1‑150 (MQ=255)
                                                                                cacaGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCtt                                                                     <  1:3983930/150‑1 (MQ=255)
                                                                                cacaGTACGAAGGCGGTCGTCATCAACAACGCGGGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATCGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCtt                                                                     >  2:3925270/1‑150 (MQ=255)
                                                                                 acaGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCttt                                                                    <  2:4737686/150‑1 (MQ=255)
                                                                                  caGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTc                                                                   >  2:4492628/1‑150 (MQ=255)
                                                                                  caGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTc                                                                   >  1:1573700/1‑150 (MQ=255)
                                                                                   aGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTcc                                                                  <  2:4052660/150‑1 (MQ=255)
                                                                                    gTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCg                                                                 >  1:1001716/1‑150 (MQ=255)
                                                                                      aCGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAAGAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCt                                                               <  2:3073626/150‑1 (MQ=255)
                                                                                       cGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTa                                                              >  2:1152345/1‑150 (MQ=255)
                                                                                          aGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTAttt                                                           <  2:760571/150‑1 (MQ=255)
                                                                                            gCGGGCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCt                                                         <  2:573541/150‑1 (MQ=255)
                                                                                              ggTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCttt                                                       <  2:3171557/150‑1 (MQ=255)
                                                                                               gtcgtcATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCtttt                                                      <  1:2024933/150‑1 (MQ=255)
                                                                                               gtcgtcATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCtttt                                                      <  2:3400042/150‑1 (MQ=255)
                                                                                                tcgtcATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTa                                                     <  1:743527/150‑1 (MQ=255)
                                                                                                tcgtcATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTa                                                     <  1:3407292/150‑1 (MQ=255)
                                                                                                  gtcATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTAtt                                                   >  1:4995446/1‑150 (MQ=255)
                                                                                                   tcatcaACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTc                                                  >  2:5057372/1‑150 (MQ=255)
                                                                                                   tcatcaACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGAGTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATTGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTc                                                  <  1:1605045/150‑1 (MQ=255)
                                                                                                         acaacGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTg                                            >  1:825781/1‑150 (MQ=255)
                                                                                                         acaacGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTg                                            >  1:657096/1‑150 (MQ=255)
                                                                                                           aacGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCt                                          <  1:911171/150‑1 (MQ=255)
                                                                                                           aacGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCt                                          <  1:3385985/150‑1 (MQ=255)
                                                                                                           aacGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCt                                          <  1:3871563/150‑1 (MQ=255)
                                                                                                            acGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTc                                         <  1:2164527/150‑1 (MQ=255)
                                                                                                             cgcgTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCa                                        >  1:637126/1‑150 (MQ=255)
                                                                                                              gcgTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCAc                                       <  1:3257402/150‑1 (MQ=255)
                                                                                                                gTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGc                                     >  1:4535687/1‑150 (MQ=255)
                                                                                                                          ttACGGCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGt                           >  1:2314220/1‑150 (MQ=255)
                                                                                                                             cGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTcc                        <  1:2805656/150‑1 (MQ=255)
                                                                                                                              ggCAAGAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCg                       <  2:2854923/150‑1 (MQ=255)
                                                                                                                               gCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCgg                      >  2:767015/1‑150 (MQ=255)
                                                                                                                               gCAATAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCgg                      >  2:5114548/1‑150 (MQ=255)
                                                                                                                                   tAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCa                  <  2:4664991/150‑1 (MQ=255)
                                                                                                                                   tAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCa                  <  1:1395080/150‑1 (MQ=255)
                                                                                                                                   tAGAGCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCa                  <  2:2469623/150‑1 (MQ=255)
                                                                                                                                    agagCAGCGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAg                 <  2:1063500/150‑1 (MQ=255)
                                                                                                                                       gcagcGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGcc              >  1:359441/1‑150 (MQ=255)
                                                                                                                                       gcagcGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGcc              <  2:1553696/150‑1 (MQ=255)
                                                                                                                                         agcGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGa            <  1:782203/150‑1 (MQ=255)
                                                                                                                                         agcGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTGTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGa            >  1:1946814/1‑150 (MQ=255)
                                                                                                                                          gcGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGat           <  1:4499832/150‑1 (MQ=255)
                                                                                                                                          gcGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGat           <  1:1074456/150‑1 (MQ=255)
                                                                                                                                           cGGAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGata          <  2:732786/150‑1 (MQ=255)
                                                                                                                                            ggAGAAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGatat         >  1:954216/1‑150 (MQ=255)
                                                                                                                                             gagaAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGAtatt        <  1:4904063/150‑1 (MQ=255)
                                                                                                                                              agaAATAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGAtatta       <  2:2317749/150‑1 (MQ=255)
                                                                                                                                                aaaTAGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATc     >  2:4077373/1‑150 (MQ=255)
                                                                                                                                                   ttGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATcgac  <  1:2277091/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CGACACCAACGTGCTGGCTTTTGGTTCACGAGTGGTTGGCCTCGAACTGGCAAAAATGATTGTGGATGCGTGGCTGGGCGCACAGTACGAAGGCGGTCGTCATCAACAACGCGTGGAGGCGATTACGGCAATAGAGCAGCGGAGAAATTGAGATTCATCCACTACTTGCATGGATGAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCACCGTTCCGGCCCAGGCCGATATTATCGAC  >  CP009273/4303466‑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: