Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,067,399 G→T 7.9% A794S (GCG→TCG)  flu → CP4‑44 prophage; antigen 43 (Ag43) phase‑variable biofilm formation autotransporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,067,3990GT 7.9% 1249.7 / 51.0 329A794S (GCG→TCG) fluCP4‑44 prophage; antigen 43 (Ag43) phase‑variable biofilm formation autotransporter
Reads supporting (aligned to +/- strand):  ref base G (139/164);  new base T (14/12);  total (153/176)
Fisher's exact test for biased strand distribution p-value = 5.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 335 total aligned reads displayed.

CCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGACAT  >  CP009273/2067251‑2067544
                                                                                                                                                    |                                                                                                                                                 
ccgtcTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACcgc                                                                                                                                                  >  1:5074642/1‑150 (MQ=255)
    ctcAGCATGCCGGGCGGTCATCTCGGTCCCGATAACCAGGGCGGTCGTGCCCGTGGGGCCACGCCGGACAGCAGCGGCAGCGATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCgggg                                                                                                                                              <  2:3110203/150‑1 (MQ=255)
       aGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGta                                                                                                                                           >  2:1386626/1‑150 (MQ=255)
        gCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGtat                                                                                                                                          <  1:4962507/150‑1 (MQ=255)
         cATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCTCGGGGGtata                                                                                                                                         >  2:2247084/1‑150 (MQ=255)
             cAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGt                                                                                                                                     <  2:577846/150‑1 (MQ=255)
              aGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGtg                                                                                                                                    <  2:3467744/150‑1 (MQ=255)
              aGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGtg                                                                                                                                    <  2:337503/150‑1 (MQ=255)
               gggCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCGCTGGATTCGTCCGGCGGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGtgc                                                                                                                                   >  2:4666713/1‑150 (MQ=255)
                 gCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGtgctg                                                                                                                                 >  1:164214/1‑150 (MQ=255)
                 gCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGtgctg                                                                                                                                 <  1:1681394/150‑1 (MQ=255)
                    gTCATCTCGGTCACGATAACCATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGctgctg                                                                                                                              <  2:4244574/150‑1 (MQ=255)
                      cATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGc                                                                                                                            >  2:2277477/1‑150 (MQ=255)
                         ctcGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCAt                                                                                                                         >  1:2468465/1‑150 (MQ=255)
                           cGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGGGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCAttc                                                                                                                       >  1:2831049/1‑150 (MQ=255)
                            ggTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCAttct                                                                                                                      >  2:3827130/1‑150 (MQ=255)
                            ggTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCAttct                                                                                                                      >  1:1095966/1‑150 (MQ=255)
                            ggTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCAttct                                                                                                                      >  1:3200600/1‑150 (MQ=255)
                             gTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCAttctt                                                                                                                     >  1:3747625/1‑150 (MQ=255)
                             gTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCAttctt                                                                                                                     >  2:2127617/1‑150 (MQ=255)
                              tCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCTCGGGGGTATATGGTGCTGCTGGCCAttcttc                                                                                                                    >  2:1070158/1‑150 (MQ=255)
                               cACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTcc                                                                                                                   >  2:4559060/1‑150 (MQ=255)
                                    tAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCTCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGa                                                                                                              <  2:3597330/150‑1 (MQ=255)
                                     aaCAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGAt                                                                                                             <  2:5101871/150‑1 (MQ=255)
                                      aCAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATg                                                                                                            <  1:1444324/150‑1 (MQ=255)
                                        aaTGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGtt                                                                                                          >  2:1644781/1‑150 (MQ=255)
                                           gtcGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATCTCTGTGACCGCGGGGGGATATGGTGCTGCTGGCCATTCTACCGATGATGTTAAg                                                                                                       >  2:4893121/3‑150 (MQ=255)
                                            gCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAgg                                                                                                      >  2:2596008/1‑150 (MQ=255)
                                             cGGTAGTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCTCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGga                                                                                                     <  1:3477764/150‑1 (MQ=255)
                                               gTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGgatg                                                                                                   <  2:1851804/150‑1 (MQ=255)
                                                tATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGgatga                                                                                                  <  1:3672124/150‑1 (MQ=255)
                                                tATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGgatga                                                                                                  <  1:4648256/150‑1 (MQ=255)
                                                   tGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCTCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGAtgatga                                                                                               >  1:3003008/1‑150 (MQ=255)
                                                    gCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGAc                                                                                              >  1:4976392/1‑150 (MQ=255)
                                                     cccGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACg                                                                                             <  1:1520554/150‑1 (MQ=255)
                                                     cccGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACg                                                                                             >  2:174091/1‑150 (MQ=255)
                                                      ccGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACgg                                                                                            >  2:2973580/1‑150 (MQ=255)
                                                       cGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGc                                                                                           <  2:684112/150‑1 (MQ=255)
                                                          ggggCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTcc                                                                                        <  2:4855413/150‑1 (MQ=255)
                                                           gggCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTccc                                                                                       <  1:457908/150‑1 (MQ=255)
                                                           gggCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTccc                                                                                       <  2:3731701/150‑1 (MQ=255)
                                                            ggCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCg                                                                                      >  2:2417037/1‑150 (MQ=255)
                                                            ggCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCg                                                                                      >  2:3779909/1‑150 (MQ=255)
                                                              ccACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTg                                                                                    <  1:3866099/150‑1 (MQ=255)
                                                                 cGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCTCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCg                                                                                 >  2:4771442/1‑150 (MQ=255)
                                                                   ccGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGc                                                                               <  1:347310/150‑1 (MQ=255)
                                                                     ggAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCAc                                                                             >  1:965838/1‑150 (MQ=255)
                                                                     ggAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGAACGCGGGGGTATATGGTGCTGCTTGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCAc                                                                             <  1:1062896/150‑1 (MQ=255)
                                                                           cagcGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCg                                                                       <  1:712687/150‑1 (MQ=255)
                                                                            agcGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCTCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCgg                                                                      <  1:707429/150‑1 (MQ=255)
                                                                             gcGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCggg                                                                     <  1:3336440/150‑1 (MQ=255)
                                                                              cGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGga                                                                    <  2:3329361/150‑1 (MQ=255)
                                                                                  aGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCTCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGgatgat                                                                >  1:3324478/1‑150 (MQ=255)
                                                                                      aTGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACCGAGGGTGCCGGTATGTCTGTGACCTCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCg                                                            <  1:1240297/150‑1 (MQ=255)
                                                                                        ggATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGc                                                          <  1:4157385/150‑1 (MQ=255)
                                                                                        ggATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGc                                                          >  2:2707069/1‑150 (MQ=255)
                                                                                        ggATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTGGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGc                                                          >  2:3096300/1‑150 (MQ=255)
                                                                                         gATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCa                                                         <  2:3148183/150‑1 (MQ=255)
                                                                                         gATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGAGGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCa                                                         >  2:2771973/1‑150 (MQ=255)
                                                                                           ttCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGc                                                       >  2:3018748/1‑150 (MQ=255)
                                                                                           ttCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGc                                                       >  1:3652642/1‑150 (MQ=255)
                                                                                             cgtccgtcTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCTCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCt                                                     >  1:2402125/1‑150 (MQ=255)
                                                                                             cgtccgtcTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCt                                                     <  2:5077339/150‑1 (MQ=255)
                                                                                             cgtccgtcTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCt                                                     >  2:189941/1‑150 (MQ=255)
                                                                                             cgtccgtcTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCt                                                     <  1:2220627/150‑1 (MQ=255)
                                                                                             cgtccgtcTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCt                                                     <  1:3577730/150‑1 (MQ=255)
                                                                                              gtccgtcTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTg                                                    >  2:1943515/1‑150 (MQ=255)
                                                                                              gtccgtcTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGACATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTg                                                    <  2:87935/150‑1 (MQ=255)
                                                                                               tccgtcTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTgg                                                   <  2:1306426/150‑1 (MQ=255)
                                                                                                  gtcTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCg                                                >  1:1406052/1‑150 (MQ=255)
                                                                                                       gAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATAc                                           >  2:825303/1‑150 (MQ=255)
                                                                                                       gAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATAc                                           >  2:3402598/1‑150 (MQ=255)
                                                                                                       gAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATAc                                           >  2:3949724/1‑150 (MQ=255)
                                                                                                        aGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGAGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATAcc                                          >  1:1210419/1‑150 (MQ=255)
                                                                                                         gggTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATACCt                                         >  2:5008543/1‑150 (MQ=255)
                                                                                                          ggTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTg                                        <  2:1252564/150‑1 (MQ=255)
                                                                                                            tGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGaa                                      <  1:3909483/150‑1 (MQ=255)
                                                                                                             gACCTAATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGAAt                                     <  2:4797832/150‑1 (MQ=255)
                                                                                                              aCCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGAATc                                    >  2:4343773/1‑150 (MQ=255)
                                                                                                              aCCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGAATc                                    >  2:3675919/1‑150 (MQ=255)
                                                                                                                 tgatgaGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGAATCTgg                                 <  2:2907358/150‑1 (MQ=255)
                                                                                                                 tgatgaGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGAATCTgg                                 <  1:3120204/150‑1 (MQ=255)
                                                                                                                 tgatgaGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGAATCTgg                                 >  2:3913551/1‑150 (MQ=255)
                                                                                                                 tgatgaGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGAATCTgg                                 <  1:1412544/150‑1 (MQ=255)
                                                                                                                   atgaGAACAGAGGTTGCCGGTTTGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGAATCTGGTa                               <  1:2588939/150‑1 (MQ=255)
                                                                                                                    tgaGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGAATCTGGTac                              <  2:3934008/150‑1 (MQ=255)
                                                                                                                    tgaGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGAATCTGGTac                              <  1:2157910/150‑1 (MQ=255)
                                                                                                                     gagaACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGAATCTGGTaca                             >  2:4074920/1‑150 (MQ=255)
                                                                                                                       gaACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGAATCTGGTacaca                           <  1:128514/150‑1 (MQ=255)
                                                                                                                         aCAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGAATCTGGTACACACg                         <  1:5116147/150‑1 (MQ=255)
                                                                                                                              ggTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGAATCTGGTACACACGtcctc                    >  1:4241560/1‑150 (MQ=255)
                                                                                                                              ggTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGAATCTGGTACACACGtcctc                    <  2:233697/150‑1 (MQ=255)
                                                                                                                               gTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGAATCTGGTACACACGtcctcc                   <  1:3828475/150‑1 (MQ=255)
                                                                                                                               gTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGAATCTGGTACACACGtcctcc                   >  1:3986521/1‑150 (MQ=255)
                                                                                                                                       tATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGAATCTGGTACACACGTCCTCCGGCCtgtg           <  1:217448/150‑1 (MQ=255)
                                                                                                                                        aTGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTgg          >  1:4810622/1‑150 (MQ=255)
                                                                                                                                         tGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTggg         >  2:2949921/1‑150 (MQ=255)
                                                                                                                                           tCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCt       >  1:1833281/1‑150 (MQ=255)
                                                                                                                                              gtgACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGAc    <  1:2328378/150‑1 (MQ=255)
                                                                                                                                              gtgACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGAc    >  1:1114538/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CCGTCTCAGCATTCAGGGCGGTCATCTCGGTCACGATAACAATGGCGGTATTGCCCGTGGGGCCACGCCGGAAAGCAGCGGCAGCTATGGATTCGTCCGTCTGGAGGGTGACCTGATGAGAACAGAGGTTGCCGGTATGTCTGTGACCGCGGGGGTATATGGTGCTGCTGGCCATTCTTCCGTTGATGTTAAGGATGATGACGGCTCCCGTGCCGGCACGGTCCGGGATGATGCCGGCAGCCTGGGCGGATACCTGAATCTGGTACACACGTCCTCCGGCCTGTGGGCTGACAT  >  CP009273/2067251‑2067544

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: