Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,191,496 C→A 37.0% 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,4960CA37.0% 344.4 / 228.6 236pseudogene (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 (76/72);  new base A (38/50);  total (114/122)
Fisher's exact test for biased strand distribution p-value = 2.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 240 total aligned reads displayed.

ACAACCAGCACTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGCCATGAAGAAAGTA  >  CP009273/2191347‑2191645
                                                                                                                                                     |                                                                                                                                                     
acaACCAGCACTGGCAGAACAGGTCATGCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTc                                                                                                                                                       <  2:3668810/150‑1 (MQ=255)
  tagcAGCTCTGGCAGAACAGGTCATTCCCTGGCTCTCGTCACCCGCAGTTGCCGTTCTGTAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCga                                                                                                                                                     <  1:64049/147‑1 (MQ=255)
  aaCCAGCACTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCga                                                                                                                                                     >  1:3656303/1‑150 (MQ=255)
    ccAGCACTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAgaaa                                                                                                                                                   <  1:2533210/150‑1 (MQ=255)
      aGCACTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCgaaaaa                                                                                                                                                 <  1:715008/150‑1 (MQ=255)
             gCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGaaaaaa                                                                                                                                          <  2:2146331/150‑1 (MQ=255)
              cAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAAt                                                                                                                                         <  2:620559/150‑1 (MQ=255)
              cAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAAt                                                                                                                                         <  1:3716506/150‑1 (MQ=255)
                gAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCg                                                                                                                                       <  2:1641958/150‑1 (MQ=255)
                gAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCg                                                                                                                                       >  2:406707/1‑150 (MQ=255)
                gAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCg                                                                                                                                       >  1:963243/1‑150 (MQ=255)
                        cATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTcc                                                                                                                               >  1:1741804/1‑150 (MQ=255)
                        cATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTcc                                                                                                                               <  1:2713596/150‑1 (MQ=255)
                         aTTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCg                                                                                                                              >  2:2639523/1‑150 (MQ=255)
                          ttCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCgg                                                                                                                             <  2:1582328/150‑1 (MQ=255)
                           tCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTAAAACCATGCCAGCGCCGATCTACTGCCAGCCGTTGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGt                                                                                                                            <  1:84739/150‑1 (MQ=255)
                               tGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCGCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTg                                                                                                                        <  1:882792/150‑1 (MQ=255)
                                 gCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTTGTCTCCCCTCCCTGGCTTTAGAAAATTAAAAAATCGCCGATTCCGGTGCTGGa                                                                                                                      <  2:4161760/150‑1 (MQ=255)
                                  ctctCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGAt                                                                                                                     <  1:3581149/150‑1 (MQ=255)
                                    ctcGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGAttt                                                                                                                   <  2:4057758/150‑1 (MQ=255)
                                    ctcGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGAttt                                                                                                                   <  2:1797731/150‑1 (MQ=255)
                                      cGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAg                                                                                                                 >  2:2452687/1‑150 (MQ=255)
                                      cGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAg                                                                                                                 >  2:3805949/1‑150 (MQ=255)
                                             cGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCCTGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAACGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCAtt                                                                                                          <  2:1466753/150‑1 (MQ=255)
                                               cAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACCGCCCTCAACCCAGGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAg                                                                                                        <  1:55745/150‑1 (MQ=255)
                                               cAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAg                                                                                                        >  1:2369903/1‑150 (MQ=255)
                                               cAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAg                                                                                                        >  1:1567681/1‑150 (MQ=255)
                                               cAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAg                                                                                                        <  2:3244166/150‑1 (MQ=255)
                                                 gTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGc                                                                                                      <  2:1234510/150‑1 (MQ=255)
                                                  ttGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCa                                                                                                     <  1:1330082/150‑1 (MQ=255)
                                                     ccGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTg                                                                                                  <  2:316047/150‑1 (MQ=255)
                                                        tACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGc                                                                                               <  2:718335/150‑1 (MQ=255)
                                                        tACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGAACTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTGTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGc                                                                                               >  1:2054694/1‑150 (MQ=255)
                                                         aCTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGcc                                                                                              >  1:1007620/1‑150 (MQ=255)
                                                             aaaTCATGCCAGCAACAACTGACACAGCCCTCAAACCCTGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATc                                                                                          <  2:2576831/150‑1 (MQ=255)
                                                             aaaTCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATc                                                                                          <  1:2421209/150‑1 (MQ=255)
                                                              aaTCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCt                                                                                         >  2:1502948/1‑150 (MQ=255)
                                                                tCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGt                                                                                       >  2:225340/1‑150 (MQ=255)
                                                                  aTGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTa                                                                                     <  2:1333031/150‑1 (MQ=255)
                                                                         caacaaCTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACa                                                                              <  2:574646/150‑1 (MQ=255)
                                                                          accaacTGACACAGCCCTCAACCCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACag                                                                             <  2:1012624/148‑1 (MQ=255)
                                                                             aacTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACagaag                                                                          >  1:4108708/1‑150 (MQ=255)
                                                                                tGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGaaa                                                                       >  1:864263/1‑150 (MQ=255)
                                                                                    acaGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAg                                                                   >  2:3590057/1‑150 (MQ=255)
                                                                                    acaGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCAGTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAg                                                                   >  1:2467341/1‑150 (MQ=255)
                                                                                      aGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTa                                                                 <  2:2077458/150‑1 (MQ=255)
                                                                                       gCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTaa                                                                >  1:2956179/1‑150 (MQ=255)
                                                                                        cccTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAAt                                                               <  1:1399327/150‑1 (MQ=255)
                                                                                            cAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAAc                                                           <  2:3626773/150‑1 (MQ=255)
                                                                                            cAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAAc                                                           >  1:1959685/1‑150 (MQ=255)
                                                                                            cAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAAc                                                           >  1:1856396/1‑150 (MQ=255)
                                                                                             aaaCCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACt                                                          <  2:2177889/150‑1 (MQ=255)
                                                                                             aaaCCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACt                                                          <  2:2973731/150‑1 (MQ=255)
                                                                                             aaaCCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACt                                                          <  2:883508/150‑1 (MQ=255)
                                                                                              aaCCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACtt                                                         >  1:2227878/1‑150 (MQ=255)
                                                                                              aaCCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACtt                                                         <  1:450017/150‑1 (MQ=255)
                                                                                                ccATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACtttt                                                       <  2:3553410/150‑1 (MQ=255)
                                                                                                ccATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACtttt                                                       >  1:450565/1‑150 (MQ=255)
                                                                                                 cATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTg                                                      >  1:3241968/1‑150 (MQ=255)
                                                                                                 cATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTg                                                      <  2:3089394/150‑1 (MQ=255)
                                                                                                       aGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGaaa                                                >  1:1246451/1‑150 (MQ=255)
                                                                                                       aGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGaaa                                                >  1:1187982/1‑150 (MQ=255)
                                                                                                        gcgcCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAat                                               <  2:1955124/150‑1 (MQ=255)
                                                                                                          gcCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAatat                                             >  2:1183500/1‑150 (MQ=255)
                                                                                                           ccGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAatata                                            >  2:3323275/1‑150 (MQ=255)
                                                                                                             gATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATc                                          <  2:338068/150‑1 (MQ=255)
                                                                                                              aTCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCt                                         >  2:1930229/1‑150 (MQ=255)
                                                                                                                   cTGCCAGCCGTAGTAGTGTGCCCTACCGGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTATGCATTGAGCCTATCTGTTATCTGACAGAAGAAATCTGTAATAAACTTTTGGCGAAATTTATCTGGTAt                                    <  2:2297642/150‑1 (MQ=255)
                                                                                                                   cTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTAt                                    <  2:2488415/150‑1 (MQ=255)
                                                                                                                     gCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTc                                  <  2:169964/150‑1 (MQ=255)
                                                                                                                       cAGCCGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCaa                                >  2:395302/1‑150 (MQ=255)
                                                                                                                         gCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCaaaa                              <  2:3459319/150‑1 (MQ=255)
                                                                                                                           cGTAGTAGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAaca                            <  2:2790697/150‑1 (MQ=255)
                                                                                                                            gtagtaGTCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAacac                           <  2:427788/150‑1 (MQ=255)
                                                                                                                             tagtagTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAAcaca                          <  2:2076322/150‑1 (MQ=255)
                                                                                                                              agtagtCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACAt                         >  2:3821484/1‑150 (MQ=255)
                                                                                                                              agtagtCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACAt                         >  1:1648026/1‑150 (MQ=255)
                                                                                                                                tagtCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCa                       <  1:2912321/150‑1 (MQ=255)
                                                                                                                                 agtCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCAc                      >  1:3337230/1‑150 (MQ=255)
                                                                                                                                 agtCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCAc                      >  1:2985036/1‑150 (MQ=255)
                                                                                                                                  gtCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACg                     >  1:791729/1‑150 (MQ=255)
                                                                                                                                  gtCTCCCCTCCCTGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACg                     <  1:2331116/150‑1 (MQ=255)
                                                                                                                                   tctcCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACgg                    <  1:741187/150‑1 (MQ=255)
                                                                                                                                    ctcccctcccTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGt                   <  1:1983426/150‑1 (MQ=255)
                                                                                                                                     tcccctcccTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGtt                  <  1:152865/150‑1 (MQ=255)
                                                                                                                                       ccctccctGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAg                <  1:380238/150‑1 (MQ=255)
                                                                                                                                       ccctccctGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAg                <  1:2582583/150‑1 (MQ=255)
                                                                                                                                       ccctccctGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAg                <  1:3062619/150‑1 (MQ=255)
                                                                                                                                         ctccctGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCAGTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGcc              >  1:501944/1‑150 (MQ=255)
                                                                                                                                         ctccctGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGcc              >  1:1450850/1‑150 (MQ=255)
                                                                                                                                           ccctGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGCCAt            <  1:4155614/150‑1 (MQ=255)
                                                                                                                                            cctGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGCCATg           >  2:2886764/1‑150 (MQ=255)
                                                                                                                                            catGGCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCAATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACAAATCACGGTAAGCCATg           >  2:510032/3‑150 (MQ=255)
                                                                                                                                             ctGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGCCATga          <  2:317574/150‑1 (MQ=255)
                                                                                                                                              tGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGCCATgaa         >  2:2110780/1‑150 (MQ=255)
                                                                                                                                               ggCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGCCATgaag        <  1:3622399/150‑1 (MQ=255)
                                                                                                                                               ggCTTTAGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGCCATgaag        >  2:2068052/1‑150 (MQ=255)
                                                                                                                                                   ttCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGCCATGAAGAAAg    <  2:3081571/150‑1 (MQ=255)
                                                                                                                                                   ttCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGCCATGAAGAAAg    >  2:2913564/1‑150 (MQ=255)
                                                                                                                                                     cGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGCCATGAAGAAAGta  >  2:1840070/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACAACCAGCACTGGCAGAACAGGTCATTCCCTGGCTCTCGACACCCGCAGTTGCCGTACTGAAATCATGCCAGCAACAACTGACACAGCCCTCAAACCATGCCAGCGCCGATCTACTGCCAGCCGTAGTAGTCTCCCCTCCCTGGCTTTCGAAAAAGAAAAAATCGCCGATTCCGGTGCTGGATTTAGCGCCATTAGGCATTGAGCCAATCTGTTATCTGACAGAAGAAATCAGTAATCAACTTTTGGCGAAATATATCTGGTATTCAAAACACATCACGGTTAGCCATGAAGAAAGTA  >  CP009273/2191347‑2191645

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