Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,978,810 G→T 16.7% I186I (ATC→ATA aslA ← putative Ser‑type periplasmic non‑aryl sulfatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,978,8100GT16.7% 1029.3 / 125.1 321I186I (ATC→ATAaslAputative Ser‑type periplasmic non‑aryl sulfatase
Reads supporting (aligned to +/- strand):  ref base G (121/146);  new base T (32/22);  total (153/168)
Fisher's exact test for biased strand distribution p-value = 7.31e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 325 total aligned reads displayed.

GAACGGTCCGGACTCAGGGCCACTTCCGGATTGACGTGAACGTCGCGCCATTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGGTGAGAATGGTGGCGCGGGTTGGGGA  >  CP009273/3978661‑3978959
                                                                                                                                                     |                                                                                                                                                     
 aaCGGTCCGGACTCAGGGCCACTTCCGGATTGACGTGAACGTCGCGCCATTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGa                                                                                                                                                      >  2:352363/1‑150 (MQ=255)
  aCGGTCCGGACTCAGGGCCACTTCCGGATTGACGTGAACGTCGCGCCATTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGAt                                                                                                                                                     <  1:3972280/150‑1 (MQ=255)
    ggTCCGGACTCAGGGCCACTTCCGGATTGACGTGAACGTCGCGCCATTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGTTgg                                                                                                                                                   <  1:1425240/150‑1 (MQ=255)
    ggTCCGGACTCAGGGCCACTTCCGGATTGACGTGAACGTCGCGCCATTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATgg                                                                                                                                                   >  1:2917565/1‑150 (MQ=255)
    ggTCCGGACTCAGGGCCACTTCCGGATTGACGAGAACGTCGCGCCATTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTATCCGATgg                                                                                                                                                   >  1:2271454/1‑150 (MQ=255)
     gTCCGGACTCAGGGCCACTTCCGGATTGACGTGAACGTCGCGCCATTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGc                                                                                                                                                  <  1:1837886/150‑1 (MQ=255)
      tCCGGACTCAGGGCCACTTCCTGATTGACGTGAACGTCGCGCGATTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCACAGCCAACGTTATGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGcc                                                                                                                                                 <  1:4391734/150‑1 (MQ=255)
      tCCGGACTCAGGGCCACTTCCGGATTGACGTGTACGTCGCGCCTTTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGcc                                                                                                                                                 <  1:4305836/150‑1 (MQ=255)
       ccGGACTCAGGGCCACTTCCGGATTGACGTGAACGTCGCGCCATTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCt                                                                                                                                                <  1:3572210/150‑1 (MQ=255)
        cGGACTCAGGGCCACTTCCGGATTGACGTGAACGTCGCGCCATTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTg                                                                                                                                               <  1:1726809/150‑1 (MQ=255)
         ggACTCAGGGCCACTTCCGGATTGACGTGAACGTCGCGCCATTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGa                                                                                                                                              >  1:901142/1‑150 (MQ=255)
           aCTCAGGGCCACTTCCGGATTGACGTGAACGTCGCGCCATTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCTATGGCCTGAGt                                                                                                                                            >  2:1824982/1‑150 (MQ=255)
             tCAGGGCCACTTCCGGATTGACGTGAACGTCGCGCCATTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGa                                                                                                                                          <  2:2233295/150‑1 (MQ=255)
              cAGGGCCACTTCCGGATTGACGTGAACGTCGCGCCATTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGAc                                                                                                                                         >  1:353008/1‑150 (MQ=255)
                gggCCACTTCCGGATTGACGTGAACGTCGCGCCATTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTATGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCTATGGCCTGAGTGACGt                                                                                                                                       <  2:3686140/150‑1 (MQ=255)
                 ggCCACTTCCGGATTGACGTGAACGTCGCGCCATTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGATATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTa                                                                                                                                      <  1:2517472/150‑1 (MQ=255)
                      cTTCCGGATTGACGTGAACGTCGCGCCATTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCt                                                                                                                                 <  2:2467668/150‑1 (MQ=255)
                      cTTCCGGATTGACGTGAACGTCGCGCCATTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCt                                                                                                                                 <  1:295224/150‑1 (MQ=255)
                       ttCCGGATTGACGTGAACGTCGCGCCATTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTg                                                                                                                                <  2:2616219/150‑1 (MQ=255)
                         ccGGATTGACGTGAACGTCGCGCCATTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCTATGGCCTGAGTGACGTAGCCCTGAt                                                                                                                              >  1:3800220/1‑150 (MQ=255)
                         ccGGATTGACGTGAACGTCGCGCCATTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGAt                                                                                                                              <  2:1104065/150‑1 (MQ=255)
                          cGGATTGACGTGAACGTCGCGCCATTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATc                                                                                                                             <  2:2068304/150‑1 (MQ=255)
                           ggATTGACGTGAACGTCGCGCCATTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCTATGGCCTGAGTGACGTAGCCCTGATCg                                                                                                                            >  1:3603011/1‑150 (MQ=255)
                                    tGAACGTCGCGCCATTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAAc                                                                                                                   <  1:4621611/150‑1 (MQ=255)
                                    tGAACGTCGCGCCATTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAAc                                                                                                                   >  1:2990119/1‑150 (MQ=255)
                                         gTCGCGCCATTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCgg                                                                                                              >  1:1550257/1‑150 (MQ=255)
                                             cgcCATTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGc                                                                                                          <  2:358694/150‑1 (MQ=255)
                                             cgcCATTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGc                                                                                                          <  2:1957577/150‑1 (MQ=255)
                                                 aTTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTgg                                                                                                      >  2:2407585/1‑150 (MQ=255)
                                                   tCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGtt                                                                                                    <  2:606503/150‑1 (MQ=255)
                                                   tCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGtt                                                                                                    <  2:3605693/150‑1 (MQ=255)
                                                     ggTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTaa                                                                                                  >  1:329710/1‑150 (MQ=255)
                                                      gtgtACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCTATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAAc                                                                                                 <  1:575225/150‑1 (MQ=255)
                                                         tACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCt                                                                                              >  1:2471839/1‑150 (MQ=255)
                                                           cATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTg                                                                                            <  2:1731373/150‑1 (MQ=255)
                                                            atatCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGc                                                                                           <  1:4541589/150‑1 (MQ=255)
                                                               tCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAgc                                                                                        >  1:1993511/1‑150 (MQ=255)
                                                                cAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAgcc                                                                                       <  1:4246695/150‑1 (MQ=255)
                                                                 aGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAgccc                                                                                      <  2:3463272/150‑1 (MQ=255)
                                                                     acCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAgcccgccc                                                                                  <  2:1324641/150‑1 (MQ=255)
                                                                      ccGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCNGCAGCGTGGTTAACCCTTGCAGcccgcccg                                                                                 <  1:4281017/150‑1 (MQ=255)
                                                                      ccGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTCGCAGcccgcccg                                                                                 >  2:1197370/1‑150 (MQ=255)
                                                                           ttAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGc                                                                            <  1:1239841/150‑1 (MQ=255)
                                                                                gCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCTATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTa                                                                       <  1:2649643/150‑1 (MQ=255)
                                                                                  cACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACa                                                                     >  2:3292555/1‑150 (MQ=255)
                                                                                   aCGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACAt                                                                    <  2:3422365/150‑1 (MQ=255)
                                                                                    cGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACAtt                                                                   <  1:2569002/150‑1 (MQ=255)
                                                                                     ggAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTg                                                                  >  1:3148853/1‑150 (MQ=255)
                                                                                     ggAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTg                                                                  <  1:1580076/150‑1 (MQ=255)
                                                                                      gAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCTATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTgg                                                                 <  1:3886507/150‑1 (MQ=255)
                                                                                      gAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTgg                                                                 >  1:4524444/1‑150 (MQ=255)
                                                                                      gAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTgg                                                                 >  1:4014204/1‑150 (MQ=255)
                                                                                        aaTCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCTATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTggcg                                                               >  2:2181489/1‑150 (MQ=255)
                                                                                          tcatcaAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTggcggc                                                             <  2:2302615/150‑1 (MQ=255)
                                                                                            atcaAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCAt                                                           >  1:4531369/1‑150 (MQ=255)
                                                                                             tcaAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATc                                                          <  2:4463729/150‑1 (MQ=255)
                                                                                                 aGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGaa                                                      >  1:464591/1‑150 (MQ=255)
                                                                                                 aGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGaa                                                      <  2:2167797/150‑1 (MQ=255)
                                                                                                 aGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGaa                                                      <  2:3885000/150‑1 (MQ=255)
                                                                                                  gCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAAt                                                     <  2:472438/150‑1 (MQ=255)
                                                                                                    cAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGc                                                   <  1:1617882/150‑1 (MQ=255)
                                                                                                     aaCGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGcc                                                  >  2:333858/1‑150 (MQ=255)
                                                                                                     aaCGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGcc                                                  >  2:452138/1‑150 (MQ=255)
                                                                                                     aaCGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGcc                                                  >  1:3241762/1‑150 (MQ=255)
                                                                                                      aCGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCTATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGGACATTGGCGGCATCAGAATGCCg                                                 >  2:4581717/1‑150 (MQ=255)
                                                                                                        gTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCgtg                                               >  1:945789/1‑150 (MQ=255)
                                                                                                         ttCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGGTGCCCGTACATTGGCGGCATCAGAATGCCgtgg                                              <  2:1672053/150‑1 (MQ=255)
                                                                                                          tCTGCGGATGCGACTCTGTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGTTTGCCCGTACATTGGCGGAATCAGAATGCCgtggt                                             <  1:1236761/150‑1 (MQ=255)
                                                                                                           cTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCTATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCgtggtg                                            <  2:826501/150‑1 (MQ=255)
                                                                                                           cTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCgtggtg                                            <  2:3016603/150‑1 (MQ=255)
                                                                                                           cTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCgtggtg                                            <  2:3058733/150‑1 (MQ=255)
                                                                                                            tGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGtggtgg                                           >  2:3216384/1‑150 (MQ=255)
                                                                                                               ggCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGtggatg                                        >  1:3371324/1‑150 (MQ=255)
                                                                                                                  tGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCTATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAg                                     <  1:87172/150‑1 (MQ=255)
                                                                                                                  tGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAg                                     >  2:2276431/1‑150 (MQ=255)
                                                                                                                  tGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAg                                     >  2:1246989/1‑150 (MQ=255)
                                                                                                                   gCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGt                                    <  1:3925136/150‑1 (MQ=255)
                                                                                                                    cGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTa                                   >  2:161221/1‑150 (MQ=255)
                                                                                                                    cGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTa                                   >  1:858261/1‑150 (MQ=255)
                                                                                                                     gACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTAt                                  >  1:1288039/1‑150 (MQ=255)
                                                                                                                      aCTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTAtt                                 >  2:746874/1‑150 (MQ=255)
                                                                                                                      aCTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTAtt                                 >  2:2919623/1‑150 (MQ=255)
                                                                                                                       ctctTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTg                                >  2:1846559/1‑150 (MQ=255)
                                                                                                                         ctTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTc                              <  2:2841278/150‑1 (MQ=255)
                                                                                                                          tttGTTTTCCCCCATATGCCATTTTCCTATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTcc                             >  2:2445983/1‑150 (MQ=255)
                                                                                                                           ttGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCg                            >  2:3816924/1‑150 (MQ=255)
                                                                                                                           ttGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCg                            >  2:1490996/1‑150 (MQ=255)
                                                                                                                               gttCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGGTga                        <  2:510523/149‑1 (MQ=255)
                                                                                                                                 tCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTCTTGTCCGGTgaga                      >  1:1026425/1‑150 (MQ=255)
                                                                                                                                  cccccATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGTAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGGTGAGaa                     >  2:437315/1‑150 (MQ=255)
                                                                                                                                     ccATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGGTGAGAAtgg                  <  2:3927523/150‑1 (MQ=255)
                                                                                                                                      cATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGGTGAGAAtggt                 <  1:279918/150‑1 (MQ=255)
                                                                                                                                         atGCCATTTTCCTATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATTCCGTGGTGGATGGAGTATTGTCCGGTGAGAATGGTGgc              >  1:1764048/1‑150 (MQ=255)
                                                                                                                                         atGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGGTGAGAATGGTGgc              >  1:2773706/1‑150 (MQ=255)
                                                                                                                                              aTTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGGTGAGAATGGTGGCGCggg         >  1:3698674/1‑150 (MQ=255)
                                                                                                                                              aTTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGGTGAGAATGGTGGCGCggg         >  2:1408950/1‑150 (MQ=255)
                                                                                                                                               ttttCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGGTGAGAATGGTGGCGCGGGt        <  1:697157/150‑1 (MQ=255)
                                                                                                                                                  tCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGGTGAGAATGGTGGCGCGGGTTgg     >  1:3215555/1‑150 (MQ=255)
                                                                                                                                                   ccGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGGTGAGAATGGTGGCGCGGGTTggg    <  1:1122239/150‑1 (MQ=255)
                                                                                                                                                     gATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGGTGAGAATGGTGGCGCGGGTTGGGGa  >  2:4025939/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAACGGTCCGGACTCAGGGCCACTTCCGGATTGACGTGAACGTCGCGCCATTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGGTGAGAATGGTGGCGCGGGTTGGGGA  >  CP009273/3978661‑3978959

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: