Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,192,748 T→A 19.9% pseudogene (2796/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,192,7480TA19.9% 659.8 / 157.5 271pseudogene (2796/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 T (102/115);  new base A (30/24);  total (132/139)
Fisher's exact test for biased strand distribution p-value = 2.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 280 total aligned reads displayed.

GCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCGGGTTAGCACCGGATCTCGGTCTGGATGATAACGGTTCGCTGCTGCTGGATTTTGGCCCACGGCAGTTCACCGTCAGCTTTGATGAAACCTTAAAACCGT  >  CP009273/2192599‑2192897
                                                                                                                                                     |                                                                                                                                                     
gcctggcctgGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAAtt                                                                                                                                                       >  1:535611/1‑150 (MQ=255)
gcctggcctgGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAAtt                                                                                                                                                       <  1:2683018/150‑1 (MQ=255)
gcctggcctgGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAAtt                                                                                                                                                       <  1:28307/150‑1 (MQ=255)
    ggcctggTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCg                                                                                                                                                   <  2:3289799/150‑1 (MQ=255)
       ctggTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACa                                                                                                                                                <  1:961644/150‑1 (MQ=255)
          gTGAATCCCAGCATAAACGCGCCACCGTTGGTTTGGATATTCTCGCTGCTATAGGTATTGATATCGCCCTTATGCAGCTTAACGGAATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTg                                                                                                                                             <  1:3332368/150‑1 (MQ=255)
          gTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTg                                                                                                                                             <  1:1252470/150‑1 (MQ=255)
            gAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGcc                                                                                                                                           <  2:2465090/150‑1 (MQ=255)
                cccAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGCCgaga                                                                                                                                       >  2:1269281/1‑150 (MQ=255)
                 ccAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGCCGagag                                                                                                                                      <  2:3760077/150‑1 (MQ=255)
                  cAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGCCGAGAGc                                                                                                                                     <  2:4019075/150‑1 (MQ=255)
                  cAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGCCGAGAGc                                                                                                                                     <  1:2349774/150‑1 (MQ=255)
                   aGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGCCGAGAGcc                                                                                                                                    >  2:4029166/1‑150 (MQ=255)
                   aGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGCCGAGAGcc                                                                                                                                    <  1:120256/150‑1 (MQ=255)
                    gCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGCCGAGAGCcg                                                                                                                                   <  2:1716491/150‑1 (MQ=255)
                     cATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGCCGAGAGCcgc                                                                                                                                  <  2:2268965/150‑1 (MQ=255)
                      aTAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGCCGAGAGCcgcg                                                                                                                                 <  1:1854458/150‑1 (MQ=255)
                       tAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGCCGAGAGCCGCGa                                                                                                                                <  2:3116893/150‑1 (MQ=255)
                            gcgcCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGCCGAGAGCCGCGAACTCa                                                                                                                           <  2:2881890/150‑1 (MQ=255)
                            gcgcCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATAGCCGACATTGCCGAGAGCCGCGAACTCa                                                                                                                           >  2:4496443/1‑150 (MQ=255)
                              gcCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGCCGAGAGCCGCGAACTCACg                                                                                                                         <  2:455496/150‑1 (MQ=255)
                                 aCCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGCCGAGAGCCGCGAACTCACGGTg                                                                                                                      <  1:2487849/150‑1 (MQ=255)
                                   cggtgggttggATAGTCTCGCTGCTATCGGGAGTGATATCGCCCTTATGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGCCATTGCCGAGAGCCGCGAACTCACGGTGGc                                                                                                                    <  1:3896938/147‑1 (MQ=255)
                                   cGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGc                                                                                                                    >  2:301067/1‑150 (MQ=255)
                                   cGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGc                                                                                                                    <  1:3896956/150‑1 (MQ=255)
                                   cGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATAGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGc                                                                                                                    <  2:1162278/150‑1 (MQ=255)
                                   cGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATAGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGc                                                                                                                    >  1:1622067/1‑150 (MQ=255)
                                   cGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCTTAACGGCATCGCCAAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTTCCGACATTGCCGAGAGCCGCGAACTCACGGTGGc                                                                                                                    <  2:902528/150‑1 (MQ=255)
                                    gttgggttggATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATAGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCg                                                                                                                   <  2:2062624/150‑1 (MQ=255)
                                      tgggttggATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGa                                                                                                                 <  2:2353411/150‑1 (MQ=255)
                                       gggttggATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAg                                                                                                                >  2:1331699/1‑150 (MQ=255)
                                         gttggATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAACGAAAAAATAGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCt                                                                                                              <  1:1750777/150‑1 (MQ=255)
                                          ttggATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATAGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCtt                                                                                                             <  1:21595/150‑1 (MQ=255)
                                               tatTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATAGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGa                                                                                                        >  1:4465353/1‑150 (MQ=255)
                                                 ttCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATc                                                                                                      <  1:3454566/150‑1 (MQ=255)
                                                  tctcGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCg                                                                                                     >  1:3085356/1‑150 (MQ=255)
                                                  tctcGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCg                                                                                                     >  2:592546/1‑150 (MQ=255)
                                                   ctcGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCgg                                                                                                    >  1:161673/1‑150 (MQ=255)
                                                    tcGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCggg                                                                                                   >  2:1951996/1‑150 (MQ=255)
                                                     cGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCGGGt                                                                                                  >  2:825573/1‑150 (MQ=255)
                                                      gctgctATCGGTAGTGATATCGCCCTTATGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCGGGtt                                                                                                 <  1:1982072/150‑1 (MQ=255)
                                                       ctgctATCGGTAGTGATATCGCCCTTATGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCGGGTTa                                                                                                <  2:1445621/150‑1 (MQ=255)
                                                           tATCGGTAGTGATATCGCCCTTATGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCGGGTTAGCAc                                                                                            >  2:3350058/1‑150 (MQ=255)
                                                            aTCGGTAGTGATATCGCCCTTATGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATAGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCGGGTTAGCAcc                                                                                           >  2:3409805/1‑150 (MQ=255)
                                                               ggTAGTGATATCGCCCTTATGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCGGGTTAGCACCGGa                                                                                        >  1:4119138/1‑150 (MQ=255)
                                                                  aGTGATATCGCCCTTATGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCGGGTTAGCACCGGAtct                                                                                     >  2:2599817/1‑150 (MQ=255)
                                                                  aGTGATATCGCCCTTATGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACAGTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCGGGTTAGCACCGGAtct                                                                                     <  1:1601422/150‑1 (MQ=255)
                                                                   gTGATATCGCCCTTATGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGTCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCGGGTTAGCACCGGAtctc                                                                                    >  2:2307253/1‑150 (MQ=255)
                                                                    tGATATCGCCCTTATGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCGGGTTAGCACCGGATCTCg                                                                                   >  2:1653230/1‑150 (MQ=255)
                                                                    tGATATCGCCCTTATGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCGGGTTAGCACCGGATCTCg                                                                                   >  2:1534304/1‑150 (MQ=255)
                                                                    tGATATCGCCCTTATGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCGGGTTAGCACCGGATCTCg                                                                                   <  2:3182190/150‑1 (MQ=255)
                                                                     gATATCGCCCTTATGCAGCTTAACGGGATCGCCCAGAAACTGAAATTCATAGCATGAAAGGAGCGGGCAAAAGCAAAAATTGGCGACATTGCCGAGAGCCGCGATCTCACGGTGGCGGAGCTTGATGATCGGGTTAGCACCGGAtctagt                                                                                  <  2:3032413/150‑4 (MQ=255)
                                                                      atatCGCCCTTATGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCGGGTTAGCACCGGATCTCGGt                                                                                 >  1:852258/1‑150 (MQ=255)
                                                                         tCGCCCTTATGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCGGGTTAGCACCGGATCTCGGTCTg                                                                              <  2:1517949/150‑1 (MQ=255)
                                                                             ccTTATGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGACGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCGGGTTAGCACCGGATCTCGGTCTGgatg                                                                          >  1:1034347/1‑150 (MQ=255)
                                                                                 aTGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCGGGTTAGCACCGGATCTCGGTCTGGATGATaa                                                                      >  2:1504512/1‑150 (MQ=255)
                                                                                 aTGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCGGGTTAGCACCGGATCTCGGTCTGGATGATaa                                                                      <  1:762201/150‑1 (MQ=255)
                                                                                 aTGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCGGGTTAGCACCGGATCTCGGTCTGGATGATaa                                                                      >  2:2036021/1‑150 (MQ=255)
                                                                                  tGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCGGGTTAGCACCGGATCTCGGTCTGGATGATAAc                                                                     <  2:2604459/150‑1 (MQ=255)
                                                                                   gCAGCTTAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCGGGTTAGCACCGGATCTCGGTCTGGATGATAACg                                                                    <  1:2178836/150‑1 (MQ=255)
                                                                                     aGCTTAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGCTGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCGGGTTAGCACCGGATCTCGGTCTGGATGATAACGGt                                                                  >  1:1004448/1‑150 (MQ=255)
                                                                                       cTTAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCGGGTTAGCACCGGATCTCGGTCTGGATGATAACGGTTc                                                                >  2:3593326/1‑150 (MQ=255)
                                                                                        ttAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCGGGTTAGCACCGGATCTCGGTCTGGATGATAACGGTTCg                                                               <  1:801004/150‑1 (MQ=255)
                                                                                        ttAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCGGGTTAGCACCGGATCTCGGTCTGGATGATAACGGTTCg                                                               <  1:2203776/150‑1 (MQ=255)
                                                                                          aaCGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATAGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCGGGTTAGCACCGGATCTCGGTCTGGATGATAACGGTTCgct                                                             >  2:4195302/1‑150 (MQ=255)
                                                                                            cGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCGGGTTAGCACCGGATCTCGGTCTGGATGATAACGGTTCgctgc                                                           >  1:3832292/1‑150 (MQ=255)
                                                                                                aTCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATAGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCGGGTTAGCACCGGATCTCGGTCTGGATGATAACGGTTCgctgctgct                                                       >  1:2525068/1‑150 (MQ=255)
                                                                                                aTCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATAGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCGGGTTAGCACCGGATCTCGGTCTGGATGATAACGGTTCgctgctgct                                                       >  2:521257/1‑150 (MQ=255)
                                                                                                 tCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCGGGTTAGCACCGGATCTCGGTCTGGATGATAACGGTTCGctgctgctg                                                      <  1:397542/150‑1 (MQ=255)
                                                                                                 tCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCGGGTTAGCACCGGATCTCGGTCTGGATGATAACGGTTCGctgctgctg                                                      <  1:1420359/150‑1 (MQ=255)
                                                                                                 tCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCGGGTTAGCACCGGATCTCGGTCTGGATGATAACGGTTCGctgctgctg                                                      <  1:2390939/150‑1 (MQ=255)
                                                                                                  cGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCGGGTTAGCACCGGATCTCGGTCTGGATGATAACGGTTCGCTGCTGCTgg                                                     >  1:1431526/1‑150 (MQ=255)
                                                                                                     ccAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATAGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCGGGTTAGCACCGGATCTCGGTCTGGATGATAACGGTTCGCTGCTGCTGGAtt                                                  <  1:698746/150‑1 (MQ=255)
                                                                                                        gAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCGGGTTAGCACCGGATCTCGGTCTGGATGATAACGGTTCGCTGCTGCTGGATTTTg                                               <  2:2771433/150‑1 (MQ=255)
                                                                                                         aaaCTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCGGGTTAGCACCGGATCTCGGTCTGGATGATAACGGTTCGCTGCTGCTGGATTTTgg                                              >  1:818679/1‑150 (MQ=255)
                                                                                                            cTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCGGGTTAGCACCGGATCTCGGTCTGGATGATAACGGTTCGCTGCTGCTGGATTTTGGccc                                           >  2:2535500/1‑150 (MQ=255)
                                                                                                              gAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCGGGTTAGCACCGGATCTCGGTCTGGATGATAACGGTTCGCTGCTGCTGGATTTTGGCCCAc                                         <  2:1766616/150‑1 (MQ=255)
                                                                                                               aaaTTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCGGGTTAGCACCGGATCTCGGTCTGGATGATAACGGTTCGCTGCTGCTGGATTTTGGCCCACg                                        <  1:4331340/150‑1 (MQ=255)
                                                                                                                aaTTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCGGGTTAGCACCGGATCTCGGTCTGGATGATAACGGTTCGCTGCTGCTGGATTTTGGCCCACgg                                       >  1:3398041/1‑150 (MQ=255)
                                                                                                                aaTTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACAGTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCGGGTTAGCACCGGATCTCGGTCTGGATGATAACGGTTCGCTGCTGCTGGATTTTGGCCCACgg                                       <  2:83084/150‑1 (MQ=255)
                                                                                                                     aaaGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCGGGTTAGCACCGGATCTCGGTCTGGATGATAACGGTTCGCTGCTGCTGGATTTTGGCCCACGGCAGtt                                  <  1:4377129/150‑1 (MQ=255)
                                                                                                                     aaaGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCGGGTTAGCACCGGATCTCGGTCTGGATGATAACGGTTCGCTGCTGCTGGATTTTGGCCCACGGCAGtt                                  <  2:3677851/150‑1 (MQ=255)
                                                                                                                         cATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCGGGTTAGCACCGGATCTCGGTCTGGATGATAACGGTTCGCTGCTGCTGGATTTTGGCCCACGGCAGTTCAcc                              >  2:2714244/1‑150 (MQ=255)
                                                                                                                         cATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCGGGTTAGCACCGGATCTCGGTCTGGATGATAACGGTTCGCTGCTGCTGGATTTTGGCCCACGGCAGTTCAcc                              >  2:3808731/1‑150 (MQ=255)
                                                                                                                           ttACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCGGGTTAGCACCGGATCTCGGTCTGGATGATAACGGTTCGCTGCTGCTGGATTTTGGCCCACGGCAGTTCACCGt                            >  1:1532885/1‑150 (MQ=255)
                                                                                                                                ggAGCGGGCAAAAGAAAAAATTGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCGGGTTAGCACCGGATCTCGGTCTGGATGATAACGGTTCGCTGCTGCTGGATTTTGGCCCACGGCAGTTCACCGTCAGCt                       >  2:732998/1‑150 (MQ=255)
                                                                                                                                  aGCGGGCAAAAGAAAAAATAGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCGGGTTAGCACCGGATCTCGGTCTGGATGATAACGGTTCGCTGCTGCTGGATTTTGGCCCACGGCAGTTCACCGTCAGCttt                     >  2:3169973/1‑150 (MQ=255)
                                                                                                                                        cAAAAGAAAAAATTGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCGGGTTAGCACCGGATCTCGGTCTGGATGATAACGGTTCGCTGCTGCTGGATTTTGGCCCACGGCAGTTCACCGTCAGCTTTGATGaa               <  1:4528531/150‑1 (MQ=255)
                                                                                                                                          acaGAAATAATAGCCGACATTGCCGAGAGCCGAGATCTCACGGTGGCGGAGCTTGAAGATCGGGTTTGCACCGGATCTCGGTCTGGATGATATCGGTTCGCTGCTGCTGGATTTTGGCCCACGGAAGTTCACCGTCAGCTTTGATGAAAc             <  2:242485/148‑1 (MQ=255)
                                                                                                                                             gAAAAAATTGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCGGGTTAGCACCGGATCTCGGTCTGGATGATAACGGTTCGCTGCTGCTGGATTTTGGCCCACGGCAGTTCACCGTCAGCTTTGATGAAACCtt          <  2:4405599/150‑1 (MQ=255)
                                                                                                                                             gAAAAAATAGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCGGGTTAGCACCGGATCTCGGTCTGGATGATAACGGTTCGCTGCTGCTGGATTTTGGCCCACGGCAGTTCACCGTCAGCTTTGATGAAACCtt          >  1:1107739/1‑150 (MQ=255)
                                                                                                                                               aaaaaTAGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCGGGTTAGCACCGGATCTCGGTCTGGATGATAACGGTTCGCTGCTGCTGGATTTTGGCCCACGGCAGTTCACCGTCAGCTTTGATGAAACCTTaa        >  1:646693/1‑150 (MQ=255)
                                                                                                                                               aaaaaTAGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCGGGTTAGCACCGGATCTCGGTCTGGATGATAACGGTTCGCTGCTGCTGGATTTTGGCCCACGGCAGTTCACCGTCAGCTTTGATGAAACCTTaa        <  2:2176946/150‑1 (MQ=255)
                                                                                                                                                aaaaTTGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCGGGTTAGCACCGGATCTCGGTCTGGATGATAACGGTTCGCTGCTGCTGGATTTTGGCCCACGGCAGTTCACCGTCAGCTTTGATGAAACCTTaaa       >  1:983077/1‑150 (MQ=255)
                                                                                                                                                aaaaTTGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCGGGTTAGCACCGGATCTCGGTCTGGATGATAACGGTTCGCTGCTGCTGGATTTTGGCCCACGGCAGTTCACCGTCAGCTTTGATGAAACCTTaaa       <  2:1467687/150‑1 (MQ=255)
                                                                                                                                                 aaaTTGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCGGGTTAGCACCGGATCTCGGTCTGGATGATAACGGTTCGCTGCTGCTGGATTTTGGCCCACGGCAGTTCACCGTCAGCTTTGATGAAACCTTaaaa      <  2:4286335/150‑1 (MQ=255)
                                                                                                                                                  aaTTGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCGGGTTAGCACCGGATCTCGGTCTGGATGATAACGGTTCGCTGCTGCTGGATTTTGGCCCACGGCAGTTCACCGTCAGCTTTGATGAAACCTTAAAAc     <  2:1412332/150‑1 (MQ=255)
                                                                                                                                                    ttGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCGGGTTAGCACCGGATCTCGGTCTGGATGATAACGGTTCGCTGCTGCTGGATTTTGGCCCACGGCAGTTCACCGTCAGCTTTGATGAAACCTTAAAACCg   >  1:4202960/1‑150 (MQ=255)
                                                                                                                                                    ttGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCGGGTTAGCACCGGATCTCGGTCTGGATGATAACGGTTCGCTGCTGCTGGATTTTGGCCCACGGCAGTTCACCGTCAGCTTTGATGAAACCTTAAAACCg   <  1:3251393/150‑1 (MQ=255)
                                                                                                                                                     tGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCGGGTTAGCACCGGATCTCGGTCTGGATGATAACGGTTCGCTGCTGCTGGATTTTGGCCCACGGCAGTTCACCGTCAGCTTTGATGAAACCTTAAAACCGt  >  1:1793434/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCCTGGCCTGGTGAATCCCAGCATAAACGCGCCACCGTTGGGTTGGATATTCTCGCTGCTATCGGTAGTGATATCGCCCTTATGCAGCTTAACGGCATCGCCCAGAAACTGAAATTCAAAGCATTACAGGAGCGGGCAAAAGAAAAAATTGCCGACATTGCCGAGAGCCGCGAACTCACGGTGGCGGAGCTTGAAGATCGGGTTAGCACCGGATCTCGGTCTGGATGATAACGGTTCGCTGCTGCTGGATTTTGGCCCACGGCAGTTCACCGTCAGCTTTGATGAAACCTTAAAACCGT  >  CP009273/2192599‑2192897

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: