Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,191,496 C→A 28.8% pseudogene (1544/2847 nt) yehH → DUF4132 family pseudogene; defective yehI paralog; molybdate metabolism regulator, first fragment; molybdate metabolism regulator, second fragment 2; yehI' molR mutations have phenotypes and the N‑terminal molR fragment of yehI' may be expressed and may have retained or evolved MolR function

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,191,4960CA28.8% 642.0 / 242.2 314pseudogene (1544/2847 nt)yehHDUF4132 family pseudogene; defective yehI paralog; molybdate metabolism regulator, first fragment; molybdate metabolism regulator, second fragment 2; yehI' molR mutations have phenotypes and the N‑terminal molR fragment of yehI' may be expressed and may have retained or evolved MolR function
Reads supporting (aligned to +/- strand):  ref base C (105/118);  new base A (48/43);  total (153/161)
Fisher's exact test for biased strand distribution p-value = 3.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 317 total aligned reads displayed.

ACAACCAGCACTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGCCATGAAGAAAGT  >  CP009273/2191347‑2191644
                                                                                                                                                     |                                                                                                                                                    
acaACCAGCACTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTc                                                                                                                                                      >  2:4125755/1‑150 (MQ=255)
acaACCAGCACTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTc                                                                                                                                                      <  1:2025444/150‑1 (MQ=255)
    ccAGCACTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAgaaa                                                                                                                                                  >  1:2217137/1‑150 (MQ=255)
      aGCACTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCgaaaaa                                                                                                                                                >  2:2165496/1‑150 (MQ=255)
       gCACTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCgaaaaag                                                                                                                                               <  2:426587/150‑1 (MQ=255)
            ggCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCgaaaaagaaaaa                                                                                                                                          <  1:2512355/150‑1 (MQ=255)
            ggCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAgaaaaagaaaaa                                                                                                                                          >  1:1015882/1‑150 (MQ=255)
               aGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATc                                                                                                                                       >  1:1606299/1‑150 (MQ=255)
               aGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATc                                                                                                                                       <  1:2498403/150‑1 (MQ=255)
                  aCAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGcc                                                                                                                                    <  1:1958709/150‑1 (MQ=255)
                    aGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGa                                                                                                                                  >  1:2340554/1‑150 (MQ=255)
                     ggTCATTCCCTGGCTCTCGACACCCGCTGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTTGTTGTCTCCCCTCCCGGGCTTTCGAAAAAGAAAAAATCGCCGAt                                                                                                                                 <  2:3834811/150‑1 (MQ=255)
                       tCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTc                                                                                                                               <  2:3582215/150‑1 (MQ=255)
                        cagtCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAACCCAGGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTcc                                                                                                                              <  1:2911455/147‑1 (MQ=255)
                         agtCCCGGGCTCTCGACACCCGCAGGTGCCGTACGGAAATCATGCCCGCAACAACTGACACAGCCCTCACCCCCTGCCAGCGCCGATCTACGGCCAGCCGTAGTAGGCGCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCg                                                                                                                             <  2:3322974/148‑1 (MQ=255)
                         aTTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCg                                                                                                                             >  1:1359685/1‑150 (MQ=255)
                         aTTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCg                                                                                                                             >  2:2421444/1‑150 (MQ=255)
                            cccTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTg                                                                                                                          <  2:4047648/150‑1 (MQ=255)
                             ccTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGc                                                                                                                         <  2:1232376/150‑1 (MQ=255)
                             ccTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGc                                                                                                                         >  2:25981/1‑150 (MQ=255)
                              cTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCt                                                                                                                        <  2:4884461/150‑1 (MQ=255)
                                 gCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACCGACACCGCCCTCAACCCAGGCCAGCGCCGAACTACTGCCAGCCGTCGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGa                                                                                                                     <  1:2835054/150‑1 (MQ=255)
                                   tctcGACACCCGCTGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACCGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGAtt                                                                                                                   <  1:2664539/150‑1 (MQ=255)
                                   tctcGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCCGCCGGAGTAGTCTCCCCTCCCTGGCGTTCGAACAAGAAAAAATCGCCGATTCCGGGGCTGGAtg                                                                                                                   >  2:3454460/1‑149 (MQ=255)
                                   tctcGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGAtt                                                                                                                   >  2:2337249/1‑150 (MQ=255)
                                      cGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCTTGCCAGCTCCTATCTACTGCCAGCCGTAGTATTCTCCCCTCCCTGGATTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAg                                                                                                                <  2:1314140/150‑1 (MQ=255)
                                       gACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAgc                                                                                                               >  2:4926911/1‑150 (MQ=255)
                                         cacCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAgcgc                                                                                                             >  1:1692338/1‑150 (MQ=255)
                                             cGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCAtt                                                                                                         >  2:4383207/1‑150 (MQ=255)
                                             cGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCAGGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCAtt                                                                                                         <  1:3559640/150‑1 (MQ=255)
                                              gCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTa                                                                                                        <  2:5004296/150‑1 (MQ=255)
                                                 gTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGc                                                                                                     >  2:1000978/1‑150 (MQ=255)
                                                  ttGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCa                                                                                                    <  2:1855825/150‑1 (MQ=255)
                                                      cGTACTGAAATCATGCCAGCAAAAACTCACACAGCCCTCAAACCATCCCAGCGCCGATCTACTGCCAGCCGTTGTAGTCACCCCTCCCTGGCTTTCCAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGa                                                                                                <  1:3911409/150‑1 (MQ=255)
                                                           tGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCTATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCaa                                                                                           <  1:51621/150‑1 (MQ=255)
                                                           tGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCaa                                                                                           >  1:2257330/1‑150 (MQ=255)
                                                            gAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAAt                                                                                          >  1:4495312/1‑150 (MQ=255)
                                                               aTCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTg                                                                                       >  1:1907806/1‑150 (MQ=255)
                                                                tCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGt                                                                                      >  2:3225127/1‑150 (MQ=255)
                                                                tCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGt                                                                                      <  2:5062113/150‑1 (MQ=255)
                                                                  aTGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTa                                                                                    <  2:2371648/150‑1 (MQ=255)
                                                                     ccAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCt                                                                                 <  2:4452163/150‑1 (MQ=255)
                                                                     ccAGCAACAACTGACACAGCCCTCAAACCATGCCACCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCCCCCTGGCTTTCGAAAAATAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCt                                                                                 <  1:477929/150‑1 (MQ=255)
                                                                        gCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGAc                                                                              <  1:2091300/150‑1 (MQ=255)
                                                                         caacaaCTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACa                                                                             <  1:154162/150‑1 (MQ=255)
                                                                          aacaacTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACag                                                                            >  1:2137731/1‑150 (MQ=255)
                                                                              acTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACagaaga                                                                        >  1:1245701/1‑150 (MQ=255)
                                                                               cTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAgaagaa                                                                       <  1:3636273/150‑1 (MQ=255)
                                                                                tGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGaaa                                                                      <  1:834511/150‑1 (MQ=255)
                                                                                tGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGaaa                                                                      >  2:3284676/1‑150 (MQ=255)
                                                                                tGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGaaa                                                                      >  2:2879886/1‑150 (MQ=255)
                                                                                tGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGaaa                                                                      >  2:819670/1‑150 (MQ=255)
                                                                                  acacAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATc                                                                    <  2:3398678/150‑1 (MQ=255)
                                                                                   cacaGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCa                                                                   <  2:372951/150‑1 (MQ=255)
                                                                                    acaGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAg                                                                  >  1:4761048/1‑150 (MQ=255)
                                                                                    acaGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGGGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAg                                                                  >  2:2535434/1‑150 (MQ=255)
                                                                                     caGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGt                                                                 >  2:1412534/1‑150 (MQ=255)
                                                                                      aGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTa                                                                <  2:1772811/150‑1 (MQ=255)
                                                                                      aGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTa                                                                >  1:1286904/1‑150 (MQ=255)
                                                                                           tCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCaa                                                           <  2:2105835/150‑1 (MQ=255)
                                                                                           tCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCaa                                                           <  2:4659611/150‑1 (MQ=255)
                                                                                             aaaCCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACt                                                         >  2:755784/1‑150 (MQ=255)
                                                                                                   tGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGc                                                   <  2:1670221/150‑1 (MQ=255)
                                                                                                    gCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCg                                                  <  2:1015885/150‑1 (MQ=255)
                                                                                                     ccAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGa                                                 <  1:1088379/150‑1 (MQ=255)
                                                                                                     ccAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGa                                                 <  1:3661227/150‑1 (MQ=255)
                                                                                                         cgcCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAata                                             <  2:2959057/150‑1 (MQ=255)
                                                                                                          gcCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAatat                                            >  2:3810799/1‑150 (MQ=255)
                                                                                                          gcCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAatat                                            >  2:1238831/1‑150 (MQ=255)
                                                                                                            cGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAatatat                                          >  2:2732355/1‑150 (MQ=255)
                                                                                                            cGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAatatat                                          >  2:4777785/1‑150 (MQ=255)
                                                                                                             gATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATc                                         >  1:4449412/1‑150 (MQ=255)
                                                                                                             gATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATc                                         <  1:4869593/150‑1 (MQ=255)
                                                                                                                cTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTgg                                      <  2:2080304/150‑1 (MQ=255)
                                                                                                                 tACAGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGt                                     >  1:465921/1‑150 (MQ=255)
                                                                                                                   cTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTAt                                   >  1:4451510/1‑150 (MQ=255)
                                                                                                                    tGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTAtt                                  <  1:5020162/150‑1 (MQ=255)
                                                                                                                      ccAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAGATATATCTGGTATTCa                                >  2:3543977/1‑150 (MQ=255)
                                                                                                                       cAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCaa                               >  1:2233298/1‑150 (MQ=255)
                                                                                                                       cAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCaa                               <  2:2824490/150‑1 (MQ=255)
                                                                                                                         gCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCaaaa                             >  1:3446783/1‑150 (MQ=255)
                                                                                                                           cGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAaca                           >  1:3471241/1‑150 (MQ=255)
                                                                                                                           cGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAaca                           <  2:192227/150‑1 (MQ=255)
                                                                                                                            gtagtaGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAacac                          <  2:1483989/150‑1 (MQ=255)
                                                                                                                              agtagtCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACAt                        >  2:716725/1‑150 (MQ=255)
                                                                                                                               gtagtCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATc                       <  2:4423675/150‑1 (MQ=255)
                                                                                                                                 agtCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCAc                     >  2:2178196/1‑150 (MQ=255)
                                                                                                                                   tctcCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACgg                   <  2:3675584/150‑1 (MQ=255)
                                                                                                                                     tcccctcccTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGtt                 >  1:4054565/1‑150 (MQ=255)
                                                                                                                                      cccctcccTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTa                <  1:636552/150‑1 (MQ=255)
                                                                                                                                       ccctcccgGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAg               <  1:1157558/150‑1 (MQ=255)
                                                                                                                                        cctccctGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGc              <  2:870181/150‑1 (MQ=255)
                                                                                                                                          tccctGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGCCa            >  2:1269294/1‑150 (MQ=255)
                                                                                                                                           ccctGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGCCAt           >  2:4550148/1‑150 (MQ=255)
                                                                                                                                            cctGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGCCATg          <  1:4354156/150‑1 (MQ=255)
                                                                                                                                              tGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGGTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGCCATgaa        >  1:2985337/1‑150 (MQ=255)
                                                                                                                                                 cTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGCCATgaagaa     <  2:2760166/150‑1 (MQ=255)
                                                                                                                                                  tttCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGCCATGAAGaaa    >  1:865410/1‑150 (MQ=255)
                                                                                                                                                   ttCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGCCATGAAGAAAg   >  2:1095842/1‑150 (MQ=255)
                                                                                                                                                    tCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGCCATGAAGAAAGt  <  2:2503434/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ACAACCAGCACTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGCCATGAAGAAAGT  >  CP009273/2191347‑2191644

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: