Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,066,133 T→G 94.9% S267R (AGC→CGC)  yihU ← gamma‑hydroxybutyrate dehydrogenase, NADH‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,066,1330TG94.9% 1276.5 / 29.7 340S267R (AGC→CGC) yihUgamma‑hydroxybutyrate dehydrogenase, NADH‑dependent
Reads supporting (aligned to +/- strand):  ref base T (10/8);  new base G (160/162);  total (170/170)
Fisher's exact test for biased strand distribution p-value = 8.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 342 total aligned reads displayed.

ATGTCGTTGATGGTGTACTTATTCATTCATTTATTCCTTTATCCAGTCGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAG  >  CP009273/4065986‑4066279
                                                                                                                                                   |                                                                                                                                                  
   tCGTTGATGGTGTACTTATTCATTCATTTATTCCTTTATCCAGTCGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATaaa                                                                                                                                               <  1:4514641/150‑1 (MQ=255)
     gTTGATGGTGTACTTATTCATTCATTTATTCCTTTATCCAGTCGTTACCATTTTAATTTTGCAGTCCTCCCGGCACAGACACGGACCTGTTCCAGAATGGCGGTCCAGTCCAGGCGACCGCGACCCCCGGCGCGCGCCAGGCGATAAAcc                                                                                                                                             >  2:2104776/1‑150 (MQ=255)
     gTTGATGGTGTACTTATTCATTCATTTATTCCTTTATCCAGTCGTTACATTTTTACTTTTGCAGTCATCCCGGCACAGACACGGACCAGTTCCAGAATTGCGGAACACTCCAGGCGACCGCGACCCGCGGCGCGCGCCTGGCGATAAAcc                                                                                                                                             >  2:2190205/1‑150 (MQ=255)
       tGATGGTGTACTTATTCATTCATTTATTCCTTTATCCAGTCGTGACAATTTTAATTTTGGAGTCATACCCGCACAGACACGGACCGGTTCCAGAAAGGCGGACCAATCCAGGCGACCCCGACCCCCCGCGCGCGCCAGGCGATAAAccac                                                                                                                                           >  1:4659922/1‑148 (MQ=255)
         aTGGTGTACTTATTCATTCATTTATTCCTTTATCCAGTCGTTACATTTTTAATTTGGCAGTCATCCCGGCACAGACACGGACCAGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCCCAGCGCGCGCCTGGCGATAAAACAccc                                                                                                                                         >  2:1470511/1‑150 (MQ=255)
         aTGGTGTACTTAGTCATTCATTTATTCCTTTATCCAGTCGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGACTGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTccc                                                                                                                                         <  2:4711176/150‑1 (MQ=255)
          tGGTGTACTTATTCATTCATTTATTCCTTTATCCAGTCGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCg                                                                                                                                        <  2:2624614/150‑1 (MQ=255)
          tGGTGTACTTATTCATTCATTTATTCCTTTATCCAGTCGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCg                                                                                                                                        <  1:171254/150‑1 (MQ=255)
            gtgtACTTATTCATTCATTTATTCCTTTATCCAGTCGTTACAATTTTAATTTTGCAGTCCACCCCGCACAGACACGGACCAGTTCCAGAAAGGCGGACCAGTACAGGCCACCGCCACCCCCAGCGCGCGCCAGGCGATAAAACAccccag                                                                                                                                      >  1:1808084/1‑147 (MQ=255)
                 cTTATTCATTCATTTATTCCTTTATCCAGTCGTTACAATTTTAAATTTGCAGTCAAACCGGCACAGACACCGACCAGTTCCAGAAAGGCGGACCACTCCAGGCGACCGCCACCCCCTGCGCGCGCCGGGCGATAAAACACCCCAGAggcg                                                                                                                                 >  1:1887110/1‑150 (MQ=255)
                  ttATTCATTCATTTATTCCTTTATCCAGTCGTTACATTTTTAAATTTGCAGTCAAACCCGCACAGACAACGACCAGTTCCAGAAAGGCGGACCACTCCAGGCGACCCCCACCCCCGGCGCGCGCCAGGCGATAAAACACCCCTGAtgcgg                                                                                                                                >  1:1564505/1‑150 (MQ=255)
                  ttATTCATTCATTTATTCCTTTATCCAGTCGTGACAATTTTAATTTTGCAGTCAAACCCGCACAGAAACCGACCAGTTACACAAATGCGGACCAATCCAGGCGAACCCGACCCCCTGCGCGCGCCAGGCGATAAACCACCCCTGAggcgg                                                                                                                                >  2:1200034/1‑150 (MQ=255)
                    attcattcATTTATTCCTTTATCCAGTCGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGcc                                                                                                                              <  2:1728239/150‑1 (MQ=255)
                         ttcattTATTCCTTTATCCAGTCGTTACATTTTTTCTTTGGCAGTCATACCCGCACAGACACCGACCTGTTCCAGAATTGCGGACCAGTCCTGGCGACCCCGACCCGCAGCGCGCGCCTGGCGATAAACCTCCCCGGAGGCGGCCCCCac                                                                                                                         >  2:3720317/1‑149 (MQ=255)
                         ttcattTATTCCTTTATCCAGTCGTTACATTTTTAATTTTGCAGTCAAACCCGCACTGACACCGACCAGTTCCAGAATTGCGGACCCGTCCAGGCGACCGCGACCCCCTGCGCGCGCCCGGCGATAAAACACCCCAGAGGCGGccccccg                                                                                                                         >  1:1085484/1‑148 (MQ=255)
                         ttcattTATTCCTTTATCCAGTCGTTACATTTTTAATTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATTGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAg                                                                                                                         >  2:2327863/1‑150 (MQ=255)
                         ttcattTATTCCTTTATCCAGTCGTGACATTTTTAAATTGGCCGTCAAACCGGCACAGACACGGACCAGTTCCAGAATTGCGGACCAGTACAGGCGACCCCCACCCCCGGCGCGCGCCAGGCGCTAAACCACCCCTGAGGCGGCCCCCAg                                                                                                                         >  2:3330095/1‑150 (MQ=255)
                          tcattTATTCCTTTATCCAGTCGTGTCATTTTTTATTTTGCAGTCAACCCCGCACAGACACCGACCGGTTACACAATGGCGGACCACTCCAGGCGACCCCGCCCCCCCGCGCCCGCCCGGCGCTAAAACAACCCGGAGGCGGcccccccc                                                                                                                        >  1:1276252/1‑147 (MQ=255)
                           cattTATTCCTTTATCCAGTCGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCg                                                                                                                       <  2:2599503/150‑1 (MQ=255)
                           cattTATTCCTTTATCCAGTCGTTACATTTTTAATTTTGCAGTCATACCCGCACAGAAACCGAACAGTTACAGAAAGGCGGGCCACTACTGGCGACCGCGACCCCCCGCGCGCGCCAGGCGATAAAACACCCCTGAGGCGGCCCCCcgcg                                                                                                                       >  2:479382/1‑150 (MQ=255)
                            attTATTCCTTTATCCAGTCGTTTCATTTTTACTTTGGCAGGACTACCCGCACAGACACCGACCAGGTCCAGAATGGCGGACCAGTCCTGGCGCCCGCGACCCCCCGCGCGCGCCAGGCGATAAAACGCCCCGGAGGCGGCCCCCACCgg                                                                                                                      >  2:4575125/1‑150 (MQ=255)
                             tttATTCCTTTATCCAGTCGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACAAGGACCGGTTCAAGAATTGCGGACCAGTCCATGCGACCCCGACCCCCAGCGCGCGCCAGGCGATAAACCACCCCTGAGGCGGCCCCCCCCGGc                                                                                                                     >  2:1336280/1‑150 (MQ=255)
                                  tCCTTTATCCAGTCGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTgg                                                                                                                <  1:1102285/150‑1 (MQ=255)
                                  tCCTTTATCCAGTCGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTgg                                                                                                                <  1:2682650/150‑1 (MQ=255)
                                          ccAGTCGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCa                                                                                                        <  2:2024213/150‑1 (MQ=255)
                                            aGTCGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGc                                                                                                      <  2:3179664/150‑1 (MQ=255)
                                              tCGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTg                                                                                                    >  2:3098345/1‑150 (MQ=255)
                                               cGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTgg                                                                                                   >  1:3416972/1‑150 (MQ=255)
                                               cGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTgg                                                                                                   <  1:953948/150‑1 (MQ=255)
                                               cGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTgg                                                                                                   >  1:101240/1‑150 (MQ=255)
                                                gTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGt                                                                                                  <  1:4281892/150‑1 (MQ=255)
                                                 ttACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGtt                                                                                                 >  2:1842232/1‑150 (MQ=255)
                                                  tACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTg                                                                                                >  1:2376065/1‑150 (MQ=255)
                                                  tACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTg                                                                                                <  1:119530/150‑1 (MQ=255)
                                                  tACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTg                                                                                                >  1:4454255/1‑150 (MQ=255)
                                                        tttACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACa                                                                                          <  2:2134733/150‑1 (MQ=255)
                                                        tttACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACa                                                                                          <  1:2432518/150‑1 (MQ=255)
                                                        tttACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACa                                                                                          >  2:681333/1‑150 (MQ=255)
                                                         ttACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACAt                                                                                         >  1:4414668/1‑150 (MQ=255)
                                                            cTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCaa                                                                                      >  1:2187943/1‑150 (MQ=255)
                                                            cTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCaa                                                                                      >  2:626458/1‑150 (MQ=255)
                                                              ttGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAgg                                                                                    >  2:797778/1‑150 (MQ=255)
                                                              ttGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAgg                                                                                    >  2:1606787/1‑150 (MQ=255)
                                                                ggCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGc                                                                                  <  1:477569/150‑1 (MQ=255)
                                                                ggCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGc                                                                                  <  1:4003098/150‑1 (MQ=255)
                                                                   aGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGAt                                                                               <  2:3840575/150‑1 (MQ=255)
                                                                      cATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGcc                                                                            <  2:2504629/150‑1 (MQ=255)
                                                                       aTCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCa                                                                           <  1:4131596/150‑1 (MQ=255)
                                                                       aTCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCa                                                                           <  1:200502/150‑1 (MQ=255)
                                                                       aTCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCa                                                                           >  1:2749990/1‑150 (MQ=255)
                                                                            ggCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATc                                                                      >  1:2646383/1‑150 (MQ=255)
                                                                            ggCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATc                                                                      >  2:1892116/1‑150 (MQ=255)
                                                                              cACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCt                                                                    <  2:3040923/150‑1 (MQ=255)
                                                                              cACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCt                                                                    <  2:2766107/150‑1 (MQ=255)
                                                                               aCTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGATGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCtt                                                                   >  1:84353/1‑150 (MQ=255)
                                                                                cTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTa                                                                  >  1:615802/1‑150 (MQ=255)
                                                                                 tGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTAt                                                                 <  2:736297/150‑1 (MQ=255)
                                                                                 tGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTAt                                                                 >  2:4074786/1‑150 (MQ=255)
                                                                                    cacGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATggg                                                              >  2:4587282/1‑150 (MQ=255)
                                                                                    cacGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATggg                                                              >  2:153897/1‑150 (MQ=255)
                                                                                     acGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGc                                                             >  1:3539367/1‑150 (MQ=255)
                                                                                      cGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCa                                                            <  1:3681940/150‑1 (MQ=255)
                                                                                         aCCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCACga                                                         >  2:1031030/1‑150 (MQ=255)
                                                                                           cTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAgatc                                                       <  2:3375819/150‑1 (MQ=255)
                                                                                            tGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAgatcg                                                      <  1:3860660/150‑1 (MQ=255)
                                                                                               tCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAgatcgatc                                                   <  2:1098699/150‑1 (MQ=255)
                                                                                                 cAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCAt                                                 >  2:4368141/1‑150 (MQ=255)
                                                                                                 cAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCAt                                                 <  1:2584887/150‑1 (MQ=255)
                                                                                                    aaTGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGaa                                              <  2:3483085/150‑1 (MQ=255)
                                                                                                     aTGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAg                                             <  1:1470511/150‑1 (MQ=255)
                                                                                                      tGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAgg                                            >  2:2540390/1‑150 (MQ=255)
                                                                                                      tGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAgg                                            <  2:184780/150‑1 (MQ=255)
                                                                                                          ggACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCACGATCGATCATGAAGGCggg                                        >  1:2254010/1‑150 (MQ=255)
                                                                                                          ggACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCggg                                        <  2:3937623/150‑1 (MQ=255)
                                                                                                           gACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGga                                       <  2:1564955/150‑1 (MQ=255)
                                                                                                           gACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGga                                       <  1:3694445/150‑1 (MQ=255)
                                                                                                           gACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGga                                       <  1:13171/150‑1 (MQ=255)
                                                                                                           gACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGga                                       <  1:1300763/150‑1 (MQ=255)
                                                                                                            aCCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGgag                                      <  1:4078844/150‑1 (MQ=255)
                                                                                                             ccAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGgaga                                     <  1:2024450/150‑1 (MQ=255)
                                                                                                              cAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCACGATCGATCATGAAGGCGGGAGaa                                    >  2:2974895/1‑150 (MQ=255)
                                                                                                                gTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTAGGGGCAAGATCGATCATGAAGGCGGGAGAAAg                                  >  2:4755211/1‑150 (MQ=255)
                                                                                                                 tCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGTAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGa                                 <  1:2819975/150‑1 (MQ=255)
                                                                                                                  ccTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGAt                                <  1:1842927/150‑1 (MQ=255)
                                                                                                                   cTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGAtc                               >  2:2816528/1‑150 (MQ=255)
                                                                                                                    tGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGAtct                              <  1:3053293/150‑1 (MQ=255)
                                                                                                                    tGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGAtct                              <  1:681333/150‑1 (MQ=255)
                                                                                                                        gaccgcgaccCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGTGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCACGATCGATCATGAAGGCGGGAGAAAGATCTCCGc                          >  2:4391674/1‑150 (MQ=255)
                                                                                                                        gaccgcgaccCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGc                          >  2:112917/1‑150 (MQ=255)
                                                                                                                         accgcgaccCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCt                         <  1:3072878/150‑1 (MQ=255)
                                                                                                                         accgcgaccCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCt                         <  2:3647285/150‑1 (MQ=255)
                                                                                                                          ccgcgaccCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTg                        <  1:1576124/150‑1 (MQ=255)
                                                                                                                           cgcgaccCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGa                       <  1:3898771/150‑1 (MQ=255)
                                                                                                                            gcgaccCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAg                      <  2:3588257/150‑1 (MQ=255)
                                                                                                                              gaccCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGa                    >  2:4264159/1‑150 (MQ=255)
                                                                                                                                cccGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACt                  >  1:1673597/1‑150 (MQ=255)
                                                                                                                                           cgcCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCa       >  2:2866260/1‑150 (MQ=255)
                                                                                                                                           cgcCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCa       >  1:4204980/1‑150 (MQ=255)
                                                                                                                                             ccTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAgg     >  1:309733/1‑150 (MQ=255)
                                                                                                                                                ggCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAg  <  1:3249423/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ATGTCGTTGATGGTGTACTTATTCATTCATTTATTCCTTTATCCAGTCGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAG  >  CP009273/4065986‑4066279

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: ‑