Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,375,260 G→T 88.7% D387Y (GAT→TAT)  degQ → serine endoprotease, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,375,2600GT88.7% 776.8 / 93.6 284D387Y (GAT→TAT) degQserine endoprotease, periplasmic
Reads supporting (aligned to +/- strand):  ref base G (20/12);  new base T (117/135);  total (137/147)
Fisher's exact test for biased strand distribution p-value = 9.42e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.77e-01

Only 100 of 290 total aligned reads displayed.

GGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCAAAAC  >  CP009273/3375113‑3375407
                                                                                                                                                   |                                                                                                                                                   
ggCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAGGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAatat                                                                                                                                                   <  1:2750951/150‑1 (MQ=255)
ggCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAatat                                                                                                                                                   <  1:525680/150‑1 (MQ=255)
ggCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAatat                                                                                                                                                   <  1:579532/150‑1 (MQ=255)
 gCACGAAAGTGAAGCTTGGCCTGCTGCGTACCGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATg                                                                                                                                                  <  2:1552870/150‑1 (MQ=255)
  cACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATgg                                                                                                                                                 >  1:1425395/1‑150 (MQ=255)
  cACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATgg                                                                                                                                                 <  2:384041/150‑1 (MQ=255)
  cACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATgg                                                                                                                                                 <  1:1753303/150‑1 (MQ=255)
   aCGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATggc                                                                                                                                                >  1:970581/1‑150 (MQ=255)
     gAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATggcgg                                                                                                                                              <  2:3266875/150‑1 (MQ=255)
          tGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAg                                                                                                                                         <  1:2591576/150‑1 (MQ=255)
           gAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAgg                                                                                                                                        <  2:1371652/150‑1 (MQ=255)
                 tGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTaa                                                                                                                                  >  1:4312471/1‑150 (MQ=255)
                   gCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTaaaa                                                                                                                                <  2:2738687/150‑1 (MQ=255)
                     ctgctgCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATc                                                                                                                              <  2:354284/150‑1 (MQ=255)
                     ctgctgCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATc                                                                                                                              <  2:1469661/150‑1 (MQ=255)
                      tgctgcGTAACGGCAAACCCCTTGAAGTAGAAGTGACGCTCGATACCAGCACCTCGACGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGTTCAGCTAAAATATTGCGGCAAAGGTATTAAAATCg                                                                                                                             <  1:885061/150‑1 (MQ=255)
                        ctgcGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGag                                                                                                                           <  2:1649989/150‑2 (MQ=255)
                        ctgcGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGAt                                                                                                                           <  1:805270/150‑1 (MQ=255)
                         tgcGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATg                                                                                                                          <  1:3515854/150‑1 (MQ=255)
                         tgcGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATg                                                                                                                          <  2:532275/150‑1 (MQ=255)
                          gcGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGa                                                                                                                         <  1:1043227/150‑1 (MQ=255)
                           cGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCAATGaa                                                                                                                        >  1:1570813/1‑150 (MQ=255)
                            gTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCCGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAg                                                                                                                       >  2:77367/1‑150 (MQ=255)
                              aaCGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGtt                                                                                                                     <  2:3282011/150‑1 (MQ=255)
                               aCGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTg                                                                                                                    <  1:404286/150‑1 (MQ=255)
                                cGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGt                                                                                                                   >  1:3465264/1‑150 (MQ=255)
                                  gCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCa                                                                                                                 >  2:2525063/1‑150 (MQ=255)
                                   cAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCaa                                                                                                                <  2:2130731/150‑1 (MQ=255)
                                    aaaCCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCaaa                                                                                                               >  1:3729900/1‑150 (MQ=255)
                                     aaCCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAg                                                                                                              >  1:1728866/1‑150 (MQ=255)
                                      aCCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAgg                                                                                                             >  1:2742641/1‑150 (MQ=255)
                                      aCCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAgg                                                                                                             <  1:3051281/150‑1 (MQ=255)
                                          cTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGc                                                                                                         >  1:3756227/1‑150 (MQ=255)
                                          cTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGc                                                                                                         >  2:1847050/1‑150 (MQ=255)
                                          cTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGc                                                                                                         <  2:1426018/150‑1 (MQ=255)
                                            ggAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGccc                                                                                                       >  2:583739/1‑150 (MQ=255)
                                              aaGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAg                                                                                                     >  2:1499901/1‑150 (MQ=255)
                                                gTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAgct                                                                                                   <  1:2257204/150‑1 (MQ=255)
                                                gTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAgct                                                                                                   >  1:2061115/1‑150 (MQ=255)
                                                 tAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAgctg                                                                                                  <  2:2483463/150‑1 (MQ=255)
                                                     aGTGACGCTCGAGACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCa                                                                                              <  2:3397303/150‑1 (MQ=255)
                                                      gTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAg                                                                                             >  2:1307763/1‑150 (MQ=255)
                                                      gTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAg                                                                                             <  1:64767/150‑1 (MQ=255)
                                                         aCGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAggct                                                                                          >  2:966240/1‑150 (MQ=255)
                                                           gCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAggctgg                                                                                        >  2:1891744/1‑150 (MQ=255)
                                                              agATACCATAACCTCTTCGTCGGCCAGCGCTGAAATTATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCACAGGAATTAAAATCGATGAAGTTGTCAAAGGAAGCACAGCTTCTCAGGCTGGCtt                                                                                     >  2:4296568/2‑150 (MQ=255)
                                                               gATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTg                                                                                    <  1:4133751/150‑1 (MQ=255)
                                                               gATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTg                                                                                    >  2:3185683/1‑150 (MQ=255)
                                                                   ccAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCaaa                                                                                >  2:872829/1‑150 (MQ=255)
                                                                    cAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCaaaa                                                                               >  2:1233823/1‑150 (MQ=255)
                                                                     aGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCaaaaa                                                                              <  1:272110/150‑1 (MQ=255)
                                                                      gCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAg                                                                             <  2:3587960/150‑1 (MQ=255)
                                                                        aCCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGAc                                                                           <  1:3173716/150‑1 (MQ=255)
                                                                        aCCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGAc                                                                           >  1:1209346/1‑150 (MQ=255)
                                                                        aCCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGAc                                                                           <  1:4287443/150‑1 (MQ=255)
                                                                           tctTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGAt                                                                        <  2:3081886/150‑1 (MQ=255)
                                                                              tcgtcgGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGAtgtg                                                                     >  1:4219444/1‑150 (MQ=255)
                                                                               cgtcgGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGa                                                                    <  2:2022069/150‑1 (MQ=255)
                                                                                 tcgGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATc                                                                  <  2:3255318/150‑1 (MQ=255)
                                                                                 tcgGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATc                                                                  >  2:3483211/1‑150 (MQ=255)
                                                                                   ggCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCAt                                                                >  2:1112942/1‑150 (MQ=255)
                                                                                    gCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCAtt                                                               <  2:1436993/150‑1 (MQ=255)
                                                                                    gCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCAtt                                                               >  1:1482638/1‑150 (MQ=255)
                                                                                     ccAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTg                                                              <  1:1290228/150‑1 (MQ=255)
                                                                                      cAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTgg                                                             >  1:3470847/1‑150 (MQ=255)
                                                                                         cgcTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGt                                                          <  2:1553824/150‑1 (MQ=255)
                                                                                         cgcTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGt                                                          <  2:3342439/150‑1 (MQ=255)
                                                                                          gcTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTc                                                         >  1:4181374/1‑150 (MQ=255)
                                                                                          gcTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTc                                                         >  1:3607389/1‑150 (MQ=255)
                                                                                            tGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCaa                                                       >  1:4008950/1‑150 (MQ=255)
                                                                                             gAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAAc                                                      <  2:3680931/150‑1 (MQ=255)
                                                                                                aTGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACcgc                                                   >  1:2476957/1‑150 (MQ=255)
                                                                                                 tGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACcgcg                                                  >  1:2515290/1‑150 (MQ=255)
                                                                                                   aTCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGAt                                                <  2:927784/150‑1 (MQ=255)
                                                                                                     cACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCg                                              <  2:2212375/150‑1 (MQ=255)
                                                                                                     cACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCg                                              <  2:3956442/150‑1 (MQ=255)
                                                                                                     cACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCg                                              <  1:990356/150‑1 (MQ=255)
                                                                                                     cACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCg                                              <  1:3226817/150‑1 (MQ=255)
                                                                                                        gCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAACTATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGt                                           <  1:220413/150‑1 (MQ=255)
                                                                                                          cAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGa                                         >  1:334881/1‑150 (MQ=255)
                                                                                                           aGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGaa                                        >  1:186798/1‑150 (MQ=255)
                                                                                                           aGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGaa                                        >  1:1093487/1‑150 (MQ=255)
                                                                                                               cTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCg                                    >  2:1586792/1‑150 (MQ=255)
                                                                                                               cTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCg                                    >  1:1921440/1‑150 (MQ=255)
                                                                                                                  gAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGAtt                                 >  1:2342668/1‑150 (MQ=255)
                                                                                                                   aaGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTg                                >  1:22325/1‑150 (MQ=255)
                                                                                                                     ggTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCt                              <  2:1916019/150‑1 (MQ=255)
                                                                                                                     ggTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCt                              <  1:748977/150‑1 (MQ=255)
                                                                                                                      gTGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTg                             >  2:4256593/1‑150 (MQ=255)
                                                                                                                       tGCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGa                            <  1:4350312/150‑1 (MQ=255)
                                                                                                                        gCAACGTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGaa                           >  2:2994559/1‑150 (MQ=255)
                                                                                                                             gTTGAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCg                      >  2:683678/1‑150 (MQ=255)
                                                                                                                                gAGCGATGGTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTaa                   <  2:803242/150‑1 (MQ=255)
                                                                                                                                       ggTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTg            >  2:4131409/1‑150 (MQ=255)
                                                                                                                                       ggTCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTg            <  2:1275357/150‑1 (MQ=255)
                                                                                                                                         tCAGCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTggc          >  1:3309278/1‑150 (MQ=255)
                                                                                                                                            gCTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTggcggc       >  2:116216/1‑150 (MQ=255)
                                                                                                                                             cTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCa      >  1:3091133/1‑150 (MQ=255)
                                                                                                                                             cTAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCa      <  2:2811374/150‑1 (MQ=255)
                                                                                                                                              tAAAATATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCaa     >  2:3749051/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GGCACGAAAGTGAAGCTTGGCCTGCTGCGTAACGGCAAACCACTGGAAGTAGAAGTGACGCTCGATACCAGCACCTCTTCGTCGGCCAGCGCTGAAATGATCACGCCAGCGCTGGAAGGTGCAACGTTGAGCGATGGTCAGCTAAAAGATGGCGGCAAAGGTATTAAAATCGATGAAGTTGTCAAAGGAAGCCCAGCTGCTCAGGCTGGCTTGCAAAAAGACGATGTGATCATTGGCGTCAACCGCGATCGGGTGAACTCGATTGCTGAAATGCGTAAAGTGCTGGCGGCAAAAC  >  CP009273/3375113‑3375407

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: