Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 709,827 G→T 43.7% D272Y (GAT→TAT)  pgm → phosphoglucomutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273709,8270GT43.7% 280.9 / inf 292D272Y (GAT→TAT) pgmphosphoglucomutase
Reads supporting (aligned to +/- strand):  ref base G (77/87);  new base T (66/61);  total (144/148)
Fisher's exact test for biased strand distribution p-value = 4.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 294 total aligned reads displayed.

CGCGATTCAGAAAGCGGGCCTGACGCTGGGCGTTGATCCGCTGGGCGGTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGACCGTCACGGTATCGTCACTCCGGCAGGTTTGATGAA  >  CP009273/709679‑709972
                                                                                                                                                    |                                                                                                                                                 
   gATTCAGAAAGCGGGCCTGACGCTGGGCGTTGATCCGCTGGGCGGTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTTATaa                                                                                                                                               <  2:3644706/150‑1 (MQ=255)
     ttCAGAAAGCGGGCCTGACGCTGGGCGTTGATCCGCTGGGCGGTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAg                                                                                                                                             <  2:2264637/150‑1 (MQ=255)
      tCAGAAAGCGGGCCTGACGCTGGGCGTTGATCCGCTGGGCGGTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGa                                                                                                                                            >  1:1829615/1‑150 (MQ=255)
       cAGAAAGCGGGCCTGACGCTGGGCGTTGATCCGCTGGGCGGTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTTATAAAGAc                                                                                                                                           <  1:3988054/150‑1 (MQ=255)
       cAGAAAGCGGGCCTGACGCTGGGCGTTGATCCGCTGGGCGGTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGAc                                                                                                                                           <  1:3837632/150‑1 (MQ=255)
        aGAAAGCGGGCCTGACGCTGGGCGTTGATCCGCTGGGCGGTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTTATAAAGACg                                                                                                                                          >  1:127933/1‑150 (MQ=255)
         gAAAGCGGGCCTGACGCTGGGCGTTGATCCGCTGGGCGGTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACgg                                                                                                                                         >  2:242621/1‑150 (MQ=255)
          aaaGCGGGCCTGACGCTGGGCGTTGATCCGCTGGGCGGTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGgc                                                                                                                                        <  2:1570690/150‑1 (MQ=255)
             gCGGGCCTGACGCTGGGCGTTGATCCGCTGGGCGGTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGcgcg                                                                                                                                     <  1:2777912/150‑1 (MQ=255)
             gCGGGCCTGACGCTGGGCGTTGATCCGCTGGGCGGTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGcgcg                                                                                                                                     <  1:3048282/150‑1 (MQ=255)
               gggCCTGACGCTGGGCGTTGATCCGCTGGGCGGTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTTATAAAGACGGCGCGAt                                                                                                                                   >  1:4553345/1‑150 (MQ=255)
               gggCCTGACGCTGGGCGTTGATCCGCTGGGCGGTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTTATAAAGACGGCGCGAt                                                                                                                                   <  1:3090641/150‑1 (MQ=255)
               gggCCTGACGCTGGGCGTTGATCCGCTGGGCGGTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGAt                                                                                                                                   <  2:4013703/150‑1 (MQ=255)
                   cTGACGCTGGGCGTTGATCCGCTGGGCGGTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGt                                                                                                                               >  2:2993695/1‑150 (MQ=255)
                     gACGCTGGGCGTTGATCCGCTGGGCGGTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTTATAAAGACGGCGCGATCCGTAt                                                                                                                             >  2:3752296/1‑150 (MQ=255)
                       cGCTGGGCGTTGATCCGCTGGGCGGTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATgg                                                                                                                           <  2:3546520/150‑1 (MQ=255)
                         cTGGGCGTTGATCCGCTGGGCGGTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGAc                                                                                                                         <  1:3155718/150‑1 (MQ=255)
                            ggCGTTGATCCGCTGGGCGGTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGc                                                                                                                      >  2:2675472/1‑150 (MQ=255)
                                 tGATCCGCTGGGCGGTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGctcctc                                                                                                                 <  1:1546599/150‑1 (MQ=255)
                                 tGATCCGCTGGGCGGTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGctcctc                                                                                                                 <  2:3809082/150‑1 (MQ=255)
                                   aTCCGCTGGGCGGTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCg                                                                                                               <  2:1844956/150‑1 (MQ=255)
                                   aTCCGCTGGGCGGTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCg                                                                                                               >  2:1667165/1‑150 (MQ=255)
                                   aTCCGCTGGGCGGTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCg                                                                                                               <  2:515229/150‑1 (MQ=255)
                                    tCCGCTGGGCGGTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTTATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGa                                                                                                              >  2:1358045/1‑150 (MQ=255)
                                      cGCTGGGCGGTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGCTCAGGTCGATCAAACCTTCCGCTTTATGCACCTTTATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAgt                                                                                                            <  2:2066165/150‑1 (MQ=255)
                                       gCTGGGCGGTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTTATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAgtg                                                                                                           <  2:2744014/150‑1 (MQ=255)
                                       gCTGGGCGGTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTTATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAgtg                                                                                                           <  2:4166690/150‑1 (MQ=255)
                                       gCTGGGCGGTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAgtg                                                                                                           <  1:3299455/150‑1 (MQ=255)
                                          gggCGGTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTTATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGc                                                                                                        <  2:2448023/150‑1 (MQ=255)
                                          gggCGGTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGc                                                                                                        <  2:1986265/150‑1 (MQ=255)
                                              ggTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATg                                                                                                    >  2:954859/1‑150 (MQ=255)
                                              ggTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTAACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTTATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATg                                                                                                    >  1:3197500/1‑150 (MQ=255)
                                               gTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATgg                                                                                                   <  1:74429/150‑1 (MQ=255)
                                                 tCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTAGCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAATGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGCTGGCg                                                                                                 >  2:3292088/1‑150 (MQ=255)
                                                  ccGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCgg                                                                                                >  2:3095498/1‑150 (MQ=255)
                                                    ggTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGc                                                                                              >  1:3514544/1‑150 (MQ=255)
                                                         cGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTTATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCctgct                                                                                         <  1:2485691/150‑1 (MQ=255)
                                                          gAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTTATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCctgctg                                                                                        >  1:362060/1‑150 (MQ=255)
                                                          gAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCctgctg                                                                                        <  2:1919908/150‑1 (MQ=255)
                                                          gAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCctgctg                                                                                        <  2:3046736/150‑1 (MQ=255)
                                                           aaTACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTTATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTgg                                                                                       >  1:4181228/1‑150 (MQ=255)
                                                             tACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCa                                                                                     <  1:703959/150‑1 (MQ=255)
                                                                tGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTTATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTg                                                                                  >  1:624131/1‑150 (MQ=255)
                                                                 ggAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTTATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGc                                                                                 >  1:579920/1‑150 (MQ=255)
                                                                      cGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGAt                                                                            <  2:732258/150‑1 (MQ=255)
                                                                       gTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTTATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATa                                                                           >  2:297990/1‑150 (MQ=255)
                                                                       gTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTTATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATa                                                                           >  1:4198916/1‑150 (MQ=255)
                                                                        tATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATaa                                                                          <  1:2653354/150‑1 (MQ=255)
                                                                         aTTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAg                                                                         <  1:1037956/150‑1 (MQ=255)
                                                                         aTTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAg                                                                         >  1:239279/1‑150 (MQ=255)
                                                                            ggCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTc                                                                      <  2:4364503/150‑1 (MQ=255)
                                                                              cGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGa                                                                    <  2:3232091/150‑1 (MQ=255)
                                                                                   aTTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAATCCTTCCGCTTTATGCACCTTTATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTgg                                                               >  2:2537913/1‑150 (MQ=255)
                                                                                     tACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTTATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCg                                                             >  1:3360745/1‑150 (MQ=255)
                                                                                       caacctcaacctGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTTATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGtt                                                           <  1:2961683/150‑1 (MQ=255)
                                                                                       caacctcaacctGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGtt                                                           <  2:1701348/150‑1 (MQ=255)
                                                                                        aacctcaacctGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTAATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGttt                                                          >  1:2498094/1‑150 (MQ=255)
                                                                                          cctcaacctGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTTATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGc                                                        >  1:562591/1‑150 (MQ=255)
                                                                                           ctcaacctGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTTATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCt                                                       <  1:349550/150‑1 (MQ=255)
                                                                                              aacctGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAAc                                                    <  1:496871/150‑1 (MQ=255)
                                                                                               acctGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTTATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACg                                                   >  2:983108/1‑150 (MQ=255)
                                                                                                  tGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGAcc                                                <  2:4623089/150‑1 (MQ=255)
                                                                                                      tATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTTATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGa                                            <  2:1529758/150‑1 (MQ=255)
                                                                                                      tATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTTATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGa                                            <  2:2605843/150‑1 (MQ=255)
                                                                                                      tATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTTATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGa                                            >  1:2028396/1‑150 (MQ=255)
                                                                                                       aTCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTTATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGAt                                           <  2:176438/150‑1 (MQ=255)
                                                                                                       aTCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTTATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGAt                                           >  1:2874524/1‑150 (MQ=255)
                                                                                                        tCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTTATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGAtt                                          <  1:3981525/150‑1 (MQ=255)
                                                                                                         cGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTa                                         >  1:4429296/1‑150 (MQ=255)
                                                                                                          gTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTAt                                        >  2:4572855/1‑150 (MQ=255)
                                                                                                           ttAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATg                                       <  1:854490/150‑1 (MQ=255)
                                                                                                            tAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGa                                      >  2:1734801/1‑150 (MQ=255)
                                                                                                                 aTCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGACCGTc                                 >  2:609890/1‑150 (MQ=255)
                                                                                                                 aTCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGACCGTc                                 >  2:1343529/1‑150 (MQ=255)
                                                                                                                   cAGGTCGATCAAACCTTCCGCTTTATGCACCTTTATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGACCGTCAc                               >  2:3759993/1‑150 (MQ=255)
                                                                                                                   cAGGTCGATCAAACCTTCCGCTTTATGCACCTTTATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGACCGTCAc                               <  2:1155301/150‑1 (MQ=255)
                                                                                                                   cAGGTCGATCAAACCTTCCGCTTTATGCACCTTTATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGACCGTCAc                               <  1:925221/150‑1 (MQ=255)
                                                                                                                        cGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGACCGTCACGGTAt                          <  1:3860738/150‑1 (MQ=255)
                                                                                                                        cGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGACCGTCACGGTAt                          <  2:3457465/150‑1 (MQ=255)
                                                                                                                         gATCAAACCTTCCGCTTTATGCACCTTTATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGACCGTCACGGTATc                         <  1:1117472/150‑1 (MQ=255)
                                                                                                                          aTCAAACCTTCCGCTTTATGCACCTTTATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGACCGTCACGGTATCg                        >  1:1610083/1‑150 (MQ=255)
                                                                                                                          aTCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGACCGTCACGGTATCg                        >  2:3147147/1‑150 (MQ=255)
                                                                                                                            cAAACCTTCCGCTTTATGCACCTTTATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGACCGTCACGGTATCGTc                      >  2:796769/1‑150 (MQ=255)
                                                                                                                             aaaCCTTCCGCTTTATGCACCTTTATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGACCGTCACGGTATCGTCa                     <  1:4592407/150‑1 (MQ=255)
                                                                                                                             aaaCCTTCCGCTTTATGCACCTTTATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGACCGTCACGGTATCGTCa                     >  2:125294/1‑150 (MQ=255)
                                                                                                                             aaaCCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGACCGTCACGGTATCGTCa                     <  1:3655230/150‑1 (MQ=255)
                                                                                                                              aaCCTTCCGCTTTATGCACCTTTATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGACCGTCACGGTATCGTCAc                    >  1:1340217/1‑150 (MQ=255)
                                                                                                                                ccTTCCGCTTTATGCACCTTTATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGACCGTCACGGTATCGTCACTc                  <  2:2073029/150‑1 (MQ=255)
                                                                                                                                 cTTCCTCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGACCGTCACGGTATCGTCACTcc                 >  2:4249654/1‑150 (MQ=255)
                                                                                                                                  ttCCGCTTTATGCACCTTTATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGACCGTCACGGTATCGTCACTCCg                <  2:3677290/150‑1 (MQ=255)
                                                                                                                                    ccGCTTTATGCACCTTTATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGACCGTCACGGTATCGTCACTCCGGc              >  2:4094639/1‑150 (MQ=255)
                                                                                                                                     cGCTTTATGCACCTTTATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGACCGTCACGGTATCGTCACTCCGGCa             <  1:2184432/150‑1 (MQ=255)
                                                                                                                                     cGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGACCGTCACGGTATCGTCACTCCGGCa             >  1:4335927/1‑150 (MQ=255)
                                                                                                                                     cGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGACCGTCACGGTATCGTCACTCCGGCa             >  1:4335873/1‑150 (MQ=255)
                                                                                                                                      gCTTTATGCACCTTTATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGACCGTCACGGTATCGTCACTCCGGCAg            <  1:1505305/150‑1 (MQ=255)
                                                                                                                                      gCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGACCGTCACGGTATCGTCACTCCGGCAg            >  1:2426516/1‑150 (MQ=255)
                                                                                                                                       cTTTATGCACCTTTATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGACCGTCACGGTATCGTCACTCCGGCAgg           >  1:1940304/1‑150 (MQ=255)
                                                                                                                                          tATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGACCGTCACGGTATCGTCACTCCGGCAGGttt        >  2:812908/1‑150 (MQ=255)
                                                                                                                                             gCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGACCGTCACGGTATCGTCACTCCGGCAGGTTtgat     <  2:348188/150‑1 (MQ=255)
                                                                                                                                               aCCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGACCGTCACGGTATCGTCACTCCGGCAGGTTtgatga   >  2:2850355/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CGCGATTCAGAAAGCGGGCCTGACGCTGGGCGTTGATCCGCTGGGCGGTTCCGGTATCGAATACTGGAAGCGTATTGGCGAGTATTACAACCTCAACCTGACTATCGTTAACGATCAGGTCGATCAAACCTTCCGCTTTATGCACCTTGATAAAGACGGCGCGATCCGTATGGACTGCTCCTCCGAGTGTGCGATGGCGGGCCTGCTGGCACTGCGTGATAAGTTCGATCTGGCGTTTGCTAACGACCCGGATTATGACCGTCACGGTATCGTCACTCCGGCAGGTTTGATGAA  >  CP009273/709679‑709972

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: ‑