Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,243,278 G→A 13.4% S431S (AGC→AGT dhaM ← putative dihydroxyacetone‑specific PTS enzymes: HPr component/EI component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,243,2780GA13.4% 936.5 / 101.7 314S431S (AGC→AGTdhaMputative dihydroxyacetone‑specific PTS enzymes: HPr component/EI component
Reads supporting (aligned to +/- strand):  ref base G (125/146);  new base A (20/22);  total (145/169)
Fisher's exact test for biased strand distribution p-value = 8.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 322 total aligned reads displayed.

AACGATGCCATCCGAACAGTGGCTTAACCCTGACGGTTGAAACGTTGCGTTTTAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAGAATAGTTGGCGAGTTAAACTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTG  >  CP009273/1243129‑1243427
                                                                                                                                                     |                                                                                                                                                     
aaCGATGCCATCCGAACAGTGGCTTAACCCTGACGGTTGAAACGTTGCGTTTTAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGgcg                                                                                                                                                       >  2:165247/1‑150 (MQ=255)
  cGATGCCATCCGAACAGTGGCTTAACCCTGACGGTTGAAACGTTGCGTTTTAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCt                                                                                                                                                     >  2:4264057/1‑150 (MQ=255)
   gATGCCATCCGAACAGTGGCTTAACCCTGACGGTTGAAACGTTGCGTTTTAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCtg                                                                                                                                                    >  1:461426/1‑150 (MQ=255)
   gATGCCATCCGAACAGTGGCTTAACCCTGACGGTTGAAACGTTGCGTTTTAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCtg                                                                                                                                                    >  1:2900179/1‑150 (MQ=255)
      gCCATCCGAACAGTGGCTTAACCCTGACGGTTGAAACGTTGCGTTTTAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCGGTgg                                                                                                                                                 >  2:3991823/1‑150 (MQ=255)
        cATCCGAACAGTGGCTTAACCCTGACGGTTGAAACGTTGCGTTTTAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGa                                                                                                                                               <  2:1889794/150‑1 (MQ=255)
                cAGTGGCTTAACCCTGACGGTTGAAACGTTGCGTTTTAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGAc                                                                                                                                       <  1:888499/150‑1 (MQ=255)
                 aGTGGCTTAACCCTGACGGTTGAAACGTTGCGTTTTAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACt                                                                                                                                      >  2:2939808/1‑150 (MQ=255)
                   tGGCTTAACCCTGACGGTTGAAACGTTGCGTTTTAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTc                                                                                                                                    >  1:1413949/1‑150 (MQ=255)
                    ggCTTAACCCTGACGGTTGAAACGTTGCGTTTTAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCACTGTGGGATACCGGACTTcc                                                                                                                                   <  2:2568955/150‑1 (MQ=255)
                         aaCCCTGACGGTTGAAACGTTGCGTTTTAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATACAGCCAATCCCCCGGTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGgcgc                                                                                                                              >  2:4582090/1‑150 (MQ=255)
                          aCCCTGACGGTTGAAACGTTGCGTTTTAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCt                                                                                                                             >  2:1083456/1‑150 (MQ=255)
                          aCCCTGACGGTTGAAACGTTGCGTTTTAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCACTGTGGGATACCGGACTTCCGGCGCt                                                                                                                             <  2:2041175/150‑1 (MQ=255)
                           cccTGACGGTTGAAACGTTGCGTTTTAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTa                                                                                                                            <  1:2809324/150‑1 (MQ=255)
                           cccTGACGGTTGAAACGTTGCGTTTTAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTa                                                                                                                            <  1:4151856/150‑1 (MQ=255)
                              tGACGGTTGAAACGTTGCGTTTTAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAgg                                                                                                                         >  1:1299735/1‑150 (MQ=255)
                               gACGGTTGAAACGTTGCGTTTTAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGc                                                                                                                        >  1:5099695/1‑150 (MQ=255)
                               gACGGTTGAAACGTTGCGTTTTAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGc                                                                                                                        >  1:1213952/1‑150 (MQ=255)
                                 cGGTTGAAACGTTGCGTTTTAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAg                                                                                                                      >  1:1926300/1‑150 (MQ=255)
                                 cGGTTGAAACGTTGCGTTTTAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCACTGTGGGATACCGGACTTCCGGCGCTAAGGCAg                                                                                                                      <  2:5162484/150‑1 (MQ=255)
                                   gTTGAAACGTTGCGTTTTAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGAt                                                                                                                    <  1:1536109/150‑1 (MQ=255)
                                    ttGAAACGTTGCGTTTTAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATa                                                                                                                   <  2:353010/150‑1 (MQ=255)
                                     tGAAACGGTGCGTTTTAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGAGGTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATAc                                                                                                                  >  2:738018/1‑150 (MQ=255)
                                      gAAACGTTGCGTTTTAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATAcc                                                                                                                 >  1:210141/1‑150 (MQ=255)
                                        aaCGTTGCGTTTTAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCtt                                                                                                               >  2:5080890/1‑150 (MQ=255)
                                            ttGCGATTTAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTAc                                                                                                           <  2:1547147/150‑1 (MQ=255)
                                             tGCGTTTTAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACa                                                                                                          >  1:4643930/1‑150 (MQ=255)
                                              gCGTTTTAACGTCCAGCGTTAGCGTTTCGTCTGGTTGTATCGCATACAGTTTCGCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACaa                                                                                                         <  2:2683837/150‑1 (MQ=255)
                                               cGTTTTAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAAc                                                                                                        <  1:3360453/150‑1 (MQ=255)
                                                    tAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCACTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAAccgccg                                                                                                   >  1:4531227/1‑150 (MQ=255)
                                                     aaCGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCgg                                                                                                  >  1:2621828/1‑150 (MQ=255)
                                                     aaCGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCgg                                                                                                  >  1:5210880/1‑150 (MQ=255)
                                                      aCGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGa                                                                                                 <  1:4032916/150‑1 (MQ=255)
                                                       cGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGAt                                                                                                <  2:2141004/150‑1 (MQ=255)
                                                        gTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATc                                                                                               <  1:2612593/150‑1 (MQ=255)
                                                         tCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATcc                                                                                              <  1:863181/150‑1 (MQ=255)
                                                          ccAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCa                                                                                             >  1:2464377/1‑150 (MQ=255)
                                                          ccAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCa                                                                                             <  2:443595/150‑1 (MQ=255)
                                                           cAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAg                                                                                            <  2:486129/150‑1 (MQ=255)
                                                               gttagcgttTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGc                                                                                        >  1:2711642/1‑150 (MQ=255)
                                                                ttagcgttTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCa                                                                                       >  2:2469327/1‑150 (MQ=255)
                                                                ttagcgttTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCa                                                                                       >  2:21777/1‑150 (MQ=255)
                                                                 tagcgttTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAg                                                                                      >  1:4883549/1‑150 (MQ=255)
                                                                  agcgttTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGt                                                                                     <  2:3541078/150‑1 (MQ=255)
                                                                  agcgttTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGt                                                                                     <  1:3577466/150‑1 (MQ=255)
                                                                  agcgttTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCACTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGt                                                                                     >  1:632233/1‑150 (MQ=255)
                                                                   gcgttTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCACTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTa                                                                                    <  2:2066488/150‑1 (MQ=255)
                                                                   gagttTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTa                                                                                    <  1:2823562/148‑1 (MQ=255)
                                                                     gttTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCACTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACt                                                                                  <  2:4982062/150‑1 (MQ=255)
                                                                      tttCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACtg                                                                                 <  2:4518659/150‑1 (MQ=255)
                                                                       ttcttcTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCACTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACtgt                                                                                <  2:2121623/150‑1 (MQ=255)
                                                                        tcttctGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACtgtg                                                                               >  1:4745771/1‑150 (MQ=255)
                                                                         cttctGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTgg                                                                              <  2:4557836/150‑1 (MQ=255)
                                                                         cttctGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTgg                                                                              <  2:4761337/150‑1 (MQ=255)
                                                                            ctGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCACTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAg                                                                           >  1:4472676/1‑150 (MQ=255)
                                                                             tGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAgg                                                                          <  1:3927162/150‑1 (MQ=255)
                                                                                    aTCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATg                                                                   <  2:4680830/150‑1 (MQ=255)
                                                                                      cGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGtt                                                                 <  1:3464685/150‑1 (MQ=255)
                                                                                       gCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTtc                                                                >  1:3613879/1‑150 (MQ=255)
                                                                                        cATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTtct                                                               >  2:5044165/1‑150 (MQ=255)
                                                                                        cATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTtct                                                               >  1:2335184/1‑150 (MQ=255)
                                                                                         aTACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTtctc                                                              >  2:1181484/1‑150 (MQ=255)
                                                                                           aCAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCACTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCg                                                            <  1:1324023/150‑1 (MQ=255)
                                                                                            cAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGTTACCGGACTTCCGGCGCTAAGTCAGATACCTTTAACAACCGCCGGATCCCGTTGCTGGACTGTGGCAGGAGAGAGGTGCTcctc                                                           >  2:3240512/1‑148 (MQ=255)
                                                                                            cAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGc                                                           >  1:305207/1‑150 (MQ=255)
                                                                                            cAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCACTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGc                                                           >  1:3301551/1‑150 (MQ=255)
                                                                                              gTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAg                                                                           <  2:4965443/132‑1 (MQ=255)
                                                                                              gTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCa                                                         <  2:4579684/150‑1 (MQ=255)
                                                                                              gTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCACTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATACATGTTCTCCGCCa                                                         >  2:3970678/1‑150 (MQ=255)
                                                                                                  ctcACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAg                                                     <  2:836521/150‑1 (MQ=255)
                                                                                                  ctcACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAg                                                     <  2:3189583/150‑1 (MQ=255)
                                                                                                   tcACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCACTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAGa                                                    >  1:430367/1‑150 (MQ=255)
                                                                                                    cACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAGaa                                                   <  1:1031432/150‑1 (MQ=255)
                                                                                                    cACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAGaa                                                   <  1:3652018/150‑1 (MQ=255)
                                                                                                    cACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCACTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAGaa                                                   >  1:2199989/1‑150 (MQ=255)
                                                                                                      cccTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAGAATa                                                 <  1:684236/150‑1 (MQ=255)
                                                                                                        ctgctgGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAGAATAGt                                               >  1:819733/1‑150 (MQ=255)
                                                                                                        ctgctgGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAGAATAGt                                               >  2:1432888/1‑150 (MQ=255)
                                                                                                         tgctgGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCCGTAGACTAGtt                                              >  2:4393187/1‑150 (MQ=255)
                                                                                                           ctgGCAAATCCAGCCAATCCCCAGTTCACGGGCGAGTAGGGCACTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTAGTGTGGAAGGATTAATGTTCTCCGCCAGTAGAATAGTTgg                                            >  2:3089858/1‑150 (MQ=255)
                                                                                                             ggCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAGAATAGTTGGCg                                          <  2:2225573/150‑1 (MQ=255)
                                                                                                                aaaTCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAGAATAGTTGGCGAGt                                       <  1:4248880/150‑1 (MQ=255)
                                                                                                                 aaTCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAGAATAGTTGGCGAGtt                                      <  1:1055397/150‑1 (MQ=255)
                                                                                                                  aTCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAGAATAGTTGGCGAGTTa                                     <  2:3667994/150‑1 (MQ=255)
                                                                                                                     cAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAGAATAGTTGGCGAGTTAAAc                                  <  1:753785/150‑1 (MQ=255)
                                                                                                                        ccAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAGCATAGTTGGCGAGTTAAACTGc                               >  1:3692010/1‑150 (MQ=255)
                                                                                                                           aTCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAGAATAGTTGGCGAGTTAAACTGCggg                            <  2:4124719/150‑1 (MQ=255)
                                                                                                                            tCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAGAATAGTTGGCGAGTTAAACTGCGGGa                           <  2:1616955/150‑1 (MQ=255)
                                                                                                                             ccccAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAGAATAGTTGGCGAGTTAAACTGCGGGAg                          <  2:1299735/150‑1 (MQ=255)
                                                                                                                                cAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAGAATAGTTGGCGAGTTAAACTGCGGGAGttc                       >  1:945606/1‑150 (MQ=255)
                                                                                                                                  gTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAGAATAGTTGGCGAGTTAAACTGCGGGAGttctt                     <  2:60278/150‑1 (MQ=255)
                                                                                                                                   ttCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAGAATAGTTGGCGAGTTAAACTGCGGGAGttcttc                    >  1:4007734/1‑150 (MQ=255)
                                                                                                                                   ttCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAGAATAGTTGGCGAGTTAAACTGCGGGAGttcttc                    >  1:2425089/1‑150 (MQ=255)
                                                                                                                                     cACGGGCGATTAGGGCACTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAGAATAGTTGGCGAGTTAAACTGCGGGAGTTcttctt                  >  2:3384438/1‑150 (MQ=255)
                                                                                                                                       cGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACACCCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAGAATAGTTGGCGAGTTAAACTGCGGGAGTTCTTCTTTc                <  1:4449950/150‑1 (MQ=255)
                                                                                                                                           cGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAGAATAGTTGGCGAGTTAAACTGCGGGAGTTCTTCTTTCGttt            <  2:146824/150‑1 (MQ=255)
                                                                                                                                            gATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAGAATAGTTGGCGAGTTAAACTGCGGGAGTTCTTCTTTCGTTTg           >  2:1712782/1‑150 (MQ=255)
                                                                                                                                               tAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAGAATAGTTGGCGAGTTAAACTGCGGGAGTTCTTCTTTCGTTTGGGt        >  1:176324/1‑150 (MQ=255)
                                                                                                                                                 gggCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAGAATAGTTGGCGAGTTAAACTGCGGGAGTTCTTCTTTCGTTTGGGTCa      <  2:4465114/150‑1 (MQ=255)
                                                                                                                                                     gcTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAGAATAGTTGGCGAGTTAAACTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTg  <  2:3153718/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AACGATGCCATCCGAACAGTGGCTTAACCCTGACGGTTGAAACGTTGCGTTTTAACGTCCAGCGTTAGCGTTTCTTCTGGTTGTATCGCATACAGTTTCTCACCCTGCTGGCAAATCCAGCCAATCCCCAGTTCACGGGCGATTAGGGCGCTGTGGGATACCGGACTTCCGGCGCTAAGGCAGATACCTTTTACAACCGCCGGATCCAGTTGCAGTACTGTGGAAGGATAAATGTTCTCCGCCAGTAGAATAGTTGGCGAGTTAAACTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTG  >  CP009273/1243129‑1243427

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: