Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,292,914 C→A 6.2% D121Y (GAT→TAT)  yhbP ← UPF0306 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,292,9140CA 6.2% 1332.7 / 39.3 340D121Y (GAT→TAT) yhbPUPF0306 family protein
Reads supporting (aligned to +/- strand):  ref base C (158/161);  new base A (11/10);  total (169/171)
Fisher's exact test for biased strand distribution p-value = 8.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 346 total aligned reads displayed.

GTTGATAAAGATCAGCTGGTTACCAGCCGGACACCGGACGATCTGCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTGCACACCGCGAATTAACGCTACCGTTTTCG  >  CP009273/3292765‑3293063
                                                                                                                                                     |                                                                                                                                                     
gTTGATAAAGATCAGCTGGTTACCAGCCGGACACCGGACGATCTGCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATc                                                                                                                                                       <  2:2350207/150‑1 (MQ=255)
gTTGATAAAGATCAGCTGGTTACCAGCCGGACACCGGACGATCTGCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATc                                                                                                                                                       <  2:2810691/150‑1 (MQ=255)
   gATAAAGATCAGCTGGTTACCAGCCGGACACCGGACGATCTGCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAg                                                                                                                                                    >  1:886841/1‑150 (MQ=255)
          aTCAGCTGGTTACCAGCCGGACACCGGACGATCTGCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGAtt                                                                                                                                             <  2:786605/150‑1 (MQ=255)
           tCAGCTGGTTACCAGCCGGACACCGGACGATCTGCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGAttt                                                                                                                                            >  1:5114788/1‑150 (MQ=255)
             aGCTGGTTACCAGCCGGACACCGGACGATCTGCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTcc                                                                                                                                          <  2:409656/150‑1 (MQ=255)
              gCTGGTTACCAGCCGGACACCGGACGATCTGCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTccc                                                                                                                                         <  2:4012111/150‑1 (MQ=255)
               cTGGTTACCAGCCGGACACCGGACGATCTGCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCca                                                                                                                                        <  1:1664059/150‑1 (MQ=255)
                 ggTTACCAGCCGGACACCGGACGATCTGCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCcaca                                                                                                                                      <  2:4838854/150‑1 (MQ=255)
                  gTTACCAGCCGGACACCGGACGATCTGCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCacac                                                                                                                                     >  1:118984/1‑150 (MQ=255)
                  gTTACCAGCCGGACACCGGACGATCTGCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCacac                                                                                                                                     >  1:1530285/1‑150 (MQ=255)
                   ttACCAGCCGGACACCGGACGATCTGCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACAcc                                                                                                                                    >  1:2976085/1‑150 (MQ=255)
                    tACCAGCCGGACACCGGACGATCTGCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCg                                                                                                                                   <  1:3876969/150‑1 (MQ=255)
                       cAGCCGGACACCGGACGATCTGCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCgggg                                                                                                                                >  2:308102/1‑150 (MQ=255)
                        aGCCGGACACCGGACGATCTGCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGc                                                                                                                               >  2:919911/1‑150 (MQ=255)
                             gACACCGGACGATCTGCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACa                                                                                                                          <  1:4760622/150‑1 (MQ=255)
                             gACACCGGACGATCTGCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACa                                                                                                                          <  2:215010/150‑1 (MQ=255)
                             gACACCGGACGATCTGCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACa                                                                                                                          >  2:2635995/1‑150 (MQ=255)
                                   ggACGATCTGCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTc                                                                                                                    >  1:614097/1‑150 (MQ=255)
                                    gACGATCTGCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCt                                                                                                                   >  1:2934646/1‑150 (MQ=255)
                                       gATCTGCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCCTCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGc                                                                                                                >  2:1048187/1‑150 (MQ=255)
                                       gATCTGCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGc                                                                                                                >  2:690551/1‑150 (MQ=255)
                                          cTGCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAAc                                                                                                             <  1:2157546/150‑1 (MQ=255)
                                          cTGCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAAc                                                                                                             <  1:761514/150‑1 (MQ=255)
                                            gCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCg                                                                                                           <  1:2293405/150‑1 (MQ=255)
                                            gCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCg                                                                                                           <  1:3250042/150‑1 (MQ=255)
                                             ccAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCgg                                                                                                          <  1:3601325/150‑1 (MQ=255)
                                               aGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGaa                                                                                                        >  1:4093362/1‑150 (MQ=255)
                                               aGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGaa                                                                                                        <  2:1158710/150‑1 (MQ=255)
                                                 cGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAg                                                                                                      <  1:2403606/150‑1 (MQ=255)
                                                  gTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGc                                                                                                     >  2:3566916/1‑150 (MQ=255)
                                                    ttAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATAGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGcga                                                                                                   >  1:5065616/1‑150 (MQ=255)
                                                     tAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGcgac                                                                                                  <  1:3628569/150‑1 (MQ=255)
                                                     tAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGcgac                                                                                                  <  1:1345192/150‑1 (MQ=255)
                                                         cgcgAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGAtt                                                                                              <  1:5032266/150‑1 (MQ=255)
                                                         cgcgAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGAtt                                                                                              <  2:5010196/150‑1 (MQ=255)
                                                          gcgAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTg                                                                                             <  1:172257/150‑1 (MQ=255)
                                                            taGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTa                                                                                           >  2:3535587/2‑150 (MQ=255)
                                                            gAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTa                                                                                           <  1:3944522/150‑1 (MQ=255)
                                                              ggCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACg                                                                                         <  1:2036294/150‑1 (MQ=255)
                                                                 gTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGGCGATTGTACGCCt                                                                                      >  2:2207172/1‑150 (MQ=255)
                                                                  ttACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCtt                                                                                     <  1:1004226/150‑1 (MQ=255)
                                                                   tACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTg                                                                                    >  2:3700471/1‑150 (MQ=255)
                                                                     cGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTgcg                                                                                  <  1:279083/150‑1 (MQ=255)
                                                                     cGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATAGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTgcg                                                                                  <  2:488407/150‑1 (MQ=255)
                                                                          gCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGa                                                                             <  2:2105639/150‑1 (MQ=255)
                                                                           cTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAg                                                                            >  2:4684098/1‑150 (MQ=255)
                                                                             cGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGt                                                                          <  2:5117971/150‑1 (MQ=255)
                                                                               gTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCg                                                                        <  1:723138/150‑1 (MQ=255)
                                                                                tGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGc                                                                       >  2:2569567/1‑150 (MQ=255)
                                                                                tGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGc                                                                       >  1:2555380/1‑150 (MQ=255)
                                                                                 gCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCt                                                                      <  1:3514624/150‑1 (MQ=255)
                                                                                 gCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCt                                                                      >  2:555905/1‑150 (MQ=255)
                                                                                 gCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCt                                                                      <  1:2205085/150‑1 (MQ=255)
                                                                                 gCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCt                                                                      >  2:1256530/1‑150 (MQ=255)
                                                                                  ccTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCtt                                                                     <  2:2601215/150‑1 (MQ=255)
                                                                                  ccTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCtt                                                                     <  1:2701350/150‑1 (MQ=255)
                                                                                    tGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTc                                                                   <  1:220868/150‑1 (MQ=255)
                                                                                     gAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTcc                                                                  >  1:3410152/1‑150 (MQ=255)
                                                                                       gTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTc                                                                <  2:2529244/150‑1 (MQ=255)
                                                                                        tCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCa                                                               >  1:3766439/1‑150 (MQ=255)
                                                                                         cgcgTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATAGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCAc                                                              <  2:4335088/150‑1 (MQ=255)
                                                                                              aGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTcc                                                         <  2:3465841/150‑1 (MQ=255)
                                                                                               gCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCa                                                        <  2:1596627/150‑1 (MQ=255)
                                                                                                ccAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGACAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAg                                                       >  2:4209684/1‑150 (MQ=255)
                                                                                                ccAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAg                                                       >  1:1168623/1‑150 (MQ=255)
                                                                                                 cAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGc                                                      >  1:4334215/1‑150 (MQ=255)
                                                                                                 aaaaTCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACAACGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAga                                                      >  2:2922855/2‑149 (MQ=255)
                                                                                                  aaaTCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGcc                                                     <  1:3170080/150‑1 (MQ=255)
                                                                                                   aaTCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCt                                                    <  2:2678792/150‑1 (MQ=255)
                                                                                                    aTCATTTTTTTACCAAAGCCCAGCGGGTTCTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCCGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTg                                                   >  2:1242455/1‑150 (MQ=255)
                                                                                                       aTTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATAGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCgg                                                >  2:1781663/1‑150 (MQ=255)
                                                                                                             ttACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCa                                          <  1:2114595/150‑1 (MQ=255)
                                                                                                             ttACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCa                                          <  1:4669273/150‑1 (MQ=255)
                                                                                                              tACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCAc                                         >  1:3123631/1‑150 (MQ=255)
                                                                                                               aCCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCAcc                                        <  2:3176994/150‑1 (MQ=255)
                                                                                                                 cAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCtt                                      >  2:1491637/1‑150 (MQ=255)
                                                                                                                  aaaGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCttt                                     <  2:3284142/150‑1 (MQ=255)
                                                                                                                  aaaGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCttt                                     <  1:1708213/150‑1 (MQ=255)
                                                                                                                   aaGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTa                                    >  1:4488614/1‑150 (MQ=255)
                                                                                                                   aaGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTa                                    >  2:200023/1‑150 (MQ=255)
                                                                                                                   aaGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTa                                    >  1:1935464/1‑150 (MQ=255)
                                                                                                                    aGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTaa                                   >  2:1697669/1‑150 (MQ=255)
                                                                                                                     gCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTaaa                                  <  1:1147159/150‑1 (MQ=255)
                                                                                                                      cccAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAAc                                 <  2:2321491/150‑1 (MQ=255)
                                                                                                                      cccAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAAc                                 >  1:4079406/1‑150 (MQ=255)
                                                                                                                         aGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTGc                              <  1:4176087/150‑1 (MQ=255)
                                                                                                                           cGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTGcac                            >  2:4053391/1‑150 (MQ=255)
                                                                                                                            gtgtTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTGcaca                           >  1:33129/1‑150 (MQ=255)
                                                                                                                             tgttgtCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTGCacac                          <  2:3779284/150‑1 (MQ=255)
                                                                                                                              gttgtCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTGCACAcc                         >  2:1400152/1‑150 (MQ=255)
                                                                                                                               ttgtCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTGCACACcg                        >  1:1600387/1‑150 (MQ=255)
                                                                                                                               ttgtCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTGCACACcg                        >  1:4892164/1‑150 (MQ=255)
                                                                                                                                    ggTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTGCACACCGCGAAt                   <  1:2489389/150‑1 (MQ=255)
                                                                                                                                       gAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTGCACACCGCGAATTaa                >  1:2606929/1‑150 (MQ=255)
                                                                                                                                       gAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTGCACACCGCGAATTaa                >  1:118493/1‑150 (MQ=255)
                                                                                                                                         aTTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTGCACACCGCGAATTAACg              <  2:1505798/150‑1 (MQ=255)
                                                                                                                                         aTTTGATTTCATAGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCTCCTTCCAGCCTGCGGATCTCACCTTTAAACTGCACACCGCGAATTAACg              >  2:221363/1‑150 (MQ=255)
                                                                                                                                          tttGATTTCATCGAGCCGGATTTCCCACACCGGGTCCGACAGCATTCTGGCAACCGGATAGCGACGATTGTACGCCTTGCGCGCGATGTCGCTTTCCTCACGTTCCAGCCTGCGGATCTCACCTTTATACTGCACACCGCGAATTAACGc             <  1:342312/150‑1 (MQ=255)
                                                                                                                                                   aTCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTGCACACCGCGAATTAACGCTACCGtttt    <  2:4578962/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTTGATAAAGATCAGCTGGTTACCAGCCGGACACCGGACGATCTGCCAGCGTTTAACCGCGAGGCGTTACGCCTGCTCGGTGCCTGAGTCGCGTAGCCAAATCATTTTTTTACCAAAGCCCAGCGTGTTGTCGGTGAATTTGATTTCATCGAGCCGGATTTCCCACACCGGGGCCGACAGCATTCTGGCAACCGGAAAGCGACGATTGTACGCCTTGCGCGCGAGGTCGCTTTCCTCACCTTCCAGCCTGCGGATCTCACCTTTAAACTGCACACCGCGAATTAACGCTACCGTTTTCG  >  CP009273/3292765‑3293063

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑