Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 745,893 G→T 6.0% S578* (TCG→TAG)  ybgQ ← putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273745,8930GT 6.0% 1200.3 / 36.4 315S578* (TCG→TAG) ybgQputative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base G (138/158);  new base T (12/7);  total (150/165)
Fisher's exact test for biased strand distribution p-value = 2.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 316 total aligned reads displayed.

TGGCGGCATAGGTACCAGAAATATCCCACTCACCCGCTGAACTCAGGTGCTGATAGTTACCGCTCACCTGGGCTCCATTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATAGCCAACCCGACCACCGTTACCGAAGGGGAGCGAAATCGACAGATAAATCTGGTTGTCGCGATCTTTATCTTCGTAGTGGGTCGTATTGAATGACGTCGAAATCGA  >  CP009273/745744‑746041
                                                                                                                                                     |                                                                                                                                                    
tGGCGGCATAGGTACCAGAAATATCCCACTCACCCGCTGAACTCAGGTGCTGATAGTTACCGCTCACCTGGGCTCCATTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACg                                                                                                                                                      >  1:1547167/1‑150 (MQ=255)
    ggcATAGGTACCAGAAATATCCCACTCACCCGCTGAACTCAGGTGCTGATAGTTACCGCTCACCTGGGCTTCATTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACAt                                                                                                                                                  >  1:2434345/1‑150 (MQ=255)
    ggcATAGGTACCAGAAATATCCCACTCACCCGCTGAACTCAGGTGCTGATAGTTACCGCTCACCTGGGCTCCATTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACAt                                                                                                                                                  >  2:3644043/1‑150 (MQ=255)
      cATAGGTACCAGAAATATCCCACTCACCCGCTGAACTCAGGTGCTGATAGTTACCGCTCACCTGGGCTCCATTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGc                                                                                                                                                >  2:2334427/1‑150 (MQ=255)
              ccAGAAATATCCCACTCACCCGCTGAACTCAGGTGCTGATAGTTACCGCTCACCTGGGCTCCATTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATg                                                                                                                                        >  2:1082050/1‑150 (MQ=255)
                aGAAATATCCCACTCACCCGCTGAACTCAGGTGCTGATAGTTACCGCTCACCTGGGCTCCATTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGt                                                                                                                                      <  2:1640507/150‑1 (MQ=255)
                 gAAATATCCCACTCACCCGCTGAACTCAGGTGCTGATAGTTACCGCTCACCTGGGCTCCATTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTg                                                                                                                                     >  2:4701955/1‑150 (MQ=255)
                 gAAATATCCCACTCACCCGCTGAACTCAGGTGCTGATAGTTACCGCTCACCTGGGCTCCATTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTg                                                                                                                                     <  1:2104638/150‑1 (MQ=255)
                  aaaTATCCCACTCACCCGCTGAACTCAGGTGCTGATAGTTACCGCTCACCTGGGCTCCATTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGc                                                                                                                                    >  1:3072624/1‑150 (MQ=255)
                   aaTATCCCACTCACCCGCTGAACTCAGGTGCTGATAGTTACCGCTCACCTGGGCTCCATTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCt                                                                                                                                   <  2:2045495/150‑1 (MQ=255)
                   aaTATCCCACTCACCCGCTGAACTCAGGTGCTGATAGTTACCGCTCACCTGGGCTCCATTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCt                                                                                                                                   <  2:898009/150‑1 (MQ=255)
                      atCCCACTCACCCGCTGAACTCAGGTGCTGATAGTTACCGCTCACCTGGGCTCCATTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCtgtg                                                                                                                                >  1:2947912/1‑150 (MQ=255)
                       tCCCACTCACCCGCTGAACTCAGGTGCTGATAGTTACCGCTCACCTGGGCTCCATTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTgg                                                                                                                               >  1:1747756/1‑150 (MQ=255)
                        cccACTCACCCGCTGAACTCAGGTGCTGATAGTTACCGCTCACCTGGGCTCCATTGTCCGGACGATCGGATTGCAGTCCAGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGc                                                                                                                              >  2:2221507/1‑150 (MQ=255)
                         ccACTCACCCGCTGAACTCAGGTGCTGATAGTTACCGCTCACCTGGGCTCCATTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCt                                                                                                                             <  2:1898055/150‑1 (MQ=255)
                         ccACTCACCCGCTGAACTCAGGTGCTGATAGTTACCGCTCACCTGGGCTCCATTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCt                                                                                                                             <  2:2721212/150‑1 (MQ=255)
                         ccACTCACCCGCTGAACTCAGGTGCTGATAGTTACCGCTCACCTGGGCTCCATTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCt                                                                                                                             >  2:2508588/1‑150 (MQ=255)
                           aCTCACCCGCTGAACTCAGGTGCTGATAGTTACCGCTCACCTGGGCTCCATTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACTACATGCGGTGTATGGTGCTGTGGCTAc                                                                                                                           >  2:1998320/1‑150 (MQ=255)
                           aCTCACCCGCTGAACTCAGGTGCTGATAGTTACCGCTCACCTGGGCTCCATTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTAc                                                                                                                           >  1:506798/1‑150 (MQ=255)
                            cTCACCCGCTGAACTCAGGTGCTGATAGTTACCGCTCACCTGGGCTCCATTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACt                                                                                                                          <  1:4342309/150‑1 (MQ=255)
                             tCACCCGCTGAACTCAGGTGCTGATAGTTACCGCTCACCTGGGCTCCATTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTg                                                                                                                         <  2:4607645/150‑1 (MQ=255)
                              cACCCGCTGAACTCAGGTGCTGATAGTTACCGCTCACCTGGGCTCCATTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGt                                                                                                                        <  1:1530804/150‑1 (MQ=255)
                                cccGCTGAACTCAGGTGCTGATAGTTACCGCTCACCTGGGCTCCATTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGttt                                                                                                                      <  2:810406/150‑1 (MQ=255)
                                 ccGCTGAACTCAGGTGCTGATAGTTACCGCTCACCTGGGCTCCATTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGtttt                                                                                                                     >  1:2084628/1‑150 (MQ=255)
                                 ccGCTGAACTCAGGTGCTGATAGTTACCGCTCACCTGGGCTCCATTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGtttt                                                                                                                     <  2:4010256/150‑1 (MQ=255)
                                  cGCTGAACTCAGGTGCTGATAGTTACCGCTCACCTGGGCTCCATTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTg                                                                                                                    <  2:2930828/150‑1 (MQ=255)
                                  cGCTGAACTCAGGTGCTGATAGTTACCGCTCACCTGGGCTCCATTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTg                                                                                                                    <  1:2912499/150‑1 (MQ=255)
                                    cTGAACTCAGGTGCTGATAGTTACCGCTCACCTGGGCTCCATTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCa                                                                                                                  >  1:44553/1‑150 (MQ=255)
                                      gAACTCAGGTGCTGATAGTTACCGCTCACCTGGGCTCCATTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATg                                                                                                                <  1:3698671/150‑1 (MQ=255)
                                      gAACTCAGGTGCTGATAGTTACCGCTCACCTGGGCTCCATTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATg                                                                                                                >  2:3622288/1‑150 (MQ=255)
                                      gAACTCAGGTGCTGATAGTTACCGCTCACCTGGGCTCCATTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATg                                                                                                                >  2:382905/1‑150 (MQ=255)
                                      gAACTCAGGTGCTGATAGTTACCGCTCACCTGGGCTCCATTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATg                                                                                                                >  1:2188889/1‑150 (MQ=255)
                                       aaCTCAGGTGCTGATAGTTACCGCTCACCTGGGCTCCATTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGt                                                                                                               >  1:1621578/1‑150 (MQ=255)
                                       aaCTCAGGTGCTGATAGTTACCGCTCACCTGGGCTCCATTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGt                                                                                                               >  2:929247/1‑150 (MQ=255)
                                       aaCTCAGGTGCTGATAGTTACCGCTCACCTGGGCTCCATTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGt                                                                                                               <  1:4610086/150‑1 (MQ=255)
                                       aaCTCAGGTGCTGATAGTTACCGCTCACCTGGGCTCCATTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGt                                                                                                               >  2:929305/1‑150 (MQ=255)
                                         cTCAGGTGCTGATAGTTACCGCTCACCTGGGCTCCATTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCa                                                                                                             <  1:3834445/150‑1 (MQ=255)
                                         cTCAGGTGCTGATAGTTACCGCTCACCTGGGCTCCATTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCa                                                                                                             >  1:1220301/1‑150 (MQ=255)
                                           cAGGTGCTGATAGTTACCGCTCACCTGGGCTCCATTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATa                                                                                                           <  1:635748/150‑1 (MQ=255)
                                           cAGGTGCTGATAGTTACCGCTCACCTGGGCTCCATTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATa                                                                                                           >  2:4345256/1‑150 (MQ=255)
                                           cAGGTGCTGATAGTTACCGCTCACCTGGGCTCCATTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATa                                                                                                           >  2:4348682/1‑150 (MQ=255)
                                           cAGGTGCTGATAGTTACCGCTCACCTGGGCTCCATTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATa                                                                                                           <  1:4404586/150‑1 (MQ=255)
                                           cAGGTGCTGATAGTTACCGCTCACCTGGGCTCCATTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATa                                                                                                           >  1:2213516/1‑150 (MQ=255)
                                            aGGTGCTGATAGTTACCGCTCACCTGGGCTCCATTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATAg                                                                                                          >  1:2485122/1‑150 (MQ=255)
                                                gCTGATAGTTACCGCTCACCTGGGCTCCATTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATAGCCaa                                                                                                      <  1:836969/150‑1 (MQ=255)
                                                 cTGATAGTTACCGCTCACCTGGGCTCCATTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATAGCCAAc                                                                                                     <  1:1422013/150‑1 (MQ=255)
                                                 cTGATAGTTACCGCTCACCTGGGCTCCATTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATAGCCAAc                                                                                                     <  2:2608782/150‑1 (MQ=255)
                                                  tGATAGTTACCGCTCACCTGGGCTCCATTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATAGCCAAcc                                                                                                    >  2:1564797/1‑150 (MQ=255)
                                                     tAGTTACCGCTCACCTGGGCTCCATTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATAGCCAACCCGa                                                                                                 >  1:3132248/1‑150 (MQ=255)
                                                       gTTACCGCTCACCTGGGCTCCATTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATAGCCAACCCGacc                                                                                               <  2:2054533/150‑1 (MQ=255)
                                                          aCCGCTCACCTGGGCTCCATTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATAGCCAACCCGaccacc                                                                                            >  2:2262077/1‑150 (MQ=255)
                                                           ccGCTCACCTGGGCTCCATTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATAGCCAACCCGACCACCg                                                                                           <  1:696026/150‑1 (MQ=255)
                                                              cTCACCTGGGCTCCATTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATAGCCAACCCGACCACCGTTa                                                                                        <  1:2260497/150‑1 (MQ=255)
                                                               tCACCTGGGCTCCATTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATAGCCAACCCGACCACCGTTAc                                                                                       >  2:3696271/1‑150 (MQ=255)
                                                                  ccTGGGCTCCATTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATAGCCAACCCGACCACCGTTACCGa                                                                                    >  1:3291056/1‑150 (MQ=255)
                                                                    tGGGCTCCATTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACTACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATAGCCAACCCGACCACCGTTACCGAAg                                                                                  >  1:4207262/1‑150 (MQ=255)
                                                                         tCCATTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATAGCCAACCCGACCACCGTTACCGAAGGGGAg                                                                             >  2:3135282/1‑150 (MQ=255)
                                                                          ccATTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATAGCCAACCCGACCACCGTTACCGAAGGGGAGc                                                                            <  1:2056496/150‑1 (MQ=255)
                                                                          ccATTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATAGCCAACCCGACCACCGTTACCGAAGGGGAGc                                                                            <  1:4073229/150‑1 (MQ=255)
                                                                            aTTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATAGCCAACCCGACCACCGTTACCGAAGGGGAGCGa                                                                          >  2:3113875/1‑150 (MQ=255)
                                                                            aTTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATAGCCAACCCGACCACCGTTACCGAAGGGGAGCGa                                                                          <  1:323839/150‑1 (MQ=255)
                                                                             ttGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATAGCCAACCCGACCACCGTTACCGAAGGGGAGCGaa                                                                         <  1:3558770/150‑1 (MQ=255)
                                                                               gTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATAGCCAACCCGACCACCGTTACCGAAGGGGAGCGAAAt                                                                       <  1:2136996/150‑1 (MQ=255)
                                                                                tCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATAGCCAACCCGACCACCGTTACCGAAGGGGAGCGAAATc                                                                      <  1:2522230/150‑1 (MQ=255)
                                                                                 ccGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATAGCCAACCCGACCACCGTTACCGAAGGGGAGCGAAATCg                                                                     <  1:3196641/150‑1 (MQ=255)
                                                                                   ggACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATAGCCAACCCGACCACCGTTACCGAAGGGGAGCGAAATCGAc                                                                   <  2:1329517/150‑1 (MQ=255)
                                                                                   ggACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATAGCCAACCCGACCACCGTTACCGAAGGGGAGCGAAATCGAc                                                                   <  1:912266/150‑1 (MQ=255)
                                                                                   ggACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATAGCCAACCCGACCACCGTTACCGAAGGGGAGCGAAATCGAc                                                                   <  2:506798/150‑1 (MQ=255)
                                                                                    gACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATAGCCAACCCGACCACCGTTACCGAAGGGGAGCGAAATCGACa                                                                  <  1:2480634/150‑1 (MQ=255)
                                                                                     aCGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATAGCCAACCCGACCACCGTTACCGAAGGGGAGCGAAATCGACAg                                                                 <  2:1219985/150‑1 (MQ=255)
                                                                                     aCGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATAGCCAACCCGACCACCGTTACCGAAGGGGAGCGAAATCGACAg                                                                 <  1:1788775/150‑1 (MQ=255)
                                                                                       gATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATAGCCAACCCGACCACCGTTACCGAAGGGGAGCGAAATCGACAGAt                                                               >  2:2193954/1‑150 (MQ=255)
                                                                                       gATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATAGCCAACCCGACCACCGTTACCGAAGGGGAGCGAAATCGACAGAt                                                               >  1:1680347/1‑150 (MQ=255)
                                                                                          cGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATAGCCAACCCGACCACCGTTACCGAAGGGGAGCGAAATCGACAGATaaa                                                            <  1:1880488/150‑1 (MQ=255)
                                                                                            gATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATAGCCAACCCGACCACCGTTACCGAAGGGGAGCGAAATCGACAGATAAATc                                                          <  2:3854067/150‑1 (MQ=255)
                                                                                               tGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATAGCCAACCCGACCACCGTTACCGAAGGGGAGCGAAATCGACAGATAAATCTgg                                                       >  2:4445969/1‑150 (MQ=255)
                                                                                                 cAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATAGCCAACCCGACCACCGTTACCGAAGGGGAGCGAAATCGACAGATAAATCTGGtt                                                     <  1:1014674/150‑1 (MQ=255)
                                                                                                  aGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACTACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATAGCCAACCCGACCACCGTTACCGAAGGGGAGCGAAATCGACAGATAAATCTGGTTg                                                    >  2:3217522/1‑150 (MQ=255)
                                                                                                     ccGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATAGCCAACCCGACCACCGTTACCGAAGGGGAGCGAAATCGACAGATAAATCTGGTTGTcg                                                 <  2:3320100/150‑1 (MQ=255)
                                                                                                       ggCAGACAGGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATAGCCAACCCGACCACCGTTACCGAAGGGGAGCGAAATCGACAGATAAATCTGGTTGTcgcg                                               <  1:1183420/150‑1 (MQ=255)
                                                                                                         cAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATAGCCAACCCGACCACCGTTACCGAAGGGGAGCGAAATCGACAGATAAATCTGGTTGTCGCGAt                                             <  1:4301531/150‑1 (MQ=255)
                                                                                                          aGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATAGCCAACCCGACCACCGTTACCGAAGGGGAGCGAAATCGACAGATAAATCTGGTTGTCGCGATc                                            <  1:3644043/150‑1 (MQ=255)
                                                                                                          aGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCGACTGTTTTGCATGTCATAGCCAACCCGACCACCGTTACCGAAGGGGAGCGAAATCGACAGATAAATCTGGGTGTCGCGATc                                            >  1:4608915/1‑150 (MQ=255)
                                                                                                              aTGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATAGCCAACCCGACCACCGTTACCGAAGGGGAGCGAAATCGACAGATAAATCTGGTTGTCGCGATCTTTa                                        <  2:2434345/150‑1 (MQ=255)
                                                                                                               tGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACTACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATAGCCAACCCGACCACCGTTACCGAAGGGGAGCGAAATCGACAGATAAATCTGGTTGTCGCGATCTTTAt                                       <  1:3799776/150‑1 (MQ=255)
                                                                                                               tGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATAGCCAACCCGACCACCGTTACCGAAGGGGAGCGAAATCGACAGATAAATCTGGTTGTCGCGATCTTTAt                                       >  2:105556/1‑150 (MQ=255)
                                                                                                               tGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATAGCCAACCCGACCACCGTTACCGAAGGGGAGCGAAATCGACAGATAAATCTGGTTGTCGCGATCTTTAt                                       <  1:4755643/150‑1 (MQ=255)
                                                                                                                    cAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATAGCCAACCCGACCACCGTTACCGAAGGGGAGCGAAATCGACAGATAAATCTGGTTGTCGCGATCTTTATCTTCg                                  <  2:4245711/150‑1 (MQ=255)
                                                                                                                    cAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATAGCCAACCCGACCACCGTTACCGAAGGGGAGCGAAATCGACAGATAAATCTGGTTGTCGCGATCTTTATCTTCg                                  <  1:3945232/150‑1 (MQ=255)
                                                                                                                     aGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATAGCCAACCCGACCACCGTTACCGAAGGGGAGCGAAATCGACAGATAAATCTGGTTGTCGCGATCTTTATCTTCGt                                 >  1:4573087/1‑150 (MQ=255)
                                                                                                                        tATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATAGCCAACCCGACCACCGTTACCGAAGGGGAGCGAAATCGACAGATAAATCTGGTTGTCGCGATCGTTATCTTCGTAGt                              >  2:2522745/1‑150 (MQ=255)
                                                                                                                         aTTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATAGCCAACCCGACCACCGTTACCGAAGGGGAGCGAAATCGACAGATAAATCTGGTTGTCGCGATCTTTATCTTCGTAGTg                             >  1:1934865/1‑150 (MQ=255)
                                                                                                                             cGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATAGCCAACCCGACCACCGTTACCGAAGGGGAGCGAAATCGACAGATAAATCTGGTTGTCGCGATCTTTATCTTCGTAGTGGGTc                         <  1:3042789/150‑1 (MQ=255)
                                                                                                                                 cATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATAGCCAACCCGACCACCGTTACCGAAGGGGAGCGAAATCGACAGATAAATCTGGTTGTCGCGATCTTTATCTTCGTAGTGGGTCGTAt                     >  1:2806766/1‑150 (MQ=255)
                                                                                                                                 cATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCGACTGTTTTGCATGTCATAGCCAACCCGACCACCGTTACCGAAGGGGAGCGAAATCGACAGATAAATCTGGTTGTCGCGATCTTTATCTTCGTAGTGGGTCGTAt                     >  2:1098921/1‑150 (MQ=255)
                                                                                                                                   tCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATAGCCAACCCGACCACCGTTACCGAAGGGGAGCGAAATCGACAGATAAATCTGGTTGTCGCGATCTTTATCTTCGTAGTGGGTCGTATTg                   <  2:2218159/150‑1 (MQ=255)
                                                                                                                                     cAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATAGCCAACCCGACCACCGTTACCGAAGGGGAGCGAAATCGACAGATAAATCTGGTTGTCGCGATCTTTATCTTCGTAGTGGGTCGTATTGaa                 >  1:4667515/1‑150 (MQ=255)
                                                                                                                                           aTCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATAGCCAACCCGACCACCGTTACCGAAGGGGAGCGAAATCGACAGATAAATCTGGTTGTCGCGATCTTTATCTTCGTAGTGGGTCGTATTGAATGACGt           >  2:623253/1‑150 (MQ=255)
                                                                                                                                             cGTTCCACTACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATAGCCAACCCGACCACCGTTACCGAAGGGGAGCGAAATCGACAGATAAATCTGGTTGTCGCGATCTTTATCTTCGTAGTGGGTCGTATTGAATGACGTCg         <  2:3415010/150‑1 (MQ=255)
                                                                                                                                               ttCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATAGCCAACCCGACCACCGTTACCGAAGGGGAGCGAAATCGACAGATAAATCTGGTTGTCGCGATCTTTATCTTCGTAGTGGGTCGTATTGAATGACGTCGaa       >  2:1970450/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGGCGGCATAGGTACCAGAAATATCCCACTCACCCGCTGAACTCAGGTGCTGATAGTTACCGCTCACCTGGGCTCCATTGTCCGGACGATCGGATTGCAGTCCGGCAGACATGCCCCAGCTATTACGTTCATCCAGCGTATCGTTCCACGACATGCGGTGTATGGTGCTGTGGCTACTGTTTTGCATGTCATAGCCAACCCGACCACCGTTACCGAAGGGGAGCGAAATCGACAGATAAATCTGGTTGTCGCGATCTTTATCTTCGTAGTGGGTCGTATTGAATGACGTCGAAATCGA  >  CP009273/745744‑746041

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: