Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,897,754 T→G 8.6% T402P (ACG→CCG)  bglF ← fused beta‑glucoside‑specific PTS enzymes: IIA component/IIB component/IIC component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,897,7540TG 8.6% 677.5 / 53.7 240T402P (ACG→CCG) bglFfused beta‑glucoside‑specific PTS enzymes: IIA component/IIB component/IIC component
Reads supporting (aligned to +/- strand):  ref base T (117/102);  new base G (8/13);  total (125/115)
Fisher's exact test for biased strand distribution p-value = 2.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 243 total aligned reads displayed.

GAAATGAAGCATCACCGTACCGACAAATGCGCAACCGATGGCAATGACACCGCCAATAACGCTGGCCCAGACGGTGAAATCAATTCCCGTTGACGGGATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATACCAAACAGACTCGTCAACGCCGCTGATCCCGCCACCACTTTTT  >  CP009273/3897605‑3897903
                                                                                                                                                     |                                                                                                                                                     
  aatgaaGCATCACCGTACCGACAAATGCGCAACCGATGGCAATGACACCGCCAATAACGCTGGCCCAGACGGTGAAATCAATTCCCGTTGACGGGATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGttt                                                                                                                                                     >  2:2068021/1‑150 (MQ=255)
     gaaGCATCACCGTACCGACAAATGCGCAACCGATGGCAATGACACCGCCAATAACGCTGGCCCAGACGGTGAAATCAATTCCCGTTGACGGGATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTgcg                                                                                                                                                  <  1:2666759/150‑1 (MQ=255)
        gCATCACCGTACCGACAAATGCGCAACCGATGGCAATGACACCGCCAATAACGCTGGCCCAGACGGTGAAATCAATTCCCGTTGACGGGATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGt                                                                                                                                               >  1:2156034/1‑150 (MQ=255)
        gCATCACCGTACCGACAAATGCGCAACCGATGGCAATGACACCGCCAATAACGCTGGCCCAGACGGTGAAATCAATTCCCGTTGACGGGATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGt                                                                                                                                               >  1:1841180/1‑150 (MQ=255)
        gCATCACCGTACCGACAAATGCGCAACCGATGGCAATGACACCGCAAATAACGCTGGCCCAGACGGTGAAATCAATTCCCGTTGACGGGATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGt                                                                                                                                               >  1:4475313/1‑150 (MQ=255)
            cACCGTACCGACAAATGCGCAACCGATGGCAATGACACCGCCAATAACGCTGGCCCAGACGGTGAAATCAATTCCCGTTGACGGGATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGcc                                                                                                                                           <  2:3826335/150‑1 (MQ=255)
            cACCGTACCGACAAATGCGCAACCGATGGCAATGACACCGCCAATAACGCTGGCCCAGACGGTGAAATCAATTCCCGTTGACGGGATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGcc                                                                                                                                           <  2:3787482/150‑1 (MQ=255)
            cACCGTACCGACAAATGCGCAACCGATGGCAATGACACCGCCAATAACGCTGGCCCAGACGGTGAAATCAATTCCCGTTGACGGGATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGcc                                                                                                                                           >  2:1653934/1‑150 (MQ=255)
             aCCGTACCGACAAATGCGCAACCGATGGCAATGACACCGCCAATAACGCTGGCCCAGACGGTGAAATCAATTCCCGTTGACGGGATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCa                                                                                                                                          >  1:2911694/1‑150 (MQ=255)
                     gACAAATGCGCAACCGATGGCAATGACACCGCCAATAACGCTGGCCCAGACGGTGAAATCAATTCCCGTTGACGGGATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAAtggt                                                                                                                                  >  1:235432/1‑150 (MQ=255)
                         aaTGCGCAACCGATGGCAATGACACCGCCAATAACGCTGGCCCAGACGGTGAAATCAATTCCCGTTGACGGGATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGcc                                                                                                                              <  2:794216/150‑1 (MQ=255)
                         aaTGCGCAACCGATGGCAATGACACCGCCAATAACGCTGGCCCAGACGGTGAAATCAATTCCCGTTGACGGGATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGcc                                                                                                                              <  2:794183/150‑1 (MQ=255)
                         aaTGCGCAACCGATGGCAATGACACCGCCAATAACGCTGGCCCAGACGGTGAAATCAATTCCCGTTGACGGGATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGcc                                                                                                                              <  2:3695654/150‑1 (MQ=255)
                         aaTGCGCAACCGATGGCAATGACACCGCCAATAACGCTGGCCCAGACGGTGAAATCAATTCCCGTTGACGGGATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGcc                                                                                                                              >  2:2856182/1‑150 (MQ=255)
                         aaTGCGCAACCGATGGCAATGACACCGCCAATAACGCTGGCCCAGACGGTGAAATCAATTCCCGTTGACGGGATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGcc                                                                                                                              >  2:2112025/1‑150 (MQ=255)
                                aaCCGATGGCAATGACACCGCCAATAACGCTGGCCCAGACGGTGAAATCAATTCCCGTTGACGGGATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAg                                                                                                                       >  1:1687201/1‑150 (MQ=255)
                                aaCCGATGGCAATGACACCGCCAATAACGCTGGCCCAGACGGTGAAATCAATTCCCGTTGACGGGATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGGTTGCGCGTAGCCAATAATGGTGGCCCCCAAAg                                                                                                                       >  1:3578435/1‑150 (MQ=255)
                                  ccGATGGCAATGACACCGCCAATAACGCTGGCCCAGACGGTGAAATCAATTCCCGTTGACGGGATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGcc                                                                                                                     >  2:2882390/1‑150 (MQ=255)
                                    gATGGCAATGACACCGCCAATAACGCTGGCCCAGACGGTGAAATCAATTCCCGTTGACGGGATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGcccc                                                                                                                   >  2:337188/1‑150 (MQ=255)
                                    gATGGCAATGACACCGCCAATAACGCTGGCCCAGACGGTGAAATCAATTCCCGTTGACGGGATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGcccc                                                                                                                   >  1:2047448/1‑150 (MQ=255)
                                    gATGGCAATGACACCGCCAATAACGCTGGCCCAGACGGTGAAATCAATTCCCGTTGACGGGATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGGTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGcccc                                                                                                                   <  2:2365710/150‑1 (MQ=255)
                                      tGGCAATGACACCGCCAATAACGCTGGCCCAGACGGTGAAATCAATTCCCGTTGACGGGATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCAc                                                                                                                 >  1:1706124/1‑150 (MQ=255)
                                         cAATGACACCGCCAATAACGCTGGCCCAGACGGTGAAATCAATTCCCGTTGACGGGATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGa                                                                                                              >  2:317056/1‑150 (MQ=255)
                                          aaTGACACCGCCAATAACGCTGGCCCAGACGGTGAAATCAATTCCCGTTGACGGGATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGAt                                                                                                             <  2:729211/150‑1 (MQ=255)
                                          aaTGACACCGCCAATAACGCTGGCCCAGACGGTGAAATCAATTCCCGTTGACGGGATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGAt                                                                                                             <  2:3046817/150‑1 (MQ=255)
                                            tGACACCGCCAATAACGCTGGCCCAGACGGTGAAATCAATTCCCGTTGACGGGATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATAc                                                                                                           <  1:2137549/150‑1 (MQ=255)
                                            tGACACCGCCAATAACGCTGGCCCAGACGGTGAAATCAATTCCCGTTGACGGGATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATAc                                                                                                           <  1:775100/150‑1 (MQ=255)
                                              acacCGCCAATAACGCTGGCCCAGACGGTGAAATCAATTCCCGTTGACGGGATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAg                                                                                                         <  1:4171408/150‑1 (MQ=255)
                                                 ccGCCAATAACGCTGGCCCAGACGGTGAAATCAATTCCCGTTGACGGGATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCg                                                                                                      <  1:1080364/150‑1 (MQ=255)
                                                  cGCCAATAACGCTGGCCCAGACGGTGAAATCAATTCCCGTTGACGGGATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGa                                                                                                     >  2:4284968/1‑150 (MQ=255)
                                                   gCCAATAACGCTGGCCCAGACGGTGAAATCAATTCCCGTTGACGGGATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGAt                                                                                                    <  2:1535889/150‑1 (MQ=255)
                                                     cAATAACGCTGGCCCAGACGGTGAAATCAATTCCCGTTGACGGGATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATaa                                                                                                  <  2:2180145/150‑1 (MQ=255)
                                                     cAATAACGCTGGCCCAGACGGTGAAATCAATTCCCGTTGACGGGATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATaa                                                                                                  <  1:1444796/150‑1 (MQ=255)
                                                      aaTAACGCTGGCCCAGACGGTGAAATCAATTCCCGTTGACGGGATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGGTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAAc                                                                                                 >  1:1001364/1‑150 (MQ=255)
                                                       aTAACGCTGGCCCAGACGGTGAAATCAATTCCCGTTGACGGGATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACa                                                                                                >  1:230672/1‑150 (MQ=255)
                                                        tAACGCTGGCCCAGACGGTGAAATCAATTCCCGTTGACGGGATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACaa                                                                                               <  2:4232524/150‑1 (MQ=255)
                                                              tGGCCCAGACGGTGAAATCAATTCCCGTTGACGGGATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGt                                                                                         <  1:1722950/150‑1 (MQ=255)
                                                              tGGCCCAGACGGTGAAATCAATTCCCGTTGACGGGATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGt                                                                                         >  1:4218430/1‑150 (MQ=255)
                                                               ggCCCAGACGGTGAAATCAATTCCCGTTGACGGGATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTa                                                                                        >  1:1168809/1‑150 (MQ=255)
                                                               ggCCCAGACGGTGAAATCAATTCCCGTTGACGGGATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTa                                                                                        >  1:3442112/1‑150 (MQ=255)
                                                                 cccAGACGGTGAAATCAATTCCCGTTGACGGGATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACt                                                                                      <  1:3201822/150‑1 (MQ=255)
                                                                  ccAGACGGTGAAATCAATTCCCGTTGACGGGATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACtt                                                                                     <  2:3955872/150‑1 (MQ=255)
                                                                  ccAGACGGTGAAATCAATTCCCGTTGACGGGATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGGTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACtt                                                                                     >  1:4341172/1‑150 (MQ=255)
                                                                   cAGACGGTGAAATCAATTCCCGTTGACGGGATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTa                                                                                    >  1:485491/1‑150 (MQ=255)
                                                                     gACGGTGAAATCAATTCCCTTTGACGGGATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGGTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTAcg                                                                                  <  2:3467381/150‑1 (MQ=255)
                                                                     gACGGTGAAATCAAGTCCCGTTGACGCGAGGGTTTGCATGAAGGTGAAAATACATGGCAAACCAAAGGAGTAGACTTTCGGTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGCGTACTTAcg                                                                                  <  2:3467372/150‑1 (MQ=255)
                                                                          tGAAATCAATTCCCGTTGACGGGATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTACGCGGCa                                                                             >  1:3388201/1‑150 (MQ=255)
                                                                            aaaTCAATTCCCGTTGACGGGATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTACGCGGCAgg                                                                           >  1:437306/1‑150 (MQ=255)
                                                                             aaTCAATTCCCGTTGACGGGATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTACGCGGCAGGt                                                                          <  1:1797318/150‑1 (MQ=255)
                                                                             aaTCAATTCCCGTTGACGGGATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTACGCGGCAGGt                                                                          >  2:1196810/1‑150 (MQ=255)
                                                                              aTCAATTCCCGTTGACGGGATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGGTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTACGCGGCAGGtt                                                                         <  2:2648177/150‑1 (MQ=255)
                                                                               tCAATTCCCGTTGACGGGATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTACGCGGCAGGTTg                                                                        <  1:11555/150‑1 (MQ=255)
                                                                               tCAATTCCCGTTGACGGGATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTACGCGGCAGGTTg                                                                        <  1:11824/150‑1 (MQ=255)
                                                                               tCAATTCCCGTTGACGGGATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTACGCGGCAGGTTg                                                                        <  1:11700/150‑1 (MQ=255)
                                                                                cAATTCCCGTTGACGGGATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTACGCGGCAGGTTGa                                                                       <  1:2404971/150‑1 (MQ=255)
                                                                                     cccGTTGACGGGATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGGTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTACGCGGCAGGTTGACGCCa                                                                  <  2:1962437/150‑1 (MQ=255)
                                                                                         ttGACGGGATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTACGCGGCAGGTTGACGCCAtata                                                              >  1:175235/1‑150 (MQ=255)
                                                                                         ttGACGGGATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTACGCGGCAGGTTGACGCCAtata                                                              >  1:2047107/1‑150 (MQ=255)
                                                                                           gACGGGATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATAcc                                                            <  1:2926673/150‑1 (MQ=255)
                                                                                              gggATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTTGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCt                                                         >  1:1643128/1‑150 (MQ=255)
                                                                                              gggATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCt                                                         >  1:3551684/1‑150 (MQ=255)
                                                                                               ggATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTg                                                        <  1:3516815/150‑1 (MQ=255)
                                                                                                   ggTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTc                                                    >  1:3208874/1‑150 (MQ=255)
                                                                                                   ggTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGGTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTc                                                    <  1:689505/150‑1 (MQ=255)
                                                                                                    gTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCg                                                   >  2:1380107/1‑150 (MQ=255)
                                                                                                    gTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCg                                                   >  2:3370955/1‑150 (MQ=255)
                                                                                                     tttGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCgg                                                  <  1:2790287/150‑1 (MQ=255)
                                                                                                       tGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTg                                                <  1:1869175/150‑1 (MQ=255)
                                                                                                       tGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTg                                                <  2:621462/150‑1 (MQ=255)
                                                                                                        gCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGa                                               <  1:57855/150‑1 (MQ=255)
                                                                                                           tGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATAc                                            >  2:3959887/1‑150 (MQ=255)
                                                                                                           tGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATAc                                            >  2:1819077/1‑150 (MQ=255)
                                                                                                           tGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATAc                                            >  1:3212774/1‑150 (MQ=255)
                                                                                                           tGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATAc                                            >  2:2666128/1‑150 (MQ=255)
                                                                                                                   aaaaTACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATACCAAACAGa                                    <  1:3761664/150‑1 (MQ=255)
                                                                                                                    aaaTACTTGGCAAACCAAAGGAGTAGACTTTCGGTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATACCAAACAGAc                                   >  2:4480367/1‑150 (MQ=255)
                                                                                                                     aaTACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATACCAAACAGACt                                  >  2:520169/1‑150 (MQ=255)
                                                                                                                      aTACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATACCAAACAGACTc                                 <  2:1708263/150‑1 (MQ=255)
                                                                                                                       tACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATACCAAACAGACTCg                                <  1:217592/150‑1 (MQ=255)
                                                                                                                         cTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATACCAAACAGACTCGTc                              <  2:1039139/150‑1 (MQ=255)
                                                                                                                          ttGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATACCAAACAGACTCGTCa                             >  1:2033085/1‑150 (MQ=255)
                                                                                                                          ttGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATACCAAACAGACTCGTCa                             >  1:3390881/1‑150 (MQ=255)
                                                                                                                           tGGCAAACCAAAGGAGTAGACTATCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGATACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATACCAAACAGACTCGTCaa                            >  2:2904090/1‑150 (MQ=255)
                                                                                                                            ggCAAACCAAAGGCGTAGACTTTCGTTTGCGCGTAGCCAATAATGGGGGCCCCCACAGCCCCACTGATACAGGCGATAACAAAGGGGAACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATACCAAACAGACTCGTCAAc                           <  1:677898/150‑1 (MQ=255)
                                                                                                                              cAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATACCAAACAGACTCGTCAAcgc                         >  1:2191716/1‑150 (MQ=255)
                                                                                                                               aaaCCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATACCAAACAGACTCGTCAAcgcc                        <  1:2318703/150‑1 (MQ=255)
                                                                                                                                aaCCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATACCAAACAGACTCGTCAAcgccg                       >  1:556521/1‑150 (MQ=255)
                                                                                                                                  ccAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATACCAAACAGACTCGTCAACGACGCt                     <  1:1644501/150‑1 (MQ=255)
                                                                                                                                     aaGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATACCAAACAGACTCGTCAACGCCGCTGAt                  >  1:2530352/1‑150 (MQ=255)
                                                                                                                                     aaGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATACCAAACAGACTCGTCAACGCCGCTGAt                  >  1:419732/1‑150 (MQ=255)
                                                                                                                                      aGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATACCAAACAGACTCGTCAACGCCGCTGATc                 >  1:535556/1‑150 (MQ=255)
                                                                                                                                        gAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATACCAAACAGACTCGTCAACGCCGCTGATccc               >  1:547980/1‑150 (MQ=255)
                                                                                                                                         aGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATACCAAACAGACTCGTCAACGCCGCTGATCCCg              >  1:413997/1‑150 (MQ=255)
                                                                                                                                         aGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATACCAAACAGACTCGTCAACGCCGCTGATCCCg              <  2:175235/150‑1 (MQ=255)
                                                                                                                                          gTAGACTTTCGGTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATACCAAACAGACTCGTCAACGCCGCTGATCCCGc             <  2:2358796/150‑1 (MQ=255)
                                                                                                                                            aGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATACCAAACAGACTCGTCAACGCCGCTGATCCCGcca           >  1:3989168/1‑150 (MQ=255)
                                                                                                                                             gACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATACCAAACAGACTCGTCAACGCCGCTGATCCCGccac          >  1:1347620/1‑150 (MQ=255)
                                                                                                                                              aCTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATACCAAACAGACTCGTCAACGCCGCTGATCCCGccacc         >  2:2484920/1‑150 (MQ=255)
                                                                                                                                               cTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATACCAAACAGACTCGTCAACGCCGCTGATCCCGccacca        <  2:4341287/150‑1 (MQ=255)
                                                                                                                                                    gTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATACCAAACAGACTCGTCAACGCCGCTGATCCCGCCACCACtttt   >  2:1280718/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAAATGAAGCATCACCGTACCGACAAATGCGCAACCGATGGCAATGACACCGCCAATAACGCTGGCCCAGACGGTGAAATCAATTCCCGTTGACGGGATGGTTTGCATGAAGGTGAAAATACTTGGCAAACCAAAGGAGTAGACTTTCGTTTGCGCGTAGCCAATAATGGTGGCCCCCAAAGCCCCACTGATACAGGCGATAACAAAGGGGTACTTACGCGGCAGGTTGACGCCATATACCGCTGGTTCGGTGATACCAAACAGACTCGTCAACGCCGCTGATCCCGCCACCACTTTTT  >  CP009273/3897605‑3897903

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: