Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,191,496 C→A 5.7% 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,4960CA 5.7% 1077.1 / 23.9 281pseudogene (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 (133/131);  new base A (4/12);  total (137/144)
Fisher's exact test for biased strand distribution p-value = 6.96e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 285 total aligned reads displayed.

ACAACCAGCACTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGCCATGAAGAAAGTA  >  CP009273/2191347‑2191645
                                                                                                                                                     |                                                                                                                                                     
acaACCAGCACTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTc                                                                                                                                                       <  1:135614/150‑1 (MQ=255)
 caACCAGCACTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCg                                                                                                                                                      >  1:1584515/1‑150 (MQ=255)
   aCCAGCACTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCgaa                                                                                                                                                    <  1:165708/150‑1 (MQ=255)
    ccAGCACTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCgaaa                                                                                                                                                   >  1:3394299/1‑150 (MQ=255)
      aGCACTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCgaaaaa                                                                                                                                                 >  2:2234453/1‑150 (MQ=255)
      aGCACTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCgaaaaa                                                                                                                                                 >  2:426807/1‑150 (MQ=255)
      aGCACTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCgaaaaa                                                                                                                                                 >  1:2551712/1‑150 (MQ=255)
        cACTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCgaaaaaga                                                                                                                                               >  2:841338/1‑150 (MQ=255)
        cACTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCgaaaaaga                                                                                                                                               >  2:70393/1‑150 (MQ=255)
         aCTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCgaaaaagaa                                                                                                                                              <  1:4177421/150‑1 (MQ=255)
          cTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCgaaaaagaaa                                                                                                                                             <  1:2696599/150‑1 (MQ=255)
           tGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATAATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCgaaaaagaaaa                                                                                                                                            >  2:1996976/1‑150 (MQ=255)
            ggCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCgaaaaagaaaaa                                                                                                                                           >  1:1021307/1‑150 (MQ=255)
                gAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCg                                                                                                                                       >  1:944496/1‑150 (MQ=255)
                gAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCg                                                                                                                                       <  1:868004/150‑1 (MQ=255)
                    aGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGa                                                                                                                                   >  2:3902770/1‑150 (MQ=255)
                       tCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTc                                                                                                                                >  1:3832648/1‑150 (MQ=255)
                       tCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTc                                                                                                                                >  1:3787376/1‑150 (MQ=255)
                       tCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCGCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTc                                                                                                                                <  1:1392423/150‑1 (MQ=255)
                        cATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTcc                                                                                                                               >  2:3679143/1‑150 (MQ=255)
                        cATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTcc                                                                                                                               >  2:3971215/1‑150 (MQ=255)
                            cccTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTg                                                                                                                           <  1:3018409/150‑1 (MQ=255)
                            cccTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTg                                                                                                                           >  1:328804/1‑150 (MQ=255)
                               tGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTg                                                                                                                        >  1:249443/1‑150 (MQ=255)
                                 gCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGa                                                                                                                      >  2:3378931/1‑150 (MQ=255)
                                     tcGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTa                                                                                                                  >  2:1638768/1‑150 (MQ=255)
                                      cGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAg                                                                                                                 >  2:879142/1‑150 (MQ=255)
                                       gACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAgc                                                                                                                <  2:2988885/150‑1 (MQ=255)
                                        acacCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAgcg                                                                                                               >  2:363637/1‑150 (MQ=255)
                                          acCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGcc                                                                                                             >  1:4297040/1‑150 (MQ=255)
                                          acCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGcc                                                                                                             >  1:2572377/1‑150 (MQ=255)
                                           cccGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCa                                                                                                            >  2:2620944/1‑150 (MQ=255)
                                           cccGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCa                                                                                                            <  1:4574983/150‑1 (MQ=255)
                                                   tGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCAt                                                                                                    >  1:2399780/1‑150 (MQ=255)
                                                      agTACTGAAATCATGCCAGCAACAACTGTCAACGCCATCAAAGCATGCCTGAGCCGAGCTAGTGCCAGCCGTAGTTGTCTCCCCTCCCTGGCTTTCGAAATACAAAAAATCGCCGATTCCTGGGCTGGATTTAGCGCCATTTGGCATTGa                                                                                                 <  2:2565525/149‑1 (MQ=255)
                                                       gTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCACACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAg                                                                                                <  1:907468/150‑1 (MQ=255)
                                                       gTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAg                                                                                                >  1:3071898/1‑150 (MQ=255)
                                                          cTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCa                                                                                             <  1:1113394/150‑1 (MQ=255)
                                                           tGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCaa                                                                                            >  2:1581011/1‑150 (MQ=255)
                                                           tGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCaa                                                                                            >  1:49692/1‑150 (MQ=255)
                                                            gAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAAt                                                                                           <  2:3784171/150‑1 (MQ=255)
                                                            gAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAAt                                                                                           <  2:2518040/150‑1 (MQ=255)
                                                              aaTCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCt                                                                                         >  2:1444888/1‑150 (MQ=255)
                                                                   tGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTAt                                                                                    <  1:1780292/150‑1 (MQ=255)
                                                                      cAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTg                                                                                 <  2:4493826/150‑1 (MQ=255)
                                                                       aGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGa                                                                                >  2:3368683/1‑150 (MQ=255)
                                                                       aGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGa                                                                                >  1:4603072/1‑150 (MQ=255)
                                                                         caacaaCTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACa                                                                              >  2:139190/1‑150 (MQ=255)
                                                                          aacaacTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACag                                                                             >  1:3711181/1‑150 (MQ=255)
                                                                               cTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAgaagaa                                                                        >  1:3102436/1‑150 (MQ=255)
                                                                                   cacaGCCCTCAAACCATGCCAGCGCCGATCTCCTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCa                                                                    <  1:1048268/150‑1 (MQ=255)
                                                                                   cacaGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCa                                                                    <  2:1516190/150‑1 (MQ=255)
                                                                                    acaGCCCTCAACCCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCTAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAg                                                                   <  1:1198238/150‑1 (MQ=255)
                                                                                    acaGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAg                                                                   <  2:486497/150‑1 (MQ=255)
                                                                                      aGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTa                                                                 >  1:1478266/1‑150 (MQ=255)
                                                                                      aGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGATTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTa                                                                 <  1:1227673/150‑1 (MQ=255)
                                                                                       gCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTaa                                                                >  1:234345/1‑150 (MQ=255)
                                                                                         ccTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATc                                                              >  1:3503678/1‑150 (MQ=255)
                                                                                         ccTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATc                                                              <  2:1698257/150‑1 (MQ=255)
                                                                                          cTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTTGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCa                                                             <  1:4101954/150‑1 (MQ=255)
                                                                                          cTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCa                                                             <  2:1174251/150‑1 (MQ=255)
                                                                                          cTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCa                                                             >  2:2057313/1‑150 (MQ=255)
                                                                                            cAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAAc                                                           <  2:2674119/150‑1 (MQ=255)
                                                                                               aCCATGCCAGCGCCGATCTACTGCCAGCCGTTGTAGTCTCCCCTCCCTTGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACttt                                                        <  2:819851/150‑1 (MQ=255)
                                                                                               aCCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACttt                                                        <  1:625465/150‑1 (MQ=255)
                                                                                                 cATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTg                                                      >  1:2842779/1‑150 (MQ=255)
                                                                                                 cATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTg                                                      <  1:1927739/150‑1 (MQ=255)
                                                                                                  aTGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTgg                                                     <  2:1383481/150‑1 (MQ=255)
                                                                                                        gcgcCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAat                                               <  2:3401000/150‑1 (MQ=255)
                                                                                                            cGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAatatat                                           >  1:3591362/1‑150 (MQ=255)
                                                                                                            cGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAatatat                                           >  2:572145/1‑150 (MQ=255)
                                                                                                              aTCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCt                                         <  2:2814544/150‑1 (MQ=255)
                                                                                                               tCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTg                                        >  1:127520/1‑150 (MQ=255)
                                                                                                               tCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTg                                        <  1:219105/150‑1 (MQ=255)
                                                                                                                cTACTGCCAGACTTAGTAGTCTCCCCTCCCTGGTTTTGTAAAAAGAAAAAATCGCGGATTCCGGTGCTGGATTTAGCGCCTTTTGTGATTGTGACAATCTGTTATCTGACAGAATAAATCAGTAATCTTCTTTTGGCGAAATATATCTgg                                       <  2:3125435/150‑1 (MQ=255)
                                                                                                                       cAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCaa                                >  2:3184120/1‑150 (MQ=255)
                                                                                                                       cAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGACATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCaa                                >  1:2031068/1‑150 (MQ=255)
                                                                                                                          ccGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAac                             >  2:130962/1‑150 (MQ=255)
                                                                                                                             tagtagTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGGTTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAAcaca                          >  1:25027/1‑150 (MQ=255)
                                                                                                                                    ctcccctcccTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGt                   <  1:1127091/150‑1 (MQ=255)
                                                                                                                                     tcccctcccTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGtt                  >  2:1079190/1‑150 (MQ=255)
                                                                                                                                        cctccctGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGc               >  1:3264551/1‑150 (MQ=255)
                                                                                                                                         ctccctGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGcc              >  1:4423052/1‑150 (MQ=255)
                                                                                                                                           ccctGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGCCAt            <  2:1627514/150‑1 (MQ=255)
                                                                                                                                           ccctGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGCCAt            <  2:235573/150‑1 (MQ=255)
                                                                                                                                           ccctGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGCCAt            <  2:4038575/150‑1 (MQ=255)
                                                                                                                                           ccctGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGCCAt            <  2:542821/150‑1 (MQ=255)
                                                                                                                                             ctGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGCCATga          >  2:733872/1‑150 (MQ=255)
                                                                                                                                             ctGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGCCATga          <  1:3320517/150‑1 (MQ=255)
                                                                                                                                               ggCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGCCATgaag        >  2:3054638/1‑150 (MQ=255)
                                                                                                                                               ggCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGCCATgaag        <  1:4561938/150‑1 (MQ=255)
                                                                                                                                               ggCTATCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGCCATgaag        <  2:4105975/150‑1 (MQ=255)
                                                                                                                                                gCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGCCATgaaga       >  2:188350/1‑150 (MQ=255)
                                                                                                                                                 cTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGCCATgaagaa      >  1:3176312/1‑150 (MQ=255)
                                                                                                                                                 cTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGCCATgaagaa      >  1:535879/1‑150 (MQ=255)
                                                                                                                                                  tttCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGCCATGAAGaaa     >  1:2257412/1‑150 (MQ=255)
                                                                                                                                                   ttCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGCCATGAAGAAAg    <  2:4590104/150‑1 (MQ=255)
                                                                                                                                                   ttCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGCCATGAAGAAAg    >  1:1704125/1‑150 (MQ=255)
                                                                                                                                                    tCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGCCATGAAGAAAGt   <  1:841351/150‑1 (MQ=255)
                                                                                                                                                     cGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGCCATGAAGAAAGta  <  1:4048956/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACAACCAGCACTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGCCATGAAGAAAGTA  >  CP009273/2191347‑2191645

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: