Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,012,953 C→A 58.8% R247S (CGT→AGT)  ubiE → bifunctional 2‑octaprenyl‑6‑methoxy‑1,4‑benzoquinone methylase/ S‑adenosylmethionine:2‑DMK methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,012,9530CA58.8% 53.0 / inf 322R247S (CGT→AGT) ubiEbifunctional 2‑octaprenyl‑6‑methoxy‑1,4‑benzoquinone methylase/ S‑adenosylmethionine:2‑DMK methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (65/67);  new base A (103/84);  total (168/154)
Fisher's exact test for biased strand distribution p-value = 3.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.33e-01

Only 100 of 325 total aligned reads displayed.

TGGTCGCGAACGACGCCGACAGCTACCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTCTGCTGGGTAAAGTATTGCGCGTG  >  CP009273/4012804‑4013100
                                                                                                                                                     |                                                                                                                                                   
 ggTCGCGAACGACGCCGACAGCTACCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCg                                                                                                                                                    >  2:1303507/1‑150 (MQ=255)
       gAACGACGCCGACAGCTACCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTta                                                                                                                                              >  1:708621/1‑150 (MQ=255)
       gAACGACGCCGACAGCTACCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTta                                                                                                                                              >  1:3907664/1‑150 (MQ=255)
       gAACGACGCCGACAGCTACCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTta                                                                                                                                              >  1:3262612/1‑150 (MQ=255)
        aaCGACGCCGACAGCTACCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTtat                                                                                                                                             >  1:4151797/1‑150 (MQ=255)
                cGACAGCTACCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTg                                                                                                                                     >  2:2661770/1‑150 (MQ=255)
                 gACAGCTACCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGa                                                                                                                                    >  1:3129576/1‑150 (MQ=255)
                     gCTACCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAgg                                                                                                                                <  1:4825450/150‑1 (MQ=255)
                      cTACCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGga                                                                                                                               >  2:108465/1‑150 (MQ=255)
                        aCCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCTGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGgaga                                                                                                                             <  2:631627/150‑1 (MQ=255)
                        aCCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGgaga                                                                                                                             <  1:4651941/150‑1 (MQ=255)
                        aCCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGgaga                                                                                                                             <  2:4071538/150‑1 (MQ=255)
                        aCCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGgaga                                                                                                                             <  2:1728462/150‑1 (MQ=255)
                          cGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGAcc                                                                                                                           >  1:3941786/1‑150 (MQ=255)
                          cGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGAcc                                                                                                                           >  2:1571586/1‑150 (MQ=255)
                           gTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCg                                                                                                                          >  2:4806938/1‑150 (MQ=255)
                             tATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGa                                                                                                                        >  1:4918393/1‑150 (MQ=255)
                             tATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGa                                                                                                                        >  1:553416/1‑150 (MQ=255)
                                cTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAAt                                                                                                                     <  1:4547704/150‑1 (MQ=255)
                                cTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAAt                                                                                                                     <  2:1225231/150‑1 (MQ=255)
                                     aGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCtt                                                                                                                >  1:215925/1‑150 (MQ=255)
                                      gAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCttt                                                                                                               >  2:4734411/1‑150 (MQ=255)
                                        atccatccGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTa                                                                                                             >  1:4931990/1‑150 (MQ=255)
                                         tccatccGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTaa                                                                                                            >  2:2214138/1‑150 (MQ=255)
                                           catccGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAAc                                                                                                          >  2:30684/1‑150 (MQ=255)
                                           catccGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAAc                                                                                                          >  2:4603505/1‑150 (MQ=255)
                                           catccGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAAc                                                                                                          <  1:2414179/150‑1 (MQ=255)
                                            atccGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAAcc                                                                                                         >  1:3656797/1‑150 (MQ=255)
                                              ccGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCtt                                                                                                       <  2:4473899/150‑1 (MQ=255)
                                               cGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCttt                                                                                                      <  2:4355436/150‑1 (MQ=255)
                                               cGTATGCATCCAGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGTCTGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGAAAGGAGACCGGAAATGCCTTTTAAACCttt                                                                                                      <  1:614740/150‑1 (MQ=255)
                                                gTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTa                                                                                                     <  1:889709/150‑1 (MQ=255)
                                                gTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTa                                                                                                     <  2:1620427/150‑1 (MQ=255)
                                                 tATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAg                                                                                                    <  2:514152/150‑1 (MQ=255)
                                                 tATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAg                                                                                                    <  2:2853069/150‑1 (MQ=255)
                                                  aTGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGt                                                                                                   <  2:5307636/150‑1 (MQ=255)
                                                  aTGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGt                                                                                                   <  2:2173908/150‑1 (MQ=255)
                                                  aTGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGt                                                                                                   >  2:1891199/1‑150 (MQ=255)
                                                  aTGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGt                                                                                                   >  2:823705/1‑150 (MQ=255)
                                                     cATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGAc                                                                                                <  1:804732/150‑1 (MQ=255)
                                                     cATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGAc                                                                                                >  2:4142151/1‑150 (MQ=255)
                                                      aTCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACg                                                                                               >  2:1026085/1‑150 (MQ=255)
                                                      aTCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGGTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACg                                                                                               >  2:750958/1‑150 (MQ=255)
                                                           gATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAgg                                                                                          >  1:1010536/1‑150 (MQ=255)
                                                           gATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAgg                                                                                          >  1:1449644/1‑150 (MQ=255)
                                                           gATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAgg                                                                                          >  1:2002533/1‑150 (MQ=255)
                                                              cAGGATACCCAGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAAt                                                                                       >  1:1122237/1‑150 (MQ=255)
                                                               aGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAAtt                                                                                      <  1:4331047/150‑1 (MQ=255)
                                                               aGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGAATACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCGTTAGTGACGGCAGGAAtt                                                                                      >  1:1190266/1‑150 (MQ=255)
                                                                ggATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTg                                                                                     >  2:3327105/1‑150 (MQ=255)
                                                                 gATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGACTTGa                                                                                    >  2:2767715/1‑150 (MQ=255)
                                                                 gATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGa                                                                                    <  1:2530391/150‑1 (MQ=255)
                                                                    aCCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAg                                                                                 <  2:3426868/150‑1 (MQ=255)
                                                                     cccTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGt                                                                                <  2:581820/150‑1 (MQ=255)
                                                                     cccTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGt                                                                                >  1:1709898/1‑150 (MQ=255)
                                                                       cTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCt                                                                              <  2:4896088/150‑1 (MQ=255)
                                                                          aaaGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCt                                                                           <  2:378709/150‑1 (MQ=255)
                                                                          aaaGCCAAGATGCAGGATGCCGGATTCGAAAGTGTCGACAACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTAAAACCATTAGTGACGGCAGGAATTGAAAGTCTGCt                                                                           >  2:1782269/1‑150 (MQ=255)
                                                                              ccATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAac                                                                       >  2:4212203/1‑150 (MQ=255)
                                                                               cATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAaca                                                                      >  1:268422/1‑150 (MQ=255)
                                                                                  gatgCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCt                                                                   <  2:2303796/150‑1 (MQ=255)
                                                                                    tgCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTc                                                                 >  1:2093653/1‑150 (MQ=255)
                                                                                    tgCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTc                                                                 <  2:1429869/150‑1 (MQ=255)
                                                                                       aGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTg                                                              >  1:1521910/1‑150 (MQ=255)
                                                                                          aTGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTAt                                                           <  1:4785100/150‑1 (MQ=255)
                                                                                          aTGCCGGATTAGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTAt                                                           <  1:4514562/150‑1 (MQ=255)
                                                                                           tGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATc                                                          >  1:3571090/1‑150 (MQ=255)
                                                                                              cGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCt                                                       <  1:4446144/150‑1 (MQ=255)
                                                                                              cGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCt                                                       >  2:2490569/1‑150 (MQ=255)
                                                                                               ggATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTc                                                      >  1:186930/1‑150 (MQ=255)
                                                                                                 aTTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCAc                                                    <  2:3031141/150‑1 (MQ=255)
                                                                                                 aTTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCAc                                                    <  2:2057878/150‑1 (MQ=255)
                                                                                                      aaaGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCcgcg                                               >  1:2981968/1‑150 (MQ=255)
                                                                                                      aaaGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCcgcg                                               >  1:1102013/1‑150 (MQ=255)
                                                                                                        aGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCt                                             >  1:4041750/1‑150 (MQ=255)
                                                                                                         gtgtCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTg                                            >  2:822999/1‑150 (MQ=255)
                                                                                                         gtgtCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTg                                            <  1:2337402/150‑1 (MQ=255)
                                                                                                             cGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGaaaa                                        >  1:3757446/1‑150 (MQ=255)
                                                                                                              gACTACTACACTCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAAc                                       <  2:5057011/150‑1 (MQ=255)
                                                                                                               actactACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACg                                      <  2:4195418/150‑1 (MQ=255)
                                                                                                                  actacAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGcc                                   <  1:3062998/150‑1 (MQ=255)
                                                                                                                   ctacAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGccc                                  >  2:4866422/1‑150 (MQ=255)
                                                                                                                   ctacAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGccc                                  >  1:5212477/1‑150 (MQ=255)
                                                                                                                   ctacAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGccc                                  >  2:4275485/1‑150 (MQ=255)
                                                                                                                   ctacAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGccc                                  >  2:4495672/1‑150 (MQ=255)
                                                                                                                    tacAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAACACGGCCCg                                 <  2:3907682/150‑1 (MQ=255)
                                                                                                                      cAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAATGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCt                               <  1:413894/150‑1 (MQ=255)
                                                                                                                       aaTCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTc                              >  1:1308470/1‑150 (MQ=255)
                                                                                                                       aaTCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTc                              >  1:2850622/1‑150 (MQ=255)
                                                                                                                         tCTGACGGCTGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTcgc                            <  2:3718695/150‑1 (MQ=255)
                                                                                                                              cGGCTGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCATTTTAAACCTTTTGTGACGGCAGTAATTTAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTctg                       <  1:4638745/150‑1 (MQ=255)
                                                                                                                                gctggggTTGTGGCGCTGAATCGGGGTTATCTGTTCTGACTGGAGACCGGAATTGCCTTTTAAACCTTTTGTGCGGGCTGGAATTTAATGTCTGCTCAACACCTTCCTTTATCGCTCACCCGCGCTTAAAACGGCCCGCTCGCGTctgct                     <  1:512951/147‑1 (MQ=255)
                                                                                                                                gCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTctgct                     <  2:2597349/150‑1 (MQ=255)
                                                                                                                                  aGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGACTTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTCTGCTgg                   >  2:728633/1‑150 (MQ=255)
                                                                                                                                   ggggTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGTAAATGCCTTTTAAACCTTTCGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTCTGCTggg                  <  1:5201245/150‑1 (MQ=255)
                                                                                                                                       ttGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTCTGCTGGGTaaa              >  2:29473/1‑150 (MQ=255)
                                                                                                                                               cTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTCTGCTGGGTAAAGTATTgcg      >  1:4316855/1‑150 (MQ=255)
                                                                                                                                                tGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTCTGCTGGGTAAAGTATTgcgc     <  2:3123172/150‑1 (MQ=255)
                                                                                                                                                tGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTCTGCTGGGTAAAGTATTgcgc     >  2:2474174/1‑150 (MQ=255)
                                                                                                                                                   aTCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTCTGCTGGGTAAAGTATTGCGCGTg  >  2:3518036/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TGGTCGCGAACGACGCCGACAGCTACCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTCTGCTGGGTAAAGTATTGCGCGTG  >  CP009273/4012804‑4013100

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: