Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,069,390 G→T 10.9% D160Y (GAT→TAT)  yihX → alpha‑D‑Glucose‑1‑P phosphatase, anomer‑specific

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,069,3900GT10.9% 919.6 / 60.5 256D160Y (GAT→TAT) yihXalpha‑D‑Glucose‑1‑P phosphatase, anomer‑specific
Reads supporting (aligned to +/- strand):  ref base G (112/116);  new base T (14/14);  total (126/130)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

Only 100 of 261 total aligned reads displayed.

CCTGCATACCACCTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTATTCTGGTGAAAGATAAAACCACCATCCCGGACTATTTCGCGAAGGTGTTATGCTAAAAACCATTCAGGACAAAGCCAGGCACC  >  CP009273/4069245‑4069539
                                                                                                                                                 |                                                                                                                                                     
ccTGCATACCACCTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATAc                                                                                                                                                   <  1:2364466/150‑1 (MQ=255)
   gCATACCACCTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGt                                                                                                                                                <  1:2627286/150‑1 (MQ=255)
    cATACCACCTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTc                                                                                                                                               <  1:3670476/150‑1 (MQ=255)
    cATACCACCTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTc                                                                                                                                               <  1:485024/150‑1 (MQ=255)
      tACCACCTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCtt                                                                                                                                             <  1:1617935/150‑1 (MQ=255)
       accaccTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCttt                                                                                                                                            >  1:1126079/1‑150 (MQ=255)
       accaccTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCttt                                                                                                                                            <  2:3174675/150‑1 (MQ=255)
       accaccTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCttt                                                                                                                                            <  1:1748360/150‑1 (MQ=255)
         caccTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCttttt                                                                                                                                          >  1:3889700/1‑150 (MQ=255)
          accTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTc                                                                                                                                         <  1:514871/150‑1 (MQ=255)
          accTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTc                                                                                                                                         <  2:1119629/150‑1 (MQ=255)
          accTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTc                                                                                                                                         <  2:1634832/150‑1 (MQ=255)
           ccTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTcg                                                                                                                                        <  1:1566532/150‑1 (MQ=255)
            cTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTcga                                                                                                                                       <  2:2298366/150‑1 (MQ=255)
             ttCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGTTCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACTGTCTTTTTcgac                                                                                                                                      <  2:3605812/150‑1 (MQ=255)
                tGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGAt                                                                                                                                   >  1:3714387/1‑150 (MQ=255)
                tGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGAt                                                                                                                                   <  2:1621775/150‑1 (MQ=255)
                tGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGAt                                                                                                                                   >  1:2879002/1‑150 (MQ=255)
                tGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGAt                                                                                                                                   >  2:3050427/1‑150 (MQ=255)
                   ccGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAAc                                                                                                                                >  1:3645623/1‑150 (MQ=255)
                    cGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACg                                                                                                                               >  1:3414384/1‑150 (MQ=255)
                        agaaTACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCTATACGGTCTTTTTCGACGATAACGCCGa                                                                                                                           >  1:633145/1‑150 (MQ=255)
                        agaaTACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGa                                                                                                                           >  2:2772462/1‑150 (MQ=255)
                          aaTACCCGGAATTTCGTGTTGCTGATGAACATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGAACGAATTTACCAGCATGTTTTGCAGGCGGAAGTTTTTTCACCCTGCGATACGGTCTTTTTCGACGATAACGCCGata                                                                                                                         <  1:1417772/150‑1 (MQ=255)
                          aaTACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGata                                                                                                                         >  1:796920/1‑150 (MQ=255)
                          aaTACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGata                                                                                                                         <  2:3425591/150‑1 (MQ=255)
                           aTACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGataa                                                                                                                        >  2:414498/1‑150 (MQ=255)
                            tACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGataat                                                                                                                       <  2:3132670/150‑1 (MQ=255)
                              cccGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAatat                                                                                                                     >  2:352105/1‑150 (MQ=255)
                                cGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAg                                                                                                                   >  1:2131293/1‑150 (MQ=255)
                                cGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAg                                                                                                                   <  2:1858300/150‑1 (MQ=255)
                                 ggAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGa                                                                                                                  >  2:1720312/1‑150 (MQ=255)
                                 ggAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGa                                                                                                                  >  1:2003144/1‑150 (MQ=255)
                                  gAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGaa                                                                                                                 >  2:3727819/1‑150 (MQ=255)
                                     aTTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGa                                                                                                              <  1:901721/150‑1 (MQ=255)
                                        cGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCTATACGTTCTTTTTCGACGTTAACGCCGATAATATAGAAGGAGcc                                                                                                           <  1:2155815/150‑1 (MQ=255)
                                        cGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGcc                                                                                                           >  2:406682/1‑150 (MQ=255)
                                         gTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCa                                                                                                          >  2:1302652/1‑150 (MQ=255)
                                         gTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCa                                                                                                          <  1:10417/150‑1 (MQ=255)
                                             tgctgcTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCTATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCa                                                                                                      >  2:2372333/1‑150 (MQ=255)
                                                tgctgACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCt                                                                                                   <  1:2728134/150‑1 (MQ=255)
                                                    gACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGc                                                                                               >  2:1774164/1‑150 (MQ=255)
                                                    gACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGc                                                                                               >  2:2264995/1‑150 (MQ=255)
                                                     aCCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCa                                                                                              >  1:472602/1‑150 (MQ=255)
                                                      ccATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCAt                                                                                             <  2:2774692/150‑1 (MQ=255)
                                                         tatctatcTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCTATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTAc                                                                                          >  1:2178783/1‑150 (MQ=255)
                                                         tatctatcTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTAc                                                                                          <  2:388981/150‑1 (MQ=255)
                                                          atctatctGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTAcc                                                                                         >  1:1180414/1‑150 (MQ=255)
                                                            ctatctGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAg                                                                                       >  1:1945249/1‑150 (MQ=255)
                                                            ctatctGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAg                                                                                       >  1:1772895/1‑150 (MQ=255)
                                                              atctGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTa                                                                                     >  2:3269968/1‑150 (MQ=255)
                                                                 tGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTATTc                                                                                  >  2:2740239/1‑150 (MQ=255)
                                                                  gTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTATTCt                                                                                 <  2:541921/150‑1 (MQ=255)
                                                                      cAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGAGTAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTATTCTGGTg                                                                             >  1:1469780/1‑150 (MQ=255)
                                                                         gATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTATTCTGGTGaaa                                                                          >  1:252552/1‑150 (MQ=255)
                                                                          aTCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCTATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTATTCTGGTGAAAg                                                                         <  1:3262542/150‑1 (MQ=255)
                                                                           tCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTATTCTGGTGAAAGa                                                                        <  1:993226/150‑1 (MQ=255)
                                                                              ggggATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCTATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTATTCTGGTGAAAGATaa                                                                     >  2:2442310/1‑150 (MQ=255)
                                                                              ggggATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTATTCTGGTGAAAGATaa                                                                     <  2:74214/150‑1 (MQ=255)
                                                                                ggATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTATTCTGGTGAAAGATaaaa                                                                   <  2:1211044/150‑1 (MQ=255)
                                                                                     cgcAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCNGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTATTCTGGNGAAAGATNAAaccacn                                                              <  2:1436516/150‑2 (MQ=255)
                                                                                     cgcAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTATTCTGGTGAAAGATAAAaccacc                                                              <  2:2562043/150‑1 (MQ=255)
                                                                                     cgcAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTATTCTGGTGAAAGATAAAaccacc                                                              >  2:3563647/1‑150 (MQ=255)
                                                                                        aaaCCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTATTCTGGTGAAAGATAAAACCACCATc                                                           >  1:3218662/1‑150 (MQ=255)
                                                                                          aCCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCTATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTATTCTGGTGAAAGATAAAACCACCATccc                                                         <  2:3317823/150‑1 (MQ=255)
                                                                                            cTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTATTCTGGTGAAAGATAAAACCACCATCCCgg                                                       >  2:268257/1‑150 (MQ=255)
                                                                                            cTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTATTCTGGTGAAAGATAAAACCACCATCCCgg                                                       <  1:728790/150‑1 (MQ=255)
                                                                                             tGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTATTCTGGTGAAAGATAAAACCACCATCCCGGa                                                      <  1:2396545/150‑1 (MQ=255)
                                                                                               aaGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTATTCTGGTGAAAGATAAAACCACCATCCCGGACt                                                    >  2:698403/1‑150 (MQ=255)
                                                                                                 gCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTATTCTGGTGAAAGATAAAACCACCATCCCGGACTAt                                                  <  2:3690617/150‑1 (MQ=255)
                                                                                                  cACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTATTCTGGTGAAAGATAAAACCACCATCCCGGACTAtt                                                 <  1:2608084/150‑1 (MQ=255)
                                                                                                    cGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTATTCTGGTGAAAGATAAAACCACCATCCCGGACTATTTc                                               >  2:1168431/1‑150 (MQ=255)
                                                                                                    cGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTATTCTGGTGAAAGATAAAACCACCATCCCGGACTATTTc                                               <  1:3386495/150‑1 (MQ=255)
                                                                                                     gAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTATTCTGGTGAAAGATAAAACCACCATCCCGGACTATTTcg                                              >  2:2061631/1‑150 (MQ=255)
                                                                                                     gAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTATTCTGGTGAAAGATAAAACCACCATCCCGGACTATTTcg                                              >  1:1265195/1‑150 (MQ=255)
                                                                                                       aTTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCTATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTATTCTGGTGAAAGATAAAACCACCATCCCGGACTATTTcgcg                                            <  2:2652721/150‑1 (MQ=255)
                                                                                                          tACCAGCCTGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTATTCTGGTGAAAGATAAAACCACCATCCCGGACTATTTCGCGAAg                                         <  2:1470079/150‑1 (MQ=255)
                                                                                                           aCCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTATTCTGGTGAAAGATAAAACCACCATCCCGGACTATTTCGCGAAgg                                        <  2:4098099/150‑1 (MQ=255)
                                                                                                              aGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTATTCTGGTGAAAGATAAAACCACCATCCCGGACTATTTCGCGAAGgtgt                                     >  1:4134636/1‑150 (MQ=255)
                                                                                                               gCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTATTCTGGTGAAAGATAAAACCACCATCCCGGACTATTTCGCGAAGGTGtt                                    <  2:1972349/150‑1 (MQ=255)
                                                                                                                cATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCTATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTATTCTGGTGAAAGATAAAACCACCATCCCGGACTATTTCGCGAAGGTGTTa                                   <  1:2804656/150‑1 (MQ=255)
                                                                                                                 aTGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTATTCTGGTGAAAGATAAAACCACCATCCCGGACTATTTCGCGAAGGTGTTAt                                  >  2:3800979/1‑150 (MQ=255)
                                                                                                                  tGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTATTCTGGTGAAAGATAAAACCACCATCCCGGACTATTTCGCGAAGGTGTTATg                                 >  2:864592/1‑150 (MQ=255)
                                                                                                                  tGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTATTCTGGTGAAAGATAAAACCACCATCCCGGACTATTTCGCGAAGGTGTTATg                                 >  2:2920711/1‑150 (MQ=255)
                                                                                                                    ttttGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTATTCTGGTGAAAGATAAAACCACCATCCCGGACTATTTCGCGAAGGTGTTATGCt                               >  1:1129455/1‑150 (MQ=255)
                                                                                                                    ttttGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATGACCAGTATTCTGGTGAAAGATAAAACCACCATCCCGGACTATTTCGCGAAGGTGTTATGCt                               <  2:1903021/150‑1 (MQ=255)
                                                                                                                         cAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTATTCTGGTGAAAGATAAAACCACCATCCCGGACTATTTCGCGAAGGTGTTATGCTaaaaa                          <  1:2091941/150‑1 (MQ=255)
                                                                                                                          aGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTATTCTGGTGAAAGATAAAACCACCATCCCGGACTATTTCGCGAAGGTGTTATGCTAAAAAc                         <  2:4086153/150‑1 (MQ=255)
                                                                                                                                aaGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTATTCTGGTGAAAGATAAAACCACCATCCCGGACTATTTCGCGAAGGTGTTATGCTAAAAACCATTCa                   >  1:3787175/1‑150 (MQ=255)
                                                                                                                                 aGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTATTCTGGTGAAAGATAAAACCACCATCCCGGACTATTTCGCGAAGGTGTTATGCTAAAAACCATTCAg                  >  1:1466816/1‑150 (MQ=255)
                                                                                                                                  ggTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTATTCTGGTGAAAGATAAAACCACCATCCCGGACTATTTCGCGAAGGTGTTATGCTAAAAACCATTCAgg                 <  2:2254682/150‑1 (MQ=255)
                                                                                                                                    tttttCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTATTCTGGTGAAAGATAAAACCACCATCCCGGACTATTTCGCGAAGGTGTTATGCTAAAAACCATTCAGGAc               >  1:3152879/1‑150 (MQ=255)
                                                                                                                                        tCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTATTCTGGTGAAAGATAAAACCACCATCCCGGACTATTTCGCGAAGGTGTTATGCTAAAAACCATTCAGGACAAAg           >  1:1078509/1‑150 (MQ=255)
                                                                                                                                         cACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTATTCTGGTGAAAGATAAAACCACCATCCCGGACTATTTCGCGAAGGTGTTATGCTAAAAACCATTCAGGACAAAGc          <  2:2100102/150‑1 (MQ=255)
                                                                                                                                          aCCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTATTCTGGTGAAAGATAAAACCACCATCCCGGACTATTTCGCGAAGGTGTTATGCTAAAAACCATTCAGGACAAAGcc         <  1:2064695/150‑1 (MQ=255)
                                                                                                                                           cccAGCTATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTATTCTGGTGAAAGATAAAACCACCATCCCGGACTATTTCGCGAAGGTGTTATGCTAAAAACCATTCAGGACAAAGCCa        <  1:2494453/150‑1 (MQ=255)
                                                                                                                                            ccAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTATTCTGGTGAAAGATAAAACCACCATCCCGGACTATTTCGCGAAGGTGTTATGCTAAAAACCATTCAGGACAAAGCCAg       >  2:774987/1‑150 (MQ=255)
                                                                                                                                              aGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTATTCTGGTGAAAGATAAAACCACCATCCCGGACTATTTCGCGAAGGTGTTATGCTAAAAACCATTCAGGACAAAGCCAGGc     >  2:3055510/1‑150 (MQ=255)
                                                                                                                                                 gATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTATTCTGGTGAAAGATAAAACCACCATCCCGGACTATTTCGCGAAGGTGTTATGCTAAAAACCATTCAGGACAAAGCCAGGCAcc  <  2:2498175/150‑1 (MQ=255)
                                                                                                                                                 gATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCGGGGCATTACCAGTATTCTGGTGAAAGATAAAACCACCATCCCGGACTATTTCGCGAAGGTGTTATGCTAAAAACCATTCAGGACAAAGCCAGGCAcc  >  1:2368095/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
CCTGCATACCACCTTCTGGCCGGAAGAATACCCGGAAATTCGTGATGCTGCTGACCATATCTATCTGTCGCAAGATCTGGGGATGCGCAAACCTGAAGCACGAATTTACCAGCATGTTTTGCAGGCGGAAGGTTTTTCACCCAGCGATACGGTCTTTTTCGACGATAACGCCGATAATATAGAAGGAGCCAATCAGCTGGGCATTACCAGTATTCTGGTGAAAGATAAAACCACCATCCCGGACTATTTCGCGAAGGTGTTATGCTAAAAACCATTCAGGACAAAGCCAGGCACC  >  CP009273/4069245‑4069539

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: