Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,012,953 C→A 73.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,9530CA73.8% 367.7 / 199.4 266R247S (CGT→AGT) ubiEbifunctional 2‑octaprenyl‑6‑methoxy‑1,4‑benzoquinone methylase/ S‑adenosylmethionine:2‑DMK methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (37/32);  new base A (92/103);  total (129/137)
Fisher's exact test for biased strand distribution p-value = 4.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 268 total aligned reads displayed.

TGGTCGCGAACGACGCCGACAGCTACCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTCTGCTGGGTAAAGTATTGCGCGTG  >  CP009273/4012804‑4013100
                                                                                                                                                     |                                                                                                                                                   
   tCGCGAACGACGCCGACAGCTACCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTg                                                                                                                                                  <  1:4641725/150‑1 (MQ=255)
     gcgAACGACGCCGACAGCTACCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGt                                                                                                                                                <  2:3757380/150‑1 (MQ=255)
          cgacgCCGACAGCTACCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATaa                                                                                                                                           <  2:1888259/150‑1 (MQ=255)
             cgCCGACAGCTACCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGtt                                                                                                                                        >  2:4427327/1‑150 (MQ=255)
               ccGACAGCTACCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCt                                                                                                                                      >  2:3055183/1‑150 (MQ=255)
               ccGACAGCTACCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCt                                                                                                                                      >  1:2385020/1‑150 (MQ=255)
                   cAGCTACCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACa                                                                                                                                  <  2:4251387/150‑1 (MQ=255)
                    aGCTACCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAg                                                                                                                                 >  1:1919512/1‑150 (MQ=255)
                     gCTACCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAgg                                                                                                                                <  2:4448837/150‑1 (MQ=255)
                       tACCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGgag                                                                                                                              >  2:971143/1‑150 (MQ=255)
                        aCCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGgaga                                                                                                                             >  2:4931483/1‑150 (MQ=255)
                          cGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGAcc                                                                                                                           <  1:1921194/150‑1 (MQ=255)
                          cGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGgagaaa                                                                                                                           <  2:4324818/150‑3 (MQ=255)
                          cGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGAcc                                                                                                                           <  1:3264568/150‑1 (MQ=255)
                             tATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGa                                                                                                                        <  2:5131978/150‑1 (MQ=255)
                               tCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGaaa                                                                                                                      <  1:541354/150‑1 (MQ=255)
                                 tGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATg                                                                                                                    >  2:4075050/1‑150 (MQ=255)
                                 tGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATg                                                                                                                    >  1:5120365/1‑150 (MQ=255)
                                 tGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATg                                                                                                                    >  1:1801912/1‑150 (MQ=255)
                                 tGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATg                                                                                                                    <  1:1047402/150‑1 (MQ=255)
                                  ggCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGc                                                                                                                   <  1:3397307/150‑1 (MQ=255)
                                    cAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCt                                                                                                                 >  2:403324/1‑150 (MQ=255)
                                    cAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCt                                                                                                                 <  1:1578182/150‑1 (MQ=255)
                                      gAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCttt                                                                                                               <  2:1745457/150‑1 (MQ=255)
                                          ccatccGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTaaa                                                                                                           >  2:3483284/1‑150 (MQ=255)
                                            atccGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAAcc                                                                                                         <  1:4726057/150‑1 (MQ=255)
                                                gTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTa                                                                                                     >  2:672909/1‑150 (MQ=255)
                                                  aTGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGt                                                                                                   <  2:1476444/150‑1 (MQ=255)
                                                  aTGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGt                                                                                                   <  2:2816434/150‑1 (MQ=255)
                                                     cATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGAc                                                                                                <  1:1991621/150‑1 (MQ=255)
                                                        cccGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGc                                                                                             <  2:3102192/150‑1 (MQ=255)
                                                        cccGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGc                                                                                             <  2:801866/150‑1 (MQ=255)
                                                        cccGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACGACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGc                                                                                             <  2:4677193/150‑1 (MQ=255)
                                                         ccGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCa                                                                                            >  2:5000958/1‑150 (MQ=255)
                                                         ccGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCa                                                                                            <  2:4752306/150‑1 (MQ=255)
                                                           gATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAgg                                                                                          >  2:1019151/1‑150 (MQ=255)
                                                            aTCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGa                                                                                         >  1:462312/1‑150 (MQ=255)
                                                                ggATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTg                                                                                     >  1:4840213/1‑150 (MQ=255)
                                                                  aTACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGaa                                                                                   <  1:2199412/150‑1 (MQ=255)
                                                                     cccTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGt                                                                                >  2:2056320/1‑150 (MQ=255)
                                                                     cccTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGt                                                                                <  2:1142678/150‑1 (MQ=255)
                                                                      ccTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTc                                                                               >  1:631552/1‑150 (MQ=255)
                                                                      ccTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTc                                                                               <  2:3653612/150‑1 (MQ=255)
                                                                       cTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCt                                                                              <  1:1172242/150‑1 (MQ=255)
                                                                        tGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTg                                                                             >  2:2750994/1‑150 (MQ=255)
                                                                          aaaGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCt                                                                           <  1:3417943/150‑1 (MQ=255)
                                                                             gCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCaa                                                                        <  2:632418/150‑1 (MQ=255)
                                                                              ccATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAac                                                                       >  1:3565863/1‑150 (MQ=255)
                                                                              ccATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAac                                                                       >  1:974035/1‑150 (MQ=255)
                                                                               cATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTTACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAaca                                                                      <  2:2552420/150‑1 (MQ=255)
                                                                               cATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAaca                                                                      >  1:1341781/1‑150 (MQ=255)
                                                                                atgatgCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAacac                                                                     >  2:2296117/1‑150 (MQ=255)
                                                                                    tgCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTc                                                                 >  2:2808302/1‑150 (MQ=255)
                                                                                     gCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTcc                                                                >  2:2521295/1‑150 (MQ=255)
                                                                                      cAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCt                                                               <  2:4147636/150‑1 (MQ=255)
                                                                                      cAGGATGCCGGATTCGAAAGTGTAGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCt                                                               <  2:502666/150‑1 (MQ=255)
                                                                                       aGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTg                                                              >  1:4125620/1‑150 (MQ=255)
                                                                                       aGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTg                                                              <  2:1623695/150‑1 (MQ=255)
                                                                                        ggATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGt                                                             <  1:1444165/150‑1 (MQ=255)
                                                                                        cgATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGt                                                             >  1:1838088/2‑150 (MQ=255)
                                                                                          aTGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTAt                                                           <  2:2099820/150‑1 (MQ=255)
                                                                                           tGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATc                                                          <  2:4498187/150‑1 (MQ=255)
                                                                                            gCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCg                                                         >  1:3357699/1‑150 (MQ=255)
                                                                                             ccGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCATGTATCGc                                                        >  1:4409201/1‑150 (MQ=255)
                                                                                              cGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCt                                                       >  2:4354629/1‑150 (MQ=255)
                                                                                               ggATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTc                                                      <  2:2757803/150‑1 (MQ=255)
                                                                                                 aTTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCAc                                                    >  1:4359082/1‑150 (MQ=255)
                                                                                                   tCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCAccc                                                  <  1:4681629/150‑1 (MQ=255)
                                                                                                    cGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGACTTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCcg                                                 >  2:3448775/1‑150 (MQ=255)
                                                                                                      aaaGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCcgcg                                               >  1:4445054/1‑150 (MQ=255)
                                                                                                        aGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCt                                             <  2:1533027/150‑1 (MQ=255)
                                                                                                           gtCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGaa                                          <  1:672908/150‑1 (MQ=255)
                                                                                                           gtCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGaa                                          <  2:1196382/150‑1 (MQ=255)
                                                                                                            tCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGGTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGaaa                                         <  2:1265303/150‑1 (MQ=255)
                                                                                                             cGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGaaaa                                        >  2:2842949/1‑150 (MQ=255)
                                                                                                                ctactaCAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACgg                                     >  1:4727543/1‑150 (MQ=255)
                                                                                                                       aaTCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTc                              >  2:395995/1‑150 (MQ=255)
                                                                                                                          cTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTcgcg                           <  1:1728170/150‑1 (MQ=255)
                                                                                                                          cTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACTGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTcgcg                           <  2:712349/150‑1 (MQ=255)
                                                                                                                           tGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAATACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGt                          <  2:847233/150‑1 (MQ=255)
                                                                                                                           tGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGt                          >  1:755554/1‑150 (MQ=255)
                                                                                                                            gACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTc                         >  2:3719792/1‑150 (MQ=255)
                                                                                                                            gACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTc                         <  1:1715465/150‑1 (MQ=255)
                                                                                                                             aCGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTct                        <  1:853808/150‑1 (MQ=255)
                                                                                                                              cGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTctg                       >  2:1695922/1‑150 (MQ=255)
                                                                                                                               ggCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTctgc                      >  1:295014/1‑150 (MQ=255)
                                                                                                                                gCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTctgct                     >  2:266552/1‑150 (MQ=255)
                                                                                                                                 cAGGGGTTGTGGCTCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAACCCTTTAGTGACGGCTGGAATTGAAAGTCTGCTCAACACCTTCCTGTCTCGCTCTCCCGCGCTGAAAACGGCCCGCTCGCGTctgctg                    <  1:3005150/150‑1 (MQ=255)
                                                                                                                                 cAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCAGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTctgctg                    <  1:3833838/150‑1 (MQ=255)
                                                                                                                                  aGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTCTGCTgg                   >  2:3248081/1‑150 (MQ=255)
                                                                                                                                           ggCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTCGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTCTGCTGGGTAAAGTAt          <  1:9839/150‑1 (MQ=255)
                                                                                                                                           ggCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTCTGCTGGGTAAAGTAt          >  1:2526108/1‑150 (MQ=255)
                                                                                                                                              gcTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTCTGCTGGGTAAAGTATTgc       >  2:3402763/1‑150 (MQ=255)
                                                                                                                                               cTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTCTGCTGGGTAAAGTATTgcg      >  2:356723/1‑150 (MQ=255)
                                                                                                                                               cTGAATCGTGGTTATAAGTTCTGACAGGAGACGGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTCTGCTGGGTAAAGTATTgcg      <  1:2431286/150‑1 (MQ=255)
                                                                                                                                                tGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTCTGCTGGGTAAAGTATTgcgc     >  2:2599324/1‑150 (MQ=255)
                                                                                                                                                 gCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTCTGCTGGGTAAAGTATTGcgcg    >  1:3806951/1‑150 (MQ=255)
                                                                                                                                                  cATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTCTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTCTGCTGGGTAAAGTATTGCGCGt   >  1:1220648/1‑150 (MQ=255)
                                                                                                                                                   agcGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTCTGCTGGGTAAAGTATTGCGCGTg  <  2:4589530/148‑1 (MQ=255)
                                                                                                                                                   aTCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTCTGCTGGGTAAAGTATTGCGCGTg  >  1:1279682/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: