Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,012,953 C→A 92.9% 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,9530CA92.9% 564.6 / 31.7 214R247S (CGT→AGT) ubiEbifunctional 2‑octaprenyl‑6‑methoxy‑1,4‑benzoquinone methylase/ S‑adenosylmethionine:2‑DMK methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (10/6);  new base A (100/96);  total (110/104)
Fisher's exact test for biased strand distribution p-value = 4.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 215 total aligned reads displayed.

CGCGAACGACGCCGACAGCTACCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTCTGCTGGGTAAAGTATTGCGCGTG  >  CP009273/4012808‑4013100
                                                                                                                                                 |                                                                                                                                                   
cgcgAACGACGCCGACAGCTACCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGt                                                                                                                                                   <  2:884825/148‑1 (MQ=255)
cgcgAACGACGCCGACAGCTACCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGt                                                                                                                                                   <  1:1308951/148‑1 (MQ=255)
cgcgAACGACGCCGACAGCTACCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGt                                                                                                                                                   <  1:3056454/148‑1 (MQ=255)
 gcgAACGACGCCGACAGCTACCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTg                                                                                                                                                  <  2:2801571/148‑1 (MQ=255)
   gAACGACGCCGACAGCTACCGGTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGCTTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGt                                                                                                                                                <  2:2542655/148‑1 (MQ=255)
    aaCGACGCCGACAGCTACCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGtt                                                                                                                                               <  1:964977/148‑1 (MQ=255)
    aaCGACGCCGACAGCTACCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGtt                                                                                                                                               >  1:3276669/1‑148 (MQ=255)
     acgacgCCGACAGCTACCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCGGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAACGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCAGAGTGGTta                                                                                                                                              >  2:1915675/1‑148 (MQ=255)
     acgacgCCGACAGCTACCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTta                                                                                                                                              >  2:858228/1‑148 (MQ=255)
     acgacgCCGACAGCTACCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTta                                                                                                                                              <  2:2473450/148‑1 (MQ=255)
     acgacgCCGACAGCTACCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTta                                                                                                                                              >  2:1910845/1‑148 (MQ=255)
      cgacgCCGACAGCTACCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTtat                                                                                                                                             <  1:2005769/148‑1 (MQ=255)
         cgCCGACAGCTACCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAg                                                                                                                                          <  2:1912312/148‑1 (MQ=255)
          gCCGACAGCTACCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGt                                                                                                                                         <  1:1717727/148‑1 (MQ=255)
           ccGACAGCTACCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGtt                                                                                                                                        <  2:662203/148‑1 (MQ=255)
            cGACAGCTACCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTc                                                                                                                                       <  1:2846139/148‑1 (MQ=255)
             gACAGCTACCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCt                                                                                                                                      >  1:572217/1‑148 (MQ=255)
              aCAGCTACCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTg                                                                                                                                     >  2:523603/1‑148 (MQ=255)
                 gCTACCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACa                                                                                                                                  >  1:2739295/1‑148 (MQ=255)
                  cTACCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAg                                                                                                                                 >  2:173714/1‑148 (MQ=255)
                      cGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGgaga                                                                                                                             >  2:3412937/1‑148 (MQ=255)
                        ttATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGAcc                                                                                                                           <  1:1367564/148‑1 (MQ=255)
                          aTCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCgg                                                                                                                         <  1:2243802/148‑1 (MQ=255)
                            cTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGaa                                                                                                                       <  1:614544/148‑1 (MQ=255)
                             tGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGaaa                                                                                                                      <  1:447101/148‑1 (MQ=255)
                               gCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATg                                                                                                                    >  1:3334702/1‑148 (MQ=255)
                                  gAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCt                                                                                                                 <  1:847699/148‑1 (MQ=255)
                                  gAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCt                                                                                                                 <  1:2652353/148‑1 (MQ=255)
                                  gAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCt                                                                                                                 >  1:1686889/1‑148 (MQ=255)
                                   aaTCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCtt                                                                                                                >  2:1855120/1‑148 (MQ=255)
                                    atccatccGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTCCAATCTGACGGCAGGGGTGGTGGCGCGGCATCGTGGGGCTTAGTTCTGACAGGAGACCGGAAGTGCGttt                                                                                                               >  2:2724245/1‑148 (MQ=255)
                                      ccatccGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTa                                                                                                             >  2:2852905/1‑148 (MQ=255)
                                       catccGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTaa                                                                                                            <  2:1216042/148‑1 (MQ=255)
                                        atccGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTaaa                                                                                                           >  2:2352263/1‑148 (MQ=255)
                                        atccGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTaaa                                                                                                           >  2:3185560/1‑148 (MQ=255)
                                          ccGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAAcc                                                                                                         >  1:92734/1‑148 (MQ=255)
                                          ccGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAAcc                                                                                                         >  1:1549711/1‑148 (MQ=255)
                                             tATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCttt                                                                                                      >  1:157042/1‑148 (MQ=255)
                                              aTGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTa                                                                                                     >  1:113560/1‑148 (MQ=255)
                                              aTGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACATTTa                                                                                                     >  2:2031325/1‑148 (MQ=255)
                                                 cATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTg                                                                                                  <  1:1158027/148‑1 (MQ=255)
                                                   tCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGAc                                                                                                >  1:3049889/1‑148 (MQ=255)
                                                      cGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGc                                                                                             <  1:1369930/148‑1 (MQ=255)
                                                       gATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCa                                                                                            <  2:1871985/148‑1 (MQ=255)
                                                       gATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCa                                                                                            <  1:247045/148‑1 (MQ=255)
                                                            ggATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAAt                                                                                       >  2:962560/1‑148 (MQ=255)
                                                            ggATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAAt                                                                                       >  1:1957059/1‑148 (MQ=255)
                                                             gATACCCTGAAAGCGTTGATGAAGGATGCGGGATTCGAAAGTGTCGAGTAATAGTCTCTGACGTCAGGGGGTGTGGCGCTGATTAGTGGTTATAAGTTCTGTCAGTAGTCCGGAAATGCCATTTAATCCTTTAGTGACGGCAGGAAtt                                                                                      <  2:55538/148‑1 (MQ=255)
                                                             gATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAAtt                                                                                      <  2:741253/148‑1 (MQ=255)
                                                                aCCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGaa                                                                                   >  1:2199473/1‑148 (MQ=255)
                                                                 cccTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGaaa                                                                                  >  1:3200897/1‑148 (MQ=255)
                                                                      aaaGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTg                                                                             >  2:241658/1‑148 (MQ=255)
                                                                      aaaGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTg                                                                             <  1:2337057/148‑1 (MQ=255)
                                                                      aaaGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTg                                                                             <  2:1420789/148‑1 (MQ=255)
                                                                        aGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCt                                                                           >  2:1549053/1‑148 (MQ=255)
                                                                            atgatgCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAac                                                                       <  1:2539866/148‑1 (MQ=255)
                                                                            atgatgCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAac                                                                       >  2:2283764/1‑148 (MQ=255)
                                                                            atgatgCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAac                                                                       >  2:2192924/1‑148 (MQ=255)
                                                                              gatgCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAacac                                                                     <  1:1738679/148‑1 (MQ=255)
                                                                               atgCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACAcc                                                                    <  2:1580019/148‑1 (MQ=255)
                                                                                 gCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCtt                                                                  >  2:3340357/1‑148 (MQ=255)
                                                                                    ggATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCt                                                               >  1:2719419/1‑148 (MQ=255)
                                                                                             aTTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTc                                                      >  1:3082284/1‑148 (MQ=255)
                                                                                             aTTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTc                                                      >  1:2852612/1‑148 (MQ=255)
                                                                                             aTTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTc                                                      >  2:823745/1‑148 (MQ=255)
                                                                                             aTTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTc                                                      >  1:1166484/1‑148 (MQ=255)
                                                                                                     gtgtCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCgcgc                                              >  1:3425482/1‑148 (MQ=255)
                                                                                                      tgtCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCt                                             >  1:740139/1‑148 (MQ=255)
                                                                                                      tgtCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCt                                             >  1:739270/1‑148 (MQ=255)
                                                                                                      tgtCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCt                                             <  2:601141/148‑1 (MQ=255)
                                                                                                        tCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGa                                           >  2:404222/1‑148 (MQ=255)
                                                                                                          gacgacTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGaaa                                         <  1:1059485/148‑1 (MQ=255)
                                                                                                          gACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGaaa                                         <  2:1963079/148‑1 (MQ=255)
                                                                                                          gACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGaaa                                         <  2:1162059/148‑1 (MQ=255)
                                                                                                           actactACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGaaaa                                        <  2:1474687/148‑1 (MQ=255)
                                                                                                           actactACAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGaaaa                                        >  1:2168897/1‑148 (MQ=255)
                                                                                                            ctactaCAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAAc                                       >  1:3141851/1‑148 (MQ=255)
                                                                                                             tactacAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACg                                      <  1:472401/148‑1 (MQ=255)
                                                                                                                 acAATCTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGccc                                  >  1:1427545/1‑148 (MQ=255)
                                                                                                                    aTCTGAGGGCTGGGTTTGTTGCGCTGCATAGTGGTTATAATTTCTGACAGTAGACCGGAAATGCCTTTTAAACCTTTATTGACGGCAGGAATTGAAATTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCt                               <  2:2852543/148‑1 (MQ=255)
                                                                                                                      cTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTcg                             >  2:2816701/1‑148 (MQ=255)
                                                                                                                      cTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTcg                             >  2:2969115/1‑148 (MQ=255)
                                                                                                                      cTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTcg                             >  2:2969638/1‑148 (MQ=255)
                                                                                                                      cTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTcg                             >  2:3271508/1‑148 (MQ=255)
                                                                                                                      cTGACGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTcg                             >  1:2098096/1‑148 (MQ=255)
                                                                                                                         aCGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGt                          <  2:3494415/148‑1 (MQ=255)
                                                                                                                          cGGCAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTc                         >  1:2564577/1‑148 (MQ=255)
                                                                                                                             cAGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCGGAAAACGGCCCGCTCGCGTctgc                      <  2:1714642/148‑1 (MQ=255)
                                                                                                                              aGGGGTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTctgct                     <  2:2342929/148‑1 (MQ=255)
                                                                                                                               ggggTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTctgctg                    >  2:1408048/1‑148 (MQ=255)
                                                                                                                                  gTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTCTGCTGGGt                 >  1:589543/1‑148 (MQ=255)
                                                                                                                                  gTTGTGGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTCTGCTGGGt                 >  1:2352297/1‑148 (MQ=255)
                                                                                                                                   ttGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTCTGCTGGGTa                >  1:302230/1‑148 (MQ=255)
                                                                                                                                    tgtgGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAATTGCCTTTTAAACCTTTAGTGACGGCTGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTCTGCTGGGTaa               <  2:2739295/148‑1 (MQ=255)
                                                                                                                                    tgtgGCGCTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTCTGCTGGGTaa               >  2:2146841/1‑148 (MQ=255)
                                                                                                                                        gcgcTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGACTTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTCTGCTGGGGAAAgga           >  2:3279605/1‑146 (MQ=255)
                                                                                                                                        gcgcTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTCTGCTGGGTAAAGTa           <  2:2905776/148‑1 (MQ=255)
                                                                                                                                          gcTGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTCTGCTGGGTAAAGTAtt         <  2:3049889/148‑1 (MQ=255)
                                                                                                                                            tGCATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTCTGCTGGGTAAAGTATTgc       >  2:800126/1‑148 (MQ=255)
                                                                                                                                              cATAGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTCTGCTGGGTAAAGTCTTgcgc     >  2:1262349/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CGCGAACGACGCCGACAGCTACCGTTATCTGGCAGAATCCATCCGTATGCATCCCGATCAGGATACCCTGAAAGCCATGATGCAGGATGCCGGATTCGAAAGTGTCGACTACTACAATCTGACGGCAGGGGTTGTGGCGCTGCATCGTGGTTATAAGTTCTGACAGGAGACCGGAAATGCCTTTTAAACCTTTAGTGACGGCAGGAATTGAAAGTCTGCTCAACACCTTCCTGTATCGCTCACCCGCGCTGAAAACGGCCCGCTCGCGTCTGCTGGGTAAAGTATTGCGCGTG  >  CP009273/4012808‑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: