Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,191,496 C→A 23.6% 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,4960CA23.6% 659.6 / 176.8 299pseudogene (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 (104/124);  new base A (35/36);  total (139/160)
Fisher's exact test for biased strand distribution p-value = 5.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.60e-01

Only 100 of 302 total aligned reads displayed.

ACAACCAGCACTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGCCATGAAGAAAGT  >  CP009273/2191347‑2191644
                                                                                                                                                     |                                                                                                                                                    
acaACCAGCACTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTc                                                                                                                                                      <  1:4785311/150‑1 (MQ=255)
acaACCAGCACTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTc                                                                                                                                                      <  1:1591398/150‑1 (MQ=255)
   aCCAGCACTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCgaa                                                                                                                                                   <  1:834679/150‑1 (MQ=255)
   aCCAGCACTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCgaa                                                                                                                                                   >  2:234276/1‑150 (MQ=255)
   aCCAGCACTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAgaa                                                                                                                                                   >  1:1377005/1‑150 (MQ=255)
    ccAGCACTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCgaaa                                                                                                                                                  <  2:2048096/150‑1 (MQ=255)
      aGCACTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCgaaaaa                                                                                                                                                <  2:3925903/150‑1 (MQ=255)
      aGCACTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAgaaaaa                                                                                                                                                <  2:177824/150‑1 (MQ=255)
       gCACTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCgaaaaag                                                                                                                                               <  2:2140238/150‑1 (MQ=255)
        cACTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCgaaaaaga                                                                                                                                              >  1:3294615/1‑150 (MQ=255)
        cACTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCgaaaaaga                                                                                                                                              >  2:4744231/1‑150 (MQ=255)
        cACTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAgaaaaaga                                                                                                                                              <  2:2956169/150‑1 (MQ=255)
         aCTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCgaaaaagaa                                                                                                                                             >  2:3680271/1‑150 (MQ=255)
          cTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCgaaaaagaaa                                                                                                                                            >  2:4271956/1‑150 (MQ=255)
            ggCAGAACGGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCgaaaaagaaaaa                                                                                                                                          >  2:3453465/1‑150 (MQ=255)
            ggCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCgaaaaagaaaaa                                                                                                                                          >  2:2535692/1‑150 (MQ=255)
            ggCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCgaaaaagaaaaa                                                                                                                                          <  1:4013055/150‑1 (MQ=255)
               aGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATc                                                                                                                                       <  1:2435095/150‑1 (MQ=255)
               aGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATc                                                                                                                                       <  2:2281165/150‑1 (MQ=255)
                gAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCg                                                                                                                                      <  1:290100/150‑1 (MQ=255)
                 aaCAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGc                                                                                                                                     >  2:1317095/1‑150 (MQ=255)
                  aCAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGcc                                                                                                                                    <  1:2865686/150‑1 (MQ=255)
                   cAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTTGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCg                                                                                                                                   <  2:168199/150‑1 (MQ=255)
                    aGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGa                                                                                                                                  <  1:2769350/150‑1 (MQ=255)
                        cATTCCCTGGCTCTCTACACCCGCAGTTGCCGTACTGAAATCATGCCTGAAACAACTTACACAGCCCTCAAACCATGCACGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAATAAAAAATCGCCGATTcc                                                                                                                              <  2:1783178/150‑1 (MQ=255)
                         aTTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCg                                                                                                                             >  2:558720/1‑150 (MQ=255)
                          ttCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCgg                                                                                                                            >  1:748353/1‑150 (MQ=255)
                           tCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGt                                                                                                                           >  2:3993045/1‑150 (MQ=255)
                            cccTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTg                                                                                                                          >  1:641793/1‑150 (MQ=255)
                             ccTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGc                                                                                                                         <  2:1830583/150‑1 (MQ=255)
                              cTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCt                                                                                                                        >  2:4339384/1‑150 (MQ=255)
                              cTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCt                                                                                                                        >  1:1873763/1‑150 (MQ=255)
                              cTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCt                                                                                                                        >  1:1333405/1‑150 (MQ=255)
                               tgtttctcTAAACCGGCTGTTGCCTTACTTAAATATTGCAAGCAACAAGTTACACAGCCCTAAAACCTTGCCAGCGCCTATATACTGCCAGCGTTTTTTGTCTCCCCTCCGTGGATTTATAAAAAGAAAAAATCGCCGATTCCGGTGCTg                                                                                                                       <  2:2814533/146‑1 (MQ=255)
                               tGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTg                                                                                                                       >  1:968806/1‑150 (MQ=255)
                               gggCTCTCGACACCCGCAGTTGCCTTACTTAAAACATGCCAGCAACAACTGACACAGCCCACAACCCCTGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCCAAAAAGAAAAAATCGCCGATTCCGGTGCTg                                                                                                                       <  1:3589927/149‑1 (MQ=255)
                                      cGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAg                                                                                                                >  2:3731005/1‑150 (MQ=255)
                                        acacCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAgcg                                                                                                              >  2:2613547/1‑150 (MQ=255)
                                              gCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTa                                                                                                        >  2:2780354/1‑150 (MQ=255)
                                                     ccttACTGAAATCATCCCACCAACAACTCACACCCCCCACAAACCATCCCAGCGCCTAACTACTGCCCGCCGTATTAGACCCCCCCCCCTGGCTTTCTAAAAACAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTg                                                                                                 <  1:1274004/147‑1 (MQ=255)
                                                     ccGTACTGAAATCATCCCCGCAACACCTGACACCGCCCTCAAACCCTGCCCGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCGGGCTTTCTAAAAATAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTg                                                                                                 <  1:1448144/150‑1 (MQ=255)
                                                      cGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGa                                                                                                >  2:3873301/1‑150 (MQ=255)
                                                      cGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGa                                                                                                >  2:672233/1‑150 (MQ=255)
                                                       gTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAg                                                                                               >  1:4591212/1‑150 (MQ=255)
                                                         aCTGAAATCATGCCAGCAACAACTGACACACCCCACAAACCATGCCAGCGCCGATCTACGGCCAGCCTTAGTTGTCTCCCCTCCCTGGCTTTCGAAAAATAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGcc                                                                                             <  1:4702718/150‑1 (MQ=255)
                                                          cTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCCGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTGTCGAAAAAGAAAAAATCGCCGAGTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCa                                                                                            >  1:91873/1‑150 (MQ=255)
                                                            taaaTCAGGCAAGCAACAACTGACACCGCCCTCAAACCATCCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCCCCCTGGCTTTAGAAAAACAAAAAACCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAAt                                                                                          <  2:1030157/149‑1 (MQ=255)
                                                                     ccAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCt                                                                                 <  2:3772629/150‑1 (MQ=255)
                                                                      cAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTg                                                                                >  2:1472939/1‑150 (MQ=255)
                                                                       aGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCTAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGa                                                                               <  1:2253450/150‑1 (MQ=255)
                                                                        gCAACAACTGACACCGCCCTCAAACCAGGCCAGCGCCGAGCTACTGCCAGCCGTAGTAGTCTCCCCCCCCTGGCTTTCGAAAAACAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGAc                                                                              <  1:1825629/150‑1 (MQ=255)
                                                                             aacTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACagaag                                                                         >  1:884686/1‑150 (MQ=255)
                                                                             aacTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACagaag                                                                         >  1:1994484/1‑150 (MQ=255)
                                                                                  acacAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATc                                                                    >  2:3033434/1‑150 (MQ=255)
                                                                                   cacaGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCa                                                                   >  1:2997778/1‑150 (MQ=255)
                                                                                   cacaGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCa                                                                   <  2:4630810/150‑1 (MQ=255)
                                                                                      aGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTa                                                                >  1:2026539/1‑150 (MQ=255)
                                                                                        cccTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAAt                                                              <  2:1673189/150‑1 (MQ=255)
                                                                                         ccTCAAACCAGGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCACCCTGGCTTTCGAAAAATAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATc                                                             <  1:4779703/150‑1 (MQ=255)
                                                                                            aaaaCCAGCCCAGCGCCGATCTCCTGCCAGCGGTAGTTGTCTCCCCTCCGGGGATTTCGAAAAACAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAAc                                                          <  1:2481021/149‑1 (MQ=255)
                                                                                               cccAGGCCCGCGCCGATCTCCTGCCAGCCGTTGTTGTCTCCCCTCCCGGGATTTCGAAAAACAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACttt                                                       <  1:4108617/149‑1 (MQ=255)
                                                                                               aCCTTGCCAGCGCCGATCTACTGCCGGCCGTTTTTGTCTCCCCTCCCTGGGTTTAGAAAAATAAATAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACttt                                                       <  1:2291077/150‑1 (MQ=255)
                                                                                               aCCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACttt                                                       >  1:1533487/1‑150 (MQ=255)
                                                                                                ccATGCCAGCGCCTATCTACTGCCAGCCGTTGTATTCTCCCCTCCGTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACtttt                                                      <  2:1193786/150‑1 (MQ=255)
                                                                                                 cctGCCAGCGCCTATATACTGCCAGCGTTATTAGTCCCCCCTCCCGGGATTTCTAAAAATAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTg                                                     <  1:2127136/148‑1 (MQ=255)
                                                                                                 caggCCCGCGCCTAGCTCCGGCCCGCCTTTGTCGTCCCCCCTCCGTGTCTTTCTAAAAATAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTg                                                     <  1:1590463/147‑1 (MQ=255)
                                                                                                    gCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCg                                                  >  1:2920193/1‑150 (MQ=255)
                                                                                                    gCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCg                                                  >  1:3834481/1‑150 (MQ=255)
                                                                                                     ccAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGa                                                 >  1:4633054/1‑150 (MQ=255)
                                                                                                       aGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGaaa                                               >  2:3564301/1‑150 (MQ=255)
                                                                                                       aGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGaaa                                               >  2:249206/1‑150 (MQ=255)
                                                                                                          gcCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAatat                                            >  2:497545/1‑150 (MQ=255)
                                                                                                           ccGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAatata                                           >  2:4608283/1‑150 (MQ=255)
                                                                                                            cGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAatatat                                          <  2:752646/150‑1 (MQ=255)
                                                                                                              aTCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCt                                        >  2:851254/1‑150 (MQ=255)
                                                                                                                cTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAATAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTgg                                      <  2:487033/150‑1 (MQ=255)
                                                                                                                   cTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTAt                                   <  2:221893/150‑1 (MQ=255)
                                                                                                                    tGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTCtt                                  >  2:4477773/1‑150 (MQ=255)
                                                                                                                     gCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTc                                 >  2:1938951/1‑150 (MQ=255)
                                                                                                                          ccttagttgTCTCCCCTCCCTGGCTTTCTAAAAATAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAac                            <  1:3351172/147‑1 (MQ=255)
                                                                                                                          ccGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAac                            >  1:2616084/1‑150 (MQ=255)
                                                                                                                           cGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAaca                           <  2:3294615/150‑1 (MQ=255)
                                                                                                                            gtagtaGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAacac                          >  2:341422/1‑150 (MQ=255)
                                                                                                                             tagtagTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGTATTTAGCGACATTAGGCATTGAGACAATCTGTTATCTGACAGAAGAAATCAGTAATAAACATTTGGCGAAATATATCTGGTATTCAAAAcaca                         >  2:589287/1‑150 (MQ=255)
                                                                                                                             tagtagTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAAcaca                         >  2:4588880/1‑150 (MQ=255)
                                                                                                                               gttgtCTCCCCTCCCGGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATc                       <  1:4111965/147‑1 (MQ=255)
                                                                                                                               gtagtCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATc                       >  2:1566088/1‑150 (MQ=255)
                                                                                                                                 agtCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCAc                     >  2:1232931/1‑150 (MQ=255)
                                                                                                                                 agtCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCAc                     >  1:2729170/1‑150 (MQ=255)
                                                                                                                                  gtCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACg                    <  1:2780361/150‑1 (MQ=255)
                                                                                                                                    ctcccctcccTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGTCAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGt                  >  2:3370461/1‑150 (MQ=255)
                                                                                                                                      cccctcccTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTa                <  1:1104428/150‑1 (MQ=255)
                                                                                                                                         ctccctGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGcc             <  1:3737562/150‑1 (MQ=255)
                                                                                                                                          tccctGGCTTTATAAAAATAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGCCa            <  1:1283228/150‑1 (MQ=255)
                                                                                                                                          tccctGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGCCa            >  1:4289775/1‑150 (MQ=255)
                                                                                                                                            cctGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGCCATg          <  1:804756/150‑1 (MQ=255)
                                                                                                                                            cctGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGCCATg          <  2:3132336/150‑1 (MQ=255)
                                                                                                                                             ctGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGCCCTga         >  1:1009555/1‑150 (MQ=255)
                                                                                                                                                 cTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGCCATgaagaa     <  1:3786975/150‑1 (MQ=255)
                                                                                                                                                   ttCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGCCATGAAGAAAg   <  2:1431167/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: