Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 165,759 G→T 57.1% 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,7590GT57.1% 61.8 / inf 328E597* (GAG→TAG) fhuAferrichrome outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base G (76/65);  new base T (95/92);  total (171/157)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 332 total aligned reads displayed.

CACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACTACCTATAAAGGCAATACGCCTGCACAGGTGCCAAAACACATG  >  CP009273/165610‑165908
                                                                                                                                                     |                                                                                                                                                     
    gTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGAt                                                                                                                                                   >  1:4094666/1‑150 (MQ=255)
     tCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATc                                                                                                                                                  >  2:377199/1‑150 (MQ=255)
     tCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATc                                                                                                                                                  <  1:1741940/150‑1 (MQ=255)
     tCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATc                                                                                                                                                  >  1:181852/1‑150 (MQ=255)
      cTAAATGTAAGCAGTATGAAGTCGGAGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTTCTGGTGCAGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGTGGGTTCCTTCTTCTCTGTTGAAGGTGGCTAGATcc                                                                                                                                                 <  1:3295289/150‑1 (MQ=255)
       tAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCcg                                                                                                                                                >  2:1002465/1‑150 (MQ=255)
        aaaGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCcgc                                                                                                                                               >  2:4799291/1‑150 (MQ=255)
        aaaGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCcgc                                                                                                                                               >  2:2187331/1‑150 (MQ=255)
         aaGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCcgcg                                                                                                                                              <  1:3107067/150‑1 (MQ=255)
           ggTAAGCAGTATGAAGTCGGCGTGAAATATGTACAGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCa                                                                                                                                            >  1:652097/1‑150 (MQ=255)
              aaGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGt                                                                                                                                         >  1:4209630/1‑150 (MQ=255)
               aGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTg                                                                                                                                        >  1:2270648/1‑150 (MQ=255)
                  aGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGATGAAGGTGGCTAGATCCGCGCACGTGGCg                                                                                                                                     >  1:1786466/1‑150 (MQ=255)
                      tGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGa                                                                                                                                 <  1:2936987/150‑1 (MQ=255)
                       gAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGAAGGCGGACCCTGAGGGTTCCTTCTACTCGGTTGAAGGTGGCTAGATCCGCTCACGTGGCGTAGaa                                                                                                                                >  2:5100148/1‑150 (MQ=255)
                        aaGTCGGCGTGAAATATGTACCGTAAGATCGTCCGATTGTAGTTACTGGTGCCGTGAATAATCTCACTAAAACCAACAACCTGATGGCGGACCCAGAGGGTTCCTACTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGaaa                                                                                                                               >  1:2668740/1‑150 (MQ=255)
                           tCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCg                                                                                                                            <  2:2367705/150‑1 (MQ=255)
                             ggCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGaa                                                                                                                          >  1:2490593/1‑150 (MQ=255)
                             ggCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGGAGAAATCGaa                                                                                                                          >  1:1010901/1‑150 (MQ=255)
                                     atatGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAgc                                                                                                                  >  2:702737/1‑150 (MQ=255)
                                     atatGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAgc                                                                                                                  >  2:5192701/1‑150 (MQ=255)
                                            ccGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTg                                                                                                           >  1:4760623/1‑150 (MQ=255)
                                              ggAAGATCGTCCGTTTGTTGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTc                                                                                                         <  2:3182960/150‑1 (MQ=255)
                                                acgATCGTCCGATTGTAGTTACTGGTGCCGTGTATACTCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGGGGCGTAGAACTCGAAGCGAAAGCGGCGCTGTCgg                                                                                                       >  2:4901385/3‑150 (MQ=255)
                                                aaGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCgg                                                                                                       >  2:84781/1‑150 (MQ=255)
                                                  gATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCTGACCCTGAGGGTTCCTTATTCTCGGTTGAAGGTGGCGAGATCAGCGCACGTGGCGTAGAAATCTAAGCGAAAGCGGCGCTGTCGGCg                                                                                                     >  2:2010228/1‑150 (MQ=255)
                                                   aTCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGa                                                                                                    >  1:1312833/1‑150 (MQ=255)
                                                    tCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGTGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGACGCGAAAGCGGCGCTGTCGGCGAg                                                                                                   <  1:4310195/150‑1 (MQ=255)
                                                      gTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAgtg                                                                                                 >  2:288768/1‑150 (MQ=255)
                                                      gTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAgtg                                                                                                 <  2:30960/150‑1 (MQ=255)
                                                       tCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAgtgt                                                                                                <  2:436378/150‑1 (MQ=255)
                                                        ccGATTGTATTTACTGGTGCCGTGTATAATCACACTAAAACCAACAACCTGTTGGCGGACCCTGTGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTACAAATCGACGCGAAAGCGGCGCTGTCGGCGAGTGtt                                                                                               <  1:2831781/150‑1 (MQ=255)
                                                           aTTGTAGTTACTGGTGCCGTGTATAAACTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTATAAATCGACGCGAAAGCGGCGCTGTCGGCGAGTGTTAAc                                                                                            <  1:1656593/150‑1 (MQ=255)
                                                            ttGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACg                                                                                           <  1:2187331/150‑1 (MQ=255)
                                                            ttGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACg                                                                                           >  2:4425684/1‑150 (MQ=255)
                                                                aGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGt                                                                                       >  1:2747741/1‑150 (MQ=255)
                                                                 gTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCTTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTc                                                                                      <  1:2152763/150‑1 (MQ=255)
                                                                 gTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTc                                                                                      <  2:3903992/150‑1 (MQ=255)
                                                                  ttACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCg                                                                                     <  2:2966810/150‑1 (MQ=255)
                                                                     cTGGTGCCGTGTATAATCTCACTAAAACCAACACCCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGtt                                                                                  <  2:3635318/150‑1 (MQ=255)
                                                                      tGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTc                                                                                 <  2:2573831/150‑1 (MQ=255)
                                                                         tGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTta                                                                              >  2:1295746/1‑150 (MQ=255)
                                                                         tGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTta                                                                              <  2:3345157/150‑1 (MQ=255)
                                                                          gCCGTGTATAATCTCACTAAAACCAACACCCTGTTGGCGGACCCTGAGGGTTCCTTCTTCTGGGTTGTTGGTGGCTTGATCCGCGCACGTGGCTTAGAAATCGACGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTtat                                                                             <  1:566100/150‑1 (MQ=255)
                                                                           ccGTTTATAATCTCACTAAAACCAACACCCTGATGGCGGACCCTGTGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGACGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTtata                                                                            <  1:51279/150‑1 (MQ=255)
                                                                           ccGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTtata                                                                            <  1:4884674/150‑1 (MQ=255)
                                                                            cGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATAc                                                                           >  2:5259367/1‑150 (MQ=255)
                                                                            cGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATAc                                                                           >  1:3196160/1‑150 (MQ=255)
                                                                            cGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATAc                                                                           >  2:3295464/1‑150 (MQ=255)
                                                                             gtgtATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACt                                                                          >  2:4604284/1‑150 (MQ=255)
                                                                              tgtATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACtt                                                                         >  2:2115524/1‑150 (MQ=255)
                                                                               gtATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTa                                                                        >  1:1378500/1‑150 (MQ=255)
                                                                               gtATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTa                                                                        >  2:3947310/1‑150 (MQ=255)
                                                                                tataATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGATGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGACGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTac                                                                       <  1:3844323/150‑1 (MQ=255)
                                                                                 ataATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTaca                                                                      >  2:4826703/1‑150 (MQ=255)
                                                                                   aaTCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACAcc                                                                    >  1:349464/1‑150 (MQ=255)
                                                                                   aaTCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACAcc                                                                    >  1:2410267/1‑150 (MQ=255)
                                                                                       tcACTAACACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATg                                                                <  2:4653391/150‑1 (MQ=255)
                                                                                       tcACTAAAACCAACAACCTGATGGCGGACCCTGTGGGTTCCTTCTTCTCGGTTGATGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATg                                                                <  1:2464581/150‑1 (MQ=255)
                                                                                        cACTAAAACCAACAACCTGATGGCGGACCCTGTGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGc                                                               <  1:1182115/150‑1 (MQ=255)
                                                                                        cACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGc                                                               >  1:3577953/1‑150 (MQ=255)
                                                                                        cACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGc                                                               <  2:650050/150‑1 (MQ=255)
                                                                                           tAAAACCAACAACCTGATGGCGGACCCTGTGGGTTCCTTCTTCTCGGTTGATGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGa                                                            <  1:863502/150‑1 (MQ=255)
                                                                                             aaaCCAACACCCTGAGGGCGGACCCTGTGGGTTCCTTCTTCTCGGTTGAAGGGGGCGAGATCCGCGCACGTGGCGTAGAAATCGACGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAAt                                                          <  1:3696161/150‑1 (MQ=255)
                                                                                             aaaCCAACAACCTGATGGCGGACCCTGTGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAAt                                                          <  2:3830355/150‑1 (MQ=255)
                                                                                             aaaCCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAAt                                                          >  1:4810409/1‑150 (MQ=255)
                                                                                              aaCCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATa                                                         <  2:2302696/150‑1 (MQ=255)
                                                                                               aCCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATac                                                        >  2:543928/1‑150 (MQ=255)
                                                                                                ccAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATaca                                                       <  2:1238455/150‑1 (MQ=255)
                                                                                                ccAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATaca                                                       >  2:1559901/1‑150 (MQ=255)
                                                                                                caaacaccCTGATGGCGGACCCTGAGGGTTCCTTCTTCGCGGTTGATGGTGGCTAGATCCGCGCACGGGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATaca                                                       <  1:3273302/148‑1 (MQ=255)
                                                                                                 caacacCCTGCTGGCGGACCCTGTGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGTGGAATacac                                                      <  1:3329111/150‑1 (MQ=255)
                                                                                                  aacaacCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATAcacc                                                     <  1:398374/150‑1 (MQ=255)
                                                                                                   acaacCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATAcacca                                                    >  2:2843028/1‑150 (MQ=255)
                                                                                                    caacCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATAcaccac                                                   <  2:2441518/150‑1 (MQ=255)
                                                                                                     cacCTGATGGCGGACCCTGAGGGTTCCTTCGTCTAGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACaccacc                                                  <  2:5219637/149‑1 (MQ=255)
                                                                                                         tGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGATGAAGGTGGCGAGATCCGCGCACGAGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTCTACTTACACCGATGCGGAATACACCACCGAta                                              >  2:1474747/1‑150 (MQ=255)
                                                                                                                ggACCCTGTGGGTTCCTTCTTCTCGTTTGATGGTGGCTAGATCCGCGCACGTGGCGTTTAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACTACCta                                       <  1:1896092/150‑1 (MQ=255)
                                                                                                                ggACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACTACCta                                       <  1:4367369/150‑1 (MQ=255)
                                                                                                                ggACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACTACCta                                       >  2:3259325/1‑150 (MQ=255)
                                                                                                                    ccTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGTCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAAAGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACAACCTATaaa                                   >  2:1519804/1‑150 (MQ=255)
                                                                                                                      tGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGACTACACCACCGATACTACCTATAAAgg                                 >  2:3832902/1‑150 (MQ=255)
                                                                                                                       gAGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGGTCTTATACTTACACCGATGCGGAATACACCACCGATACTACCTATAAAGGc                                >  1:4186887/1‑150 (MQ=255)
                                                                                                                        aGGGTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACTACCTATAAAGGCa                               >  2:487799/1‑150 (MQ=255)
                                                                                                                          ggTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACTACCTATAAAGGCAAt                             <  2:4144597/150‑1 (MQ=255)
                                                                                                                          ggTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACTACCTATAAAGGCAAt                             >  2:3160440/1‑150 (MQ=255)
                                                                                                                           gTTCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACTACCTATAAAGGCAATa                            >  2:4294591/1‑150 (MQ=255)
                                                                                                                             tCCTTCTTCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACTACCTATAAAGGCAATACg                          >  2:2808363/1‑150 (MQ=255)
                                                                                                                               cttcttCTCGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACTACCTATAAAGGCAATACGcc                        >  2:3493697/1‑150 (MQ=255)
                                                                                                                                 tcttctCGGTTGTAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACTACCTATAAAGGCAATACGCCTg                      <  2:3715207/150‑1 (MQ=255)
                                                                                                                                    tctcGGTTGATGGTGGCTAGATCCGCGCAGGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACTACCTATAAAGGCAATACGCCTGcac                   <  2:1827831/150‑1 (MQ=255)
                                                                                                                                     ctcGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACTCCCTATAAAGGCAATACGCCTGcaca                  >  2:408197/1‑150 (MQ=255)
                                                                                                                                       cGGTTGAAGGTGGCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACTACCTATAAAGGCAATACGCCTGCACAgg                <  2:3683288/150‑1 (MQ=255)
                                                                                                                                        ggTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACTACCTATAAAGGCAATACGCCTGCACAGGt               >  1:5287501/1‑150 (MQ=255)
                                                                                                                                        ggTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACTACCTATAAAGGCAATACGCCTGCACAGGt               >  1:3636300/1‑150 (MQ=255)
                                                                                                                                             aaGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACTACCTATAAAGGCAATACGCCTGCACTGGTGCCaa          >  2:2471151/1‑150 (MQ=255)
                                                                                                                                                  ggCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACTACCTATAAAGGCAATACGCCTGCACAGGTGCCAAAacac     <  2:1683665/150‑1 (MQ=255)
                                                                                                                                                  ggCTAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACTACCTATAAAGGCAATACGCCTGCACAGGTGCCAAAacac     <  2:2997721/150‑1 (MQ=255)
                                                                                                                                                  ggCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACTACCTATAAAGGCAATACGCCTGCACAGGTGCCAAAacac     <  1:1546948/150‑1 (MQ=255)
                                                                                                                                                     gagaTCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACTACCTATAAAGGCAATACGCCTGCACAGGTGCCAAAACACATg  <  2:91375/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CACCGTCTAAAGGTAAGCAGTATGAAGTCGGCGTGAAATATGTACCGGAAGATCGTCCGATTGTAGTTACTGGTGCCGTGTATAATCTCACTAAAACCAACAACCTGATGGCGGACCCTGAGGGTTCCTTCTTCTCGGTTGAAGGTGGCGAGATCCGCGCACGTGGCGTAGAAATCGAAGCGAAAGCGGCGCTGTCGGCGAGTGTTAACGTAGTCGGTTCTTATACTTACACCGATGCGGAATACACCACCGATACTACCTATAAAGGCAATACGCCTGCACAGGTGCCAAAACACATG  >  CP009273/165610‑165908

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