Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,189,888 C→A 17.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,8880CA17.6% 825.4 / 116.3 275intergenic (+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 (140/86);  new base A (20/29);  total (160/115)
Fisher's exact test for biased strand distribution p-value = 1.01e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 282 total aligned reads displayed.

CCGGTCCTGGCGGGAAAGATATTTTCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCATATCAGTTG  >  CP009273/2189739‑2190035
                                                                                                                                                     |                                                                                                                                                   
ccGGTCCTGGCGGGAAAGATATTTTCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGtctc                                                                                                                                                     >  2:1750105/1‑150 (MQ=255)
ccGGTCCTGGCGGGAAAGATATTTTCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGtctc                                                                                                                                                     <  2:325786/150‑1 (MQ=255)
      cTGGCGGGAAAGATATTTTCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTa                                                                                                                                               <  1:4045251/150‑1 (MQ=255)
        ggCGGGAAAGATATTTTCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAAc                                                                                                                                             <  2:566761/150‑1 (MQ=255)
          cGGGAAAGATATTTTCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCTAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAAccc                                                                                                                                           >  2:4710375/1‑150 (MQ=255)
          cGGGAAAGATATTTTCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAAccc                                                                                                                                           >  1:5139276/1‑150 (MQ=255)
          cGGGAAAGATATTTTCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAAccc                                                                                                                                           >  2:713768/1‑150 (MQ=255)
           gggAAAGATATTTTCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTATCTCTAACCCg                                                                                                                                          >  2:3726068/1‑150 (MQ=255)
            ggAAAGATATTTTCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGc                                                                                                                                         >  1:2263649/1‑150 (MQ=255)
            ggAAAGATATTTTCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGc                                                                                                                                         >  1:2789369/1‑150 (MQ=255)
               aaGATATTTTCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATa                                                                                                                                      >  2:3925549/1‑150 (MQ=255)
               aaGATATTTTCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATa                                                                                                                                      >  1:2922431/1‑150 (MQ=255)
                     ttttCCTGCTAAGCCCGGATGCCGGTTCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAAGTCCTAAAGTTGTTTTCTTGCGATTTTGTCTATCTCTAACCCGCATAAATACt                                                                                                                                <  1:571298/150‑1 (MQ=255)
                     ttttCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTATCTCTAACCCGCATAAATACt                                                                                                                                <  1:2157720/150‑1 (MQ=255)
                      tttCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTGTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTACCCCGCATAAATACTg                                                                                                                               <  2:3083076/150‑1 (MQ=255)
                        tCCTGCTTAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGt                                                                                                                             <  2:2532363/150‑1 (MQ=255)
                         ccTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTa                                                                                                                            <  2:4573778/150‑1 (MQ=255)
                         ccTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTa                                                                                                                            <  2:401639/150‑1 (MQ=255)
                            gCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCa                                                                                                                         >  2:735824/1‑150 (MQ=255)
                             cTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCAt                                                                                                                        >  2:424186/1‑150 (MQ=255)
                             cTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGGTGTTGTCTTGCGTTTTTGTCTCTCTCTAACCCGCATCACTACGGGTAGctt                                                                                                                        >  2:4408699/1‑148 (MQ=255)
                                aGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTg                                                                                                                     >  1:1410679/1‑150 (MQ=255)
                                aGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTg                                                                                                                     >  1:54378/1‑150 (MQ=255)
                                 gCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGc                                                                                                                    <  2:1600762/150‑1 (MQ=255)
                                    cGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCagt                                                                                                                 >  2:1124416/1‑148 (MQ=255)
                                     ggATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTc                                                                                                                >  2:4134743/1‑150 (MQ=255)
                                                 tAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATaaaa                                                                                                    >  1:3663089/1‑150 (MQ=255)
                                                  aaaCCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATACAAt                                                                                                   >  2:2973952/1‑150 (MQ=255)
                                                  aaaCCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAAt                                                                                                   <  1:4158169/150‑1 (MQ=255)
                                                  aaaCCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAAt                                                                                                   <  1:5074995/150‑1 (MQ=255)
                                                    aCCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTa                                                                                                 >  1:317030/1‑150 (MQ=255)
                                                    aCCGGGTCATCAGGAGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTa                                                                                                 <  2:3412288/150‑1 (MQ=255)
                                                     ccGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGAGAAAATTAc                                                                                                >  1:1593988/1‑150 (MQ=255)
                                                      cGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACa                                                                                               <  2:4077311/150‑1 (MQ=255)
                                                      cGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACa                                                                                               >  1:1469141/1‑150 (MQ=255)
                                                       gggTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTATCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAg                                                                                              >  1:3042283/1‑150 (MQ=255)
                                                       gggTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTATCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAg                                                                                              <  2:2198966/150‑1 (MQ=255)
                                                        ggTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAgg                                                                                             >  2:1371376/1‑150 (MQ=255)
                                                        ggTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAgg                                                                                             >  2:1544299/1‑150 (MQ=255)
                                                         gTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAggg                                                                                            >  1:3309325/1‑150 (MQ=255)
                                                         gTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAggg                                                                                            >  2:630283/1‑150 (MQ=255)
                                                         gTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTATCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAggg                                                                                            >  2:4941780/1‑150 (MQ=255)
                                                           catcaGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGAt                                                                                          <  2:2738670/150‑1 (MQ=255)
                                                           catcaGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGAt                                                                                          <  2:3582725/150‑1 (MQ=255)
                                                           catcaGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGAt                                                                                          <  1:1922856/150‑1 (MQ=255)
                                                             tcaGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGc                                                                                        <  1:846645/150‑1 (MQ=255)
                                                               aGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAg                                                                                      <  1:2588048/150‑1 (MQ=255)
                                                               aGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAg                                                                                      >  1:3661350/1‑150 (MQ=255)
                                                                ggTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGa                                                                                     <  2:2081406/150‑1 (MQ=255)
                                                                   gAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTATCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATg                                                                                  <  1:3168185/150‑1 (MQ=255)
                                                                      ataatCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATgaga                                                                               >  2:4751580/1‑150 (MQ=255)
                                                                           cccccTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACtt                                                                          <  2:560646/150‑1 (MQ=255)
                                                                            ccccTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGAAGAATCTGCAATCAACTGGATAATATTAAAGGGATGCAGAATCAGACACttt                                                                         >  1:4910680/1‑150 (MQ=255)
                                                                              ccTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTtat                                                                       <  2:1359498/150‑1 (MQ=255)
                                                                                ttCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTtatct                                                                     >  1:2287207/1‑150 (MQ=255)
                                                                                ttCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTtatct                                                                     <  2:4464684/150‑1 (MQ=255)
                                                                                ttCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTATCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTtatct                                                                     <  2:3715298/150‑1 (MQ=255)
                                                                                 tCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTtatcta                                                                    <  2:1236590/150‑1 (MQ=255)
                                                                                      gcgcTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAgg                                                               <  1:779455/150‑1 (MQ=255)
                                                                                        gcTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGAc                                                             >  2:4261749/1‑150 (MQ=255)
                                                                                         cTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTATCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACg                                                            >  2:2441622/1‑150 (MQ=255)
                                                                                          tGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGa                                                           >  2:3634200/1‑150 (MQ=255)
                                                                                           gCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGaa                                                          <  2:2355723/150‑1 (MQ=255)
                                                                                           gCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGaa                                                          <  2:2355597/150‑1 (MQ=255)
                                                                                           gCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGaa                                                          >  1:3998803/1‑150 (MQ=255)
                                                                                             aTCGACAGCGCCATTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGaaaa                                                        <  1:4235251/150‑1 (MQ=255)
                                                                                              tCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTATCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGaaaaa                                                       <  2:3810855/150‑1 (MQ=255)
                                                                                                 aCAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATca                                                    >  2:873585/1‑150 (MQ=255)
                                                                                                  cAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGTATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATcac                                                   >  2:3221590/1‑150 (MQ=255)
                                                                                                   aGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATcaca                                                  >  2:4766040/1‑150 (MQ=255)
                                                                                                   aGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATATATCAGGACGAAAAATcaca                                                  >  2:1715732/1‑150 (MQ=255)
                                                                                                    gcgcCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACAt                                                 <  1:1683480/150‑1 (MQ=255)
                                                                                                     cgcCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTATCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATa                                                <  1:4732222/150‑1 (MQ=255)
                                                                                                       ccTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATaaa                                              <  1:3528076/150‑1 (MQ=255)
                                                                                                       ccTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTATCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATaaa                                              >  2:932011/1‑150 (MQ=255)
                                                                                                         ttttCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAAtt                                            >  1:4413734/1‑150 (MQ=255)
                                                                                                         ttttCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTATCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAAtt                                            <  1:1008651/150‑1 (MQ=255)
                                                                                                          tttCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTc                                           >  2:1981322/1‑150 (MQ=255)
                                                                                                            tCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTATCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAg                                         >  1:1142732/1‑150 (MQ=255)
                                                                                                                      ttCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAg                               <  1:1373370/150‑1 (MQ=255)
                                                                                                                       tcctccaGTTGTTTTCGTGCGATTTTGTCTCTCTCTCACCCGCATAAATACTGGTAGCATCTGCCTTCAACTGGATAACATTACAGGGATGCAGAATGAGACACTTTATCTATCTGGACGAAAAATCACATAAATTCAGGGCAGTTGAGc                              <  2:3358366/144‑1 (MQ=255)
                                                                                                                            aaGTTGTTTTCTTGCGATTTTGTCTATCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAg                         >  1:2165334/1‑150 (MQ=255)
                                                                                                                            aaGTTGTTTTCTTGCGATTGTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAg                         <  1:1971766/150‑1 (MQ=255)
                                                                                                                             aGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAgg                        >  1:4109421/1‑150 (MQ=255)
                                                                                                                             aGTTGTTTTCTTGCGATTTTGTCTATCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAgg                        >  2:2696223/1‑150 (MQ=255)
                                                                                                                                tgttTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGaa                     >  1:2991177/1‑150 (MQ=255)
                                                                                                                                 gttTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGaaa                    >  1:983153/1‑150 (MQ=255)
                                                                                                                                   tttCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCACCAGGGAAACg                  >  2:3554696/1‑150 (MQ=255)
                                                                                                                                   tttCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACg                  >  1:4495149/1‑150 (MQ=255)
                                                                                                                                   tttCTTGCGATTTTGTCTATCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACg                  <  1:995252/150‑1 (MQ=255)
                                                                                                                                     tCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAg                >  1:3829201/1‑150 (MQ=255)
                                                                                                                                     tCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAg                >  1:3838051/1‑150 (MQ=255)
                                                                                                                                     tCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAg                >  1:1976971/1‑150 (MQ=255)
                                                                                                                                       ttGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGtt              <  1:2929769/150‑1 (MQ=255)
                                                                                                                                       ttGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGtt              <  2:1427130/150‑1 (MQ=255)
                                                                                                                                          cGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCa           <  1:163922/150‑1 (MQ=255)
                                                                                                                                            aTTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCata         >  1:3575024/1‑150 (MQ=255)
                                                                                                                                            aTTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCata         >  1:468581/1‑150 (MQ=255)
                                                                                                                                             ttttGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCatat        <  1:4320616/150‑1 (MQ=255)
                                                                                                                                                tGTCTATCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCATATCAg     >  1:2960052/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CCGGTCCTGGCGGGAAAGATATTTTCCTGCTAAGCCCGGATGCCGGTGCTAAACCGGGTCATCAGGTGAAATAATCCCCCTTCAAGGCGCTGCATCGACAGCGCCTTTTCTTTATAAATTCCTAAAGTTGTTTTCTTGCGATTTTGTCTCTCTCTAACCCGCATAAATACTGGTAGCATCTGCATTCAACTGGATAAAATTACAGGGATGCAGAATGAGACACTTTATCTATCAGGACGAAAAATCACATAAATTCAGGGCAGTTGAGCAACAGGGAAACGAGTTGCATATCAGTTG  >  CP009273/2189739‑2190035

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: