Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,189,888 C→A 5.6% intergenic (+76/‑65) metG → / → yehH methionyl‑tRNA synthetase/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,189,8880CA 5.6% 1045.7 / 24.4 277intergenic (+76/‑65)metG/yehHmethionyl‑tRNA synthetase/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
Reads supporting (aligned to +/- strand):  ref base C (135/126);  new base A (11/5);  total (146/131)
Fisher's exact test for biased strand distribution p-value = 2.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 293 total aligned reads displayed.

CGGTCCTGGCGGGAAAGATATTTTCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCATATCAGTTGGG  >  CP009273/2189740‑2190037
                                                                                                                                                    |                                                                                                                                                     
cGGTCCTGGCGGGAAAGATATTTTCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGtctct                                                                                                                                                      <  2:2286319/150‑1 (MQ=255)
   tCCTGGCGGGAAAGATATTTTCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGtctctctc                                                                                                                                                   >  1:3409631/1‑150 (MQ=255)
    ccTGGCGGGAAAGATATTTTCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTctctctct                                                                                                                                                  <  1:1067812/150‑1 (MQ=255)
     cTGGCGGGAAAGATATTTTCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTa                                                                                                                                                 >  1:3035175/1‑150 (MQ=255)
      tGGCGGGAAAGATATTTTCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTaa                                                                                                                                                >  1:1493888/1‑150 (MQ=255)
          gggAAAGATATTTTCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCg                                                                                                                                            >  2:2113844/1‑150 (MQ=255)
          gggAAAGATATTTTCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTATCTCTAACCCg                                                                                                                                            <  1:622187/150‑1 (MQ=255)
          gggAAAGATATTTTCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCGTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCg                                                                                                                                            >  2:1002846/1‑150 (MQ=255)
            gAAAGATATTTTCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCa                                                                                                                                          <  2:4478433/150‑1 (MQ=255)
              aaGATATTTTCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATa                                                                                                                                        >  1:391809/1‑150 (MQ=255)
              aaGATATTTTCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATa                                                                                                                                        >  1:3725129/1‑150 (MQ=255)
              aaGATATTTTCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATa                                                                                                                                        <  1:3002349/150‑1 (MQ=255)
               aGATATTTTCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTATCTCTAACCCGCATaa                                                                                                                                       >  1:4310198/1‑150 (MQ=255)
                  tatTTTCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTATCTCTAACCCGCATAAATa                                                                                                                                    <  2:536000/150‑1 (MQ=255)
                    ttttCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACt                                                                                                                                  <  2:2733200/150‑1 (MQ=255)
                        ccTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTGTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTa                                                                                                                              >  2:1592632/1‑150 (MQ=255)
                         cTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAg                                                                                                                             >  1:2646148/1‑150 (MQ=255)
                           gCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCa                                                                                                                           >  2:2340741/1‑150 (MQ=255)
                            cTAAGCCCGGATGCCGTTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTAATAAATTCATAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCAt                                                                                                                          >  1:1937250/1‑150 (MQ=255)
                             tAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATc                                                                                                                         <  1:3223972/150‑1 (MQ=255)
                              aaGCCCGGATGCCGGTGCTAAACCGGGTCATCATGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTATTGCGATGTTGTCTCTCTCTAACCCGCATAAATACTCGTAGCATCt                                                                                                                        >  1:1011568/1‑150 (MQ=255)
                                 cccGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTATCTCTAACCCGCATAAATACTGGTAGCATCTGCa                                                                                                                     <  1:667472/150‑1 (MQ=255)
                                    ggATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTc                                                                                                                  <  2:2308039/150‑1 (MQ=255)
                                       tGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAAc                                                                                                               >  1:4527395/1‑150 (MQ=255)
                                       tGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAAc                                                                                                               >  1:1110617/1‑150 (MQ=255)
                                           ggTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGa                                                                                                           >  2:4317295/1‑150 (MQ=255)
                                              gCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATaa                                                                                                        <  2:4287770/150‑1 (MQ=255)
                                               cTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATaaa                                                                                                       <  2:3046073/150‑1 (MQ=255)
                                                      gggTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAg                                                                                                >  1:3093060/1‑150 (MQ=255)
                                                      gggTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAg                                                                                                <  2:4222812/150‑1 (MQ=255)
                                                      gggTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAg                                                                                                <  1:2438495/150‑1 (MQ=255)
                                                       ggTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAgg                                                                                               >  2:3277111/1‑150 (MQ=255)
                                                       ggTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAgg                                                                                               <  1:1226238/150‑1 (MQ=255)
                                                           atcaGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATg                                                                                           <  2:2599941/150‑1 (MQ=255)
                                                                gTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGaa                                                                                      <  2:1818741/150‑1 (MQ=255)
                                                                gTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGaa                                                                                      >  2:829956/1‑150 (MQ=255)
                                                                gTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGaa                                                                                      >  2:830012/1‑150 (MQ=255)
                                                                gTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGaa                                                                                      <  2:1409063/150‑1 (MQ=255)
                                                                 tGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAAt                                                                                     >  1:2984756/1‑150 (MQ=255)
                                                                   aaaTAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATga                                                                                   >  1:4071531/1‑150 (MQ=255)
                                                                      taatCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGac                                                                                <  1:1252831/150‑1 (MQ=255)
                                                                       aatCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGaca                                                                               <  1:2338840/150‑1 (MQ=255)
                                                                       aatCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGaca                                                                               <  1:1260875/150‑1 (MQ=255)
                                                                       aatCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGaca                                                                               >  2:1357347/1‑150 (MQ=255)
                                                                          cccccTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACtt                                                                            >  1:4339285/1‑150 (MQ=255)
                                                                            cccTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTta                                                                          <  1:1092516/150‑1 (MQ=255)
                                                                             ccTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTtat                                                                         >  1:877814/1‑150 (MQ=255)
                                                                               ttCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTtatct                                                                       >  2:3275746/1‑150 (MQ=255)
                                                                               ttCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTtatct                                                                       >  1:395915/1‑150 (MQ=255)
                                                                               ttCAAGGCGCTGCATCGAAAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAATTACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTtatct                                                                       >  2:999055/1‑150 (MQ=255)
                                                                                 cAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTtatctat                                                                     >  2:622829/1‑150 (MQ=255)
                                                                                 cAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTtatctat                                                                     >  1:653017/1‑150 (MQ=255)
                                                                                  aaGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTtatctatc                                                                    <  2:2709441/150‑1 (MQ=255)
                                                                                     gcgcTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAgg                                                                 >  1:99705/1‑150 (MQ=255)
                                                                                       gcTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGAc                                                               >  2:1755799/1‑150 (MQ=255)
                                                                                       gcTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGAc                                                               >  1:1315361/1‑150 (MQ=255)
                                                                                          gCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGaa                                                            <  2:3035500/150‑1 (MQ=255)
                                                                                          gCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGaa                                                            <  2:3035174/150‑1 (MQ=255)
                                                                                           cATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGaaa                                                           >  2:1608694/1‑150 (MQ=255)
                                                                                            aTCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGaaaa                                                          <  1:2129165/150‑1 (MQ=255)
                                                                                            aTCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGaaaa                                                          <  2:928476/150‑1 (MQ=255)
                                                                                              cGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAAt                                                        <  1:2454195/150‑1 (MQ=255)
                                                                                               gACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATc                                                       >  1:3454846/1‑150 (MQ=255)
                                                                                               gACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGCTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATc                                                       >  1:1558910/1‑150 (MQ=255)
                                                                                                 cAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATcac                                                     <  2:2517043/150‑1 (MQ=255)
                                                                                                 cAGCGCCTTTTCTTTATAAAGTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATcac                                                     <  1:3583331/150‑1 (MQ=255)
                                                                                                  aGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATcaca                                                    <  2:4107558/150‑1 (MQ=255)
                                                                                                  aGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATcaca                                                    <  2:1856496/150‑1 (MQ=255)
                                                                                                   gcgcCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACAt                                                   >  1:272774/1‑150 (MQ=255)
                                                                                                        ttttCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAAtt                                              <  2:3335443/150‑1 (MQ=255)
                                                                                                        ttttCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAAtt                                              >  1:1396443/1‑150 (MQ=255)
                                                                                                        ttttCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATATGCATTAAACTGGATAAAATTAAAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAAtt                                              >  2:4223152/1‑150 (MQ=255)
                                                                                                               tataAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCa                                       >  2:22144/1‑150 (MQ=255)
                                                                                                               tataAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCa                                       >  1:4212043/1‑150 (MQ=255)
                                                                                                                ataAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAg                                      <  1:2226829/150‑1 (MQ=255)
                                                                                                                  aaaTTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGtt                                    <  2:1890093/150‑1 (MQ=255)
                                                                                                                  aaaTTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGtt                                    <  2:3725129/150‑1 (MQ=255)
                                                                                                                    aTTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGa                                  >  2:96103/1‑150 (MQ=255)
                                                                                                                     ttCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAg                                 <  2:4397816/150‑1 (MQ=255)
                                                                                                                       ccTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCa                               >  2:2041770/1‑150 (MQ=255)
                                                                                                                             gttgttTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATTCTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAggg                         <  2:360334/150‑1 (MQ=255)
                                                                                                                               tgttTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGaa                       <  2:1054168/150‑1 (MQ=255)
                                                                                                                                gttTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGaaa                      <  2:145702/150‑1 (MQ=255)
                                                                                                                                 ttttCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAAc                     >  1:2871112/1‑150 (MQ=255)
                                                                                                                                     cTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGt                 <  2:1982908/150‑1 (MQ=255)
                                                                                                                                      ttGCGATTTTGTCTCTCTCTAACCCGCATATATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGtt                >  1:4461038/1‑150 (MQ=255)
                                                                                                                                         cGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCa             <  1:1682530/150‑1 (MQ=255)
                                                                                                                                         cGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCa             <  2:3203662/150‑1 (MQ=255)
                                                                                                                                          gATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCat            >  2:3454854/1‑150 (MQ=255)
                                                                                                                                           aTTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCata           >  2:3553410/1‑150 (MQ=255)
                                                                                                                                           aTTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCata           <  2:3409631/150‑1 (MQ=255)
                                                                                                                                            ttttGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCatat          <  2:818148/150‑1 (MQ=255)
                                                                                                                                             tttGTCTCTCTCTAACCCTCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCATATc         <  2:1499585/150‑1 (MQ=255)
                                                                                                                                               tGTCTATCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCATATCAg       >  1:3821973/1‑150 (MQ=255)
                                                                                                                                               tGTCTATCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTAAAGTGATGCAGAATGAGACACTTTATATATCAGTACTAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCATATCAg       >  2:4009722/1‑150 (MQ=255)
                                                                                                                                                gTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCATATCAGt      >  1:3643354/1‑150 (MQ=255)
                                                                                                                                                gTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCATATCAGt      >  2:706125/1‑150 (MQ=255)
                                                                                                                                                 tctctctcTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCATATCAGtt     <  2:3354995/150‑1 (MQ=255)
                                                                                                                                                  ctctctctAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCATATCAGTTg    >  1:4470694/1‑150 (MQ=255)
                                                                                                                                                  ctctctctAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCATATCAGTTg    >  1:2035920/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGGTCCTGGCGGGAAAGATATTTTCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCATATCAGTTGGG  >  CP009273/2189740‑2190037

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: