Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 165,759 G→T 84.0% E597* (GAG→TAG)  fhuA → ferrichrome outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273165,7590GT84.0% 599.2 / 120.5 279E597* (GAG→TAG) fhuAferrichrome outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base G (25/20);  new base T (118/116);  total (143/136)
Fisher's exact test for biased strand distribution p-value = 6.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 284 total aligned reads displayed.

CCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACTACCTATAAAGGCAATACGCCTGCACAGGTGCCAAAACAC  >  CP009273/165612‑165905
                                                                                                                                                   |                                                                                                                                                  
ccGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCtag                                                                                                                                                  >  1:4840394/1‑150 (MQ=255)
 cGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCtaga                                                                                                                                                 <  2:3912383/150‑1 (MQ=255)
   tCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATc                                                                                                                                               <  1:4123383/150‑1 (MQ=255)
    cTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATcc                                                                                                                                              <  2:2294341/150‑1 (MQ=255)
      aaaGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCcgc                                                                                                                                            <  2:3828017/150‑1 (MQ=255)
         ggTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCa                                                                                                                                         <  2:2773281/150‑1 (MQ=255)
             aGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTg                                                                                                                                     <  1:4338217/150‑1 (MQ=255)
             aGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTg                                                                                                                                     <  2:1709324/150‑1 (MQ=255)
             aGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTg                                                                                                                                     <  2:4840394/150‑1 (MQ=255)
                 gTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGATGAAGGTGGCTAGATCCGCGCACGTGGCGt                                                                                                                                 >  2:2677646/1‑150 (MQ=255)
                  tATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTa                                                                                                                                <  2:324374/150‑1 (MQ=255)
                   aTGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAg                                                                                                                               <  2:4669195/150‑1 (MQ=255)
                   aTGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAg                                                                                                                               >  2:5090827/1‑150 (MQ=255)
                   aTGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAg                                                                                                                               <  1:738233/150‑1 (MQ=255)
                      aaGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGaaa                                                                                                                            >  1:2706799/1‑150 (MQ=255)
                      aaGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCGTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGaaa                                                                                                                            <  2:2912618/150‑1 (MQ=255)
                       aGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAAt                                                                                                                           <  1:119782/150‑1 (MQ=255)
                          cGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGa                                                                                                                        >  2:136554/1‑150 (MQ=255)
                          cGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGa                                                                                                                        >  1:3105983/1‑150 (MQ=255)
                             cGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGc                                                                                                                     <  1:1580213/150‑1 (MQ=255)
                             cGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGc                                                                                                                     <  2:60645/150‑1 (MQ=255)
                              gTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCg                                                                                                                    >  2:4040270/1‑150 (MQ=255)
                               tGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGa                                                                                                                   <  2:804971/150‑1 (MQ=255)
                                  aaTATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAg                                                                                                                >  2:4601889/1‑150 (MQ=255)
                                    tatGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAgcg                                                                                                              <  1:3725471/150‑1 (MQ=255)
                                     atGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAgcgg                                                                                                             <  1:2091711/150‑1 (MQ=255)
                                      tttACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAgcggc                                                                                                            >  1:1481042/3‑150 (MQ=255)
                                      tGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAgcggc                                                                                                            <  2:5153924/150‑1 (MQ=255)
                                         aCCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCt                                                                                                         >  1:5045289/1‑150 (MQ=255)
                                          ccGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTg                                                                                                        <  2:379561/150‑1 (MQ=255)
                                          ccGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTg                                                                                                        <  2:1647758/150‑1 (MQ=255)
                                           cGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGt                                                                                                       >  1:1092799/1‑150 (MQ=255)
                                           cGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGt                                                                                                       >  1:928494/1‑150 (MQ=255)
                                               aGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGc                                                                                                   <  2:1491803/150‑1 (MQ=255)
                                                 aTCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGa                                                                                                 <  2:4507645/150‑1 (MQ=255)
                                                    gTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAgtg                                                                                              >  2:923886/1‑150 (MQ=255)
                                                    gTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAgtg                                                                                              >  1:421940/1‑150 (MQ=255)
                                                     tCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGGGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAgtgt                                                                                             <  2:3459688/150‑1 (MQ=255)
                                                      ccGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGtt                                                                                            <  1:4768848/150‑1 (MQ=255)
                                                        gATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTaa                                                                                          >  1:1259080/1‑150 (MQ=255)
                                                        gATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTaa                                                                                          >  2:2722290/1‑150 (MQ=255)
                                                        gATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGGAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTaa                                                                                          >  2:3062289/1‑150 (MQ=255)
                                                             tAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAg                                                                                     <  2:312531/150‑1 (MQ=255)
                                                              aGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGt                                                                                    >  2:3677946/1‑150 (MQ=255)
                                                                 tACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCgg                                                                                 <  1:3673802/150‑1 (MQ=255)
                                                                  aCTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGt                                                                                <  2:686185/150‑1 (MQ=255)
                                                                  aCTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGt                                                                                >  1:3427645/1‑150 (MQ=255)
                                                                    tGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTc                                                                              >  1:562658/1‑150 (MQ=255)
                                                                     ggTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCt                                                                             >  1:2407704/1‑150 (MQ=255)
                                                                     ggTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCt                                                                             >  2:3509629/1‑150 (MQ=255)
                                                                      gTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCtt                                                                            <  1:100348/150‑1 (MQ=255)
                                                                       tGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTta                                                                           >  2:639356/1‑150 (MQ=255)
                                                                         ccGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTtata                                                                         <  1:3953957/150‑1 (MQ=255)
                                                                          cGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATAc                                                                        <  2:2557077/150‑1 (MQ=255)
                                                                          cGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATAc                                                                        >  1:3129431/1‑150 (MQ=255)
                                                                             gtATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTa                                                                     <  1:2859416/150‑1 (MQ=255)
                                                                                taATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTacac                                                                  >  2:3406319/1‑150 (MQ=255)
                                                                                  aTCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCg                                                                >  1:370082/1‑150 (MQ=255)
                                                                                   tctcACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTATTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGa                                                               <  1:171187/150‑1 (MQ=255)
                                                                                      cccTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGc                                                            <  1:317423/148‑1 (MQ=255)
                                                                                       aCTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCg                                                           >  1:383651/1‑150 (MQ=255)
                                                                                       aCTAAAACAAACAACCTGTTGGCGGACCCTGTGGGTTCCTTCTTCTCTGTTGAAGGTGGCTAGATCCGCGCAGGTGGCGTAGAAATCTAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCg                                                           <  1:958213/150‑1 (MQ=255)
                                                                                        cTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCgg                                                          >  1:3634100/1‑150 (MQ=255)
                                                                                        cTAAAACCAACAACCTGAGGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCgg                                                          <  1:98080/150‑1 (MQ=255)
                                                                                         tAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGa                                                         >  2:3955908/1‑150 (MQ=255)
                                                                                          aaaaCCAACACCCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGaa                                                        <  1:5062443/150‑1 (MQ=255)
                                                                                          aaaaCCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGaa                                                        >  2:880247/1‑150 (MQ=255)
                                                                                           aaaCCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAAt                                                       >  1:1615038/1‑150 (MQ=255)
                                                                                           aaaCCAACAACCTGATGGCGGACCCGGCGGGTTCCTTCTTCGCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAACTCGAAGCGAACGCGGCGCGGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGact                                                       >  2:322211/1‑148 (MQ=255)
                                                                                                aacaacCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATAcacc                                                  >  2:2366126/1‑150 (MQ=255)
                                                                                                 acaacCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATAcacca                                                 <  1:2819751/150‑1 (MQ=255)
                                                                                                 acaacCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATAcacca                                                 >  2:1987055/1‑150 (MQ=255)
                                                                                                  caacCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATAcaccac                                                >  2:2800546/1‑150 (MQ=255)
                                                                                                   aacCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACaccacc                                               >  2:555700/1‑150 (MQ=255)
                                                                                                   aacCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACaccacc                                               >  1:776410/1‑150 (MQ=255)
                                                                                                   aacCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACaccacc                                               >  2:5000180/1‑150 (MQ=255)
                                                                                                      cTGATGGCTGACCCTGAGGGTTCATTCTTCTCGGTTGAATGTGGCTAGATACGCGAACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATAATTAAACCGATGCGGAATACACCACCGAt                                            >  2:1929103/1‑150 (MQ=255)
                                                                                                      cTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGAt                                            >  1:4572724/1‑150 (MQ=255)
                                                                                                         aTGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGAtact                                         <  1:1062165/150‑1 (MQ=255)
                                                                                                         aTGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGAtact                                         <  1:3285977/150‑1 (MQ=255)
                                                                                                         aTGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGAtact                                         >  2:1148189/1‑150 (MQ=255)
                                                                                                              ggACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACTACCta                                    <  1:3228397/150‑1 (MQ=255)
                                                                                                                   cTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACTACCTATAAAg                               >  2:773676/1‑150 (MQ=255)
                                                                                                                     gAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACTACCTATAAAGGc                             >  1:840118/1‑150 (MQ=255)
                                                                                                                       gggTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACTACCTATAAAGGCaa                           <  1:3009444/150‑1 (MQ=255)
                                                                                                                        ggTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACTACCTATAAAGGCAAt                          <  2:4157250/150‑1 (MQ=255)
                                                                                                                           tCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACTACCTATAAAGGCAATACg                       <  2:2143454/150‑1 (MQ=255)
                                                                                                                             cttcttCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACTACCTATAAAGGCAATACGcc                     >  1:2897489/1‑150 (MQ=255)
                                                                                                                             cttcttCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACTACCTATAAAGGCAATACGcc                     >  2:4008047/1‑150 (MQ=255)
                                                                                                                               tcttctCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACTACCTATAAAGGCAATACGCCTg                   >  1:2187543/1‑150 (MQ=255)
                                                                                                                               tcttctCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACTACCTATAAAGGCAATACGCCTg                   <  2:1397075/150‑1 (MQ=255)
                                                                                                                                    tcGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACTACCTATAAAGGCAATACGCCTGCACAg              >  2:4289925/1‑150 (MQ=255)
                                                                                                                                       gTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACTACCTATAAAGGCAATACGCCTGCACAGGTg           >  2:4795609/1‑150 (MQ=255)
                                                                                                                                       gTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACTACCTATAAAGGCAATACGCCTGCACAGGTg           <  2:4372026/150‑1 (MQ=255)
                                                                                                                                         tGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACTACCTATAAAGGCAATACGCCTGCACAGGTGcc         >  1:1848493/1‑150 (MQ=255)
                                                                                                                                          gAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACTACCTATAAAGGCAATACGCCTGCACAGGTGCCa        <  2:238864/150‑1 (MQ=255)
                                                                                                                                            aGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACTACCTATAAAGGCAATACGCCTGCACAGGTGCCaaa      >  1:1452373/1‑150 (MQ=255)
                                                                                                                                            aGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACTACCTATAAAGGCAATACGCCTGCACAGGTGCCaaa      >  1:4030791/1‑150 (MQ=255)
                                                                                                                                             ggTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACTACCTATAAAGGCAATACGCCTGCACAGGTGCCaaaa     >  1:447338/1‑150 (MQ=255)
                                                                                                                                              gTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACTACCTATAAAGGCAATACGCCTGCACAGGTGCCAAAac    >  1:155709/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACTACCTATAAAGGCAATACGCCTGCACAGGTGCCAAAACAC  >  CP009273/165612‑165905

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: