Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,066,133 T→G 100% 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,1330TG97.0% 1021.6 / 7.7 262S267R (AGC→CGC) yihUgamma‑hydroxybutyrate dehydrogenase, NADH‑dependent
Reads supporting (aligned to +/- strand):  ref base T (1/7);  new base G (151/103);  total (152/110)
Fisher's exact test for biased strand distribution p-value = 1.06e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 268 total aligned reads displayed.

ATGTCGTTGATGGTGTACTTATTCATTCATTTATTCCTTTATCCAGTCGTT‑ACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAA  >  CP009273/4065986‑4066278
                                                                                                                                                    |                                                                                                                                                 
aTGTCGTTGATGGTGTACTTATTCATTCATTTATTCCTTTATCCAGTCGTT‑ACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCgat                                                                                                                                                 <  2:2260218/150‑1 (MQ=255)
 tGTCGTTGATGGTGTACTTATTCATTCATTTATTCCTTTATCCAGTCGTT‑ACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCgata                                                                                                                                                >  2:3205814/1‑150 (MQ=255)
 tGTCGTTGATGGTGTACTTATTCATTCATTTATTCCTTTATCCAGTCGTT‑ACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCgata                                                                                                                                                >  1:1535027/1‑150 (MQ=255)
    cGTTGATGGTGTACTTATTCATTCATTTATTCCTTTATCCAGTCGTT‑ACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAAc                                                                                                                                             >  2:4122127/1‑150 (MQ=255)
     gTTGATGGTGTACTTATTCATTCATTTATTCCTTTATCCAGTCGAT‑ACAATTTAAATTTGTAATACAACCCGGCAATGAAACGGACCAGTTCCAAAATGGCGGACCAGTACTGGCGACCCCGAACCCCAGCGCGCGCCTGGCGATAAAcc                                                                                                                                            >  2:4375182/1‑150 (MQ=255)
      ttGATGGTGTACTTATTCATTCATTTATTCCTTTATCCAGTCGTTAAAATTTTAAATTTGGC‑CTAAAAACGCCACTGAAACAGAACTGTTACAGAAAGGCGGAACAATACTGGAGACCACGACCCCCAGAGCGCGCCATGCGATAAAcca                                                                                                                                           >  1:1667235/1‑149 (MQ=255)
      ttGATGGTGTACTTATTCATTCATTTATTCCTTTATCCAGTCGTT‑ACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCt                                                                                                                                           <  1:2878299/150‑1 (MQ=255)
           ggTGTACTTATTCATTCATTTATTCCTTTATCAAGTCGTT‑ACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGt                                                                                                                                      <  1:642614/150‑1 (MQ=255)
            gtgtACTTATTCATTCATTTATTCCTTTATCCAGTCGTT‑ACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGTGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTg                                                                                                                                     >  2:2275473/1‑150 (MQ=255)
              gtACTTATTCATTCATTTATTCCTTTATCCAGTCGTT‑TCATTTTTACATTTGCAGTAAACCCCGCACAGACACCGACCAGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCCCGGCGCGCGCCAGGCGATAAACCACCCGTGAg                                                                                                                                   >  1:629603/1‑150 (MQ=255)
               tACTTATTCATTCATTTATTCCTTTATCCAGTCGTT‑ACAATTTTACTTTTGCAGTCATACCCGCACAGACACCGACCAGTTCCAGAAATGCGGACCACTACAGGCGACCCCGACCCCCCGCGCGCCCCAGGCGAAAAACCACCCCTGAgg                                                                                                                                  >  1:796558/1‑150 (MQ=255)
                aCTTATTCATTCATTTATTCCTTTATCCAGTCGTT‑ACATTTTTAAATTTGCAGTCATCCCCGCACAGACACCGACCGGTTCCAGAAATGCGGACCAGTCCAGGCGACCGCGACCCCCTGCGCGCGCCAGGCGTTAAAACACCCCTGAggc                                                                                                                                 >  1:4446388/1‑150 (MQ=255)
                  ttATTCATTCATTTATTCCTTTATCCAGTCGTT‑ACATTTTTACTTTGGCAGGCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAggcgg                                                                                                                               <  2:1177798/150‑1 (MQ=255)
                    attcattcATTTATTCCTTTATCCAGTCGTT‑ACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGcc                                                                                                                             >  2:4182901/1‑150 (MQ=255)
                      tcattcatTTATTCCTTTATCCAGTCGTT‑ACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGcccc                                                                                                                           >  2:716160/1‑150 (MQ=255)
                        attcattTATTCCTTTATCCAGTCGTT‑ACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCa                                                                                                                         <  2:23967/150‑1 (MQ=255)
                         ttcattTATTCCTTTATCCAGTCGTT‑ACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAg                                                                                                                        >  2:3441235/1‑150 (MQ=255)
                         ttcattTATTCCTTTATCCAGTCGTT‑ACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCGGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAg                                                                                                                        >  2:2612182/1‑150 (MQ=255)
                          tcattTAGTCCTTTATCCAGTCGTT‑ACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGc                                                                                                                       <  1:3923891/150‑1 (MQ=255)
                             tttATTCCTTTATCCAGTCGTT‑ACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGc                                                                                                                    <  1:1147068/150‑1 (MQ=255)
                             tttATTCCTTTATCCAGTCGTT‑ACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGc                                                                                                                    <  1:67493/150‑1 (MQ=255)
                                       tATCCAGTCGTT‑ACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACAt                                                                                                          <  2:3372083/150‑1 (MQ=255)
                                         tCCAGTCGTT‑ACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGc                                                                                                        >  1:1888416/1‑150 (MQ=255)
                                         tCCAGTCGTT‑ACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGc                                                                                                        >  2:4555669/1‑150 (MQ=255)
                                             gTCGTT‑ACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCt                                                                                                    <  1:2995157/150‑1 (MQ=255)
                                             gTCGTT‑ACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCt                                                                                                    <  2:1568188/150‑1 (MQ=255)
                                               cGTT‑ACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGACTGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTgg                                                                                                  >  2:3812614/1‑150 (MQ=255)
                                                gTT‑ACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCAGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGt                                                                                                 >  1:3447469/1‑150 (MQ=255)
                                                    aCATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTgg                                                                                              >  1:452462/1‑150 (MQ=255)
                                                     cATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGc                                                                                             >  1:3596948/1‑150 (MQ=255)
                                                     cATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGc                                                                                             >  1:1007311/1‑150 (MQ=255)
                                                      aTTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCg                                                                                            <  1:2534405/150‑1 (MQ=255)
                                                        ttttACTTTGTCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGAc                                                                                          <  1:2630097/150‑1 (MQ=255)
                                                           tACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATc                                                                                       >  2:1504011/1‑150 (MQ=255)
                                                           tACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATc                                                                                       >  1:1169008/1‑150 (MQ=255)
                                                            aCTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCa                                                                                      >  1:3028812/1‑150 (MQ=255)
                                                            aCTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCa                                                                                      >  1:3364579/1‑150 (MQ=255)
                                                            aCTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCa                                                                                      <  1:1743659/150‑1 (MQ=255)
                                                            aCTATGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGAGAGGCGGCCCCCAGCGGCCTTGGCCCATGCAGCTGGTTGGCGACATCa                                                                                      >  2:367173/1‑150 (MQ=255)
                                                              tttGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAg                                                                                    <  2:3522153/150‑1 (MQ=255)
                                                                tGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAggg                                                                                  <  1:1271600/150‑1 (MQ=255)
                                                                     gTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATg                                                                             >  1:331564/1‑150 (MQ=255)
                                                                       cATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGcc                                                                           >  2:323931/1‑150 (MQ=255)
                                                                         tCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCaa                                                                         <  1:102514/150‑1 (MQ=255)
                                                                          cccGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAg                                                                        >  2:4395869/1‑150 (MQ=255)
                                                                           ccGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGa                                                                       >  2:1349741/1‑150 (MQ=255)
                                                                            cGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGAt                                                                      >  2:1779098/1‑150 (MQ=255)
                                                                             ggCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATc                                                                     <  2:3484177/150‑1 (MQ=255)
                                                                              gCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATcc                                                                    >  1:964866/1‑150 (MQ=255)
                                                                                aCTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCtt                                                                  >  1:4334496/1‑150 (MQ=255)
                                                                                 cTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTa                                                                 <  2:4621551/150‑1 (MQ=255)
                                                                                   gACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATg                                                               >  2:2655733/1‑150 (MQ=255)
                                                                                   gACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATg                                                               <  1:4083447/150‑1 (MQ=255)
                                                                                     cacGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATggg                                                             >  2:464320/1‑150 (MQ=255)
                                                                                       cGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCa                                                           >  2:888394/1‑150 (MQ=255)
                                                                                        ggACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCaa                                                          >  1:2324226/1‑150 (MQ=255)
                                                                                        ggACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCaa                                                          >  1:1056160/1‑150 (MQ=255)
                                                                                        ggACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCaa                                                          <  1:3833635/150‑1 (MQ=255)
                                                                                         gACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAg                                                         <  1:4378011/150‑1 (MQ=255)
                                                                                         gACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAg                                                         <  1:1800696/150‑1 (MQ=255)
                                                                                          aCCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTGATGGGCAAga                                                        >  1:298260/1‑150 (MQ=255)
                                                                                           ccTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCACgat                                                       >  1:1703634/1‑150 (MQ=255)
                                                                                           ccTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAgat                                                       <  1:200065/150‑1 (MQ=255)
                                                                                            cTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTGATGGGCAAgatc                                                      >  2:275140/1‑150 (MQ=255)
                                                                                            cTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCGTATGGGCACgatc                                                      >  2:1044994/1‑150 (MQ=255)
                                                                                                 ccAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCa                                                 >  2:3921008/1‑150 (MQ=255)
                                                                                                  cAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCAt                                                <  2:4498169/150‑1 (MQ=255)
                                                                                                     aaTGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGaa                                             <  1:1865647/150‑1 (MQ=255)
                                                                                                       tGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAgg                                           >  1:3617986/1‑150 (MQ=255)
                                                                                                           ggACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCggg                                       <  2:2742584/150‑1 (MQ=255)
                                                                                                               cAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCCCGATCGATCATGAAGGCGGGAGaa                                   >  2:3084374/1‑150 (MQ=255)
                                                                                                               cAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCACGATCGATCATGAAGGCGGGAGaa                                   >  2:2515310/1‑150 (MQ=255)
                                                                                                                  tCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGa                                <  2:1360824/150‑1 (MQ=255)
                                                                                                                     tGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGAtct                             <  2:3604771/150‑1 (MQ=255)
                                                                                                                     tGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGAtct                             <  2:3604442/150‑1 (MQ=255)
                                                                                                                     tGGCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATCCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGAtct                             <  2:2903912/150‑1 (MQ=255)
                                                                                                                      ggCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGAtctc                            <  2:2976056/150‑1 (MQ=255)
                                                                                                                      ggCGACCGCGACCCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGAtctc                            <  1:2435558/150‑1 (MQ=255)
                                                                                                                       gcgaccgcgaccCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTcc                           <  2:1551045/150‑1 (MQ=255)
                                                                                                                       gcgaccgcgaccCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTcc                           <  2:777107/150‑1 (MQ=255)
                                                                                                                        cgaccgcgaccCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCg                          <  1:878516/150‑1 (MQ=255)
                                                                                                                         gaccgcgaccCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCCAGATCGATCATGAAGGCGGGAGAAAGATCTCCGc                         >  2:3170984/1‑150 (MQ=255)
                                                                                                                         gaccgcgaccCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCACGATCGATCATGAAGGCGGGAGAAAGATCTCCGc                         >  2:1629905/1‑150 (MQ=255)
                                                                                                                         gaccgcgaccCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGc                         >  1:4174671/1‑150 (MQ=255)
                                                                                                                            cgcgaccCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTCTGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGa                      >  2:3952670/1‑150 (MQ=255)
                                                                                                                            cgcgaccCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGa                      <  2:4013473/150‑1 (MQ=255)
                                                                                                                               gaccCGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGa                   >  2:529334/1‑150 (MQ=255)
                                                                                                                               gaccCGCTGCGCGCGCCTGGCGATAAACCTCCAGTGTGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGTAAGATCTCCGTTGAGGa                   <  2:2171893/150‑1 (MQ=255)
                                                                                                                                  ccGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTGATGGGCACGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACtt                >  2:2837100/1‑150 (MQ=255)
                                                                                                                                   cGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACttt               <  1:2285392/150‑1 (MQ=255)
                                                                                                                                   cGCTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACttt               <  1:3423404/150‑1 (MQ=255)
                                                                                                                                     cTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCACGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACtttgt             >  1:4033056/1‑150 (MQ=255)
                                                                                                                                     cTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCACGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACtttgt             >  1:3980156/1‑150 (MQ=255)
                                                                                                                                     cTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACtttgt             >  1:4467297/1‑150 (MQ=255)
                                                                                                                                     cTGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACtttgt             >  2:2804557/1‑150 (MQ=255)
                                                                                                                                      tGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCACGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACtttgtt            >  2:2000848/1‑150 (MQ=255)
                                                                                                                                      tGCGCGCGCCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACtttgtt            >  1:3724243/1‑150 (MQ=255)
                                                                                                                                            cgcCTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCa      <  2:1840874/150‑1 (MQ=255)
                                                                                                                                              ccTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAgg    <  2:655209/150‑1 (MQ=255)
                                                                                                                                               cTGGCGATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCACGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGa   >  1:1454983/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
ATGTCGTTGATGGTGTACTTATTCATTCATTTATTCCTTTATCCAGTCGTT‑ACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAA  >  CP009273/4065986‑4066278

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